The contribution of photosynthesis traits and plant height components to plant height in wheat at the individual quantitative trait locus level

Plant height is an important agronomic trait for morphogenesis and grain yield formation in wheat. In this study, we performed both normal and multivariate conditional quantitative trait locus (QTL) analyses for plant height with spike length, internode number, length of the first internode to the sixth internode from the top during harvest, and photosynthesis traits at the seedling stage and heading stage based on a recombinant inbred line population. A total of 49 normal QTLs were detected, as well as 312 conditional QTLs. The genetic region Xbcd1970-Xbcd262 on chromosome 2D harbored the most QTLs, with 6 normal QTLs and 39 conditional QTLs. A comparison between the normal and conditional QTL mapping analyses suggested that the length of the third internode, fourth internode, and fifth internode from the top showed a high genetic association with plant height, whereas all photosynthesis traits showed weaker associations. This comparative analysis could serve as a platform for dissecting the genetic relation between objective traits and other phenotypic traits before manipulation of genes collocated with QTL clusters.

Scientific RepoRtS | (2020) 10:12261 | https://doi.org/10.1038/s41598-020-69138-0 www.nature.com/scientificreports/ recombinant inbred line (RIL) population to design a QTL mapping analysis of plant height and its component traits, including spike length, length of the first internode to the sixth internode from the top, and internode number. Additionally, we also performed a conditional analysis to quantify the genetic relationship between plant height and plant height components, photosynthesis traits, and transpiration rate.

Material and methods
plant materials and experimental design. The experimental design used in the present study has been described previously 17 . Briefly, the International Triticeae Mapping Initiative (ITMI) population comprising 112 RILs was used 18 . This population was derived from the synthetic hexaploid wheat cultivar W7984 and the winter wheat cultivar Opata 85. The two parents and all RILs were grown in Mianyang, Sichuan Province (31°32′ N and 104°42′ W), from 2015 to 2017. Each line was planted in five 1.5 m rows. The distance between two rows was 30 cm, and the seeding rate was 10 cm per plant. Three repetitions were included in accordance with a randomized complete block design in each year.
phenotypic evaluations. For each line, all phenotypic traits were evaluated using five representative leading tillers. Plant height and its component traits, including spike length, internode number, and the length of the first internode to the sixth internode from the top, were measured during harvest. Five photosynthesis traits, chlorophyll content, the net photosynthetic rate, intercellular CO 2 concentration, the transpiration rate, and stomatal conductance, were investigated at both the seedling and heading stages. The chlorophyll content was measured using a SPAD-502 Plus chlorophyll meter (Minolta, Osaka, Japan), and the other four photosynthesis traits were estimated using an LI-6400XT Portable Photosynthesis system (LI-COR, Lincoln, USA). Data from two consecutive years (E1 of 2015-2016 and E2 of 2016-2017) were combined and then used as phenotypic traits in E3.

Statistical analysis.
Analyses of variance and analysis of covariance (ANCOVA) were based on the minimum norm quadratic unbiased estimation proposed by Zhu 19 . The phenotypic and genetic correlations among all the traits were investigated via the multivariate restricted maximum likelihood (REML) by PROC MIXED of SAS software version 9.4 20,21 . A network analysis was carried out with the R package 'qgraph' 22,23 , and heritability was estimated following the approach suggested by Ma et al. 24 . Conditional analysis was performed using the QGAStation software package 1.0 25 , employing the 'QTLData' menu settings as described previously 15 . The output file of QGAStation provides the conditional phenotypic values y (A|B) , which indicate a new value of A without the influence of B. In the conditional analysis, plant height was considered the final complex trait, whereas the plant height components and photosynthesis traits were considered causal traits.
Genetic mapping and QTL analysis. The genetic linkage map of the ITMI population used in our study was generated by Song and coauthors 26 , was 2541 cM in size, contained 1,410 distinct loci and covered 2,541 cM of genetic distance. Inclusive composite interval mapping was performed using IciMapping 4.1 27 . The walking speed and LOD threshold parameters were as previously described 17 . Both the normal phenotypic values and the new conditional phenotypic values were used in the QTL mapping analyses. These analyses were initially performed separately for each year and then carried out using data combined across both years. Conditional analysis at the QTL level. The comparative analysis of normal and conditional QTL mapping had four possible outcomes: (1) a normal QTL for trait A with no identification for the conditional trait of y(A|B), suggesting that this QTL for trait A is derived fully from trait B; (2) a normal QTL for trait A can also be detected for y(A|B) with a phenotypic variation range greater than 10%, indicating that this QTL for trait A is partly affected by trait B; (3) a normal QTL for trait A was also detected for y(A|B) with similar phenotypic variation, indicating that this QTL for trait A is totally independent of trait B; (4) only one additional QTL for y(A|B) was identified, suggesting that the expression of the QTL for trait A is suppressed by trait B, and this expression could be identified by elimination of trait B influence.

Results phenotypic evaluations.
The data for all phenotypic traits assessed here are presented in Table 1. On average, the synthetic hexaploid parent of W7984 has more than six internodes, whereas Opata 85 has only five internodes.
We first assessed plant height and its related component traits; the W7984 parent was 22.07 cm taller on average than was the Opata 85 parent (P < 0.05). The internode length was 8.65% longer for the top node of Opata 85 but was 3.73%, 15.90%, 26.95%, and 40.42% longer for that of W7984 from the 2nd node to the 5th node, respectively (P < 0.05). W7984 had a higher chlorophyll content than Opata 85 at the seedling stage (P < 0.05), but the difference decreased at the heading stage. Similar trends were observed for photosynthetic capacity, and transgressive segregation was observed for all traits in the population.
The highest heritability was observed for the net photosynthetic rate at the heading stage, followed by the chlorophyll content, while the intercellular CO 2 concentration at the heading stage exhibited the lowest heritability value of 36.04% during the two years. Genetic correlation analysis showed that plant height was positively correlated with the spike length (34%), first internode length from the top (48%), second internode length from the top (62%), third internode length from the top (69%), fourth internode length from the top (77%), fifth internode length from the top (69%), sixth internode length from the top (47%), net photosynthetic rate at the heading stage at the heading stage (21%), and intercellular CO 2 concentration at the heading stage (20%) (  (1st); length of second internode from the top (2nd); length of third internode from the top (3rd); length of fourth internode from the top (4th); length of fifth internode from the top (5th); length of sixth internode from the top (6th); chlorophyll content at seedling stage (CNS) and heading stage (CNH); net photosynthetic rate at seedling stage (PNS) and heading stage (PNH); stomatal conductance at seedling stage (GsS) and heading stage (GsH); intercellular CO2 concentration at seedling stage (CiS) and heading stage (CiH); transpiration rate at seedling stage (TRS) and heading stage (TRH); respectively. www.nature.com/scientificreports/ also contained stable QTLs for the length of the second internode from the top, the length of the fifth internode from the top, and the total internode number.

QTL analysis for plant height conditioned upon plant height components and photosynthesis traits.
When the influences of the third internode length from the top, the fourth internode length from the top, and the fifth internode length from the top on plant height were removed, no QTLs for plant height were detected in any environment (Table 3), suggesting a high genetic association between these three component traits and plant height. Similar and weaker associations were observed for the influences of the first internode  www.nature.com/scientificreports/ length and the second internode length from the top, respectively, on plant height. At least one normal QTL for plant height was detected using conditional analysis for all photosynthesis traits, with the exception of the chlorophyll content at the seedling stage and net photosynthetic rate at both seedling and heading stages for the combined yearly data. Conditional analysis also allowed the identification of additive QTLs not detected using traditional QTL analysis. Removing the influence of the second internode length from the top allowed the identification of the greatest number of extra QTLs contributing to plant height, with three, four, and five identified in environments E1, E2, and E3, respectively. The same approach allowed for the identification of 8 QTLs contributing to the fifth internode length from the top, the sixth internode length from the top, stomatal conductance at the heading stage, and the intercellular CO 2 concentration at the seedling stage.
The combination of normal QTLs that disappeared and extraconditional QTLs indicated a higher genetic relationship between plant height and its components than between plant height and photosynthesis traits, with the exception of spike length and internode number. Similar tendencies were observed for the phenotypic and genetic correlation analyses.

Discussion
QTLs for plant height with respect to previous studies. Previous genetic studies have identified numerous major QTLs for plant height on the majority of the 21 wheat chromosomes 15 . The semidominant GAresponsive dwarfing gene Rht18 is linked to Xbarc3 on the short arm of chromosome 6A 28 and is a candidate gene for inclusion in breeding programs to improve wheat performance 29 . This gene was identified in the present study between Xcdo29 and Xtam36. The QTL on 1A is located in approximately the same region as QHt.fra-1A is 30 . A major QTL on 5A in the genetic region of Xmwg522-Xbarc330 was mapped to a position where QHt. inra-5A was found, but it was identified only in the environment of E1. Additionally, a major QTL in the genetic region of Xgwm293-Xcdo785 on 5A was identified in E2 and E3, confirming the findings of previous reports 31 . Genetic relationship between plant height and its components and photosynthesis traits at specific QTLs. Comparative analyses of normal and conditional QTL mapping can dissect genetic relationships between plant height and related traits. A normal QTL for plant height was identified on 1A in the interval between Xgwm136 and XksuD14 (Supplementary information Table S1, Fig. 2). However, this QTL was no longer detected once the influences of the first internode length to the sixth internode length were removed. Therefore, we concluded that at this locus, variation in plant height was attributable entirely to internode length (the first to the sixth internodes). While no normal QTLs for plant height components were detected at this locus, conditional QTLs for the fifth internode length were identified when the influence of intercellular CO 2 concentration at the seedling and heading stages, the net photosynthetic rate at the seedling stage, and the transpiration rate at both stages were removed. This indicated that the fifth internode length was suppressed by these photosynthesis traits. However, the relationship between these photosynthesis traits and plant height was nonsignificant during the direct conditional analysis.
For the locus found at Xmwg522-Xbarc330 on 5A, plant height was entirely attributable to the variation in spike length, internode number, and the third internode length to the sixth internode length. We also identified a normal QTL for the fourth internode length and the fifth internode length, but their genetic correlations with photosynthesis traits were insignificant. A conditional QTL was also identified at this location, but its expression was suppressed by the chlorophyll content at the heading stage. Previous work has described two singlenucleotide polymorphism (SNP) markers (IWA6949 and IWA6573) that are associated with the MIKC-type MADS-box gene(s) of WM10A (TraesCS5A02G293000 in the Ensembl plant database) in the same genomic region 32,33 . MIKC-type MADS-box genes are involved in virtually all aspects of plant development in a variety of important crop species, including wheat, rice, and banana [34][35][36] . Manipulation of this gene might be an avenue to improve important agronomic traits, such as plant height.
For the locus in the region of Xcdo29-Xtam36 on 6A, plant height was entirely controlled by the chlorophyll content at the seedling stage, the transpiration rate at the seedling stage, and all plant height components, with the exception of the first internode length. A normal QTL for the fourth internode length was identified at this locus. Its expression entirely contributed to chlorophyll content at the seedling stage and heading stage and both stomatal conductance and net photosynthetic rate at the heading stage. The contribution of normal QTLs to the transpiration rate at the heading stage to the fourth internode length was insignificant. Previous studies have shown that the Xcdo29-Xtam36 region contains a SNP marker (IWA1520) associated with the Figure 2. QTL clusters carried normal QTL or conditional QTL for plant height. Traits are plant height (PH); spike length (SL); internode number (IN); length of first internode from the top (1st); length of second internode from the top (2nd); length of third internode from the top (3rd); length of fourth internode from the top (4th); length of fifth internode from the top (5th); length of sixth internode from the top (6th); chlorophyll content (CN); net photosynthetic rate (PN); stomatal conductance (GS); intercellular CO2 concentration (Ci); transpiration rate (TR); respectively. S and H indicating seedling stage and heading stage respectively. The networks are comparison of conditional QTL and normal QTL analysis. Yellow arrow indicated suppression between plant height components or photosynthetic traits with plant height, or photosynthetic traits with plant height components, whereas the red arrow indicated positive contribution between two traits. Dotted arrow indicated independent relationship between two traits. Traits with green background or green internode indicated normal QTLs were identified, whereas traits with yellow background or yellow internode indicated only conditional QTLs were identified. Utilization of conditional analysis for plant breeding and pathway network construction. In plant breeding, manipulation of a single gene typically results in changes to several phenotypic traits; for example, overexpression of H + -ATPase results in enhanced light-induced stomatal opening, photosynthesis, and plant growth in Arabidopsis 38 . Engineering plants to carry favorable traits such as those that promote disease resistance, hardiness, and yield quality can also result in the manifestation of unwanted characteristics such as excessive height, late heading and maturity, hardened glumes, and increased threshing difficulty 39 . Thus, the improvement of a single trait requires examination of the genetic relationships between the target trait and other phenotypic traits prior to gene manipulation. Taking the first internode length as an example, this trait previously showed a positive correlation with lodging resistance and kernel weight 40,41 . We identified four QTLs for this trait in the present study, all of which showed insignificant effects on plant height, especially those traits associated with the environment-independent locus within Xbarc117-Xgwm129 on 5A. Introgression of segments containing the W7984 allele in common wheat might lengthen the uppermost internode but has no effect on plant height. In classic genetic analysis, verification of gene-trait relationships need to be confirmed by fine mapping, map-based cloning, and transgenic analysis 42 , each step of which is laborious and time consuming. Improved annotation rates/quality of wheat genome reference assemblies will facilitate easier identification of orthologous genes of genes of known function among syntenic regions, and QTLs identified in model genomes can then be sought in crop genomes for potential inclusion in plant breeding strategies 43,44 . Moreover, the expression levels of candidate genes and/or enzyme activity and metabolite concentrations can also be used as phenotypic traits. Conditional analysis therefore provides a platform for the construction of pathway networks to bridge the gap between genes and traits of interest.

conclusions
In this study, we detected 4 normal QTLs and 312 conditional QTLs for plant height. In 40 QTL clusters, the genetic region of Xbcd1970-Xbcd262 on chromosome 2D carried the most QTLs, with 6 normal QTLs and 39 conditional QTLs. The combination of normal and conditional QTL mapping analyses suggested that the length