Distinct Geographical Distribution of the Miscanthus Accessions with Varied Biomass Enzymatic Saccharification

Miscanthus is a leading bioenergy candidate for biofuels, and it thus becomes essential to characterize the desire natural Miscanthus germplasm accessions with high biomass saccharification. In this study, total 171 natural Miscanthus accessions were geographically mapped using public database. According to the equation [P(H/L| East) = P(H/L∩East)/P(East)], the probability (P) parameters were calculated on relationships between geographical distributions of Miscanthus accessions in the East of China, and related factors with high(H) or low(L) values including biomass saccahrification under 1% NaOH and 1% H2SO4 pretreatments, lignocellulose features and climate conditions. Based on the maximum P value, a golden cutting line was generated from 42°25’ N, 108°22’ E to 22°58’ N, 116°28’ E on the original locations of Miscanthus accessions with high P(H|East) values (0.800–0.813), indicating that more than 90% Miscanthus accessions were originally located in the East with high biomass saccharification. Furthermore, the averaged insolation showed high P (H|East) and P(East|H) values at 0.782 and 0.754, whereas other climate factors had low P(East|H) values, suggesting that the averaged insolation is unique factor on Miscanthus distributions for biomass saccharification. In terms of cell wall compositions and wall polymer features, both hemicelluloses level and cellulose crystallinity (CrI) of Miscanthus accessions exhibited relative high P values, suggesting that they should be the major factors accounting for geographic distributions of Miscanthus accessions with high biomass digestibility.


Introduction
Miscanthus, a typical C4 perennial grass with extremely high biomass yield, is currently considered a leading bioenergy crop for the production of biofuel feedstock [1][2][3][4]. Miscanthus has been historically utilized in the industries of forage, textile, and shelter. As a fast-growing grass, it requires little water and fertilization, and can thus grow in poor soil with minimal inputs [5,6]. Miscanthus is a genus of about 20 species and distributes widely under various climate conditions across Asia and the nearby Pacific islands [5]. In China, at least seven species have been found, and each species has a unique geographical distribution [7][8][9][10]. Because of its selfincompatible property, Miscanthus has a diverse genetic variation. Recent reports indicated that Chinese Miscanthus accessions exhibited diverse cell wall compositions and wall polymer features for largely varied biomass saccharification [10,11].
Plant cell walls represent a large biomass resource for biofuel and other chemical products. However, lignocellulose recalcitrance greatly affects cost-effectiveness of biomass processing that involves in three major steps: physical and chemical pretreatment, enzymatic hydrolysis, and yeast fermentation [12,13]. Recalcitrance is determined mainly by the cell wall composition and wall polymer features [14][15][16][17]. Genetic improvement of energy crops is proposed as a promising solution for reducing recalcitrance [10,17,18].
Cellulose is a major cell wall polysaccharide that accounts for 20% to 35% of the dry matter of grasses, and its crystallinity is reportedly a negative factor affecting biomass hydrolysis in plants [14,19,20]. Xylan is a major hemicellulose that positively affects biomass enzymatic saccharification by reducing lignocellulose crystallinity (CrI) in Miscanthus [15,16]. Furthermore, the arabinose (Ara) substitution degree (reverse Xylose/Ara) of xylan could enhance biomass saccharification under various chemical pretreatments in Miscanthus [16]. Lignin is an amorphous heteropolymer containing three major phenolic components: p-coumaryl alcohol (H), coniferyl alcohol (G), and sinapyl alcohol (S) [21,22]. Recent reports have indicated that lignin is a negative factor on biomass digestibility under various physical and chemistry pretreatments in Miscanthus, the S/G ratio determines the negative effect of lignin on biomass digestibility in Miscanthus [15,23,24].
In this work, we selected 171 Miscanthus accessions that have been reported to represent four major species from a large natural Miscanthus collection (n = 1400) around China [10,11], and mapped the geographical distribution of each species. Also, we performed a correlation analysis of cell wall composition/wall polymer features, biomass saccharification, and meteorological factors at the geographical location of each species. Furthermore, we characterized and selected the Miscanthus accessions for high biomass enzymatic digestion. These informations will be useful forMiscanthus breeding toward the development of desirable bioenergy crop.  (Table 1, Fig 1 and S1 File). The samples harvested from Hunan experimental field in the 2009 and 2010 seasons were dried at 50°C after treated at 105°C for 5 min. The dried tissues were ground through a 40 meshes screen and stored in a dry container until use.

Wall polysaccharide extraction
The plant cell wall fractionation method was used for extraction of cellulose and hemicelluloses as described by Peng LC et al. [25] and Xu N et al. [15]. The crude cell wall material was suspended in 0.5% (w/v) ammonium oxalate and heated for 1 h in a boiling water bath, and the supernatants were combined as total pectin. The remaining pellet was suspended in 4 M KOH containing 1.0 mg /mL sodium borohydride for 1 h at 25°C, and the combined supernatant was neutralized, dialyzed and lyophilized as hemicelluloses. The KOH non-extractable residue was further extracted with acetic-nitric acids for 1 h at 100°C and the remaining pellet was defined as crystalline cellulose. All experiments were carried out in biological triplicate.

Colorimetric assay for hexoses and pentoses
UV/VIS Spectrometer (Shanghai MAPADA Instruments Co., Ltd. V-1100D) was used for total hexose and pentose measurement. Hexose content was determined by the anthrone (Sigma-Aldrich Co. LLC) /H 2 SO 4 method [26]. Cellulose was dissolved in 67% (v/v) H 2 SO 4 (1.0 mL) with shaking at 25°C for 1 h, and then 10.0 0μL aliquot was used to determine cellulose content and hexoses yield released from pretreatment and enzymatic hydrolysis using the anthrone/ H 2 SO 4 method. Pentose content was measured by the orcinol/HCl method [26], and the assay was used to determine pentose yield released from pretreatment and enzymatic hydrolysis.
Because the high pentose level in the sample can affect the absorbance reading at 620 nm for hexose content by anthrone/H 2 SO 4 method, the deduction from pentose reading at 660 nm was carried out for final calculation of hexose level. All experiments were carried out in triplicate.

Detections of hemicellulosic monosaccharides and cellulose crystallinity
Trifluoroacetic acid (TFA) and myo-inositol were purchased from Aladdin Reagent Inc. Acetic anhydride and acetic acid were obtained from Sinopharm Chemical Reagent Co., Ltd. 1-Methylimidazole was purchased from Sigma-Aldrich Co. LLC. Monosaccharide standards including l-rhamnose, l-arabinose, l-fucose, d-xylose, d-galactose, d-glucose and d-mannose, were obtained from Sinopham Chemical Reagent Co., Ltd. The combined supernatants from 4 M KOH fraction were dialyzed for 36 h after neutralization with acetic acid. The polysaccharides dissolved in 2.5 mL TFA (2 M) were heated in a sealed tube at 121°C in an autoclave (15 psi) for 1 h. Myo-inositol (200 μg) was added as the internal standard. The supernatant was dried under vacuum at 38°C to remove TFA. Distilled water (800 μL) and a freshly prepared solution of NaBH 4 (400 μL, 100 mg/mL in 6.5 M aqueous NH 3 ) were added to each sample. Sample was capped, mixed well and incubated at 40°C for 30 min. Excess NaBH 4 was decomposed by adding acetic acid (800 μL). 400 μL Sample was then moved into a 25 mL glass tube. Acetic anhydride (4 mL) was added to the tube and the solution mixed again. Then 1-methylimidazole (600 μL) was added. After mixing, the sample was allowed to stand for 10 min. Excess acetic anhydride was decomposed by adding distilled water (10 mL). Then dichloromethane (3 ml) was added, mixed gently, centrifuged (2,000 g, 10 seconds) for phase separation. After removing the upper phase, the sample was washed with distilled water (3 × 20.0 mL). The collected lower phase was dehydrated by adding with anhydrous sodium sulfate and stored at −20°C until analyzed by GC-MS. The sample preparation and GC-MS (SHIMADZU GCMS-QP2010 Plus) running were previously described by Li FC et al. [16]. X-ray diffraction (XRD) method was used to determine cellulose crystallinity index (CrI) using Rigaku-D/MAX instrument (Uitima III, Japan) as described by Zhang W et al. [14]. The raw biomass powder was laid on the glass sample holder (35 × 50 × 5 mm) and detected under plateau conditions. Ni-filtered Cu Kα radiation (λ = 0.154056 nm) generated at voltage of 40 kV and current of 18 mA, and scanned at speed of 0.0197°/s from 10-45°. The crystallinity index (CrI) was estimated using the intensity of the 200 peak (I 200 , θ = 22.5°) and the intensity at the minimum between the 200 and 110 peaks (I am , θ = 18.5°) as the follow: CrI = 100 × (I 200 −I am )/I 200 . I 200 represents both crystalline and amorphous materials while I am represents amorphous material. Standard error of the CrI method was determined at ±0.05~0.15 using five representative samples in triplicate.

Total lignin and monolignol assay
Total lignin was determined by the two-step acid hydrolysis method according to Laboratory Analytical Procedure of the National Renewable Energy Laboratory. The lignin includes acidinsoluble and -soluble lignin. The acid-insoluble lignin was calculated gravimetrically after correction for ash, and the acid-soluble lignin was measured by UV spectroscopy. For acid-insoluble lignin determination, a 0.5 g sample recorded as W1. Each sample was run in triplicate. The sample was extracted with benzene-ethanol (2:1, v/v) in a Soxhlet for 4 h, and then air-dried in hood overnight. The sample was hydrolyzed with 10 mL 72% H 2 SO 4 (v/v) in shaker at 30°C for 1.5 h. After hydrolysis, the acid was diluted to a concentration of 2.88%, and then placed in the autoclave for 1 h at 121°C (15 psi). The autoclaved hydrolysis solution was vacuum-filtered through the previously weighed filtering crucible. The filtrate was captured in a filtering flask for acid-soluble lignin. The lignin was washed free of acid with hot distilled water and the crucible and acid-insoluble residue was dried in an oven at 80°C until constant weight was achieved. Then, the samples were removed from the oven and cool in a dry-container. The weight of the crucible and dry residue was recorded to the nearest 0.1 mg (W2). The dried residue was ashed in the muffle furnace at 200°C for 30 min and 575°C for 4 h. The crucibles and ash were weighed to the nearest 0.1 mg and recorded the weight (W3). Acid-insoluble lignin (AIL) on original sample was calculated as the following: AIL (%) = (W2-W3) × 100/W1. For acid-soluble lignin determination, it was solubilized during the hydrolysis process, and was measured by UV spectroscopy. The hydrolysis liquor obtained previously was transfer into 250 mL volumetric flask and brought up to 250 mL with 2.88% sulfuric acid. The absorbance of the sample was read at 205 nm on a UV-vis spectroscopy (Beckman Coulter Inc., Du800), and 2.88% sulfuric acid was used as blank. The method of calculation about the amount of acid soluble lignin was as follows: ASL (%) = (A × D × V/1000 × K × W1) × 100%. A (absorption value), D (Dilution ratio of the sample), K (absorptivity constant) = 110 L/g/cm. Total lignin (%) = ASL% + AIL%. All experiments were carried out in triplicate.
Three monolignols were measured using HPLC as described by Xu N et al. [10]. The standard chemicals, p-Hydroxybenzaldehyde (H), vanillin (G) and syringaldehyde (S) were purchased from Sinopharm Chemical Reagent Co., Ltd. The sample was extracted with benzeneethanol (2:1, v/v) in a Soxhlet for 4 h, and the remaining pellet was collected as cell-wall residue (CWR). The procedure for nitrobenzene oxidation of lignin was conducted as follows: 0.05 g CWR was added with 5 mL 2 M NaOH and 0.5 mL nitrobenzene, and a stir bar was put into a 25-mL Teflon gasket in a stainless steel bomb. The bomb was sealed tightly and heated at 170°C (oil bath) for 3.5 h and stirred at 20 rpm. Then, the bomb was cooled with cold water. The chromatographic internal standard (ethyl vanillin) was added to the oxidation mixture. This alkaline oxidation mixture was washed three times with 30 mL CH 2 C1 2 /ethyl acetate mixture (1:1, v/v) to remove nitrobenzene and its reduction by-products. The alkaline solution was acidified to pH 3.0 to 4.0 with 6 M HCl, and then extracted with CH 2 C1 2 /ethyl acetate (3 × 30 mL) to obtain the lignin oxidation products, which were in the organic phase. The organic extracts were evaporated to dryness under reduced pressure at 40°C. The oxidation products were dissolved in 10 mL chromatographic pure methanol. For HPLC analysis the solution was filtered with a membrane filter (0.22 μm). Then, 20 μL solution was injected into the HPLC (Waters 1525 HPLC) column Kromat Universil C18 (4.6 mm × 250 mm, 5 μm) operating at 28°C with CH 3 OH:H 2 O:HAc (25:74:1, v/v/v) carrier liquid (flow rate: 1.1 mL/minute). Calibration curves of all analytes routinely yielded correlation coefficients 0.999 or higher, and the detection of the compounds was carried out with a UV-detector at 280 nm.

Biomass pretreatment and enzymatic hydrolysis
Biomass pretreatments with NaOH and H 2 SO 4 and sequential enzymatic hydrolysis were performed as described by Huang JF et al. [11] and Xu N et al. [15]. For NaOH pretreatment, the well-mixed powder of the biomass sample (0.5 g) was added with 10 mL NaOH at 1% (w/v). The tube was shaken at 150 rpm for 2 h at 50°C, and centrifuged at 3,000 g for 5 minutes. The pellet was washed three times with 10 mL distilled water, and stored at −20°C for enzymatic hydrolysis. All supernatants were collected for determination of total sugars (pentoses and hexoses) released from alkali pretreatment, and samples with 10 mL distilled water were shaken for 2 h at 50°C as the control. All samples were carried out in biological triplicate. For H 2 SO 4 pretreatment, the well-mixed powder of the biomass sample (0.5 g) was added with 10 mL H 2 SO 4 at 1% (v/v). The tube was sealed and heated at 121°C for 20 minutes in an autoclave (15 psi) after mixing well. Then, the tube was shaken at 150 rpm for 2 h at 50°C, and centrifuged at 3,000 g for 5 minutes. The pellet was washed three times with 10 mL distilled water, and stored at −20°C for enzymatic hydrolysis. All supernatants were collected for determination of total sugars released from acid pretreatment, and samples with 10 mL distilled water were shaken for 2 h at 50°C as the control. All samples were carried out in biological triplicate.
The remaining residues from various pretreatments were washed 2 times with 10 mL distilled water, and once with 10 mL mixed-cellulases reaction buffer (0.2 M acetic acid-sodium acetate, pH 4.8). The washed residues were added with 10 mL(2 g/L) mixed-cellulases (containing β-glucanase 6 × 10 4 U) and cellulase 600 U and xylanase 10 × 10 4 U from Imperial Jade Bio-technology Co., Ltd). During the enzymatic hydrolysis, the samples were shaking under 150 rpm at 50°C for 48 h. After centrifugation at 3,000 g for 10 min, the supernatants were collected for determining amounts of pentoses and hexoses released from enzymatic hydrolysis. The samples with 10 mL reaction buffer were shaken for 48 h at 50°C as the control. All samples were carried out in triplicate.

Acquisition and analysis of climate data
The annual data for average and minimum earth's surface temperature, solar radiation, and rainfall for 171 locations in China were obtained from the National Aeronautics and Space Administration (NASA) Surface Meteorology and Solar Energy (SSE 6.0) website [27] (https:// eosweb.larc.nasa.gov/sse/). The meteorological factors include average and minimum temperatures (°C), rainfall data over the 22-year period from Jan 1983 to Dec 2004; and solar radiation from July 1983 to Jun 2005.

Draw the geographical distribution of the Miscanthus on map of China
The geographical distribution of the Miscanthus accessions were drawn in statistical computing and graphics language R (https://www.r-project.org/) using "maps" and "mapdata" packages.

Determination of the golden cutting line and Bayes probabilities
Based on the biomass enzymatic digestibility under pretreatment with 1% NaOH or 1% H 2 SO 4 , the 171 Miscanthus accessions were equally divided into three groups with each containing 57 accessions. We mapped the geographic locations of Group I and Group III accessions, which had relatively high (H) and low (L) biomass digestibility, respectively. The boundary of high and low biomass digestibility samples on Chinese map are calculated using the Bayes theorem, and the method is known as the (empirical) Bayes method. It is a proper and standard statistical method. Suppose that the boundary divide China into two parts. Each part contains high and low biomass digestibility samples. To ensure the P(High | Area), a conditional probability of High given Area, is the maximum value by calculating the P parameter based on the equation: P(High | Area) = P(High\Area) / P(Area). Then we got the cutting line

Distinct geographic locations of four Miscanthus species
A total of 171 representative Miscanthus accessions of four major species were selected from the collection of over 1400 natural Miscanthus accessions, including 78 M. sinensis, 36 M. floridulus, 32 M. sacchariflorus, and 25 M.lutarioriparius (Fig 1A). In general, the Miscanthus  Table 1 and S1 File). However, each species exhibited a typical geographic location. M. sacchariflorus mainly distributed in the north, whereas almost all M. floridulus accessions were located in the south. In addition, M. sinensis had a broad distribution and M. lutarioriparius showed a narrow location. Due to the distinct geographic location, the Miscanthus accessions of each species were originally grown at typical climate conditions including average temperature, minimum temperature, average insolation, and average precipitation (Table 1). For instance, M. floridulus mainly grew under a relatively high average temperature (18.44°C) and precipitation (4.09 mm/day), whereas most of the M. sacchariflorus accessions exposed to high average insolation (4.05 kWh/m2/day) and low average precipitation (2.05 mm/day). By comparison, M. sinensis exhibited relatively large variation in average and minimum temperatures (CV%: 26.85%, 46.15%), while M. lutarioriparius had small variation (3.55%, 5.86%). Hence, the four major Miscanthus species had distinct geographic distributions under typical climate conditions, suggesting that each Miscanthus species should grow under its typical ecosystem.

Diverse cell wall compositions of four Miscanthus species
The Miscanthus accessions of the four species exhibited different cell wall compositions in addition to distinct geographic locations (Fig 1A, Table 2 and S1 File). For instance,

Biomass digestibility in the four Miscanthus species
Biomass digestibility (saccharification) was measured by either the hexose yield (hexoses/cellulose) released from the hydrolysis by a crude cellulase mixture of lignocellulose after pretreatment or the total sugar yield (hexoses and pentoses/dry weight) from both pretreatment and enzymatic hydrolysis (Table 3 and S1 File). In this study, we determined the sugar yield of 171 Miscanthus samples released after pretreatment with 1% NaOH or 1% H 2 SO 4 . As a result, these four Miscanthus species showed varied hexose yields with a CV% ranged from18.10% to 33.06%, and total sugar yields with a CV% from 11.87% to 46.05%. M. sacchariflorus had the highest average hexose yield, whereas M. lutarioriparius showed the lowest average hexose yield upon pretreatment with 1% NaOH or 1% H 2 SO 4 . Furthermore, M. sinensis exhibited the highest total sugar yield at 51% under pretreatment with 1% NaOH or 1% H 2 SO 4 among the four species. These results suggest that diverse cell wall compositions lead to varied biomass digestibility in four Miscanthus species, consistent with the previous reports [14-16, 23, 28].

Golden cutting line of Miscanthus distribution for high biomass saccharification
Based on the biomass enzymatic digestibility under pretreatment with 1% NaOH or 1% H 2 SO 4 , the 171 Miscanthus accessions were equally divided into three groups with each containing 57 accessions. We mapped the geographic locations of Group I and Group III accessions, which had relatively high (H) and low (L) biomass digestibility, respectively (Fig 2A and  2B). Surprisingly, we found that almost all Miscanthus accessions in Group I located in the . Based on the cutting line with P (H|East) values at 0.800 (1% NaOH) and 0.813 (1% H 2 SO 4 ), we identified that 52 Miscanthus accessions of Group I and 12-13 accessions of Group III were located in the east (Fig 2A and  2B). As a comparison, 91% of the Miscanthus accessions with high biomass saccharification distributed in the east, whereas only 21-23% of the Miscanthus accessions with low biomass saccharification distributed in the east. Hence, the results have demonstrated that the golden cutting line generated from the maximum P value, should be applicable for screening out of the desire Miscanthus accessions with high biomass saccharificaiton in the suitable growing regions.

Typical climate factors at the locations of the Miscanthus accessions with high biomass saccharification
In this study, we initially collected data about 26 meteorological factors for the original locations of the Miscanthus accessions. Four of them, including average insolation (kWh/m2/day), average and minimum temperatures (°C), and average precipitation (mm/day) were selected for probability (P) parameter analysis of the Miscanthus accessions with high biomass saccharification. The 171 Miscanthus accessions were equally divided into three groups based on the Using the cutting line of the Miscanthus locations described above, we calculated the P parameter for the Miscanthus accessions (n = 114) (Group I and Group III) exposed to high (H) average insolation in the east with P(H|East) and P(East|H) values at 0.782 and 0.754, respectively (Fig 3). These two P values indicated that the Miscanthus accessions located in the east of China have 78% probability of exposing to high average insolation or about 75% of the Miscanthus accessions exposing to high average insolation are located in the east. Furthermore, we observed that 77% (1% NaOH pretreatment) and 70% (1% H 2 SO 4 pretreatment) of the Miscanthus accessions showing high biomass saccharifcaition had been historically exposed to high average insolation in the east. In addition, we calculated P(H|East) and P(East|H) values on the average and minimum temperatures and average precipitation (Additional files 3, 4 and 5). These three meteorological factors exhibited low P(East|L) and P(East|H) values ranging from 0.316 to 0.667, indicating that they did not play a major role in biomass enzymatic digestibility of the  Table 4). As the levels of the two cell wall polymers (cellulose and lignin) and three major polymer features (CrI, Xyl/Ara, and S/G) have been determined as the negative factors on biomass saccahrification in Miscanthus, we calculated their P(East|L) values, but P(East|H) was generated for hemicellulose due to its positive effect. As a result, relatively high P values (0.684 and 0.789) were obtained for the cellulose and hemicellulose, indicating that either low (L) cellulose level or high (H) hemicellulose contents of the Miscanthus accessions were located mainly in the east. By comparison, a much lower P value (0.456) was observed for lignin, suggesting that lignin might not correlate with the Miscanthus distribution. Among the three cell wall polymer features (Table 4), cellulose CrI showed the highest P(East|L) value (0.739), whereas Xyl/Ara of hemicellulose and S/G of lignin had P values at 0.609 and 0.696, respectively. Hence, hemicellulose level and cellulose CrI should be the major factors associated with the original locations of the Miscanthus accessions with high biomass enzymatic digestibility.
Miscanthus biomass is increasingly regarded as one of leading biofuel feedstocks. Genetic improvement of the Miscanthus species will enhance their production and utilization as bioenergy crops. This study characterized four major Miscanthus species for their diverse geographical distribution, distinct cell wall compositions/wall polymer features, and biomass enzymatic saccharification. This study at first time determined the relationship of those three cell wall characteristics of the Miscanthus species with the geographical distribution and related climatic factors. These information will be very useful for the collection of Miscanthus germplasms and breeding desired Miscanthus varieties for biofuel production Notably, the golden cutting line (42°25' N, 108°22' E to 22°58' N, 116°28' E) of the Miscanthus distribution defined in this study will be helpful for selecting Miscanthus accessions with high biomass enzymatic digestibility in the east of China, which could be applied to make national policy for potentially developing bioenergy Miscanthus crops in China. Furthermore, as the average insolation is the unique positive factor on biomass saccharification among the four meteorological factors for the Miscanthus distribution, it should be considered as an important parameter for growing Miscanthus crops over the world.
Plant cell wall composition and wall polymer feature have significant effects on biomass saccharification. This study revealed that high hemicellulose content and low cellulose CrI level are major factors for enhancing biomass enzymatic digestion of Miscanthus. However, as hemicelluloses positively affect biomass saccharification by reducing cellulose CrI in Miscanthus [15], the cellulose crystallinity should be considered a primary target for genetic modification of plant cell walls in Miscanthus. In terms of the lignin level without distinct effect on Miscanthus distributions, we inferred that it is due to lignin playing important roles in plant resistance to various environmental stresses. This study provides significant information on the germplasm accessions under four major Miscanthus species, including biomass enzymatic saccharification, cell wall compositions, and their correlation with geographical distribution and climatic factors (Table 5, Additional file 2). These information will be useful for the plant geneticists and breeders to develop superior bioenergy crops from Miscanthus species. Also, they will be helpful guidelines for policy makers to develop sustainable strategy for biofuel production in China and worldwide.

Conclusions
Based on the probability parameter calculation, a golden cutting-line has been generated to identify Miscanthus accessions located in the east of China with high biomass saccharification. It has indicated that either the averaged insolation or the hemicelluloses level and lignocelluose crystallinity are the main factors affecting biomass digestibility of Miscanthus accessions. This study can suggest the potential strategy for developing dedicated bioenergy Miscanthus crops in China and over the world.