The experiment was conducted at the experimental station of Shandong Agricultural University, Tai’an, China (36° 9' N, 117° 9' E; altitude 128m) during the 2018–2019 and 2019–2020 growing seasons. The experimental land was sandy loam with organic matter of 15.8g kg⁻¹, total nitrogen of 1.34g kg⁻¹, nitrate nitrogen of 14.21mg kg⁻¹, ammonium nitrogen of 6.58mg kg⁻¹, available phosphorus of 8.87mg kg⁻¹, and available potassium of 84.62mg kg⁻¹ in the 0–20cm soil layer.

Two typical cultivars of winter wheat currently used in local production, Shannong16 (SN16) and Shannong23 (SN23), were utilized in this study. SN16 is a lodging-sensitive cultivar, and SN23 is a lodging-resistant cultivar. Three light quality treatments, normal light (CK), RL and, FRL, were set up for the experiment, and the basal internodes of stems were given RL (660nm, 1.3>R/FR>1.6) for 7h and far-red light (720nm, 0.3>R/FR>1.0) for 7h daily from 9:00 to 16:00 from the beginning of jointing (GS 31; Zadoks et al., 1974) to flowering (GS 65), respectively. The light device was installed at 5cm from the ground surface with the light source perpendicular to the stalk, as shown in Figure 1. The plot size was 9m² with three replications. The treatments were planted at a density of 2.25 million plants ha⁻¹. In addition, during the growing seasons, 120kg ha⁻¹ N, 75kg ha⁻¹ P₂O₅, and 150kg ha⁻¹ K₂O were applied as the basal fertilizer before planting, and 120kg ha⁻¹ N was top-dressed at the jointing stage (GS 31). Pests, diseases, and weeds were well controlled to avoid yield losses.

The breaking strength of wheat was measured at flowering (GS 65), milking (GS 75), and dough (GS 85) stages using a digital stem strength meter (Hangzhou Tepe Instrument Co., Ltd., Hangzhou, China; Hussain et al., 2020). The internodes without stem sheath were placed 5cm away from the groove of the support post. The center of the internode was aligned with the pressure probe. The breaking strength was the peak value when the internode was broken. It was displayed on the instrument screen and its unit was Newton (N).

Fifteen representative main stems (removing stem sheath) were taken at flowering (GS 65), milking (GS 75), and dough (GS 85) stages, and the diameter, wall thickness, and dry weight per unit length of the second internode at the base of the stem were prepared by adapting the procedure used by Wei et al. (2008). The diameter and wall thickness of the stem were measured with a vernier caliper with an accuracy of 0.001mm; dry weight per unit length was scored with the formula (the dry weight of basal second internode/the length of basal second internode).

The gibberellin content in the second internode of wheat stems was measured from the beginning of the formation of the second internode to day 35. Gibberellin acid4 (GA4) and gibberellin acid7 (GA7) were extracted according to the method described by Dave et al. (2011), with the extraction buffer (isopropanol: chromatographic acetic acid, 99:1) being used as the extractant. The qualitative and quantitative analyses of the GA4 and GA7 in the wheat stem were performed by ultra-performance liquid chromatography-electrospray tandem mass spectrometry (UPLC-MS/MS) system (ACQUITY UPLC I-Class/Xevo TQ-S, Waters Co., United States).

The lignin content of the second internode of wheat stems was determined from the beginning of the formation of the second internode to day 49. Lignin determination was done as previously reported by Peng et al. (2014). Around 0.3g of fresh stem samples with leaf sheaths removed were ground in liquid nitrogen, rinsed five times with 80% ethanol to remove soluble metabolites, rinsed with acetone, and dried in a drying oven. Samples were exposed to a 4:1 (v/v) mixture of acetic acid: acetyl bromide and incubated at 70°C for 2h. Then, the samples were cooled to room temperature and transferred to 50ml volumetric flasks to which 2M NaOH, acetic acid, and 7.5M hydroxylamine hydrochloride were added to terminate the reaction. The fixation of each sample was completed with acetic acid. The absorbance of the samples was read at 280nm using a spectrophotometer (Beijing Pu-analysis General Instrument Co., Ltd.). Lignin content was expressed as OD280 ml⁻¹ g⁻¹ fresh weight (FW).

At flowering (GS 65), milking (GS 75), and dough (GS 85) stages, the content of the lignin subunits (S, G, and H) in the second internode of the wheat stem was determined. The content of lignin subunits was prepared according to the procedure used by Zheng et al. (2017b) with slight modifications. Lignin in wheat stems was oxidatively depolymerized to lignin subunits by alkaline nitrobenzene oxidation. The qualitative and quantitative analysis of the three lignin subunits in wheat straw was performed by UPLC-MS/MS system (ACQUITY UPLC I-Class/Xevo TQ-S, Waters Co., United States). A dry sample of 20mg of the protein-free stem was accurately weighed in a 50ml white digestion tube; then 3ml of 2M NaOH and 0.5ml of nitrobenzene were added, and the sample was thoroughly mixed and heated in a Master 40 microwave digestion instrument (Shanghai, Xin Yi Microwave Chemical Technology Co., Ltd.) at 170°C for 1h. The reaction mixture was then centrifuged and extracted with 4ml of ethyl acetate. The supernatant was extracted twice. Ethyl acetate was evaporated using a Concentrator plus vacuum centrifuge concentrator (Eppendorf, Germany). The sample was re-dissolved with 6ml of 50% acetonitrile in water, filtered through 0.22μm membrane, and injected into 4μl for detection on the machine.

We searched for genes in the NCBI GenBank nucleotide database¹ and the Ensembl plant’s database,² which gave homologous gene information for a given gene, then designed primers based on the comparison results of the wheat genome sequence published by IWGSC.³ The primers were designed in the conserved regions of genes. Instead of amplifying a specific gene, we amplified that gene and all its homologs at the same time as a way to test the expression of all of them, since their expression products were functionally identical. Gene-specific primers with Tm ranging from 50 to 60°C and amplifying PCR products ranging from 75 to 200bp in length were designed using Primer 5.0 (Table 1). All the designed primers were tested for full efficiency before use for expression profiling of target genes. The gene expression encompassed three independent biological replicates of each treatment. For each biological replicate, three technical replicates of each PCR reaction were performed. Briefly, total RNA from the second basal internode was isolated using a modified Trizol extraction method and treated with DNase I to remove any contaminant genomic DNA. First-strand cDNA was synthesized from 1μg of total RNA using the Prime Script RT Reagent Kit (QUANSHIJIN, China).

The quantitative PCR (qPCR) assay was performed on the 7,500 Real-Time PCR instrument (ABI, United Kingdom). The reaction consists of 2μl of the diluted cDNA as a template, 10μl of Top/Tip Green qPCR Supermix (QUANSHIJIN, China), 0.4μl of 10μM forward primer (200nM final concentration), 0.4μl of 10μM reverse primer (200nM final concentration), 0.4μl of passive reference dye (50×; QUANSHIJIN, China), and 6.8μl nuclease-free water, total 20μl. The samples were subjected to the following thermal cycling conditions: DNA polymerase activation at 94°C for 30s followed by 45cycles of denaturation at 95°C for 5s, annealing at 50–60°C (depending on the melting temperature of the primer set) for 15s and extension at 72°C for 10s in duplicate in 96-well optical reaction plates (ABI, United Kingdom). Relative expression was calculated using equation 2^-∆∆Ct (Livak and Schmittgen, 2001).

Tables and figures were processed with Microsoft Excel 2007 and SigmaPlot 14.0 Software. Data management and analysis were performed using Data Processing System software 7.05 (DPS) and IBM SPSS statistics 22 statistical software. The means and significant differences between treatments were separated using the least significant difference (LSD) test at 5% probability.

The analysis of variance showed that wheat variety and light quality significantly (p<0.01) affected wheat stem breaking strength, dry weight per length, stem diameter, wall thickness, and S, G, and H subunit contents as well as total lignin content (Table 2). What stands out in Table 2 is the effect of variety×light quality (C×L) on breaking strength, S, G, and H subunit contents were significant (p<0.01). However, the effect of variety×light quality (C×L) on dry weight per length, stem diameter, wall thickness was not significant. In addition, variety×light quality (C×L) had a significant effect on total lignin content (p<0.05).

The results, as shown in Table 3, indicated that both varieties showed a decreasing trend in folding resistance from the flowering to the waxing stage. Further analysis showed that compared to the CK treatment, the second internode stem breaking strength was significantly increased by 54.45, 46.69, and 34.73% in the three periods of the RL treatment in 2018–2019 and by 47.44, 49.03, and 36.74% in 2019–2020 (mean of both varieties), and the three stages of the FRL treatment in 2018–2019 times had a significant reduction of 21.55, 31.51, and 43.89% in the second internode breaking strength of the stem; and 26.84, 31.24, and 44.01% in 2019–2020 (mean of both varieties).

Dry weight per length in the wheat stem of both varieties showed a trend of increasing and then decreasing from anthesis to dough stage. It was apparent from Table 4 that compared to the CK treatment, the dry weight per length in wheat stem under the RL treatment increased by 13.10, 7.23, and 18.73% in 2018–2019; 26.15, 26.44, and 26.09% in 2019–2020 (mean of both varieties), respectively, and the dry weight per length in wheat stem under the FRL treatment in 2018–2019 decreased by 8.81, 3.44, and 9.92%, respectively; and 8.92, 17.07, and 10.51% (mean of both varieties) from 2019 to 2020, respectively. In addition, dry weight per length in the wheat stem was significantly higher in SN23 than in SN16 (Table 4).

We could see that R/FR had significant effects on the second internode diameter and wall thickness at the base of the stem, and both increased with the increase in the proportion of R/FR (Figure 2). As can be seen in Figure 2, compared to the CK treatment, the second internode diameter of the RL treatment increased by 4.06% and the FRL treatment decreased by 3.07% in 2018–2019, and in 2019–2020, the second internode diameter of the RL treatment increased by 4.19% and the FRL treatment decreased by 1.70% (mean of both varieties). Compared to the CK treatment, the wall thickness of the second internode of the plants increased by 6.44% in the RL treatment and decreased by 11.59% in the FRL treatment in 2018–2019, and increased by 12.45% in the RL treatment and decreased by 6.74% in the FRL treatment in 2019–2020 (mean of both varieties). Additionally, diameter and wall thickness were significantly higher in SN23 than in SN16.

It could be seen from the data in Figure 3 that the light quality within the population significantly affected the gibberellin content (GA4, GA7) of wheat stems, with GA7 content significantly higher compared to GA4 (Figure 3). Compared with the CK treatment, the RL treatment significantly reduced the GA4 and GA7 content and the FRL treatment significantly increased the GA4 and GA7 content. Further analysis showed that GA4 and GA7 contents showed a trend of increasing and then decreasing with time and reached the maximum value 7days after the formation of the second node at the base of the wheat stem. In addition, the gibberellin content (GA4, GA7) of SN23 was always lower than that of SN16.

Lignin accumulation during wheat stem development was determined according to the effect of different R/FR conditions on the degree of lignification of the second internode of wheat stems (Figure 4). In 2018–2019, the total lignin accumulation under the RL treatment increased by 13.34% compared with the CK treatment; the FRL treatment decreased by 10.91% compared with the CK treatment (CK; mean of both varieties); in 2019–2020, the total lignin accumulation under the RL treatment increased by 13.95%; the FRL treatment decreased 13.44% compared to the CK treatment (CK; mean of both varieties). Figure 4 showed that there were differences in lignin accumulation during stem development in winter wheat, with faster accumulation rates from 0 to 35 days after the formation of the second internode. These results suggested that the jointing and anthesis stages were critical periods for lignin accumulation. In addition, the trend of lignin accumulation was the same in both varieties, but the lignin accumulation of SN23 was higher than that of SN16, under different treatments.

As shown in Figure 5, wheat stem lignin consisted of S, G, and H subunits, with S and G subunits being the most abundant and H subunits accounting for the smallest percentage. The total lignin content and subunit content increased with fertility, and the total lignin content and subunit content of SN23 were higher than those of SN16, and both responded to RL and FRL with similar trends. In 2018–2019, compared with the CK treatment, the total lignin subunit content of the RL treatment significantly increased by 35.28, 36.54, and 31.95% (mean of both varieties); the FRL treatment significantly decreased by 27.38, 30.88 and 18.84% (mean of both varieties).

As could be seen from Figure 6, in 2018–2019, compared with the CK treatment, the RL treatment significantly increased the content of S and G subunits and the proportion of S+G subunits, and significantly decreased the proportion of H subunits; the FRL treatment significantly decreased the content of S and G subunits and the proportion of S+G subunits, and significantly increased the proportion of H subunits. Overall, these results indicated that S+G and H subunit proportions play a key role in the performance of wheat in terms of resistance to lodging.

The RL treatment significantly increased the relative expression of TaPhy B, which decreased gradually with the growth and development process, and the RL treatment slowed down this decrease, while the FRL treatment accelerated it (Figure 7). The RL treatment significantly decreased the relative expression of TaPhy A, which gradually decreased with the growth and development process, and the RL treatment accelerated this rate of decrease, while the FRL treatment slowed down this rate of decrease.

As shown in Figure 7, the transcript abundance of most of the key lignin synthesis enzymes decreased gradually with the progress of growth and development. The effect of different light qualities on the relative expression of TaPAL and TaC3H gradually increased with increasing days of light exposure, and high RL and FRL proportion conditions retarded the decreasing rate. The relative expression of TaC4H and Ta4CL genes reached the maximum at 14days after illumination, and the effect of different light quality on the relative expression of TaC4H and Ta4CL genes followed the trend: RL>CK>FRL. The relative expression of TaHCT and TaCSE genes increased significantly under high RL and FRL proportion conditions, but the differences with other treatments gradually decreased over time. TaF5H and TaCOMT were key enzymes in the lignin S subunit production step, and RL treatment significantly increased the relative expression of both genes, while the differences were not significant under the CK treatment and the FRL treatment, and with the progression of fertility. TaCCR and TaCAD are the last two steps of lignin subunit synthesis, and the relative expression of TaCCR gene reached the maximum after 14days of light treatment, while TaCAD gradually decreased with the fertility process. Compared with the CK treatment, the RL treatment significantly increased the relative expression of both genes, the FRL treatment had no significant effect on the relative expression of TaCAD gene, and the FRL treatment significantly decreased the relative expression of TaCCR gene. The relative expression of all genes was higher in SN23 than in SN16 during all stages of stem development, and the lowest relative expression of all genes was achieved at maturity.

Stem morphological characteristics such as filling, diameter, and wall thickness of the second internode at the base of the stem were highly significantly and positively correlated with folding resistance (Figure 8). The correlation analysis also showed that the lignin content of the second internode at the base of the stem was highly significantly and positively correlated with folding resistance, and the proportion and distribution of lignin subunits affected the breaking strength of the stem, with the proportion of H subunits significantly and negatively correlated with the breaking strength of the stem, and the proportion of S+G subunits significantly and positively correlated with the breaking strength of the stem (Table 5). The results indicated that increasing the second internode filling, stem thickness, and wall thickness at the base of the wheat stem was beneficial to enhance its fracture resistance and increasing the lignin content and S+G subunit proportion at the base of the stem was beneficial to enhance the resistance of wheat to lodging.

Plants sense changes in the light environment through photopigments, the intracellular distribution of which is influenced by light intensity and light quality. It is involved in various photomorphogenesis by participating in the R/FR reaction and influencing the growth process of plants. It regulates various photomorphogenic responses to optimize plant growth and morphogenesis by participating in the R/FR response (Smith, 2000; Rockwell et al., 2006; Fraser et al., 2021).

The breaking strength of the internode at the base of the stem is closely related to the light environment, and this study showed that high R/FR conditions (Figure 9) significantly increased the diameter and wall thickness of the second internode at the base of the wheat stem (Figure 2) and increased stem fullness (Table 4), thereby significantly enhancing stem breaking resistance and improving the lodging resistance. This is due to the existence of two types of reactions in the photomorphogenesis process, namely the R/FR reaction and the blue light/near-ultraviolet light reaction. Among them, the proportion of R/FR has the greatest effect on plant photomorphogenesis (Smith, 1982; Park and Runkle, 2017). Among the many phytochromes, Phy A is the only FRL receptor in plants, and at low R/FR proportions, Phy A regulates the inhibition of shade avoidance responses leading to stem elongation, whereas, at high R/FR proportions, the inhibition for shade avoidance responses is mainly regulated by Phy B (Shlumukov et al., 2001; Zhou and Li, 2019). These findings confirmed that the reasons for the different expressions of Phy A and Phy B genes under different R/FR conditions (Figure 7). After the light environment changed, the distribution of phytochromes is altered to enter the nucleus to bind to some transcriptional regulators (Franklin and Quail, 2010; Wang and Zhang, 2017). In Arabidopsis, under low R/FR conditions, Phy B converts to Pr type, it transfers from the nucleus to the cytoplasm, thus enhancing the stability of PIF 4 growth regulators, at the same time, the cycling of Phy A between Pr and Pfr types compensates to some extent for the release of PIF 4 and prevents its degradation, which promotes the synthesis of GA (Tamotsu et al., 2005; Lorrain et al., 2007). Similar findings have been made recently in wheat (Pham et al., 2018; Sibbett, 2018). At the same time, in the condition of high R/FR, phytochrome inhibits the biosynthesis of GAs, and Phy A cooperates with Phy B to regulate the inactivation of GAs (Franklin, 2008). Stem elongation is one of the functions of GAs. The decrease of GAs concentration enhances the stability of DELLA protein, promotes the degradation of PIF, and finally inhibits stem elongation (Fumiaki et al., 2012; Balcerowicz, 2020). Analysis of the results with our studies confirmed that high R/FR conditions significantly reduced GA content (Figure 3).

Several reports show that increasing planting density results in lower R/FR at the base of the wheat stem, elongate growth at the base internodes, significantly reducing lodging resistance of wheat (Syros et al., 2005; Han et al., 2015), which is consistent with the results of Table 3. However, these findings were contrary to previous studies, which have suggested that the height of cucumber (Cucumis sativus Linn.), tomato (Lycopersicon esculentum Miller.), and ginger (Zingiber officinale Roscoe) plants were significantly increased under supplemental RL conditions, increasing the risk of lodging (Zhang et al., 2008; Cui et al., 2009). The difference in the results of this experimental study may be related to the differences between species.

Gibberellin regulates lignin biosynthesis by affecting lignin monomer polymerization (Biemelt et al., 2004). Light regulates gibberellin biosynthesis and metabolism (Liu et al., 2021). RL inhibits the elongation of Arabidopsis stems by controlling the synthesis of GAs (Tamotsu et al., 2005; Liu et al., 2021). The effect of GAs on xylem lignin is less, but significantly reduces the lignin content of phloem. One study shows that stem lignin content is significantly reduced by spraying exogenous GA3 onto wheat (Peng et al., 2014). GA1, GA3, GA4, and GA7 exist commonly in plants as the primary bioactive forms of GAs (Hedden and Thomas, 2012; Tavares et al., 2018). GAs have different activities in different tissues/organs, e.g., embryo in germinating seeds of rice, young leaves in developing seedlings, the base of internode elongation (Silverstone and Sun, 2000). GA1 is the active GAs in soybean seedlings, GA4 is the active GAs in Arabidopsis hypocotyl, and GA3 plays a major role in seed germination (Kaneko et al., 2003; Falcioni et al., 2018). The GA7 content in wheat stems was remarkably high compared to GA4 in Figure 3, so we guessed that GA7 might be the active GAs in wheat stems. Some studies find that the GAs content and lignin content of dwarf plants are lower than those of wild-type plants, and the GAs content is positively correlated with the lignin content of stems (Biemelt et al., 2004; Okuno et al., 2014). Overexpression gibberellin 2-oxidase (GA2oxs) in switchgrass develop plants with low level of lignin (Wuddineh et al., 2015), whereas the lignin content of stems decreased when exogenous GA3 sprays on wheat (Peng et al., 2014). That could be that low levels of GAs increased lignin content and high levels of GAs inhibited lignin synthesis. In addition, our study also found that wheat varieties with different lodging resistance were consistent in their response to R/FR, but the GAs and lignin contents varied greatly. Interestingly, the lignin and GAs contents of the lodging-sensitive wheat under high R/FR conditions were maintained at the same level as those of the lodging-resistant wheat under CK. It indicates that lignin and GA contents are closely related to genotypes and regulated by environmental factors.

The analysis of the breaking strength of stems at different periods in wheat revealed that the breaking strength of stems decreased with the growth and development process (Table 3), while total lignin increased with the growth and development process (Figure 4). Surprisingly, that was contrary to a significant positive correlation between breaking strength and lignin content (Figure 8D). This result may be explained by the fact that stem breaking strength is not only correlated with lignin content, but also with other traits of the stem. Wheat straw is composed of 40% cellulose, 35% hemicellulose, and 25% lignin (Ruiz et al., 2013). In wheat, cellulose is an important factor in regulating the mechanical strength of the stem, which has implications in crop lodging (Kaur et al., 2017). Hemicellulosic arabinose reduces the ability of rice to lodging resistance by reducing cellulose crystallinity (Li et al., 2015). The low lignin content but higher breaking strength of stems in wheat during early growth and development may be due to the high pectin and water content inside the stem. Findings have shown that downregulation of switchgrass pectin biosynthesis gene GAUT4 leads to reduced ferulic acid synthesis, lower cross-linking of lignin hydrate, and lower rigidity and fracture resistance of the stem, and these nonstructural substances are also important factors in determining the magnitude of the stem breaking strength (Handakumbura and Hazen, 2012; Li et al., 2019). Another important finding is that differences in the structure and proportion of lignin subunits determine their contribution to wheat stem resistance performance (Zheng et al., 2017a). The analysis of lignin subunits revealed that the proportion of H subunits also increased with fertility (Figure 6), and the proportion of H subunits was significantly and negatively correlated with the breaking strength of the stem (Table 5), which could also explain this phenomenon.

Phy A and Phy B promote the synthesis of lignin (Hernando et al., 2021). This study showed that high R/FR conditions significantly increased TaPhy B gene expression and significantly decreased TaPhy A gene expression (Figure 7). The activated state of photosensitized pigments activates the expression of related genes through G-protein, calmodulin, and cGMP pathways, which induce oxidative lignin deposition by regulating the flow of photosynthetic carbon sources to the lignin synthesis pathway (Anderson et al., 2015). The R/FR proportion at the base of the wheat stem may be a key factor in altering lignin synthesis and ultimately lignin composition, and a high R/FR proportion significantly increases lignin accumulation in wheat (Liu et al., 1996; Luo et al., 2019). Analysis of stem lignin subunit data revealed that high R/FR proportions significantly increased lignin subunit content and significantly reduced H subunit content and proportions (Figures 5, 6), with most of the reduced H subunit present as S and G subunits (Figure 5), and that R/FR significantly affected stem breaking strength (Table 3).

It has been suggested that the S/G proportion at the basal second internode of wheat stems was significantly and positively correlated with the stem breaking strength (Zheng et al., 2017a). However, some studies have shown that there were multiple possibilities for the relationship between the content and proportion of lignin subunits and the breaking strength of the stem (Luo et al., 2019). This study found that the content of G+S+H lignin subunits, the proportion of G+S lignin subunits to H subunits were the main factors affecting the breaking strength of wheat stems, and the stem breaking strength was significantly correlated with the S/G proportion only at the dough stage (Table 5).

Light affects lignin biosynthetic capacity at the transcriptional level, and the R/FR proportion ultimately affects lignin synthesis, subunit proportion, and distribution (Rogers et al., 2005; Mah et al., 2018), thus affecting stem lodging resistance performance. The current study found that a high R/FR proportion significantly increased the expression of lignin synthesis-related enzyme genes (Figure 7). The increase in lignin content was achieved by increasing the activity and gene expression of lignin biosynthetic enzymes in the lignin synthesis process. Genes involved in the early steps of the lignin biosynthesis pathway affect lignin content, while genes involved in the final steps alter lignin composition (Li et al., 2008; Yu et al., 2020). It has been shown that PAL, CAD, HCT, COMT, and 4CL are significantly and positively correlated with lignin content in wheat (Boudet et al., 2003; Bi et al., 2010). Under shade conditions, the R/FR proportion decreased, thereby decreasing the expression of CAD, 4CL, and PAL genes. The expression of lignin synthesis-related enzyme genes is significantly reduced in soybean and Arabidopsis, under low R/FR proportion (Rogers and Campbell, 2004; Hussain et al., 2020). Under lower R/FR conditions, C3H, C4H, CCR, CCoAOMT, and POD expression are downregulated in soybean, leading to a decrease in the content of lignin synthesis precursors such as caffeic acid, erucic acid, and ferulic acid, and alter lignin subunit content and proportion (Ren, 2017; Liu et al., 2018b).

In addition to the above genes, this study showed that the expression of CSE and F5H genes in the lignin synthesis pathway increased with the increase of R/FR proportion. CSE is a key enzyme for carbon flow partitioning in the lignin synthesis pathway, and F5H is a key enzyme for the synthesis of lignin S subunits and their free radicals. It has an important effect on lignin biosynthesis and molecular structure polymerization (Day et al., 2009; Stewart et al., 2009). Another important finding was that the expression of each gene reached its lowest level at the late stage of wheat growth and development, and the difference between treatments decreased due to the programmed cell death stage of wheat, where the stem lost its biological activity and lignin synthesis stopped, but high R/FR conditions could delay this decreasing trend (Figure 7). The most interesting finding was that most of the gene expression peaks occurred in the early stage of stem development and significantly decreased in the late stage of development, but lignin content continued to increase. A possible explanation for this might be that higher levels of gene expression occur in the early stages of stem development when the rate of lignin accumulation is high, and that lignin accumulates later in growth and development due to the stability of lignin that is difficult to break down, thus allowing lignification to be maintained.

We can have a bold guess that the light environment at the base of dense populations becomes poor (R/FR↓). Low R/FR condition makes phytochrome inactive (Franklin, 2008). Inactivated phytochrome leads to an increase in GAs content and a decrease in lignin content, resulting in slender stems and deteriorating stem quality (Tamotsu et al., 2005). After RL irradiation (R/FR↑), phytochrome transforms into the activation state (Pfr). The activation state (Pfr), which inhibits the shade avoidance response, leads to a decrease in GA content and an increase in the expression of keys enzyme genes for lignin synthesis, contributing to improved stem quality (Fumiaki et al., 2012; Balcerowicz, 2020). The hypothetical mechanism by which high R/FR promotes lignin biosynthesis is shown in Figure 10. Future genomic, transcriptomic, and proteomic analyses should reveal more details and deepen our understanding of the role of R/FR in crop growth and development.