Identification of colored wheat genotypes with suitable quality and yield traits in response to low nitrogen input

Colored wheat is a valuable resource that is rich in anthocyanins and minerals and thus contributes additional nutritional value to a healthy human diet. However, the effects of nitrogen fertilization on anthocyanin content (AC) and the balance between quality and yield still merit discussion. In this study, blue, purple, and common-colored wheat genotypes were used to investigate three nutrient quality traits, seven processing quality traits, three yield traits and seven grain morphology traits at three nitrogen levels in two years to excavate their possible plasticity under low-nitrogen stress and the tradeoffs among these traits. The highest AC was found in the blue genotypes followed by the purple genotypes. Analysis of variance (ANOVA) showed that AC could be significantly increased by reducing N application, especially in the purple genotypes. Therefore, growing colored wheat with low nitrogen input could allow efficient harvesting of grain with higher AC. However, the other nutrient quality traits and most processing quality traits were observed to decrease under low-nitrogen (LN) stress. Additionally, a correlation analysis indicated that the nutrient quality traits had stable tradeoffs with thousand kernel weight at all N levels because of the significantly negative correlations among them. Therefore, the additive main effect and multiplicative interaction (AMMI) model was used to further identify the most suitable colored genotypes with the best yield potential and also nutrient quality relative characteristics under LN stress. The blue lines Lanmai2999 and purple varieties Zhongkezinuomai 168 were found to be specifically adapted to LN stress with the highest AC values and showed stable performance in the other nutrient quality- and yield-related features. To further investigate the possible mechanism of anthocyanin accumulation in response to reduced N application, the expression of four genes (TaCHS, TaFDR, TaCHI and TaANS) involved in the anthocyanin synthesis pathway was evaluated. All four genes were downregulated under high nitrogen fertilizer application, indicating that anthocyanin synthesis in colored wheat might be inhibited by nitrogen fertilizer. Therefore, this research provided information for optimizing nitrogen fertilizer management in producing colored wheat and also demonstrated that it is efficient and economical to plant colored wheat genotypes in nitrogen-poor areas for use in a healthy human diet, improving the benefits of wheat planting and facilitating nitrogen pollution control.

Introduction Wheat (Triticum aestivum L.) is one of the most produced food crop on earth and plays an essential role in the human diet. Approximately 750 million tons are produced annually, together with maize and rice, forming the main staple food crop for 35% of the world's population (FAO, 2018). Thus, quality and yield improvement are always the main targets of wheat breeders. In particular, traits related to processing and end-use quality, such as kernel hardness, test weight, water absorption, Zeleny sedimentation score, and wet gluten content, which directly or indirectly determine the commodity value of wheat, are the crucial optimization objectives in the wheat breeding process [1,2].
Additionally, with the recent popularization of the concept of a healthy diet, in addition to food prices, which are mainly decided by production amounts and processing costs, people have become increasingly concerned about health and nutrition. Thus, the desire to "eat healthy" has become one of the most important factors influencing consumer food choices, and the nutrient quality traits in wheat, such as nutritive element content, have also gradually attracted breeders' attention [3,4,5]. Anthocyanin is a plant flavonoid pigment commonly found in traditional foods such as fruits, vegetables and grains, giving them their natural colors [6]. The primary ecophysiological function of anthocyanin is to attract seed dispersal agents and to protect plant tissues against abiotic and biotic stressors such as high radiation or pathogens [7]. Dietary intake of anthocyanins has profound health-promoting effects [8,9], especially considering the plentiful evidence showing their contributions to human health as antioxidant, anti-inflammatory, antidiabetic, cardiovascular diseases, and anticancer agents, among other benefits [10]. Thus, food rich in natural anthocyanins is increasingly popular for a healthy diet.
Colored wheat, which is rich in anthocyanin content (AC) in the pericarp (purple wheat) or aleurone layer (blue wheat) compared with that of common-colored wheat (white or red wheat) [11,12], can play a significant role in the prevention of various diseases associated with oxidative stress [13,14,15]. However, colored wheat germplasm resources, especially those with both optimal quality-and yield-related traits and availability for large-scale planting, are still lacking [5]. Thus, breeding colored wheat with both superior quality and superior yield features has become a popular objective to increase the value of processed wheat products.
In addition to genetic factors, abiotic environmental factors such as extreme photothermal characteristics, oxygen concentration and even fertilizer levels can affect anthocyanin metabolism [16]. For example, nitrogen (N) limitation can significantly promote pigmentation [17].
N is a major element critical for plant growth and affects both yield and quality formation [18]. In general, an increased N supply drives higher productivity and higher grain protein content [19], also improves some common processing and end-use relative features, such as wet gluten content, Zeleny sedimentation value, test weight and so on [20,21,22,23]. However, excessive N fertilizer input can cause environmental pollution through N leaching and runoff, causing eutrophication of freshwater and marine ecosystems and N 2 O (a greenhouse gas) emission associated with denitrification by soil bacteria [19]. Thus, wheat with high tolerance to low N (LN) stress is an environmentally friendly germplasm resource.
However, colored wheat with a high commodity value that is adapted to LN stress and their corresponding N management still needs to be developed. Therefore, the present study quantified the genetic variation of 39 colored wheat germplasm lines with respect to nutrient quality traits, processing quality traits, grain morphology traits and yield conmponenets at different N levels, with the aim of (1) identifying colored wheat lines adapted to planting under LN stress; (2) investigating the effects of N input on grain color and other quality traits; and (3) evaluating the relationship between grain color and yield related traits at different N input levels.

Plant materials and treatments
Thirty-nine colored wheat lines (S1 Table), including color cultivars widely used in breeding program in China (such as Lanlimai and Luozzhen No.1), as well as varieties released in Southwest China in recent years (such as Chuanmai 42, Chuanyu 16, Mianmai 45 and Zhongkemai 138 and so on) and some intermediate lines bred by the Chengdu Institute of Biology, Chinese Academy of Sciences (CIBCAS), were selected as representatives of colored wheat: 10 blue, 22 purple and 7 common color (Fig 1).
The field experiments were conducted at Chengdu Plain Agricultural Ecology Test Station of Chinese Academy of Sciences, China (N 30˚34 0 , E 103˚52 0 ) in 2016-2017 and 2017-2018. This test station is specifically provided to construct crop cultivation and breeding experiments by Chengdu Institute of Biology, Chinese Academy of Sciences. And we confirm that these field research in this study did not involve any endangered or protected planting species. Each year, soil samples were taken randomly at depths of 0-30 cm using a soil corer before sowing and submitted for analysis of available N (S2 Table). The wheat lines were planted in a randomized complete block using a split-plot design with three replications; the main plots consisted of three levels of N fertilizer application (low N, LN; medium N, MN; and high N, HN), and the subplot treatment was the 39 colored wheat cultivars. The subplot size was 2.0×1.0 m, containing 4 rows with 2 m long and 0.25 m apart. In the HN field, 120 kg N ha -1 was applied after sowing in both years. In MN, 90 kg N ha -1 was applied after sowing in both years. In the LN treatment, no nitrogen was applied during the entire growing period. Crop management followed local conventional agricultural practices, including disease and pest management, chemical weed control and other agronomic practices.

Phenotypic evaluation and measurement
Three nutrient quality traits of the wheat grains, including AC, N concentration (NC) and protein content (GPC), were measured. Among them, the total AC was determined according to a previously published method [15]: a mixture of methanol acidified with 1 N HCl (85:15; v/v, 1 mL) was added to the samples (~100 mg) and then incubated at 4˚C for 12 h in the dark. After centrifugation at 12 000 rpm for 10 min at 4˚C, the supernatants were transferred into an ELISA plate. The absorbance of each sample was assessed at 510 nm using a Varioskan LUX Multimode Microplate Reader. A solution of methanol acidified with 1 N HCl was used as a blank. The total AC was corrected for dry matter and expressed as cyanidin 3-O-glucoside equivalent. The NC of the plant components was analyzed according to the following method: grains were collected at harvest and oven-dried to a constant weight at 60˚C. The dried samples were first milled (<0.20 mm) using a mixer-grinder mill with mechanical modification for small-quantity samples for the measurement of N content [24]. The milled samples were then redried at 80˚C for 48 h and weighed to 0.00001 g precision, encapsulated in tin capsules and analyzed for N%. The total N content was measured using an automatic N analyzer (BUCHI AutoKjeldahl Unit K-370; BUCHI Laboratory Equipment, Flawil, Switzerland). The protein content of the colored grains was determined according to the method of Hurkman and Tanaka [25].
Eight processing quality traits, including wet gluten content (WGC), Zeleny sedimentation value (ZEL), grain hardness (GH), test weight (TW), water absorption (WA), formation time (FT), dough stabilization time (DST) and maximum tensile resistance (MTR), were measured by near-infrared reflectance spectroscopy (NIRS) with a Perten DA-7250 instrument (Perten Instruments, Huddinge, Sweden) and expressed on the basis of 14% moisture. The measurements were calibrated using calibration samples according to the manufacturer's instructions. Previous studies have confirmed that NIRS is a powerful method for measuring the quality traits of wheat; that it is rapid, effective, and convenient compared with traditional methods; and that it is appropriate for use in wheat breeding programs designed to improve qualityrelated traits [23].
Three yield components were investigated at maturity, including thousand-kernel weight (TKW), kernels per spike (KPS) and tiller number (TN). For each line in different replicates, six plants were randonmly chosen to measure TKW, KPS and TN. TKW were investigated using the Analysis System of WSeen SC-G Instrument (Zhejiang, China) by weighing more than 200 random kernels with two technical repeats; KPS were calculated from main and integrated tillers and TN were determined by the spike number harboring mature grains.

RNA extraction and qRT-PCR
Total RNA was isolated at 25 days after pollination from the grains of the purple wheat line Zhongkezinuomai168 (ZKZNM168), which is famous because of its high yield and was released by CIBCAS. TRIzol was used for RNA extraction in this study. After treatment with RNase-free DNase I (Promega), the total RNA (1 μg) was reverse transcribed using the Trans-Script II One-Step gDNA Removal and cDNA Synthesis SuperMix kits (TransGen Biotech). The reverse transcription products were used as templates for quantitative real-time PCR performed on a CFX96 Real-Time PCR system (Bio-Rad) using SYBR Premix Taq (TransGen Biotech) according to the manufacturer's protocol. GAPDH was used as an internal control gene, and the relative expression levels were measured using the 2 -ΔΔCt analysis method. The primers used in RT-qPCR are listed in S3 Table.

Statistical data analysis
One-way analysis of variance (ANOVA) was conducted to evaluate the variance, and Spearman's rank correlation was performed to analyze phenotypic correlation coefficients between traits and significance between groups using SPSS 20.0 (SPSS, Chicago, IL, USA). The additive main effect and multiplicative interaction (AMMI) analysis was conducted by Genstat 19.0 (VSN International Ltd, UK).

Correlation analysis at different N levels
The correlation analysis among three wheat nutrient quality traits, eight processing quality traits, three yield-related traits and seven grain morphology traits were evaluated under LN, MN and HN conditions (Table 1, S4 Table, S5 Table) Overall, all three nutrient quality traits (AC, NC and GPC) were significantly positively correlated with each other, regardless of N treatment, indicating a possible common genetic background controlling these traits (Table 1, S4 Table, S5 Table).
Among the processing quality traits, only WGC, ZEL, FT and DST had significantly positive correlations with three nutrient quality traits (AC, NC and GPC), not affected by N treatment. However, apart from them, the significance of the correlation coefficients between the other processing quality traits and nutrient quality traits was not stable under the different N treatments. For example, TW was obviously correlated with AC under LN, with AC and NC under MN, and with all four nutrient traits under HN (Table 1, S4 Table, S5 Table).
For the relationships between the three yield-related traits and the nutrient traits, only TKW showed a consistently significant correlation with all four nutrient traits at all N levels, that is, negative relationships with AC, NC and GPC. These results indicated that grain weight components such as grain size and morphology traits might also be related to grain nutrient traits (Table 1, S4 Table, S5 Table). For the correlation of yield components with eight processing traits, only TKW-GH and TN-TW show stale significance under all N inputs, indicating the yield-quality tradeoffs were mainly existed between TKW and three nutrient traits.
To further specialize the contribution of the grain-related traits to these nutrient quality traits, seven grain size and shape traits were investigated. In detail, GW was significantly related to all three nutrient quality traits at all N levels, similar to TKW, showing negative relationships with AC, NC and GPC (Table 1, S4 Table, S5 Table). Thus, the relationship between TKW and nutrient traits might be mainly affected by GW. However, GL was always significantly positively related to AC at all three N levels, providing a possibility to simultaneously improve grain size and AC by finding their common genetic basis.

Performance comparison of colored wheat lines
An ANOVA to analyze the different performance traits of blue, purple and common-colored wheat was conducted across all N levels. In total, the three kinds of wheat with different colored grains could be distinguished in regard to most grain features ( Table 2).
In detail, for the nutrient quality traits, AC, NC and GPC were all higher in blue wheat than in the other wheat colors. In addition, purple wheat also had obviously higher AC, NC and GPC than common-colored wheat ( Table 2).
For the processing quality traits, blue wheat also showed significantly higher values than purple and common-colored wheat for the measured processing quality traits, except for GH ( Table 2). The numbers which were highlighted by blue color indicates significance at the level of 0.05.
The numbers which were highlighted by purple color indicates significance at the level of 0.01. AC anthocyanin content, NC nitrogen concentration, For grain size/shape, it was not easy to distinguish colored wheat from common-colored wheat based on only their similar grain size parameters, such as GA and GP. However, blue wheat had the longest grains, and purple wheat had the shortest grains, and overall, the grain shapes were obviously different among blue, purple and common-colored wheat, that is, blue wheat had the largest LWR and purple wheat had the roundest grains (GR) ( Table 2).
For the three yield components, only TKW showed differences among the three kinds of wheat, with common-colored wheat > purple wheat > blue wheat, possibly due to differences in their grain shapes. However, their KPS were nearly the same ( Table 2).

Effects of different nitrogen fertilizer treatments on colored wheat lines
Overall, nitrogen application remarkably affected wheat nutrient quality ( Table 2). Under LN stress, AC was higher than that in MN and HN conditions. However, NC and GPC were positively related to N fertilizer application levels, being highest under HN and lowest under LN. Among the processing quality traits, WGC, ZEL, GH, FT and DST were also sensitive to LN stress and showed increasing values with increasing N application. However, for grain size/ shape traits and yield-related traits, only GA and TN were significantly affected by different N levels. These results indicated that quality traits, especially nutrient quality traits, might have more plasticity with respect to N levels than do grain weight and morphology traits, which have been proven to have higher heritability. And thus it provided an idea to select sepcialized and optimized N managements to improve quality of wheat with stably superior grain weight perfomance, with the aim to balance yield potential and quality characteristics. Based on the results of investigating these three different colored wheat groups, all showed trends similar to the overall trends for all the nutrient quality traits and most of the other traits ( Table 3). It is interesting that the AC in blue and purple wheat did not significantly decrease when N application was increased beyond MN, while their NC and GPC values both obviously increased, indicating that in colored wheat, the MN level might be enough to synthesize not only AC but also NC and GPC, leading to more economical and efficient N use.
To further investigate the tolerance of different colored wheat lines to low N stress to futher facilitate N-input restriction and reduce N pollution, the amplitude of the variation was calculated as follows: the difference between each measured value from LN and HN was estimated, and these differences were used to obtain the percentage of each original measured value based on the HN condition, defined as (LN-HN)/HN×100%. The results showed that AC in purple wheat was more sensitive to LN stress (increased 58.18%) and accumulated more under LN conditions than blue wheat (increased 47.11%), although its GPC also decreased slightly but not significantly ( Table 4).

Identification of colored wheat lines suitable for LN stress
As shown in Tables 2 and 3, AC appears to accumulate when N application is reduced; therefore, colored wheat, is more suitable for planting under LN conditions to maximize its inherently high AC. To further identify specific genotypes adapted to LN stress with stable quality and yield potential, the AMMI model was used to evaluate the stability and adaptation of the blue and purple wheat lines for different N applications.
For the AC in blue wheat genotypes (Fig 2, Table 5), the lowest AMMI stability values (ASVs) were observed for Lanmai3707, Lanmai2909 and Lanmai3471, which indicated a higher stability of these genotypes under all N treatments; however, their AC values were relatively low. Higher AC values were observed in Lanmai2999, Lanlimai-3 and Lanlimai-2, in decreasing order. The highest ASV belonged to Lanmai2999, followed by Lanlimai and Lanlimai-2, indicating that these genotypes were more specifically adapted to certain N levels. Fig 2  depicts genotype-environment affinity as orthogonal projections of the genotypes on the environmental vectors to identify the best genotypes with respect to N levels. Lanmai2999 and Lanlimai-2 were better adapted to both LN and MN than the other genotypes, and they also had higher AC values.
The other two nutrient quality traits (NC and GPC) were positively correlated with AC, while three yield components (TKW, TN and KPS) were negatively correlated with AC. All five traits were generally depressed by reducing N application. Thus, their stability and adaptation to LN stress were further considered to identify the most suitable genotypes with not only high AC but also satisfactory quality and yield features.
Between Lanmai2999 and Lanlimai-2, which were adapted to LN stress with high AC, Lan-mai2999 showed lower ASVs for NC, GPC, TKW, TN and KPS than did Lanlimai-2, indicating that these traits in Lanmai2999 are more stable under all N levels (Table 5). Similarly, the other processing traits are also more stable in Lanmai2999, compared with Lanlimai-2 (S6   (Fig 3, Table 6), ZKZNM168, Mianzimai 828 and Zimai1769 had higher AC across all N levels. The lowest ASV was observed in Zimai1769, which indicated that its high AC could be stable across all N levels. The highest ASV belonged to ZKZNM168, followed by Zimai1487 and Mianzimai 828, indicating that they are unstable with respect to N treatment. According to Fig 3, ZKZNM168, showing the highest AC, was the best adapted to environmental LN stress. In combination with other quality-and yield-related traits, ZKZNM168 also showed relatively low ASVs and moderate measured values (Table 6, S7  Table, S2 Fig), indicating that ZKZNM168 might be the most suitable genotype of the purple genotypes for planting under LN stress due to its high AC and stable performance for other quality and yield traits.

Expression of key genes for anthocyanin synthesis in wheat
To study the influence of nitrogen levels on AC in colored wheat, we detected the expression of four key genes involved in anthocyanin synthesis in ZKZNM168 (a purple wheat variety released by CIBCAS), by qRT-PCR. When tested at approximately 25 days after wheat pollination, the results showed that these four key genes (TaCHS, TaFDR, TaCHI and TaANS) were all downregulated when more nitrogen fertilizer was applied (Fig 4). This result was consistent with the previous result that AC in wheat decreased with increasing nitrogen fertilizer application.

Discussion
Wheat is the most widely grown crop and a major source of protein for humans. Commoncolored wheat grains are generally white or red, with very little color in the seed coat. With the improvement of living standards, an increasing pursuit of diversification and nutrition has been observed among breeders [26,27,28]. On the other hand, the benefits of planting wheat in China, especially in Southwest China, have not significantly increased in recent years, probably because of the increasing labor and land costs, which have caused fewer farmers to be willing to plant wheat on a large scale. Therefore, breeding highly nutritious wheat to increase its commercial value is an efficient method to improve wheat planting benefits. Blue and purple wheat are valuable germplasm resources with high levels of anthocyanins [12], which are responsible for oxidation resistance and might increase the nutritional value of wheat flour and processed foods [5,15].
In previous studies, wheat quality traits, such as GPC, have been found to negatively correlate with yield-related traits because of the dilution effect [23,29]. Hence, simultaneously improving quality and yield is always a difficult goal for wheat breeders. Anthocyanin is an unconventional but intriguing nutrient trait in wheat that differs significantly among varieties, highly correlated with grain color [30]. The relationship of AC with familiar quality and yield traits, such as GPC, WGC, MTR, TKW and KPS, deserves investigation, in order to facilitate to breed special usage color wheat varieties with satisfying yield. In this study, the correlation of a total of 20 grain quality-and yield-related traits with AC were investigated, including two nutrient quality traits, eight processing quality traits, three yield-related traits and seven grain morphology traits. The correlation analysis (Table 1, S4 and S5 Tables) demonstrated that AC was positively correlated with most processing quality traits, such as ZEL, WGC, FT and DST. In our opinion, the relation between AC and these processing quality traits might be indirectly controlled by its stable and remarkable positive correlation with NC (correlation coefficient: 0.617~0.681). N is a major element required for most biological compounds, such as protein, gluten and anthocyanins. NC had extremely positive relationships with GPC (0.823~0.947) and WGC (0.806~0.907) and further affected other related processing traits, such as ZEL (influenced by GPC), FT and DST (controlled by GPC and WGC) ( Table 1, S4 and S5 Tables),  in accordance with previous results [20,21,22,23]. The three yield components TKW, TN and KPS are the most important contributors to yield formation. Among these three components, TKW is mainly controlled by genetic factors and has the highest heritability, so it is the most stable trait of these three [31]. It has been found to have a balancing effect between yield traits and quality traits [32]. In this study, AC also had a negative correlation with TKW, TN and KPS and showed an especially stable, significant correlation with TKW (Table 1, S4 and  S5 Tables). Therefore, the relationship of AC with grain size and shape traits, as the major traits affecting TKW, was further investigated. The results indicated that smaller and longer grains might have higher AC, which would need further explanation at the molecular level.
Additionally, N is the crucial element affecting AC and the other quality, besides genetic factors [20,21,22,23] [23,24]. Thus, to high efficiently plant speicial usage color wheat, the optimized N managements could be used to fine tune these processing traits. However, the N pollution caused by excessive N fertilizer application has pushed breeders to select highly N-efficient varieties with optimal tolerance to LN management [33]. Therefore, in consideration of the higher nutrient and commercial value of colored wheat, breeding colored wheat varieties well adapted to LN stress might have both economic and ecological importance.
Previous studies have demonstrated that wheat quality traits such as GPC, ZEL, WGC, GH, FT and DST are affected by N application [20,21,22,23]. In this study, the results were consistent with those of these previous studies (Table 2), that is, most quality traits increased with more N fertilizer application. For yield-related traits, TN and KPS also increased when provided with a sufficient N supply, while TKW and its related morphological traits decreased, possibly due to a tradeoff with TN and KPS.
A previous study found that restricted N fertilization can lead to increased coloration and antioxidant capacity of colored wheat by shifting the plant metabolism towards increased anthocyanin synthesis. Moreover, excessive N fertilization enhances anthocyanin degradation.

PLOS ONE
Screening the suitable colored wheat varieties adapting to nitrogen-poor environment Some studies have also reported that increasing N supply decreases both anthocyanin synthesis and chlorophyll degradation but accelerates starch degradation in fruit [34]. The anthocyanin biosynthesis pathway includes phenylalanine ammonia lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3'-hydroxylase (F3'H), flavonoid 3'5'-hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR), and anthocyanidin synthase/leucoanthocyanidin dioxygenase (ANS/ LDOX). Among them, ANS catalyzes the conversion of leucoanthocyanidin into anthocyanidin, which requires glycosylation, methylation, and acylation to form stable anthocyanin. These modifications are mediated by glucosyltransferase, methyltransferase and acyltransferase [35]. These enzymes are encoded by the corresponding upstream, early, late, and modification genes. The upstream structural genes include PAL and 4CL; the early structural genes include CHS, CHI, F3H, F3'H and F3'5'H; the late structural genes include DFR and ANS; and the modification genes include GT, MT, and AT. The regulatory genes for plant anthocyanin biosynthesis include genes encoding MYB, bHLH and WD40 transcription factors, which control anthocyanin accumulation in the tissues by regulating the expression of structural genes [36,37]. The expression of CHS and DFR in purple wheat showed tissue and cultivar specificity, and DFR had 3 copies; CHS in different cultivars contained a single nucleotide polymorphism [38,39]. The results of this study, together with those of other recent publications ( Fig  4) [40,41], indicate that anthocyanins are produced in a highly regulated, fine-tuned, and nutrient-specific manner as a result of metabolic adaptation to nutrient stress. In particular, anthocyanin-specific genes, namely, TaCHS, TaFDR, TaCHI and TaANS, which are involved in anthocyanin glycosylation and sequestration in the vacuole, are highly expressed in response to nutrient deficiency (Fig 4), indicating that the observed anthocyanin accumulation is not simply the result of an upregulation of the general phenylpropanoid pathway, which would also yield other secondary metabolites, such as flavones [15]. Additionally, anthocyanin-rich fractions may be obtained from appropriate procedures prior to milling [42]. These fractions are also rich in fiber and minerals and could represent suitable ingredients for the production of staple foods such as pasta or bread. Processing whole grains to bake products such as bread, cookies or muffins improves phenolic acid bioaccessibility and bioavailability by releasing some of the insoluble bound phenolic acids and increasing the level of soluble or free phenolic acids [43]. Similarly, we found that AC could accumulate when N application was reduced in blue, purple and common-colored wheat (Table 3). In particular, the AC in purple wheat increased more when the plants were under LN stress (Table 4), but the yield and quality traits were negatively affected. Thus, overall, colored wheat is well adapted to LN stress for anthocyanin biosynthesis but better adapted to MN conditions for the purposes of yield or other quality characteristics.
To further identify colored wheat varieties suitable for LN planting, the AMMI model was used in this study. Finally, the wheat genotypes Lanmai2999 and ZKZNM168 were found to be the most suitable for LN stress, with the highest AC values and good performance in other quality and yield features (Tables 5 and 6, S6 and S7 Tables). However, the yield related traits

PLOS ONE
need further improvement. On the other hand, if regardless of maximizing the AC, this study also could provide some information to select specialized N management for each wheat line to fine tune its producing or processing features. For example, Lanlimai-3 has the highest ZEL but sensitive to N levels (S6 Table). From S1 Fig, it is could deducible that moderate N can be enough for its higher ZEL; and Mianzimai No.1 has the highest hardness (S7 Table) and was adapted for use in low-input management systems to get higher GH (S2 Fig).

Conclusion
This study presented a detailed investigation of the quality and yield-related traits of colored wheat at different N levels to provide a basis for breeding colored wheat with high tolerance to LN stress. Overall, blue wheat had the highest AC across all N levels, while purple wheat accumulated greater AC under LN stress. Additionally, although AC was found to increase when N application was reduced, some other related quality and yield traits were decreased. Thus, the colored wheat is well adapted to LN stress for anthocyanin accumalation but better adapted to MN conditions for the purposes of moderate yield or other quality characteristics. Two varieties each of blue wheat and purple wheat were identified as having the highest comprehensive