Correlation between Genetic Characteristics, Cell Structure and Material Properties of Moso Bamboo (Phyllostachys edulis (Carriere) J. Houzeau) in Different Areas of China

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau), native to China, is one of the most economically and ecologically important bamboo species. Since the economic interests and the strong clonality, it has been widely cultivated in southern China, which inevitably reduces the natural stands and leads to gene loss in this species. In this study, three natural populations of Moso bamboo distributed in Anhui, Guangxi, and Zhejiang province, were used to analyze the correlation between phenotypic traits, cell structure, and material properties from the perspective of phenotypic, genetic, and environmental. Among those traits and properties, fiber width was correlated with wall thickness at breast height and average nodes length under branch positively. Leaf length was correlated positively with fiber lumen diameter and parenchyma lumen diameter. Furthermore, it showed a very close correlation between moisture content, bending strength, modulus of elasticity, and diameter at breast height, clear height, and leaf length. The lumen diameter of fiber cell wall thickness is positively correlated with bending strength and modulus of elasticity. Density is positively correlated with parenchyma cell wall thickness. The experimental design is relatively detailed and representative, and the workload is huge. This study reflects the research objectives with scientific and rational experiments and data. This study will analyze the differences of various indicators from the perspective of genetic to build a bridge between micro-structure and macro-structure for rational utilization of the whole area of Moso bamboo resources in China.


Introduction
Bamboo is an important non-woody forest resource rich in natural lignocellulose. The global annual output of bamboo is about 1.5 billion, and continuously regenerated bamboo is an important alternative resource for about 1/2 of the annual wood gap (250 million m 3 ) [1,2]. China is one of the distribution centers of bamboo plants, ranking first in the world in terms of bamboo species, storage, and area [3][4][5]. There are more than 800 bamboo species in China, which contains more than 1/3 of the bamboo species of the world. Moso bamboo (Phyllostachys edulis (Carrière) J. Houz.) is one of the core bamboo species in the development of the bamboo industry, which has the widest distribution area, the highest output value, and the best industrial utilization in China. Under the background of modern forestry advocating "developing bamboo in place of plastic, timber and steel substitute", it is particularly important to accelerate the utilization and research Materials were collected from the nature forest of Moso bamboo. A total of 70 samples were taken from each area as a population. To avoid sampling the same clone plant, we collected the individual plant in the same direction and from more than 50 m apart. The collected materials should be 3-or 4-years-old and comply with the Chinese national standard GB/T 2690-2000 "Mao Bamboo plant". One intact and undamaged internode of each plant from diameter at breast height was cut and bought back to the lab as the experimental material.

Determination of Phenotypic Traits
The diameter at breast height (DBH), node length at DBH, clear height, and Nodes under branch were measured before cutting using the DBH ruler. The average of nodes length under branch was calculated. When the internodes brought back to the lab were cut, the wall thickness at DBH in the 4 directions of bamboo rings was measured for averaging.
A total of 50 leaves were selected randomly from the upper, middle, and lower parts of the bamboo culm. The length of every leaf was measured from the tip of the leaf to the point where the leaf joins the petiole but not including the petiole, and the widest part of the leaf was measured as the leaf width. The average length and width of 50 leaves in every plant and their ratio of length to width were calculated [32,33].

Frequency and Size of Vascular Bundle
A sample piece of the central internode was cut, 10 mm × 10 mm × t mm (bamboo wall thickness), as test materials. The materials were softened for 2-3 h, cut into 30-40 microns thin, dry gas, and made a permanent microscope section. Under the system of microscopic imaging, the bamboo wall sections were divided into 3 uniform divisions along the radial direction of the bamboo wall using image analysis software. The 3 divisions were bamboo green (BG) in the outer zone, bamboo timber (BT) in the middle zone, and bamboo yellow (BY) in the inner zone, respectively. The Frequency of the vascular bundle was determined by checking the number of the vascular bundle in BT of bamboo, the size of the vascular bundle by the ratio of the length to width of the vascular bundle [34].

Length and Width of Fiber and Parenchyma
The fiber length and width of bamboo samples were determined by the segregation process [32,35]. The BT of samples were cut into about 2 cm length as a matchstick and put into the test tube. After the complete segregation of the sample, the segregation liquid was poured out, and then the samples were washed with distilled water to make temporary sections. The sections were amplified and measured by a digital projector 50 times. Each bamboo was randomly measured with 30 fibroblasts and 30 parenchyma cells for their length and width, as shown in Figure 1.

Wall Thickness and Lumen Diameter of Fiber Cell and Parenchyma
In this study, the wall thickness of fiber cells and parenchyma refers to the wall thickness of the single-cell [33]. The prepared permanent Section 2.3.1 was imaged 400 times based on the system of microscopic imaging, and the thickness and lumen diameter were measured by the wood analysis software. The 30 fibroblasts and 30 parenchyma cells Figure 1. Schematic diagram of measure method about the wall thickness and lumen diameter of fiber cell and parenchyma. In the figure of the permanent section, a vascular bundle in BT was selected. According to the direction of the yellow arrow, the fiber cell was selected in turn to obtain the double-wall thickness and lumen diameter, and according to the direction of the red arrow, the parenchyma was also selected.

Wall Thickness and Lumen Diameter of Fiber Cell and Parenchyma
In this study, the wall thickness of fiber cells and parenchyma refers to the wall thickness of the single-cell [33]. The prepared permanent Section 2.3.1 was imaged 400 times based on the system of microscopic imaging, and the thickness and lumen diameter were measured by the wood analysis software. The 30 fibroblasts and 30 parenchyma cells of each bamboo were measured randomly for their wall thickness and lumen diameter, as shown in Figure 1.

Crystallinity
After being air-dried, the culms of mature bamboo were fractionated manually with a knife into 3 parts: bamboo green (BG), bamboo timber (BT), and bamboo yellow (BY). The BT parts were ground, and the fraction passing 80-mesh but not 100-mesh was collected for determining crystallinity of cellulose by X-ray diffractometer [36]. The copper target was used for X-ray tube and nickel plate for eliminating CuKa radiation. The control scanning voltage was 40 kV, and the scanning current was 40 mA. The scanning speed was 0.02°/s, and the sample scanning range was from 5° to 45°. According to the intensity of the diffraction pattern, the Segal method was adopted to calculate the relative crystallinity. There was a maximum diffraction peak (I200) near the scanning curve (2θ = 22°), and a minimum peak (Iam) near the scanning curve (2θ = 18°), then the calculation formula of crystallinity is the following formula: where Cr is the percentage of relative crystallinity, I002 is the maximum strength of lattice diffraction Angle (002), namely the diffraction intensity in the crystal region. Iam is the scattering intensity of non-crystal background diffraction at 2θ = 18°.
100μm Figure 1. Schematic diagram of measure method about the wall thickness and lumen diameter of fiber cell and parenchyma. In the figure of the permanent section, a vascular bundle in BT was selected. According to the direction of the yellow arrow, the fiber cell was selected in turn to obtain the double-wall thickness and lumen diameter, and according to the direction of the red arrow, the parenchyma was also selected.

Crystallinity
After being air-dried, the culms of mature bamboo were fractionated manually with a knife into 3 parts: bamboo green (BG), bamboo timber (BT), and bamboo yellow (BY). The BT parts were ground, and the fraction passing 80-mesh but not 100-mesh was collected for determining crystallinity of cellulose by X-ray diffractometer [36]. The copper target was used for X-ray tube and nickel plate for eliminating CuKa radiation. The control scanning voltage was 40 kV, and the scanning current was 40 mA. The scanning speed was 0.02 • /s, and the sample scanning range was from 5 • to 45 • . According to the intensity of the diffraction pattern, the Segal method was adopted to calculate the relative crystallinity. There was a maximum diffraction peak (I 200 ) near the scanning curve (2θ = 22 • ), and a minimum peak (I am ) near the scanning curve (2θ = 18 • ), then the calculation formula of crystallinity is the following formula: where Cr is the percentage of relative crystallinity, I 002 is the maximum strength of lattice diffraction Angle (002), namely the diffraction intensity in the crystal region. I am is the scattering intensity of non-crystal background diffraction at 2θ = 18 • .

Physical Properties
The physical properties of bamboo were carried out according to GB/T 157805-1995GB/T 157805- (1995 and the method of [37]. The measurement traits included oven-dried volume shrinking ratio, air-dried volume shrinking ratio, moisture content, oven-dried density, and air-dried density. A total of 10 samples of each bamboo were measured.

Microfibril Angle
With a sharp single-sided blade, BY and BG of the bamboo block (about 10 mm × 20 mm × t mm) were carefully removed. After the fibers were exposed, the stringy thin slices with a thickness of about 1 mm were cut out as the measurement samples of microfibril Angle (MFA). The test equipment was an X-ray diffractometer. The X-ray source was CuK, the wavelength was 0.154 and 178 nm, and the transmission diffraction pattern was adopted. The control scanning voltage was 40 kV, and the scanning current was 40 mA. The diffraction image was processed with Origin software to realize the fast calculation of the microfilament angle ( Figure 2). Each bamboo was measured 5 times by 5 samples.

Physical Properties
The physical properties of bamboo were carried out according to GB/T 157805-1995GB/T 157805- (1995 and the method of [37]. The measurement traits included oven-dried volume shrinking ratio, air-dried volume shrinking ratio, moisture content, oven-dried density, and air-dried density. A total of 10 samples of each bamboo were measured.

Microfibril Angle
With a sharp single-sided blade, BY and BG of the bamboo block (about 10 mm × 20 mm × t mm) were carefully removed. After the fibers were exposed, the stringy thin slices with a thickness of about 1 mm were cut out as the measurement samples of microfibril Angle (MFA). The test equipment was an X-ray diffractometer. The X-ray source was CuK, the wavelength was 0.154 and 178 nm, and the transmission diffraction pattern was adopted. The control scanning voltage was 40 kV, and the scanning current was 40 mA. The diffraction image was processed with Origin software to realize the fast calculation of the microfilament angle ( Figure 2). Each bamboo was measured 5 times by 5 samples. The specific calculation method was as follows: Phi scanning intensity data obtained by X-ray was imported into Origin8.0, and the scanning curve was fitted by Gaussian function. MFA was calculated using 0.6T method, and its fitting function is: where σ1 and σ2 are half peak width in degrees. T is for σ1 + σ2 in degrees. The unit of MFA is degree [38].

Mechanical Properties
The mechanical properties of bamboo were carried out according to GB/T 157805-1995 (1995) and the method of Tran [39] and Zhang [40]. The measurement traits included bending strength (MOR) and Modulus of Elasticity (MOE). The sample size was 160 mm The specific calculation method was as follows: Phi scanning intensity data obtained by X-ray was imported into Origin8.0, and the scanning curve was fitted by Gaussian function. MFA was calculated using 0.6T method, and its fitting function is: where σ 1 and σ 2 are half peak width in degrees. T is for σ 1 + σ 2 in degrees. The unit of MFA is degree [38].

Mechanical Properties
The mechanical properties of bamboo were carried out according to GB/T 157805-1995 (1995) and the method of Tran [39] and Zhang [40]. The measurement traits included bending strength (MOR) and Modulus of Elasticity (MOE). The sample size was 160 mm (length) × 10 mm (width) × t mm (bamboo wall thickness). The span of the bending test was 120 mm, and the compression speed was 8 mm/min along the string direction of the bamboo wall. Each set of tests was repeated 5 times by 5 samples of each bamboo.

Data Verification
Data were analyzed by R statistical software (R 3.5.3). The normality was tested, and the fitting graph of the normal distribution was drawn by shapiro.test function in the stats package. The homogeneity of variance was tested by the bartlett.test. The regression equation was fitted using the lm. The AOV function was used for a one-way ANOVA. The glht function in the multcomp packet was used for multiple comparative analyses. Principal component analysis (PCA) was carried out by using prcomp in the stats package to obtain eigenvalue and contribution rate. The gravel diagram of PCA was drawn by fa.parallel function in the psych package.

Correlation Analysis
Genetic correlation refers to the correlation between two traits in the same population or between the same trait in different years or at different locations due to genetic factors. The phenotypic variance could be divided into genetic variance and environmental variance. In the same way, phenotypic covariance can also be divided into genetic covariance and environmental covariance. Therefore, the corresponding phenotypic correlation, genetic correlation, and environmental correlation can be calculated. In forestry, it is one of the important contents of genetic determination to develop the correlation among different traits, also an important basis for indirect selection of forest trees [41].
The phenotypic correlation coefficient, genetic correlation coefficient, and environmental correlation coefficient were analyzed using batch and pin commands in the ASReml statistical analysis software based on R language [41,42]. The three kinds of correlation coefficients were estimated according to the variance and covariance. The formula is as follows: Phenotypic correlation coefficient : Genetic correlation coefficient : Environmental correlation coefficient : In the above formula, Cov p (x, y), Cov g (x, y), and Cov e (x, y) represented phenotype covariance, genotype covariance, and environment covariance, respectively. σ 2 p , σ 2 g , and σ 2 e represented phenotype variance, genotype variance, and environmental variance, respectively. σ 2 px and σ 2 py represented the phenotypic variance of trait X and trait Y, σ 2 gx and σ 2 gy represented the genetic variance of trait X and trait Y, and σ 2 ex and σ 2 ey represented the environmental variance of trait X and trait Y.
Generally speaking, the absolute value of the correlation coefficient in 0~0.09, 0.1~0.3, and 0.5~1.0 mean no correlation, weak correlation, and strong correlation, respectively. However, a significant difference test, namely t-test, was necessary to test whether the 2 groups of data were significantly correlated. In this study, if the correlation coefficient of the 2 groups of data was greater than 0.3 and the difference was significant, both of the data was seen to be correlated. The significance level of correlation was analyzed by the corr.test function in psych package and the correlation chart among traits was drawn by corrgram function in corrgram package.

Data Analysis of Phenotypic Traits
In order to investigate the variation of phenotypic traits of Phyllostachys edulis distributed in different subtropical provenance regions, 209 samples of Moso bamboo distributed in Anhui, Guangxi, and Zhejiang province were obtained. Nine phenotypic traits were analyzed using a one-way ANOVA and principal component analysis based on R. These traits were diameter at breast height (DBH), node length at DBH (NL_DBH), Wall thickness at DBH (WT_DBH), Clear height (CH), Average of nodes length under branch (NL_CH), Nodes under branch (N_CH), Leaf length (LL), Leaf width (LW), and Leaf length/leaf width (LL.LW), respectively. The data has been published and showed that the phenotypic traits of Moso bamboo were significantly different from 3 populations [43]. The main research results were as follows, the peaks of three DBH traits (DBH, NL_DBH, WT_DBH) and NL_CH were all from Anhui Province. The maximum of CH, N_CH, and LL were all in Guangxi Province. Within 3 populations, the variable coefficient of CH and NL_CH were the largest. The peaks of the variable coefficient of the three leaf traits (LL, LW, and LL.LW) were all from Guangxi Province.

Correlation Analysis of Phenotypic Traits
Based on the phenotypic traits of Phyllostachys edulis from Anhui, Guangxi, and Zhejiang, the analysis of phenotypic correlation and genetic correlation were carried out on the genetic variation of nine phenotypic traits ( Figure 3). The results showed that there were different degrees of correlation among the traits, both in phenotype and heredity. The one above the diagonal line is the pie chart of the phenotypic correlation coefficient, the one below the diagonal line is the phenotypic correlation coefficient (up) and genetic correlation coefficient (down). *** stands for 0.1% level prominent, ** stands for 1% level prominent, * stands for 5% level prominent. LL: leaf length, LW: leaf width, LL.LW: leaf length/leaf width. The one above the diagonal line is the pie chart of the phenotypic correlation coefficient, the one below the diagonal line is the phenotypic correlation coefficient (up) and genetic correlation coefficient (down). *** stands for 0.1% level prominent, ** stands for 1% level prominent, * stands for 5% level prominent.

Analysis of Cell Structure Traits and Correlation
The clear height (CH) was highly positively correlated with NL_CH and N_CH. The phenotypic correlation (genetic correlation) was the highest among all traits and reached the extremely significant level of 0.1%, which was 0.618 (0.201) with NL_CH and 0.744 (0.967) with N_CH, respectively. It indicated that the larger CH is, the longer NL_CH is and the more N_CH is. NL_CH was positively correlated with the DBH, NL_DBH, and WT_DBH. The phenotypic correlation (genetic correlation) were 0.483 (0.519), 0.525 (0.994) and 0.477 (0.784), respectively. In the genetic correlation, the three pairs of indicators all reached a significant positive correlation, indicating that they were controlled by the genetic mechanism.
Whether phenotypic correlation or genetic correlation, the three DBH traits (DBH, NL_DBH, WT_DBH) were positively correlated and reached a very significant level. The phenotypic correlation (genetic correlation) was 0.606(0.793), 0.393(0.797) and 0.346(0.88). It showed that the larger DBH is, the longer NL_DBH is, and the thicker WT_DBH is. In the selection of bamboo material, the DBH and NL_DBH were used to predict the WT_DBH for rational classification and utilization of bamboo.
There was a very significant correlation between the three leaf traits (LL, LW, and LL.LW). The phenotypic correlation coefficients between LL and LW, LL, and LL.LW, LW, and LL.LW were 0.526, 0.443, and −0.502, respectively. As LL.LW was the ratio of leaf length to width, it was a general rule that the LL.LW was positively correlated with LL and negatively correlated with LW. The phenotypic correlation was highly positive between LL and LW. The genetic correlation between them showed a very high significance, which the correlation coefficient was 0.666, indicating that the correlation between leaf length and width was mainly controlled by genetic mechanism. There was a high significance and medium positive correlation between LW and NL_CH. It indicated that leaf width could further reflect the trait of NL_CH, and indirectly measure the three traits of DBH.

Correlation Analysis of Cell Structure Traits
Based on the cell structure data of Phyllostachys edulis from three areas, the analysis of phenotypic correlation and genetic correlation were conducted on the 11 cell structure traits ( Figure 4). The results showed that the overall correlation among the 11 traits was low, and only the correlation of the three groups reached more than 0.3. The one above the diagonal line is the pie chart of the phenotypic correlation coefficient, the one below the diagonal line is the phenotypic correlation coefficient (up) and genetic correlation coefficient (down). *** stands for 0.1% level prominent, ** stands for 1% level prominent, * stands for 5% level prominent. The correlation between L_PAREN and W_PAREN was the highest, and the phenotypic correlation (genetic correlation) was 0.747 (0.89) with a highly positive correlation and extremely significant level of 1%. It indicated that the elongation of parenchyma cells could promote its lateral broadening. There was not the same correlation between L_FIBER and W_FIBER of fiber cells with a phenotypic correlation of approximately 0 and a genetic correlation of 0.579, unlike the L_PAREN and W_PAREN with a highly positive correlation.
There was a negative correlation between WT_PAREN and LD_PAREN of parenchyma cells with the phenotypic correlation (genetic correlation) of −0.422 (−0.734). It indicated that LD_PAREN decreased with the increase of WT_PAREN. There was a positive correlation between WT_PAREN and LD_FIBER with the phenotypic correlation (genetic correlation) of 0.36 (0.564). It indicated that LD_ FIBER increased with the increase of WT_PAREN.

Correlation Analysis of Physical Mechanic Properties
According to the physical mechanic properties of Phyllostachys edulis from the three areas, the analysis of phenotypic correlation and genetic correlation was conducted on 8 properties ( Figure 5). The results showed that the overall correlation among the 8 properties was high, and the correlation of several groups reached 0.3 or even more than 0.5.
MOR was highly positively correlated with MOE. The phenotypic correlation (genetic correlation) was the highest among the 8 traits and reached the extremely significant level of 0.1%, which was 0.974 (0.999). O_DENSITY and A_DENSITY were highly positively correlated with the phenotypic correlation of 0.89 and genetic correlation of 0.929. Both of the densities were highly positively correlated with MOE, respectively. The variation coefficients were above 0.5, which reached an extremely significant level. However, the phenotypic correlation between density traits and MOR did not reach a significant level of 5%, which was not statistically significant. The results showed that only MOE of Phyllostachys edulis would be grown with the increase of density, and MOR needs to be further studied. The results of this study that the density was highly positively correlated with MOE, consistent with general research [46].
There was not only a highly and significantly positive correlation between O_SHRINK and that of A_SHRINK with the phenotypic correlation (genetic correlation) of 0.679 (0.71) but also a medium positive correlation between O_SHRINK and MC or O_SHRINK and MFA with 0.377 (0.314) or 0.4 (0.877), respectively. The results showed that the shrinking ratio of volume was closely related to the moisture content of Moso bamboo, and the change affected its MFA. The main reason for bamboo shrinkage is that the vessels in the vascular bundle lose water and shrink [47], which also explains the close relationship between the shrinkage ratio and water content.
The moisture content of bamboo is an important parameter to affect bamboo properties. Their properties of bending and compression differ significantly with the change of water content. In this study, although MC was moderately negatively correlated with MOE and MOR, the phenotypic correlation did not reach the significant level of 5%, and there was no statistical significance. The one above the diagonal line is the pie chart of the phenotypic correlation coefficient, the one below the diagonal line is the phenotypic correlation coefficient (up) and genetic correlation coefficient (down). *** stands for 0.1% level prominent, ** stands for 1% level prominent, * stands for 5% level prominent.

Correlation Analysis between Phenotypic Traits and Cell Structure Properties
The correlation was conducted on 9 phenotypic traits and 11 cell structure properties of Moso bamboo from Anhui, Guangxi, and Zhejiang. In the total 99 groups of traits (9 × 11 = 99), the absolute value of the phenotypic correlation of 43 groups was less than 0.1, Figure 5. Correlation analysis of 8 material properties in Moso bamboo (Phyllostachys edulis). Notes: O_SHRINK: oven-dried volume shrinking ratio, A_SHRINK: air-dried volume shrinking ratio, MC: moisture content, O_DENSITY: oven-dried density, A_DENSITY: air-dried density, MFA: microfibril angle, MOR: bending strength, MOE: modulus of elasticity. The one above the diagonal line is the pie chart of the phenotypic correlation coefficient, the one below the diagonal line is the phenotypic correlation coefficient (up) and genetic correlation coefficient (down). *** stands for 0.1% level prominent, ** stands for 1% level prominent, * stands for 5% level prominent.

Correlation Analysis between Phenotypic Traits and Cell Structure Properties
The correlation was conducted on 9 phenotypic traits and 11 cell structure properties of Moso bamboo from Anhui, Guangxi, and Zhejiang. In the total 99 groups of traits (9 × 11 = 99), the absolute value of the phenotypic correlation of 43 groups was less than 0.1, which could be regarded as no correlation, accounting for 43.43% of the total groups (Table 2). That indicated that there was not much correlation between phenotypic traits and cell structure properties on the whole. There were four groups of traits to obtain medium and high correlation that showed a close correlation. The phenotypic correlation (genetic correlation) between W_FIBER and WT_DBH was 0.592 (0.953), which was very significant at the 1% level, indicating that the thickness of the bamboo wall at DBH might increase the width of the bamboo fiber at the same height. On the other side, the phenotypic correlation between W_FIBER and NL_CH was 0.377, and the genetic correlation was 0.854, which showed that the larger NL_CH was, the wider bamboo fiber was.
There was a moderate positive correlation between LL and LD_FIBER or LL and LD_PAREN. The phenotypic correlation (genetic correlation) was 0.361 (0.567) and 0.361 (0.687), respectively. They all reached a very significant level of 1%. It indicated that LD_FIBER and LD_PAREN would increase with the longer leaf length. Curiously, there was no significant relationship between LD_FIBER and LD_PAREN. Their phenotypic coefficient was 0.074, showing that LL is significantly related to LD_FIBER or LD_PAREN.

Correlation Analysis between Phenotypic Traits and Material Properties
The correlation was analyzed on 9 phenotypic traits and 8 material properties of Moso bamboo from Anhui, Guangxi, and Zhejiang. In the whole 72 groups of traits (9 × 8 = 72), there were relatively rich groups in the traits to achieve a medium-high correlation ( Table 3). The absolute value of the phenotypic correlation of 20 groups was less than 0.1, which could be regarded as no correlation, accounting for 27.78% of the total groups.  The analysis shows that MC, MOR, and MOE were highly correlated with DBH, CH, and LL, respectively, indicating that these six traits are closely related. Among them, MOR and MOE were negatively correlated with DBH, respectively. Moreover, the phenotypic correlation (genetic correlation) was −0.484 (−0.879) and −0.466 (−0.857), respectively. At the same time, MOR and MOE were significantly positively correlated with CH and LL, indicating that the smaller DBH, the larger CH, or the longer leaf might all increase MOR and MOE.
The MFA of bamboo showed a moderate negative correlation with NL_DBH and NL_CH with phenotypic correlation (genetic correlation) of −0.324 (−0.963) and −0.316 (−0.831), respectively. That showed a close relationship between MFA and the node length of bamboo. Moreover, because NL_CH and NL_DBH are highly positively correlated (3.1 above), the increase of any of the two traits may decrease the bamboo MFA.
A_SHRINK was significantly negatively correlated with LW, and the phenotypic correlation (genetic correlation) was −0.323 (−0.714), indicating that the wider the leaf was, the smaller the A_SHRINK was. The absolute value of correlation between LL and A_SHRINK was less than 0.1, which could be seen as no correlation. However, there was a significant and positive correlation between A_density and LL with the phenotypic correlation (genetic correlation) of 0.334 (0.973), indicating that the longer the leaf was, the higher the A_density was. Therefore, the A_DENSITY and A_SHRINK of bamboo can be measured indirectly by measuring leaf length and leaf width during field investigation and sampling.

Correlation Analysis between Cell Structure and Material Properties
The correlation was conducted on 11 cell structure properties and 8 material properties of Moso bamboo from Anhui, Guangxi, and Zhejiang. In the total 88 groups of traits (11 × 8 = 88), the absolute value of the phenotypic correlation of 30 groups was less than 0.1, which could be regarded as no correlation, accounting for 34.09% of the total groups (Table 4). Some physical-mechanical properties, such as SHRINK, DENSITY, MOR, and MOE) were closely correlated with wall thickness and lumen diameter, while MC and MFA were not correlated with all cell structure properties. The increase of wall thickness of fiber is the result of the thickening of the secondary wall, and the secondary wall has a multi-layer structure with alternating width. The content of lignin in the thin layer is higher than that in the thick layer, and the cell orientation in each layer is different (Liese, 1998;1996;Alvin, 1988;Murphy, 1992;Murphy, 1997;Walter, 2001;Parameswaran, 1976). In the study of Abd (1990), the wall thickness of fiber has a positive correlation with the modulus of elasticity and a negative correlation with bending strength. The lumen diameter of fiber has a positive correlation with all mechanical properties except compressive strength. In this study, the phenotypic correlation (genetic correlation) between WT_FIBER and MOE was 0.44 (0.614) to show a significant positive correlation, which was consistent with the results of Abd. However, there was still a significant positive correlation between WT_FIBER and MOR with the phenotypic correlation (genetic correlation) of 0.423 (0.597), which was contrary to the results of ABD's study. LD_FIBER has a significant positive correlation with A_DENSITY, MOR, and MOE, with the phenotypic correlation (genetic correlation) of 0.347 (0.724), 0.3 (0.559), and 0.474 (0.546), respectively. It indicated that the thickening of LD_FIBER can significantly enhance MOR and increase the air-dry density.
O_DENSITY and WT_PAREN were significantly positively correlated with the phenotypic correlation (genetic correlation) of 0.395 (0.651), and A_DENSITY and WT_PAREN were significantly positively correlated with 0.367 (0.654). That indicated that the thickening of the cell wall of parenchyma could significantly increase the density of bamboo. O_SHRINK and A_SHRINK was significantly positively correlated to LD_PAREN with the phenotypic correlation (genetic correlation) of 0.484 (0.901) and 0.417 (0.628), respectively, which showed that increasing the lumen diameter of parenchyma may increase the shrinking ratio of bamboo.

Discussion
The inheritance of traits is often controlled by multiple genes, and the interaction between genes makes the traits to be interrelated and interrelated [48]. Only by understanding the degree of correlation between traits can we balance the selection of traits in breeding improvement to improve the efficiency of bamboo improvement. As a basic genetic parameter, genetic correlation plays an important role in analyzing the correlation between important traits, estimating the progress of genetic correlation and evaluating clones, etc. [49]. At the same time, the mutual relationship among various traits affects the consequences of synchronous improvement among various traits.
From the correlation analysis of phenotypic traits, the phenotypic correlation and genetic correlation of the three DBH traits (DBH, NL_DBH, and WT_DBH) were positively correlated with each other and reached a very significant level ( Figure 6). It showed that the larger DBH is, the longer NL_DBH is, and the thicker WT_DBH is [50]. In the selection of bamboo material, the DBH and NL_DBH were used to predict the WT_DBH for rational classification and utilization of bamboo. Therefore, the DBH traits could be indirectly evaluated by NL_CH, during field investigation and sampling. Furthermore, there was a high significance and medium positive correlation between LW and NL_CH ( Figure 6). It indicated that leaf width could further reflect the trait of NL_CH, and indirectly measure the three DBH traits.  From the correlation analysis of cell structure traits, it was interesting that there was no the same correlation between L_FIBER and W_FIBER of fiber cell with a phenotypic correlation of approximately 0 and a genetic correlation of 0.579, unlike the L_PAREN and W_PAREN with highly positive correlation [51]. It showed that although the length and width of fiber cells were controlled by the genetic mechanism and the genes of L_FIBER may change the expression of W_FIBER genes, but there was no obvious phenotypic correlation to be formed.
Relevant studies have shown that bamboo density is an important index to measure the properties of bamboo. In general, the higher the bamboo density is, the higher its mechanical properties will be [52], which is consistent with the results of this study that the density is highly positively correlated with MOE. The shrinking ratio of volume was closely related to the moisture content of Moso bamboo, and the change affected its MFA. The results showed that the shrinking ratio of volume was closely related to the moisture content of Moso bamboo and changed affected its MFA. The main reason for bamboo shrinkage is that the vessels in the vascular bundle lose the water and shrink [47,53], which also explains the close relationship between the shrinkage ratio and water content. The MFA of wood can directly affect the mechanical properties and dimensional stability of wood, which is positively correlated with wood density, longitudinal shrinkage ratio, and elastic modulus, and negatively correlated with the tensile strength, tear strength, and hardness of wood [54][55][56][57]. However, in the study of previous studies [58][59][60] and this study on bamboo, the MFA did not show an obvious change pattern similar to that of wood.
From the correlation analysis of 9 phenotypic traits and 11 cell structure properties, W_FIBER was highly positively correlated with WT_DBH and moderately positively correlated with NL_CH. In combination with the previous results about a moderate positive correlation between NL_CH and WT_DBH, W_FIBER, NL_CH, and WT_DBH may be controlled by linked genes or by the same pleiotropic genes. Therefore, WT_DBH and W_FIBER can be indirectly evaluated by NL_CH during field sampling. On the other hand, it can be seen that all of DBH, NL_DBH, CH, NL_CH, LW, and LL.LW has no significant and moderate or above correlation with any cell structure properties. This may be because the correlation coefficient is affected by the sample size. That is, the larger the sample size, the smaller the correlation coefficient required to achieve a significant correlation. In this study, the sample size was more than 200 plants, which might make the correlation coefficient slightly lower and make the correlation of some groups with a moderate correlation greater than 0.3. However, it is also based on the large sample size to reach the statistical significance, the reference of the results will be greatly improved, thus that the significantly related traits can be better screened out.
From the correlation analysis of 9 phenotypic traits and 8 material properties, MC, MOR and MOE are highly correlated with DBH, CH, and LL, respectively. In combination with the previous results about the significant negative correlation between MOR and MC or MOE and MC can we infer that MC may be positively correlated with DBH and negatively correlated with CH and LL? This conclusion was confirmed in Table 2, and the correlations were all at a significant level. That is, the phenotypic correlation (genetic correlation) between MC and DBH was 0.402 (0.787), that between MC and CH was −0.405 (−0.932), and that between MC and LL was −0.429 (−0.973), respectively. In addition, NL_CH was significantly positively correlated with MOR and MOE, and the phenotypic correlation (genetic correlation) was 0.313 (0.933) and 0.309 (0.945).
The tensile strength, modulus of elasticity, hardness, density, and dimensional stability of the fiber bundle generally increase with the increase of crystallinity [61]. In this study, the crystallinity was significantly positively correlated with MOE with the phenotypic correlation (genetic correlation) of 0.305 (0.992), which was consistent with previous studies. However, the density did not show a significant correlation with the crystallinity in this study.
Depending on the breeding goal, breeders tend to pay attention to the performance of a single character, and the evaluation and analysis of germplasm resources are onesided. On the other hand, the determination and sampling of cell structure and material properties of bamboo can cause irreversible damage to the plant, and the measurement process is complex. If the relationship between those traits has been known, cell structure and material properties can be indirectly predicted and assessed using easily measurable phenotypic traits. This will greatly improve the utilization of bamboo and reduce the cost of bamboo material improvement.