Taylor’s Power Law for Leaf Bilateral Symmetry

Leaf shape and symmetry is of interest because of the importance of leaves in photosynthesis. Recently, a novel method was proposed to measure the extent of bilateral symmetry in leaves in which a leaf was divided into left and right sides by a straight line through the leaf apex and base, and a number of equidistant strips were drawn perpendicular to the straight line to generate an equivalent number of differences in area between the left and right parts. These areal differences are the basis for a measure of leaf bilateral symmetry, which was then examined to see how well it follows Taylor’s power law (TPL) using three classes of plants, namely, 10 geographical populations of Parrotia subaequalis (H.T. Chang) R.M. Hao et H.T. Wei, 10 species of Bambusoideae, and 10 species of Rosaceae. The measure of bilateral symmetry followed TPL for a single species or for a class of closely related species. The estimate of the exponent of TPL for bamboo plants was significantly larger than for the dicotyledonous trees, but its goodness of fit was the best among the three classes of plants. The heterogeneity of light falling on branches and leaves due to above-ground architectural patterns is an important contributor to leaf asymmetry.


Introduction
The leaf is the most important organ of plant photosynthesis.The scaling relationship between leaf weight and leaf area (W = aA b , where W denotes leaf weight, A denotes leaf area, a and b are constants to be estimated) is of wide interest because the leaf dry mass per unit area (LMA = W/A = aA b−1 ) has been shown to be proportional to photosynthetic efficiency [1][2][3][4].Leaf shape can affect the scaling relationship between leaf weight and leaf area [5], so that closely related species with a large variation in leaf shape, for example, Rosa, will have a very weak scaling relationship between leaf weight and leaf area for pooled data from those species.However, if the closely related species are also similar in leaf shape, for example, bamboo plants, there can be a strong scaling relationship both for a single species and for pooled leaf weight and leaf area data from different species [6].Since LMA can affect the photosynthetic rate, leaf shape can also affect it [5].Thus, leaf shape appears to be important for the investigation of leaf functional traits.Seed plants are considered to have evolved independently from ancient large-leaf vascular plants with diverse leaf shapes [5], and the leaves of many plants, especially the angiosperms, exhibit bilateral symmetry [7].Previous studies on leaf shape paid much attention to leaf teeth and leaf anatomic structure, but largely neglected leaf bilateral symmetry [5,8].Although bilateral symmetry has attracted attention in biotic and abiotic areas, there are still few quantitative studies about leaf symmetry [9][10][11][12].
Forests 2018, 9, 500 2 of 11 Taylor's power law (TPL) first described the variance-mean relationship in ecology [13]: V = αM β , where V and M represent the variance and mean, respectively, of a series of specific ecological measures, for example, the population densities of insects in different-size quadrats, and α and β are constants to be estimated.This power law has been found to be valid in many areas of study [14].Controversy surrounding recent relevant studies is about whether the exponent β of TPL is determined by an inherent biological mechanism or is only a purely statistical feature [15][16][17][18].Cohen and Xu demonstrated that the estimate of the exponent of TPL is proportional to the skewness of the distribution generated by random sampling in blocks [15].Using a feasible set approach, Xiao et al. showed that TPL arises from the constraining influences of two primary variables, namely, the number of individuals and the number of sites [16].TPL with an exponent ranging from 1.0 to 2.0 could be generated by the most possible configurations of individuals among sites.Shi et al. [17] used two clustering point process models to examine the effect of dispersal distance of plant seeds or rhizomes on the estimate of the exponent of TPL, and found that with increasing dispersal distance the exponent described a complex response curve.In addition, dry weights of crops at different investigation times and the development and growth rates of insects at different temperatures have also been shown to follow TPL [18].These measures associated with TPL appear to be related to the energy distribution among different statistical units.To use a more direct approach to demonstrate the link between TPL and energy distribution, Li et al. [19] checked the influence of different time interval divisions of the same seismic time series on the estimate of the exponent of TPL between the variance and mean of the released energies among different earthquakes within the same interval.The estimates of the exponent of TPL for the seismic time series approximated 2.0.Then, fixing the exponent of TPL to 2.0, the estimate of the logarithm of the constant α in TPL was demonstrated to be a logarithmic function of the time interval division.These studies related to TPL have added to our understanding of the variation of a large number of ecological and non-ecological measures, for example, animal and plant population density, biomass, poikilothermic developmental rate, crime, precipitation, released energies of aftershocks, trading activity of stock, and so on [13,[18][19][20][21][22].In our recent studies, weight is a typical representative of energy for biological measures, so the variance-mean relationship based on weight measures should reflect TPL well [23].
In this paper, we examine whether TPL holds for a measure of leaf bilateral symmetry.Because a relationship between leaf weight and leaf area has been demonstrated to hold for many individual species [1,2,6], we believe that a measure of bilateral symmetry based on leaf area instead of leaf weight could fit the TPL.A strong scaling relationship between leaf weight and leaf area should result in a strong TPL for a measure of bilateral symmetry based on leaf area.In addition, the light heterogeneity caused by differences of the above-ground architectural patterns among different families [24] should result in different estimates of the exponent of TPL.The scaling relationship between leaf weight and leaf area actually only affects the goodness of fit for TPL of the leaf bilateral symmetry measure.However, it might not determine the size of the estimate of the exponent of TPL, which could be closely related to the light heterogeneity.The topology of branches apparently has an important influence on the light interception for the leaves at different layers of the crown.A regular and sparse branching pattern can render leaves to capture light with two sides equally that consequently reduces the degree of leaf bilateral asymmetry.In this case, the estimate of the exponent of TPL would be expected to be smaller than that from an irregular and dense branching pattern.

Materials
Three families of plants were used: 10 geographical populations of Parrotia subaequalis (H.T. Chang) R.M. Hao & H.T. Wei (Hamamelidaceae) with 150 leaves for each population, 10 species of Bambusoideae (Poaceae) with 100 to 500 leaves for every species, and 10 species of Rosaceae with more than 300 leaves for each species.Tables S1-S3 on the online Supplementary Materials show the Forests 2018, 9, 500 3 of 11 details of these plants, including their collection sites, Latin names, sample sizes.For P. subaequalis, the populations were collected in South-Eastern China.The plants of Bambusoideae and Rosaceae were all collected from the Nanjing Forestry University campus in Nanjing, Jiangsu Province, China (see Tables S1-S3 for details).To carry out the interspecific comparison, the sample size (number of leaves collected) of each species is more than 100.To test whether the power-law relationship between the variance and mean of a measure of leaf bilateral symmetry for a single species can be representative of whole bamboo clusters, we collected 500 leaves of Pleioblastus chino (Franchet & Savatier) Makino.We collected these leaves when their shapes were basically or completely stable in order to reduce the measurement error during the calculation of the extent of leaf bilateral symmetry.

Methods
Shi et al. proposed a method to calculate the extent of leaf bilateral symmetry (or asymmetry) [7].A leaf was divided into left and right (or upper and lower) sides by a straight line passing through leaf apex and leaf base (Figure 1).Using n strips perpendicular to the straight line to intersect the leaf, the areas of the intersected parts between the leaf and strips were obtained.Let A 1i (i ranging from 1 to n) represent the areas of the intersected parts between the left (upper) side of the leaf with the i-th strip, and A 2i (i ranging from 1 to n) represent the areas of those parts between the right (lower) side of the leaf with the i-th strip.Then the following differences in a leaf were obtained: Savatier) Makino.We collected these leaves when their shapes were basically or completely stable in order to reduce the measurement error during the calculation of the extent of leaf bilateral symmetry.

Methods
Shi et al. proposed a method to calculate the extent of leaf bilateral symmetry (or asymmetry) [7].A leaf was divided into left and right (or upper and lower) sides by a straight line passing through leaf apex and leaf base (Figure 1).Using n strips perpendicular to the straight line to intersect the leaf, the areas of the intersected parts between the leaf and strips were obtained.Let A1i (i ranging from 1 to n) represent the areas of the intersected parts between the left (upper) side of the leaf with the i-th strip, and A2i (i ranging from 1 to n) represent the areas of those parts between the right (lower) side of the leaf with the i-th strip.Then the following differences in a leaf were obtained: Here,   is used to quantify the extent of the subregion on the leaf segmented by the i-th strip deviating from a perfect bilateral symmetry.Clearly,   ≥ 0. In actuality, because of the existence of measurement errors and left-right natural asymmetry,   is usually greater than zero.From these measurements on a leaf, paired values of the mean (M) and variance (V) of these differences were obtained.Assuming that there are N leaves for a species or a cluster of species, there would be N pairs of data for the mean and variance: (  ,   ), where j represents the j-th leaf and j ranges from 1 to N. According to the definition of TPL, Taking the natural logarithm of both sides of the above equation, we have where  =  is the response variable,  =  is the independent variable, and  = α and β are parameters to be estimated.Here, D i is used to quantify the extent of the subregion on the leaf segmented by the i-th strip deviating from a perfect bilateral symmetry.Clearly, D i ≥ 0. In actuality, because of the existence of measurement errors and left-right natural asymmetry, D i is usually greater than zero.From these measurements on a leaf, paired values of the mean (M) and variance (V) of these differences were obtained.Assuming that there are N leaves for a species or a cluster of species, there would be N pairs of data for the mean and variance: M j , V j , where j represents the j-th leaf and j ranges from 1 to N. According to the definition of TPL, Forests 2018, 9, 500 4 of 11 Taking the natural logarithm of both sides of the above equation, we have where y = lnV is the response variable, x = lnM is the independent variable, and γ = lnα and β are parameters to be estimated.Assume that N leaves for a species have been collected (see Tables S1-S3 for the detailed sample size information for every species).For any leaf, we can calculate the mean and variance of the absolute values of the areal differences between the left and right subregions segmented by the strips on the leaf (Figure 1).Thus, for each species, there are in total N pairs of data for mean and variance.Linear regression was then used to fit the log-transformed data of mean and variance to estimate the parameters β and γ of Equation ( 3).The bootstrap percentile method was used to test whether there was a significant difference between any of the estimates of the exponent of TPL for any two species within the same family, and the 95% confidence intervals (CIs) of the exponent of TPL were also calculated using the bootstrap percentile to indicate the reliability of the parameter estimates [25,26].The coefficient of determination (i.e., R 2 ) was calculated to measure the goodness of fit of the linear regression and the degree of the predicted values deviating from the observations.
where ŷj denotes the predicted value of y j , and y denotes the mean of y j (j ranging from 1 to N).For the comparison of the estimates of the exponent of TPL among three families (which means that the data for the mean and variance of species within the same family were pooled), the analysis of covariance was used.R statistical software was used to carry out all calculations and to draw figures [27].The extraction of the planar coordinates of the leaf edge and the measure of the extent of bilateral symmetry used the procedures in references [7,28,29].

Results and Discussion
For P. subaequalis, TPL holds for the 10 geographical populations, with the coefficients of determination ranging from 0.81 to 0.91 (Figure 2).The estimates of the exponent of TPL for the measure of leaf bilateral symmetry for the 10 populations range from 1.5 to 1.7 except for population 'TT' (Code 9) which had an estimate of 1.79, significantly larger than those of the other populations except for population 'TX' (Code 10) (Figure 3a).The estimate of the exponent of TPL of population 'HZ' (Code 4) is significantly smaller than those of populations 'HN', 'TT' and 'TX' (Codes 2, 9 and 10).There are eight populations in total whose estimates of the exponent of TPL are not significantly different (see the letter 'B' in Figure 3a).The estimate of the exponent of TPL for the pooled data for these 10 populations is 1.649 with 95% CI (1.618, 1.681).
For the plants of Bambusoideae, there are large differences in the estimates of the exponent of TPL between the pooled data and the individual species (Figure 3b).The estimate based on the pooled data is 1.925 with 95% CI (1.904, 1.946) (Figure 4), apparently larger than the estimates of the exponent of TPL, which for most species range from 1.4 to 1.9, except for species 2, 4, and 5 whose estimates are significantly larger than the estimates of the remaining seven species (Figures 3b and 4).These three species with large estimates of the exponent of TPL belong to the same genus Indocalamus.We note that these estimates of the exponent of TPL are the largest among the three classes of plant (Figure 3).
For the plants of Rosaceae, the estimates of the exponent of TPL range from 1.3 to 1.7, with the estimate for the pooled data being 1.621 with 95% CI (1.602, 1.641) (Figure 5).Although there are significant differences in the estimates of the exponent of TPL among different species, the range of the estimates for this family is smaller than that for the bamboo plants (Figure 3).
From the comparison of the exponents of TPL among three pooled data sets using analysis of covariance, we found that they were significantly different, although the difference between the estimates of the exponent for the pooled data of P. subaequalis and for the pooled data of Rosaceae was very small (Table 1 and Figure 6).However, the estimate of the exponent for the pooled data of Bambusoideae was significantly higher than for the others.After pooling the data of ln(variance) and ln(mean) for all three classes of plants, TPL still holds and the estimate of the exponent of TPL is 1.773 with 95% CI (1.760, 1.787).It appears that the measure of bilateral symmetry for all leaves of the plants investigated follow a general TPL (R 2 = 0.907) if we neglect differences among the species and among the classes (Figure 6).For the plants of Bambusoideae, there are large differences in the estimates of the exponent of TPL between the pooled data and the individual species (Figure 3b).The estimate based on the pooled data is 1.925 with 95% CI (1.904, 1.946) (Figure 4), apparently larger than the estimates of the exponent of TPL, which for most species range from 1.4 to 1.9, except for species 2, 4, and 5 whose estimates are significantly larger than the estimates of the remaining seven species (Figures 3b and Figure 4).These three species with large estimates of the exponent of TPL belong to the same genus Indocalamus.We note that these estimates of the exponent of TPL are the largest among the three classes of plant (Figure 3).S2.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.
For the plants of Rosaceae, the estimates of the exponent of TPL range from 1.3 to 1.7, with the estimate for the pooled data being 1.621 with 95% CI (1.602, 1.641) (Figure 5).Although there are significant differences in the estimates of the exponent of TPL among different species, the range of the estimates for this family is smaller than that for the bamboo plants (Figure 3).S2.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.The detailed Latin name and collecting information can be found in Table S3.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.
From the comparison of the exponents of TPL among three pooled data sets using analysis of covariance, we found that they were significantly different, although the difference between the estimates of the exponent for the pooled data of P. subaequalis and for the pooled data of Rosaceae was very small (Table 1 and Figure 6).However, the estimate of the exponent for the pooled data of Bambusoideae was significantly higher than for the others.After pooling the data of ln(variance) and ln(mean) for all three classes of plants, TPL still holds and the estimate of the exponent of TPL is 1.773 with 95% CI (1.760, 1.787).It appears that the measure of bilateral symmetry for all leaves of the plants investigated follow a general TPL (R 2 = 0.907) if we neglect differences among the species and among The detailed Latin name and collecting information can be found in Table S3.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.Our previous study demonstrated that there was a strong scaling relationship between leaf weight and leaf area using 11 bamboo species, with the estimate for the scaling exponent being 1.147, with 95% CI (1.143, 1.152) and R 2 = 0.99 (see Figure 5 in reference [6]).However, our recent studies show that the estimate of the scaling exponent between leaf dry weight and leaf area using the pooled data of the 10 populations of P. subaequalis is 1.245, with 95% CI (1.214, 1.277) and R 2 = 0.70 [30], and the estimate for the scaling exponent between leaf dry weight and leaf area using the pooled data of the 10 species of Rosaceae is 1.462, with 95% CI (1.439, 1.486) and R 2 = 0.70 (unpublished data of P. Shi).Therefore, the scaling relationship of leaf weight and leaf area of bamboo plants is significantly stronger than those of the other two classes of plants based on the comparison of the coefficients of determination (i.e., R 2 ).Because the biomass (weight) is a better measurement than the area or length in describing TPL [18,23], plants with a better scaling relationship between leaf weight and leaf area will have a better TPL for describing leaf bilateral symmetry.If we could cut a leaf into different parts with equidistant intervals from the leaf base to leaf apex and measure their weights (see Figure 1), we believe that we would obtain a better goodness of fit for TPL by using the measures of leaf weight than using the measures of leaf area.However, the cost of measurement makes this impracticable.
The estimate of the exponent of TPL for bamboo plants is higher than those for P. subaequalis and Rosaceae, and the estimate for P. subaequalis is slightly higher than that for Rosaceae (Table 1).Leaf asymmetry is closely related to light heterogeneity caused by the crown structure [8,24,31].The above-ground architecture of plants could significantly affect phyllotaxy and leaf shape [24], and neighbourhood competition among plants can also have an influence on the spatial arrangement of branches and twigs, further affecting phyllotaxy and leaf shape [32,33].P. subaequalis and Rosaceae have similar crown patterns and leaf shapes, so their estimates of the exponent of TPL are very similar (Figure 6).However, for bamboo plants, there is a large variation in crown patterns and spatial arrangement of branches, and the branching system is more complicated than for other grasses [34].For example, a branch of Phyllostachys has two divarications, whereas a branch of Pleioblastus has three divarications.The extent of variation of light heterogeneity is rather large for some bamboo species, for example, plants of Indocalamus.In a natural population, approximately half the leaves face the sun, and another half are shaded, leading to large differences in light interception.The spatial densities of the bamboos have a significant influence on leaf size [35].For a leaf, if one side can receive light but another side is shaded, the former will be invested by photosynthesis with more biomass [8] and the leaf will become more bilaterally asymmetric.Thus, plants with a more regular spatial arrangement of branches for light capture might have a smaller variation in their leaf bilateral symmetry.Apparently, bamboo plants do not have such a regular pattern.In the habitats of bamboos Our previous study demonstrated that there was a strong scaling relationship between leaf weight and leaf area using 11 bamboo species, with the estimate for the scaling exponent being 1.147, with 95% CI (1.143, 1.152) and R 2 = 0.99 (see Figure 5 in reference [6]).However, our recent studies show that the estimate of the scaling exponent between leaf dry weight and leaf area using the pooled data of the 10 populations of P. subaequalis is 1.245, with 95% CI (1.214, 1.277) and R 2 = 0.70 [30], and the estimate for the scaling exponent between leaf dry weight and leaf area using the pooled data of the 10 species of Rosaceae is 1.462, with 95% CI (1.439, 1.486) and R 2 = 0.70 (unpublished data of P. Shi).Therefore, the scaling relationship of leaf weight and leaf area of bamboo plants is significantly stronger than those of the other two classes of plants based on the comparison of the coefficients of determination (i.e., R 2 ).Because the biomass (weight) is a better measurement than the area or length in describing TPL [18,23], plants with a better scaling relationship between leaf weight and leaf area will have a better TPL for describing leaf bilateral symmetry.If we could cut a leaf into different parts with equidistant intervals from the leaf base to leaf apex and measure their weights (see Figure 1), we believe that we would obtain a better goodness of fit for TPL by using the measures of leaf weight than using the measures of leaf area.However, the cost of measurement makes this impracticable.
The estimate of the exponent of TPL for bamboo plants is higher than those for P. subaequalis and Rosaceae, and the estimate for P. subaequalis is slightly higher than that for Rosaceae (Table 1).Leaf asymmetry is closely related to light heterogeneity caused by the crown structure [8,24,31].The above-ground architecture of plants could significantly affect phyllotaxy and leaf shape [24], and neighbourhood competition among plants can also have an influence on the spatial arrangement of branches and twigs, further affecting phyllotaxy and leaf shape [32,33].P. subaequalis and Rosaceae have similar crown patterns and leaf shapes, so their estimates of the exponent of TPL are very similar (Figure 6).However, for bamboo plants, there is a large variation in crown patterns and spatial arrangement of branches, and the branching system is more complicated than for other grasses [34].For example, a branch of Phyllostachys has two divarications, whereas a branch of Pleioblastus has three divarications.The extent of variation of light heterogeneity is rather large for some bamboo species, for example, plants of Indocalamus.In a natural population, approximately half the leaves face the sun, and another half are shaded, leading to large differences in light interception.The spatial densities of the bamboos have a significant influence on leaf size [35].For a leaf, if one side can receive light but another side is shaded, the former will be invested by photosynthesis with more biomass [8] and the leaf will become more bilaterally asymmetric.Thus, plants with a more regular Forests 2018, 9, 500 9 of 11 spatial arrangement of branches for light capture might have a smaller variation in their leaf bilateral symmetry.Apparently, bamboo plants do not have such a regular pattern.In the habitats of bamboos without human management, they grow fast and reach their maximal environmental capacity in a short period.For many herbaceous bamboos, this period does not exceed 10 years and the average biomass of bamboos decreases with increasing density [36,37].The leaves usually cannot extend to their natural state in crowed environments, and the possibility of bilateral asymmetry increases.Comparing the exponent of TPL for the measure of leaf bilateral symmetry for bamboo populations with different densities merits investigation.It is also worthwhile checking whether there are significant differences in the estimates of the exponent of TPL among woody plants with different densities that have the same or similar crowns.

Conclusions
Taylor's power law describing the variance-mean scaling relationship has been demonstrated to hold for a measure of leaf bilateral symmetry in three classes of plants investigated, with estimates of the exponent of TPL in the range of 1.0 to 2.0, in accordance with previous reports in the relevant ecological literature.TPL was found to be applicable to any single species or to the pooled data of closely related different species.However, even for different populations of the same species, the estimates of the exponent of TPL are not identical.The estimate of the exponent for the pooled data for different species within the same taxon was more robust than that obtained using the data of a single species.The estimate of the exponent of TPL for the pooled data of bamboo plants is significantly higher than those for the pooled data of P. subaequalis and Rosaceae.This means that there is a larger variation in the measure of leaf bilateral symmetry for bamboos, which implies that there might be larger light heterogeneity for the leaves at different positions.The goodness of fit of the TPL relationship for bamboos is the best among three classes of plants, which might be related to a good scaling relationship between leaf weight and leaf area.Considering that biomass (or weight) is the best measure for describing TPL, a plant or a class of plants whose leaves have a better scaling relationship between weight and area will have a better TPL for measuring leaf bilateral symmetry.Thus, plants that have the same or similar spatial arrangements of the crown and branches and have the same or similar leaf shape should give a more precise estimate of the exponent of TPL and the goodness of fit for TPL than plants lacking the above characteristics in the above-ground architectural patterns.

Figure 1 .
Figure 1.Illustration of the measure of bilateral symmetry for a leaf of Prunus cerasifera f. atropurpurea.The leaf base is on the left, and the leaf apex is on the right.For convenience, only five equidistant strips are shown.In practice, we used 1000 equidistant strips intersecting the leaf to calculate the mean and variance of the absolute differences between Ai1 and Ai2, where i ranges from 1 to 1000.

Figure 1 .
Figure 1.Illustration of the measure of bilateral symmetry for a leaf of Prunus cerasifera f. atropurpurea.The leaf base is on the left, and the leaf apex is on the right.For convenience, only five equidistant strips are shown.In practice, we used 1000 equidistant strips intersecting the leaf to calculate the mean and variance of the absolute differences between A i1 and A i2 , where i ranges from 1 to 1000.

Figure 2 .
Figure 2. Taylor's power law for the measure of leaf bilateral symmetry for 10 populations of P. subaequalis (Hamamelidaceae).The two-letter title of every panel represents the population code (see Table S1 for details).The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 2 .
Figure 2. Taylor's power law for the measure of leaf bilateral symmetry for 10 populations of P. subaequalis (Hamamelidaceae).The two-letter title of every panel represents the population code (see Table S1 for details).The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 3 .
Figure 3. Boxplots of the estimates of the exponent of Taylor's power law (TPL) for three classes of plants.The values forming the box were calculated using 4000 bootstrap replications.The letters on the top of the boxes are used to indicate significant differences between pairs of plants.Codes sharing the same letter(s) indicate no significant difference; otherwise, there is a significant difference.The significance tests were obtained using the bootstrap percentile method.

Figure 3 .
Figure 3. Boxplots of the estimates of the exponent of Taylor's power law (TPL) for three classes of plants.The values forming the box were calculated using 4000 bootstrap replications.The letters on the top of the boxes are used to indicate significant differences between pairs of plants.Codes sharing the same letter(s) indicate no significant difference; otherwise, there is a significant difference.The significance tests were obtained using the bootstrap percentile method.

Figure 4 .
Figure 4. Taylor's power law for the measure of leaf bilateral symmetry for 10 species of Bambusoideae (Poaceae).The detailed Latin name and collecting information can be found in TableS2.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 4 .
Figure 4. Taylor's power law for the measure of leaf bilateral symmetry for 10 species of Bambusoideae (Poaceae).The detailed Latin name and collecting information can be found in TableS2.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 5 .
Figure 5. Taylor's power law for the measure of leaf bilateral symmetry for 10 species of Rosaceae.The detailed Latin name and collecting information can be found in TableS3.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 5 .
Figure 5. Taylor's power law for the measure of leaf bilateral symmetry for 10 species of Rosaceae.The detailed Latin name and collecting information can be found in TableS3.The small open circles represent the observations, and the straight line represents the fitted regression line based on the ordinary least-squares method.'CI' represents the 95% confidence interval of the estimate of the slope of the fitted linear equation.

Figure 6 .
Figure 6.Comparison of the estimates of the exponent of TPL among the three classes of plants.The dashed line represents the regression line that used the pooled data of the three classes of plants.The fitted results shown at the top-left corner are for the dashed line with the pooled data.

Figure 6 .
Figure 6.Comparison of the estimates of the exponent of TPL among the three classes of plants.The dashed line represents the regression line that used the pooled data of the three classes of plants.The fitted results shown at the top-left corner are for the dashed line with the pooled data.

Table 1 .
Analysis of covariance on the significance of difference of the exponent of TPL among three classes of plants.