Diversity in Phenological and Agronomic Traits of Miscanthus sinensis Collected in Korea and Eastern Asia

: Four-year ﬁeld experiments were conducted to investigate phenotypic traits associated with the biomass yield of 173 Miscanthus sinensis accessions collected from Korea and neighboring East Asian countries. Nine phenological and agronomic traits associated with biomass yield were assessed to investigate their phenotypic diversity and relationships with biomass yield as well as the latitudes of the M. sinensis accessions collection sites. Correlation analyses among phenological and agronomic traits, biomass yield, and collection site revealed that heading date, vegetative growth duration, leaf area, and stem growth traits (stem height, stem diameter, and stem dry weight) were closely related to biomass yield. The latitude of collection site exhibited a signiﬁcant negative correlation with heading date, and heading date showed a signiﬁcant positive correlation with biomass yield, indicating the high biomass potential of the accessions originating from lower latitude due to longer vegetative growth. The best biomass yield was mainly observed in M. sinensis accessions from the southern parts of Korea, such as Jeolla and Jeju provinces, with over 20 Mg DM ha − 1 . Agronomic traits measured in the second year after planting also showed a high correlation with biomass yield measured in the fourth year after planting. In particular, vegetative growth duration, leaf area, stem diameter, and stem dry weight measured in the second year were signiﬁcantly related to the fourth-year biomass yield. Therefore, these ﬁndings suggest that agronomic traits measured in the second year can be used for screening M. sinensis genetic resources and breeding lines with high biomass yield potential. study investigated the relationships between the second-year phenotypic traits and the fourth-year biomass yield. Our results suggest that potential biomass yield can be predicted in an early year of miscanthus cultivation—for example, the second year after planting. The early traits showing high correlations with the fourth-year biomass yield may be utilized as selection traits for screening breeding lines with high biomass yield potential in a miscanthus breeding program, saving time and effort in phenotyping and screening.


Introduction
Miscanthus species are exceptional among C 4 grasses for their biomass productivity in temperate climates [1] and have been considered as a dedicated bioenergy crop due to their ideal characteristics, such as high biomass yield potential, broad environmental adaptation, low inputs after plant establishment, and high carbon sequestration [2]. For commercial biomass production from miscanthus, Miscanthus × giganteus has mainly been cultivated due to its high biomass yield and environmental adaptation to cold temperate regions. Many long-term field trials conducted in Europe revealed that the maximum biomass yield from parental populations usually takes 7-10 years [25][26][27][28]. Miscanthus is a perennial outcrossing species, and even when propagated by rhizome, it requires 4 years to reach its maximum potential biomass yield after planting [2,7,18]. This indicates that a miscanthus breeding and yield evaluation may take four times longer than that of annual crops. Such a long period for phenotyping and screening is a major bottleneck in a miscanthus breeding program [21]. Thus, reducing the evaluation period to estimate miscanthus' performance in biomass production is key in reducing the breeding period. It can be hypothesized that earlier phenotypic traits assessed in the first or second year after a miscanthus planting may be related to later phenotypic traits assessed in the third or fourth year and potential biomass yield. The earlier phenotypic traits closely associated with potential biomass yield can be effectively used for selection traits for miscanthus breeding, and thus, time and effort can be significantly saved.
Therefore, this study was conducted to determine the phenotypic diversity of M. sinensis by evaluating its phenotypic traits during the first four years after transplanting, to investigate the relationship between phenotypic traits and biomass yield, and to identify key traits at an early year after planting that can be utilized to estimate biomass yield potential of mature stands.

Collection of M. sinensis Accessions
Rhizomes of M. sinensis were collected from various locations in Korea and the neighboring countries China, Japan, and Russia in autumn or early spring from 2008 to 2011. In total, 144 accessions were collected in Korea, 17 accessions in Northeastern China and Far Eastern Russia, and 12 accessions in Southern Japan ( Figure 1). Miscanthus × giganteus was used as a reference species. GPS information of collection sites was recorded using GPSMAP 60CS (Garmin Ltd., Lenexa, KS, USA).
Plant breeding is one of the main methods to improve productivity, quality, and other traits [11]. In the case of outcrossing an annual plant, breeding a new crop cultivar from parental populations usually takes 7-10 years [25][26][27][28]. Miscanthus is a perennial outcrossing species, and even when propagated by rhizome, it requires 4 years to reach its maximum potential biomass yield after planting [2,7,18]. This indicates that a miscanthus breeding and yield evaluation may take four times longer than that of annual crops. Such a long period for phenotyping and screening is a major bottleneck in a miscanthus breeding program [21]. Thus, reducing the evaluation period to estimate miscanthus' performance in biomass production is key in reducing the breeding period. It can be hypothesized that earlier phenotypic traits assessed in the first or second year after a miscanthus planting may be related to later phenotypic traits assessed in the third or fourth year and potential biomass yield. The earlier phenotypic traits closely associated with potential biomass yield can be effectively used for selection traits for miscanthus breeding, and thus, time and effort can be significantly saved.
Therefore, this study was conducted to determine the phenotypic diversity of M. sinensis by evaluating its phenotypic traits during the first four years after transplanting, to investigate the relationship between phenotypic traits and biomass yield, and to identify key traits at an early year after planting that can be utilized to estimate biomass yield potential of mature stands.

Collection of M. sinensis Accessions
Rhizomes of M. sinensis were collected from various locations in Korea and the neighboring countries China, Japan, and Russia in autumn or early spring from 2008 to 2011. In total, 144 accessions were collected in Korea, 17 accessions in Northeastern China and Far Eastern Russia, and 12 accessions in Southern Japan ( Figure 1). Miscanthus × giganteus was used as a reference species. GPS information of collection sites was recorded using GPS-MAP 60CS (Garmin Ltd., Lenexa, KS, USA).    (Figure 2). As both sites are located in a similar latitude, no significant difference in temperature and precipitation was observed during the cultivation period. All collected miscanthus rhizomes were initially planted in plastic pots and placed in a glasshouse until they were transplanted in a field evaluation nursery. In April 2011, all the accessions were planted in the field located in Suwon and Yeoju at a density of 1 plant m −2 accession −1 by completely randomized planting with three replications (2 replications in Suwon and 1 replication in Yeoju). The soils field-tested from Suwon and Yeoju were a sandy loam with a pH of 5.60 and 4.90, electrical conductivity of 0.23 and 0.36 dS/m, total organic carbon of 1.26 and 2.60%, total nitrogen of 0.18 and 0.10%, available phosphorus (P) concentration of 100.33 and 134.84 mg kg −1 , and soil cation exchange capacity (CEC) of 15.83 and 5.76 cmol kg −1 , respectively. Nitrogen fertilizer (urea granule) was applied at 60 kg N ha −1 year −1 in early June, while other fertilizers such as phosphorus and potassium were not applied. Weed management was conducted manually and by sequential treatment of the following herbicides: S-metolachlor (Dual Gold ® ; Syngenta Korea, Korea) as a pre-emergence herbicide application at 750 g a.i. ha −1 and dicamba (Banvel ® ; Sungbo Chemical Co., Ltd., Seoul, Korea) as a post-emergence herbicide application at 482 g a.i. ha −1 [29].

Diversity in Phenological and Agronomic Traits
To evaluate the phenotypic diversity of M. sinensis accessions collected in wide geographic areas from a latitude of N 32° for the southernmost to N 43° for the northernmost in association with miscanthus biomass production, phenological and agronomic traits were measured in field conditions over four years. All the traits showed large variations,

Assessment of Phenological and Agronomic Traits
From the first year after planting, 9 phenological and agronomic traits were assessed every year: first shoot emergence date, heading date, vegetative growth duration, flag leaf/second leaf/third leaf length and width for calculating leaf area, stem height, number of stems, stem dry weight, stem diameter, and biomass yield. The first shoot emergence date was measured when the first shoot emerged from the ground. The heading date was measured when the panicle emerged 1 cm from the tip of the stem. The vegetative growth duration was calculated as the difference in the first shoot emergence date and heading date. The leaf length and width of the second and third leaf from the flag leaf were measured when the panicle was fully emerged. Leaf area was calculated as follows [20]: At maturity in late November, stem height, the number of stems, and stem diameter were measured prior to harvest. Stem height was measured from the soil surface to the tip of the panicle. The number of stems was counted as the number of stems per plant in the space of 0.01 m 2 through the middle of the bunch contributing to the biomass yield (small stems were excluded according to Clifton-Brown and Lewandowski's method [18]). Stem diameter was determined at the midpoint between the second or the third basal nodes using Vernier calipers. All the stems with leaf sheaths were harvested from the whole bunch of the plant and oven-dried at 80 • C for 48 h, and the stem dry weight for the bunch was recorded. To estimate biomass yield (Mg ha −1 ), the dry biomass weight of the bunch was multiplied by the plant density of 13,300 plants ha −1 .

Statistical Analysis
All the data were initially subjected to analysis of variance (ANOVA). Genetic variation and correlation analyses among phenological traits, agronomical traits, and biomass yield of M. sinensis accessions were conducted. Genetic variation among miscanthus accession was estimated using the mean data of agronomical traits recorded during the cultivation period. Correlation efficiency was calculated using Pearson's correlation coefficient. To examine the relationships among selected traits, linear regression analyses were conducted. All statistical analyses were performed using JMP 13 (SAS Institute Inc., Cary, NC, USA).

Diversity in Phenological and Agronomic Traits
To evaluate the phenotypic diversity of M. sinensis accessions collected in wide geographic areas from a latitude of N 32 • for the southernmost to N 43 • for the northernmost in association with miscanthus biomass production, phenological and agronomic traits were measured in field conditions over four years. All the traits showed large variations, as summarized in Tables 1 and 2, Tables S1-S3 and Figure 3, which shows selected major phenological and agronomic traits. From the first year after planting to the fourth year after planting, agronomic traits related to biomass yield gradually increased every year ( Table 2 and Tables S1-S3). Particularly, more notable increases in agronomic traits were observed between the second and the third year.   [28]. ** Leaf area is the sum of the flag leaf area, the second leaf area, and the third leaf area.
Phenological traits are not biomass yield components but indirectly associated with plant biomass, which is mainly determined during vegetative growth duration, starting from the first shoot emergence and ending at the heading [22]. Therefore, variation in the first shoot emergence date and heading date is essential for evaluating phenotypic diversity in the biomass productivity of M. sinensis. In the fourth year after planting, the first shoot emergence date ranged from day 92 (Julian date) to day 100 with an average of day 96.7, showing a narrow variation of 8 days irrespective of collection site (Table 2, Figure 3A). In contrast, the first heading date ranged from day 170 to day 273 with an average of day 228.5, showing a much wider variation of 103 days and having a close association with collection site (Table 2, Figure 3B). The first heading date of Korean accessions collected at latitudes from N 33 • to N 38 • ranged from day 184 to 273, showing a maximum difference of 90 days between the earliest and the latest, suggesting high genetic diversity in heading and flowering time. Accessions from Southern Japan located between N 32 • and N 33 • showed day 205-244 for the first heading date. Accessions from Northeastern China and Far Eastern Russia located between N 42 • and N 43 • showed a much earlier first heading date, from day 170 to day 196, than those collected from lower latitudes, while accessions from Southern Korea and Southern Japan showed a very delayed first heading date ( Table 2). These findings clearly demonstrate that the first heading date is closely associated with the latitude of the collection site. The narrow variation in the first shoot emergence date and the much wider variation in the first heading date result in wide variation in vegetative growth duration. A delayed vegetative growth stage may allow miscanthus to produce more biomass [22,23]. Previous studies also reported that the latitudinal origin of miscanthus accession affected their phenological traits, such as the first heading date and flowering date [7,[30][31][32], resulting in variation in vegetative growth durations and biomass productivity. A similar result was also seen in switchgrass. Casler et al. [33] and Jensen et al. [32] revealed that biomass yields of switchgrass accessions were closely associated with the latitudes of their collection sites, which were also correlated with their phenological and agronomic traits.
The leaf area of all M. sinensis accessions gradually increased each year and showed wide variation depending on the collection site. In the fourth year, leaf areas of the second and the third leaves of M. sinensis accessions from Korea and Japan ranged from 138 to 681 cm 2 with an average of 286.0 cm 2 and from 180 to 434 cm 2 with an average of 280.2 cm 2 , respectively. Those from China and Russia ranged from 115 to 217 cm 2 with an average of 181.7 cm 2 (Table 2, Figure 3C). Miscanthus sinensis accessions from higher latitudes showed a smaller leaf area, while those from lower latitudes showed a larger leaf area, suggesting that leaf growth is negatively correlated with collection site latitude (Tables S4 and S5).  In determining miscanthus biomass yield, stem traits such as number of stems and stem dry weight are important traits as a yield component. Number of stems showed wide ranges of variation but was not correlated with collection site (Tables S4 and S5). In contrast, other stem traits, such as stem height, stem thickness, and stem dry weight, showed correlations with collection site. We found that the higher the latitude is, the smaller are the stem height, stem thickness, and stem dry weight (Tables S4 and S5). Stem height of M. sinensis accessions from higher latitudes (China and Russia) ranged from 116 to 219 cm (average = 178.2 cm), while those from lower latitudes (Korea and Japan) ranged from 82 to 363 cm (average = 237.2 cm) and from 150 to 341 cm (average = 277.1 cm), respectively, with M. × giganteus showing the tallest stem height of 380.9 cm (Table 2, Figure 3E). Stem diameter also showed broad ranges of distribution along the collection sites of miscanthus, ranging from 2.62 to 5.20 mm (average = 4.09 mm), 3.24 to 9.38 mm (average = 6.00 mm), and 4.82 to 7.98 mm (average = 6.78 mm) for China-Russia, Korea, and Japan, respectively (Table 2, Figure 3F), revealing that M. sinensis accessions from higher latitudes produce thinner stems than those from lower latitudes. Stem dry weight is determined by stem height and stem thickness.
The stem dry weight of miscanthus depends on individual stem growth, including stem height and stem thickness, and thus is an important trait in determining biomass yield [16,19,34,35]. Like stem height and stem thickness, stem dry weight also showed a clear relationship with origin latitude (Figure 4). Miscanthus sinensis accessions from lower latitudes showed greater stem dry weight than those from higher latitudes. In the fourth year, the stem dry weight of M. sinensis accessions from Korea and Japan ranged from 6.50 to 81.49 g stem −1 (average = 29.10 g stem −1 ) and 10.85 to 62.13 g stem −1 (average = 37.81 g stem −1 ), respectively Interestingly, the stem dry weight of those from China-Russia was much smaller, from 2.25 to 14.68 g stem −1 (average = 8.46 g stem −1 ), than those from the lower latitudes of Korea and Japan, suggesting that M. sinensis accessions from lower latitudes have greater potential in stem growth and consequently biomass yield (Table 2, Figure 3G Figure 3H). Some of the M. sinensis accessions from Korea and Japan showed a greater biomass yield than M. × giganteus, of which the biomass yield was 18.50 Mg ha −1 .
Overall, M. sinensis accessions from Korea and its neighboring countries showed high phenotypic diversity in their phenological and agronomic traits, except the first shoot emergence date (Tables 1 and 2, Tables S1-S3). Wide variation in heading date suggests that heading date determines the vegetative growth duration of miscanthus, resulting in wide variation in agronomic traits such as leaf area, stem height, stem diameter, and stem dry weight related to vegetative growth. Stem growth traits associated with stem and biomass also showed high phenotypic diversity, with a particular correlation between stem growth traits except the number of stem and the latitude of the collection site, as reported in previous studies [30,32]. The high diversity of our accessions in their phenotypic traits provides an opportunity to select an elite germplasm for a future breeding program.

Relationships between Latitude of Collection Site and Phenotypic Traits
For a better understanding of how the latitude of the collection site is related to phenological and agronomic traits, correlation and linear regression analyses were conducted. Correlation analyses between traits assessed in the fourth year after planting and latitude revealed that most phenotypic traits were related to the latitude of the collection site, except the first shoot emergence date (Table S5). All phenological and agronomic traits except number of stems were negatively correlated with latitude. In particular, heading date, vegetative growth duration, leaf area, stem diameter, and stem dry weight showed greater negative correlations with r approximately −0.6 than other traits such as the first shoot emergence date, suggesting that accessions collected from lower latitudes may display greater value and possess greater potential in these traits. In contrast, the number of stems showed a positive correlation with latitude, indicating that the higher the latitude is, the greater is the number of stems. However, unlike our results, other studies reported that the number of stems was highly correlated with biomass yield and a good selection trait for estimating biomass yield of miscanthus [18,19].
Linear regression analyses between latitude and phenological and agronomic traits reconfirmed that heading is delayed with decreasing latitude ( Figure 4B), resulting in increased vegetative growth duration, as the first shoot emergence date was not related to latitude ( Figure 4A). Leaf area ( Figure 4C), stem height ( Figure 4E), stem diameter ( Figure 4F), and stem dry weight ( Figure 4G) also increased with decreasing latitude, resulting in increased biomass yield ( Figure 4H). These findings suggest that accessions collected at lower latitudes have greater potential in biomass yield, although variation is wide among accessions collected in the same latitude. Studies conducted in Korea and China also showed similar results: accessions collected in lower latitudes produce greater biomass than those collected in higher latitudes [17,30,34]. Accessions collected at the higher latitudes of China-Russia started seedling emergence in the spring at a similar time to other accessions from lower latitudes ( Figure 4A) but started heading much earlier, with no great variation, almost 40-50 days earlier than accessions from lower latitudes ( Figure 4B), resulting in greatly shortened vegetative growth duration. Variation in phenological traits such as heading date and vegetative growth duration along latitude may affect the quantitative growth of miscanthus and lead to different performance in biomass components and biomass yield, although other factors may also be involved. Robson et al. [21][22][23] commented that the biomass yield of miscanthus was affected by complex associations among environmental factors, morphological traits, and phenological traits. A significant association between heading date and latitude of collection site suggests photoperiodic and thermal influences on heading and flowering date [7,[30][31][32]. Similar cases were also found in other bioenergy grasses such as switchgrass, whose phenotypic traits were also affected by latitude of collection site, a result of latitudinal adaptation [33,36]. The significant relationship between latitude and agronomic traits (Figure 4) suggests that accessions collected from different geographical latitudes will provide more phenotypically diverse germplasms, which can be used for breeding new miscanthus varieties adapting to various latitudinal locations. However, prior to using them for breeding, it is necessary to investigate their genetic relationship in association with phenotypic variation.

Relationships between Phenological and Agronomic Traits in Association with Biomass Yield
To evaluate relationships between phenological and agronomic traits in association with biomass yield, we conducted correlation and linear regression analyses between phenological traits, agronomic traits, and biomass yields of M. sinensis accessions measured in the second and fourth year after miscanthus planting ( Figure 5 and Tables S4 and S5).  Tables S4 and S5 show correlations between phenological traits, agronomic traits, and biomass yield assessed in the second year and fourth year after miscanthus planting, respectively. With the exception of the first emergence date and the number of stems, the other traits showed significant correlations with biomass yield consistently each year. In particular, stem dry weight (r = 0.619 *** and 0.467 ***), leaf area (r = 0.513 *** and 0.406 **), stem height (r = 0.467 *** and 0.373 **), heading date (r = 0.378 ** and 0.373 **), vegetative growth duration (r = 0.441 ** and 0.367 **), and stem diameter (r = 0.260 ** and 0.393 **) were significantly correlated with biomass yield assessed in the second and fourth year after miscanthus planting, respectively. These results suggest that these traits are important for biomass formation in miscanthus. In determining miscanthus biomass yield, the number of stems and stem dry weight are considered important traits as yield components [16,18,30,37]. However, in our study, the number of stems showed no significant correlation with other agronomic traits and biomass yield. This difference might be due to differences in counting the number of stems. We counted the number of stems in the space of 0.01 m 2 through the middle of the bunch, while the other studies measured the total number of stems. In contrast, stem dry weight showed a strong correlation with the main phenological and agronomic traits ( Figure 5). In the second year after miscanthus planting, heading date, vegetative growth duration, leaf area, and stem diameter were strongly correlated with stem dry weight, which was strongly correlated with biomass yield (r = 0.619 ***) but not much correlated with number of stems ( Figure 5A). These correlations were maintained similarly in the fourth year after miscanthus planting ( Figure 5B), suggesting that six phenotypic traits (except number of stem) are key traits determining miscanthus biomass yield.
Linear regression analyses revealed that stem diameter is the most significant contributor to stem dry weight (r 2 = 0.830), followed by leaf area (r 2 = 0.613), heading (r 2 = 0.321), vegetative growth duration (r 2 = 0.305), and stem height (r 2 = 0.279) ( Figure 6). Our findings that leaf area and stem dry weight were most strongly related to biomass suggest that leaf area plays an important role in stem growth, particularly stem diameter, and biomass accumulation. Previous studies reported that vegetative growth durations, stem height, stem diameter, and stem dry weight were strongly correlated with biomass yield potential [8,11,[16][17][18][19][20][21][22][23]30,34,38]. However, the linear relationships between biomass yield and phenotypic traits were not strongly supported by the regression coefficients ( Figure 6C,F,I,L,O). It can be assumed that biomass yield is affected by complex interactions of environmental factors and various phenotypic traits. Clifton-Brown et al. [7] and Robson et al. [21] reported that miscanthus biomass yield is a combined result of complex associations with environmental factors and phenotypic traits. Although our results suggested that miscanthus biomass is more attributed to leaf area, stem diameter, and stem weight than the other traits, these findings do not mean that the other traits should not be taken into account. Our results imply that the miscanthus accessions tested in this study are genetically very diverse.

Early Traits Determining Biomass Yield Potential of Miscanthus
Producing high biomass is the main goal for miscanthus cultivation. Moreover, miscanthus requires at least three to four years after planting for commercial biomass production, and phenotyping of agronomic traits associated with biomass yield requires more than three years in the field, costing significant time, effort, labor, and money. Additionally, breeding miscanthus has faced several obstacles, such as a lack of available germplasm and complex interactions of phenotypic traits in determining biomass. Thus, finding key traits determining miscanthus biomass yield is important in miscanthus breeding. Many studies have focused on those traits that are associated with biomass yield productivity [16,17,19,38], but they investigated relationships between traits and biomass yields assessed in the same year. However, if certain traits measured in the early stage of establishment are related to biomass yield potential, which can usually only be assessed in the fourth year after miscanthus planting in the field, we can significantly shorten the period for screening candidate lines.    To examine early phenotypic traits that are closely associated with the fourth-year biomass yield, we compared phenotypic traits assessed in the second year and the fourth year after miscanthus planting. Significant correlations exist between the same traits assessed in the second and fourth year after planting: heading date (r = 0.943 ***), vegetative growth duration (r = 0.909 ***), leaf area (r = 0.470 ***), stem height (r = 0.270 **), stem diameter (r = 0.869 ***), number of stem (r = 0.995 ***), and stem dry weight (r = 0.786 ***) ( Table 3). Correlation analyses between the main phenotypic traits assessed in the second year and biomass yield assessed in the fourth year showed significant correlations, except number of stems (Table 3, Figure 7). As expected, the second-year biomass yield showed the greatest positive correlation with the fourth-year biomass yield (r = 0.785 ***). Among individual phenotypic traits, stem dry weight showed the next greatest positive correlation with the fourth-year biomass yield (r = 0.525 ***), suggesting that stem dry weight is one of the most significant contributors determining miscanthus biomass yield. Other traits such as vegetative growth duration, leaf area, stem height, and stem diameter showed similar correlations with the fourth-year biomass yield (Figure 7). Anzoua et al. [16] investigated the relationship between agronomic traits and biomass yield potential and revealed that biomass production was strongly correlated with stem height and leaf length. The other studies reported that stem diameter, stem height, and heading date showed significant correlations with biomass yield [19,34,35]. Clifton-Brown and Lewandowski [18] suggested that at least 2 years of field trials are required to predict the biomass yield of miscanthus and that first-year data are insufficient to predict biomass yield. These previous studies investigated the relationships between phenotypic traits and biomass yield assessed in the same year, while our study investigated the relationships between the second-year phenotypic traits and the fourth-year biomass yield. Our results suggest that potential biomass yield can be predicted in an early year of miscanthus cultivation-for example, the second year after planting. The early traits showing high correlations with the fourthyear biomass yield may be utilized as selection traits for screening breeding lines with high biomass yield potential in a miscanthus breeding program, saving time and effort in phenotyping and screening.

Conclusions
In this study, we found that our miscanthus accessions showed high diversity in phenological and agronomic traits, and these traits were strongly correlated with latitude of Figure 7. Schematic representation of the correlation between phenological (heading date) and main agronomic traits assessed in the 2nd year after miscanthus planting, and biomass yield assessed in the 4th year after planting. *, **, and *** indicate statistical significance at 0.05, 0.01, and 0.001 levels, respectively.

Conclusions
In this study, we found that our miscanthus accessions showed high diversity in phenological and agronomic traits, and these traits were strongly correlated with latitude of collection site. The best biomass yield was mainly observed in M. sinensis accessions from the southern parts of Korea, such as Jeolla and Jeju provinces, with over 20 Mg DM ha −1 , which may be mainly attributed to longer vegetative growth. Our study clearly demonstrated that latitude of collection site affects phenological traits, particularly heading date, and consequently vegetative growth, which determine the biomass yield of miscanthus in a combined manner. Relationship analyses between earlier traits assessed in the second year and later traits assessed in the fourth year suggest that earlier traits may be useful parameters to estimate relevant parameters and potential biomass yield in perennial crops including miscanthus. Among traits assessed in the second year after miscanthus planting, leaf area of the second and third leaves from the flag leaf was closely associated with stem diameter and stem dry weight. Stem diameter and stem dry weight assessed in the second year were highly correlated with those traits and biomass assessed in the fourth year, suggesting that leaf area, stem diameter, and stem dry weight are key traits in estimating potential biomass yield. These traits can be used as a target trait for screening miscanthus breeding lines with high biomass yield potential. Considering our results, it can be concluded that the selected traits, such as vegetative growth duration, leaf area, stem diameter, and stem dry weight, observed in the second year after miscanthus planting may help us to predict the potential biomass yield of miscanthus and screen candidate lines with our target traits in a miscanthus breeding program. It is necessary to validate our findings and the phenotypic traits proposed in other independent field studies using different plant materials.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/agronomy11050900/s1, Table S1: Summary of phenological and agronomic traits of miscanthus accessions. Phenological and agronomic traits were assessed in the 1st year after miscanthus planting, Table S2: Summary of phenological and agronomic traits of miscanthus accessions. Phenological and agronomic traits were assessed in the 2nd year after miscanthus planting, Table S3: Summary of phenological and agronomic traits of miscanthus accessions. Phenological and agronomic traits were assessed in the 3rd year after miscanthus planting, Table S4: Correlation coefficients between latitude of miscanthus collection site, phenological traits, agronomic traits and biomass yield assessed in the 2nd year after miscanthus planting, Table S5: Correlation coefficients between latitude of miscanthus collection site, phenological traits, agronomic traits and biomass yield assessed in the 4th year after miscanthus planting.