Aphid Resistance Evaluation and Constitutive Resistance Analysis of Eighteen Lilies

Simple Summary Aphis gossypii is an important pest that harms lilies and limits the development of the lily industry. Improving host plant resistance is an effective and environmentally friendly method for aphid control. We studied the resistance of 16 lily cultivars and 2 wild lily species to A. gossypii and the biological characteristics of their leaves. Among the 18 tested lily plants, ‘Palazzo’, ‘Nymph’, ‘Cameleon’, and Lilium lancifolium showed strong resistance. The analysis of the correlation between the thickness of leaf palisade tissue and the number of greenhouse aphids reveals a significant negative correlation. This indicates that the thicker palisade tissue may be responsible for the observed resistance. Identifying these lilies is important for managing aphid populations and provides additional solutions for lily-integrated pest management. Abstract Lilies (Lilium spp.) are famous bulb flowers worldwide, with high ornamental value. Aphid damage has seriously constrained the development of the lily industry. In this study, the aphid resistance of 16 lily cultivars and 2 wild lily species was characterized in the field and greenhouse. Leaf color parameters, stomatal density and size, thickness of leaf layers, leaf waxy content, and leaf water content were determined to explore the constitutive resistance of lilies. The results show that there was a significant positive correlation between the number of aphids in the field and in the greenhouse (p ≤ 0.05, r = 0.47). This indicated that the level of aphid infestation in both the field and the greenhouse is generally consistent across different types of lily plants. Among these 18 lilies, ‘Palazzo’, ‘Nymph’, ‘Cameleon’ and L. lancifolium were resistant to A. gossypii, while ‘Black Beauty’ and ‘Magnefique’ had poor resistance. The correlation analysis results showed that the number of aphids was negatively correlated with leaf abaxial surface a*, stomatal size, water content, and thickness of leaf palisade tissue and positively correlated with leaf distal axial surface b*, C*, and waxy content. Among them, the correlation between the number of aphids and the thickness of leaf palisade tissue reached a significant level (p ≤ 0.05, r = −0.521). This indicated that the thickness of the palisade tissue of lily leaves might be an important factor influencing the proliferation of aphids. This study not only screened out aphid-resistant lilies but also established a crucial research foundation for the targeted breeding and molecular breeding of lilies with aphid resistance.


Introduction
Insect infestation has always been an important limiting factor for green crop production.With the large-scale development of crop production, the threat of pests has rapidly developed from a few species causing a small amount harm within a small area to large-scale damage being inflicted by multiple species with a high degree of harm [1].Among them, aphids have become a severe pest, affecting the development of production and the improvement of efficiency.Aphids often cluster on the young leaves and buds of plants, causing serious damage through nutrient feeding, honeydew deposition, and viral disease transmission [2].This infestation leads to plants shriveling, affecting flowering and fruiting, and reducing ornamental and application value, resulting in serious economic losses [3].The main method of aphid management is the use of pesticides, which not only increases production costs but also causes harm to the environment [4][5][6][7].Therefore, people are beginning to search for alternative biological control methods for pesticides.Enhancing host plant defenses through breeding is the most sustainable, effective, and eco-friendly method [8].
Plants have their own defense systems that are developed during long-term interactions with insects.These defense mechanisms can be classified into two categories: constitutive resistance and induced resistance [4,9,10].Constitutive resistance is inherent in plants and is related to their genotypes [11,12].It acts throughout their lifespan and consists of physical barriers and chemical barriers [4,10,13].Research on physical barriers in plants currently focuses on factors such as leaf color, leaf trichomes, the wax content of the leaf surface, leaf structure, and so on [14,15].For example, many phytophagous insects show a preference for a yellow-green color and have no tendency for a red color [16][17][18].Resistant cultivars usually have longer and denser trichomes, as well as increased leaf cuticle thickness, palisade tissue thickness, and spongy tissue thickness [15,19,20].
Lilies (Lilium spp.) are important perennial herbaceous bulbous plants that are primarily distributed in the temperate regions of the Northern Hemisphere, including Eastern Asia, Europe, and North America [21].Lilies are cultivated as ornamental plants worldwide, serving various purposes such as cut flowers, potted flowers, and urban landscaping [22].Due to their beautiful flowers and pleasant fragrance, lilies occupy a significant portion of the market in the global cut flower industry and are highly appreciated by customers [23].Additionally, the steroidal saponins, flavonoids, and polysaccharides contained in the Lilium are the main active substances with pharmacological activities, such as anti-tumor, anti-inflammatory, and antioxidant activities [21].Thus, lilies are also widely cultivated as edible and medicinal plants in several Asian countries [22,24] Biotic stress caused by pests and pathogens seriously hinders the sustainable development of the lily industry.Among the pests that harm lilies, cotton aphids are particularly problematic.They tend to gather on the tender leaves and buds of lilies to suck their juices, which poses a significant threat to the plant's health and ornamental value.In addition, cotton aphids can cause sooty blotch and spread viruses such as lily mosaic virus (LMV) and lily ring spot virus (LRSV) [3,25,26].Despite the threat posed by cotton aphids, there are very few studies related to aphid resistance in lilies.Currently, only a few lily cultivars have been studied for their resistance to aphids, and little is known about the mechanisms related to aphid resistance in lilies [3].In this study, the aphid resistance of eighteen lilies was evaluated, then the aphid-resistant mechanisms were systematically explored by means of morphology, physiology, and metabolite analysis.This study aimed to screen out aphid-resistant lily germplasms directly and lays an important research foundation for targeted breeding and molecular breeding.Ultimately, the cultivation of new aphid-resistant lily cultivars is expected to reduce the use of pesticides, lower production costs, and contribute to the green ecological sustainability of the lily industry.

Plant Preparation
The 18 lily plant materials for testing were planted in the National Engineering Research Center for Floriculture (China, Beijing, Changping District, 116.446 • E, 40.151 • N) (Table 1).Before planting, the lily bulbs were soaked in carbendazim solution at 500 mg/L for 40 min, then the rotting and diseased scales were peeled off, and the basal roots were pruned to 1 cm.The greenhouse lilies were grown in pots with single bulbs (diameter of pots: 17 cm, height: 15 cm).The cultivation substrate was a mixture of peat/perlite/vermiculite at a ratio of 5:1:1 with carbendazim added for sterilization.When planting, we kept the bud tip of the lily bulbs upward, and the substrate was 5~8 cm above the top of the lily bulbs.Field planting was carried out in early April.The treated lily bulbs were planted in the experimental field according to the completely randomized zone group method, with the burial depth consistent with that of potted plants.Finally, we watered the bulbs thoroughly to ensure that they were tightly bonded to the soil.

Aphid Rear
The aphids, identified as A. gossypii, were harvested from lilies at the National Engineering Research Center, in Beijing, China (Figure S2).The morphological identification was based on the Journal of Economic Insects of China and the Ecological Atlas of Aphids in Beijing [27,28]. A. gossypii used in the aphid resistance experiment were reared on the lily cultivar 'Magnefique' at 22 ± 2 • C and 65 ± 5% relative humidity with a photoperiod of 16:8 (light: dark) in the breeding greenhouse.

Greenhouse Test
Five fourth-instar A. gossypi apterous nymphs were placed on the young leaves at the top of each lily stem and confined with a voile bag (20 × 30 cm, mesh = 80) when the lilies grew to a height of 10~15 cm [29] (Figure S1).The number of A. gossypii was counted after 10 days.The environmental conditions were kept the same as in Section 2.2.Five biological replicates were set.No insecticides were used during the experiment.

Field Test
When the lilies in the field grew to a height of 10-15 cm, they were manually inoculated with aphids for resistance evaluation tests.Each plant was inoculated with five fourthinstar A. gossypi apterous nymphs, and the number of A. gossypi was counted after 10 days, with 5 biological replicates for each plant material.It was ensured that no insecticides were applied around the test field during the test period.

Determination of Biological Parameters of Lily Leaf 2.4.1. Leaf Preparation
The fifth and sixth leaves under the lily buds were collected to test when the lilies in the greenhouse grew to a height of 10-15 cm.

Leaf Color
To determine L* (the lightness), a* (red to green axis), and b* (yellow to blue axis) of the leaf's adaxial and abaxial surfaces, we used a spectrophotometer (NF555, Nippon Denshoku, Japan) to measure the middle part.We followed a previous study's measuring method [18].Five biological replicates were set.Leaf color was analyzed using two color spaces, CIELAB (L*a*b*) and CIE L*C*h*, selected based on their wide acceptance by the industry and scientific community.C* (Chroma) and h* (Hue angle) were calculated using the following formula: Leaf stomata were investigated by using the imprinting method referenced from a previous study [30].Colorless nail polish was applied to the middle of the lower epidermis of the lily fresh leaves, and gently torn off after the nail polish dried to make a clinical slide.These slides were then observed and photographed with an optical microscope (Sdptop CX40P, Sungrant, Suzhou, China).For each cultivar or species, three samples were selected, and for each sample, five random visual fields were observed.The number of stomata, stomatal length, and stomatal width of each field were counted with Image-Pro plus 6.0 (Media cybernetics, Rockville, MD, USA).Stomatal density, which is the number of stomata per unit leaf area (square millimeters), was calculated.

Leaf Anatomical Structure
The lily fresh leaves were cut crosswise into 0.5 cm-wide slices at the center and fixed with FAA fixative (Solarbio G2350) for 70 h.The fixed material was dehydrated in ethanol, made transparent in xylene, embedded in paraffin, and then sliced using a microtome (Leica RM 2016, Wetzlar, Germany) with a section thickness of 12 µm.Safranine O-fast green staining was performed after the specimen was dried.Then, it was sealed with neutral resin and photographed with an optical microscope (Sdptop CX40P, Sungrant, Suzhou, China).Three samples were selected for each plant material, and each sample was randomly observed in three fields of view.Data were measured by Image-Pro plus 6.0 (Media cybernetics, Rockville, MD, USA).The measured parameters included the thickness of the leaf, upper epidermis, palisade tissue, spongy tissue, and lower epidermis.The method of making paraffin sections and the measurement standard of data can be found in a previous study [31].

Leaf Waxy Content
Leaf waxy content measurement was referenced and improved from a previous study [32].Fresh lily leaves weighing 2 g were cut and placed into a beaker, where they were soaked in 30 mL of chloroform for a minute.Next, the liquid was filtered and transferred into a clean beaker that had been weighed beforehand.The beaker was then weighed again until all the chloroform had evaporated.Three replicates of the experiment were carried out.The waxy content of the leaves was calculated using the following formulas: Leaf waxy content (mg/g) = (W2 − W1) × 1000/2.0g (W1 is the mass of the clean beaker, W2 is the mass of the beaker after chloroform evaporation).

Leaf Water Content
Once the lily leaves were clearly separated and tagged, they were quickly weighed to avoid water loss.Then, the fresh leaves were placed in an oven at a temperature of 110 • C for 10 min and then kept at 80 • C until a constant dry weight was obtained [33].The determination was repeated three times for each plant material.The leaf water content was calculated by using the following formulas: Leaf water content (%) = (fresh weight − dry weight)/fresh weight × 100%

Statistical Analysis
All data were checked for normality and the homogeneity of variance before statistical analysis.The normality of data distributions was tested using the Shapiro-Wilk normality test.The homogeneity of variance between groups was tested using the Levene test for homogeneity of variance.The aphid number, leaf color, leaf thickness, waxy content, and other biological parameters of each lily were tested via one-way ANOVA followed by Duncan's multiple comparisons test with a 95% confidence interval of the difference.Correlation analysis was performed by means of Pearson correlation coefficient calculation.All statistical analyses were performed on SPSS20.0 (IBM SPSS, Somers, NY, USA).Images were created using Origin software (2023b, MicroCal, East Northampton, MA, USA) and Excel 2016 (Microsoft, Redmond, WA, USA).

Aphid Resistance
After 10 days of inoculation, aphids were found to have significantly increased on 'Magnefique' and 'Black Beauty' in the greenhouse, while 'Cameleon' had the lowest aphid proliferation (Figure 1).These results indicate that 'Magnefique' and 'Black Beauty' displayed weaker aphid resistance compared with the other 16 plant materials, while cultivars like 'Cameleon' demonstrated stronger aphid resistance.The field resistance evaluation results also support this finding.However, there are some inconsistencies regarding the aphid resistance of some plants in the greenhouse and field experiments.For example, 'Cameleon' displayed the strongest aphid resistance in the greenhouse, while 'Apricot Fudge' had the lowest aphid proliferation in the field (Figure 1).Moreover, greenhouse aphid populations were higher than field aphid populations for the same lily materials.As shown in the figure, many aphids in the greenhouse gather on lilies' tender leaves and flower buds to absorb the juices (Figure 2).These results are probably because greenhouse conditions are more favorable for aphid growth and reproduction and there is less threat from their natural enemies present in the field [34][35][36].
Although greenhouse cut flower production is their main form of application, lilies are still used in many applications as important flower beds and flower sea plants.Therefore, this study combines the aphid resistance evaluation results of 18 lily plants obtained from greenhouse and field evaluations, which provide higher practical application value for the pest control and cultivation management of these lilies.Correlation analysis showed a significant positive correlation between greenhouse aphid numbers and field aphid numbers (p ≤ 0.05, r = 0.47), indicating that the degree of aphid infestation in the field and greenhouse is generally consistent among different lily plants (Figure 3).This laid the foundation for the further exploration of the aphid resistance of lilies.
x FOR PEER REVIEW 6 of 16  Although greenhouse cut flower production is their main form of application, lilies are still used in many applications as important flower beds and flower sea plants.Therefore, this study combines the aphid resistance evaluation results of 18 lily plants obtained from greenhouse and field evaluations, which provide higher practical application value for the pest control and cultivation management of these lilies.Correlation analysis showed a significant positive correlation between greenhouse aphid numbers and field aphid numbers (p ≤ 0.05, r = 0.47), indicating that the degree of aphid infestation in the   Although greenhouse cut flower production is their main form of application, lilies are still used in many applications as important flower beds and flower sea plants.Therefore, this study combines the aphid resistance evaluation results of 18 lily plants obtained from greenhouse and field evaluations, which provide higher practical application value for the pest control and cultivation management of these lilies.Correlation analysis showed a significant positive correlation between greenhouse aphid numbers and field aphid numbers (p ≤ 0.05, r = 0.47), indicating that the degree of aphid infestation in the field and greenhouse is generally consistent among different lily plants (Figure 3).This laid the foundation for the further exploration of the aphid resistance of lilies.Different plants have various insect resistance strategies such as the waxiness of the epidermis, stomatal density, and deeper vascular bundle burial depth [10].However, it is not yet clear why there are significant differences in resistance among various lily plants.After ten days of inoculation, there were notable differences in aphid populations among different lilies.We hypothesize that certain biological properties of lily leaves play a crucial role in hindering the aphids (Figure 1).

Leaf Color
After analysis, significant differences (p < 0.05) were found among different lilies in L*, a*, b*, C*, and h* on the adaxial and the abaxial surfaces of leaves.It is well-documented that many phytophagous insects, such as aphids, prefer yellow-green light [37,38].Consequently, green leaves are more likely to be infested by various insect species as compared to red leaves [16,17].In this study, 'Secret Kiss' had the smallest distal a* value, indicating the most intense green color, and the largest distal b* value, indicating the yellowest color, and it was the least aphid-resistant (Table 2; Figures 1 and S3).In addition, the highly aphid sensitive 'Magnefique' also had a very low a* value and high b* value on both the adaxial and abaxial surfaces.The aphid-resistant lilies such as 'Palazzo', 'White Triumph', and 'Conca D'or' exhibited a lower value of brightness (L*), yellowness (b*), and chromaticity (C*) and higher values of redness (a*) on both the adaxial and abaxial surfaces.These results align with the theoretical understanding of aphid visual ecology, which suggests that phytophagous insects may favor specific colors or intensities in their preferred plants [39].Different plants have various insect resistance strategies such as the waxiness of the epidermis, stomatal density, and deeper vascular bundle burial depth [10].However, it is not yet clear why there are significant differences in resistance among various lily plants.After ten days of inoculation, there were notable differences in aphid populations among different lilies.We hypothesize that certain biological properties of lily leaves play a crucial role in hindering the aphids (Figure 1).

Leaf Color
After analysis, significant differences (p < 0.05) were found among different lilies in L*, a*, b*, C*, and h* on the adaxial and the abaxial surfaces of leaves.It is well-documented that many phytophagous insects, such as aphids, prefer yellow-green light [37,38].Consequently, green leaves are more likely to be infested by various insect species as compared to red leaves [16,17].In this study, 'Secret Kiss' had the smallest distal a* value, indicating the most intense green color, and the largest distal b* value, indicating the yellowest color, and it was the least aphid-resistant (Table 2; Figures 1 and S3).In addition, the highly aphid sensitive 'Magnefique' also had a very low a* value and high b* value on both the adaxial and abaxial surfaces.The aphid-resistant lilies such as 'Palazzo', 'White Triumph', and 'Conca D'or' exhibited a lower value of brightness (L*), yellowness (b*), and chromaticity (C*) and higher values of redness (a*) on both the adaxial and abaxial surfaces.These results align with the theoretical understanding of aphid visual ecology, which suggests that phytophagous insects may favor specific colors or intensities in their preferred plants [39].Values expressed as mean ± SE (n = 5); means in the same column followed by different lowercase letters are significantly different at the 0.05 level (Duncan' test).

Leaf Stomata
Lilies' stomatal length, stomatal width, and stomatal density of lilies were found to be significantly different (p < 0.05).Previous studies have shown that larger stomatal densities can alter leaf surface topography, leading to behavioral responses from insects [15].Plant aphid-resistant varieties tend to have a smaller stomatal size and higher stomatal density [30].The stomatal length and stomatal width of L. leucanthum, which is weakly resistant to aphids, were significantly smaller than those of other lilies, and the stomatal density was significantly higher than that of other lilies.The more aphid-resistant 'Conca D'or', 'Palazzo' had a larger stomatal length and stomatal width, as well as lower stomatal density (Table 3; Figures 1 and S4).This finding is contrary to previous studies, which have suggested that plants with increased stomatal sizes and greater stomatal densities enable the conversion of additional photoassimilates, resulting in strengthened cell walls and enhanced tolerance to aphids [40].The stomata play a complex role in influencing plant aphid resistance.

Leaf Anatomical Structure
Plants have evolved various constitutive defense systems to protect themselves from insect damage.One such defense mechanism is the presence of palisade tissue and spongy tissue, which act as physical barriers to resist insect attacks [41].Researchers have previously compared the resistance of cucumbers to aphids and found that the deeper the vascular bundles are buried, the stronger the cucumber's resistance to aphids.This may be due to the increased difficulty aphids face in feeding on the vascular sap [31].In this study, significant differences (p < 0.05) in the thickness were observed among all the lilies tested.'Cameleon' exhibited the highest leaf thickness, measuring 570.41 µm, and demonstrated greater resistance to aphids.The weakly aphid-resistant 'Secret Kiss' had the thinnest palisade tissue at 32.87 µm (Table 4; Figures 1 and S5).There are similar cases in previous studies, like a pear tree study which revealed that leaf thickness and palisade tissue thickness are the key factors influencing aphid feeding [31].

Leaf Waxy Content
Leaf surface wax is a layer of lipophilic compounds that covers the plant surface [42].It serves to prevent water loss and resist insect attacks.Leaves with a thicker waxy structure have smoother surfaces, which reduces the attachment ability of phytophagous insects and thus leads to lower egg drop and hatchability, resulting in reduced pest damage.A previous study found that removing the wax from the surface of the leaves promoted cochineal insect feeding on Brassicaceae [43].Additionally, the wax content of cabbage was significantly and negatively correlated with green peach aphid preference [36].
However, in this study, the most severely aphid-infested lilies, 'Magnefique' and 'Black Beauty', had the highest wax content at 1.45 mg/g and 1.44 mg/g, respectively.On the other hand, aphid-resistant lilies, such as 'Conca D or' and 'Palazzo', had the lowest wax content of 0.71 mg/g and 0.79 mg/g, respectively (Figures 1 and 4).This anomaly is not exceptional.Another previous study found that epidermal waxes can also benefit plants by influencing natural enemy insects [44].Lower wax content increases the attachment ability of natural enemy insects such as Hippodamia convergens Linnaeus (Ladybirds, Coleoptera), thereby promoting their predation on phytophagous insects [45,46].These findings highlight the complexity of the role of surface waxes in the tertiary trophic relationship between plants, phytophagous insects, and natural enemy insects.

Leaf Water Content
The water content of lily leaves ranged from 86.02% to 91.67%, with significant differences (p < 0.05).Previous studies have shown a positive correlation between host feeding preference and leaf water content [41].While L. leucanthum is susceptible to aphids, its water content was significantly lower than that of the other lilies, at 86.02% (Figures 1 and 5).The highest water content was found in 'Eyeliner', 'Armandale', at 91.67% and 91.62%, respectively, but they are not aphid-resistant.Previous studies have shown that a decrease in leaf water content and leaf water potential increases the difficulty of aphids to feed in the xylem and shortens the feeding time [47].The decrease in water potential elevated the total amino acid content of the plant, including asparagine and valine, which are critical for aphid performance.However, aphids did not benefit from improved phloem sap quality [48].At a lower water content, the expression of JA-and SA-related genes rose in plants due to drought stress, which in turn acted as a suppressor of aphids [49].The water content of leaves has a more complex relationship with aphid populations.

Leaf Water Content
The water content of lily leaves ranged from 86.02% to 91.67%, with significant differences (p < 0.05).Previous studies have shown a positive correlation between host feeding preference and leaf water content [41].While L. leucanthum is susceptible to aphids, its water content was significantly lower than that of the other lilies, at 86.02% (Figures 1 and  5).The highest water content was found in 'Eyeliner', 'Armandale', at 91.67% and 91.62%, respectively, but they are not aphid-resistant.Previous studies have shown that a decrease in leaf water content and leaf water potential increases the difficulty of aphids to feed in the xylem and shortens the feeding time [47].The decrease in water potential elevated the total amino acid content of the plant, including asparagine and valine, which are critical for aphid performance.However, aphids did not benefit from improved phloem sap quality [48].At a lower water content, the expression of JA-and SA-related genes rose in plants due to drought stress, which in turn acted as a suppressor of aphids [49].The water content of leaves has a more complex relationship with aphid populations.

Correlation between the Aphid Population and Biological Parameters of Lily Leaves
In this study, Pearson's correlation coefficient was utilized to investigate the relationship between various biological parameters of lilies and their ability to resist aphids.The results indicate that the number of aphids in the greenhouse was significantly negatively

Correlation between the Aphid Population and Biological Parameters of Lily Leaves
In this study, Pearson's correlation coefficient was utilized to investigate the relationship between various biological parameters of lilies and their ability to resist aphids.The results indicate that the number of aphids in the greenhouse was significantly negatively correlated with the thickness of the leaf palisade tissue (p ≤ 0.05, r = −0.521),suggesting that thicker leaf tissue made it harder for aphids to feed and resulted in a lower aphid population (Figure 6).This finding is consistent with a previous study which showed that the thickness of palisade tissue was one of the main factors influencing pear's resistance to pear psylla [31].Additionally, the findings indicate that parameters such as stomatal length, stomatal width, water content, and leaf abaxial surface a* were negatively correlated with the aphid population, while wax content, leaf abaxial surface b*, and C* were positively correlated.Although these correlations did not reach statistical significance, they still suggest that these leaf traits have some influence on aphid resistance in lilies.Previous studies have also reported similar findings [10,37], indicating that a combination of factors may determine the resistance of lilies to aphids.

Conclusions
As a limitation to the development of the lily industry, A. gossypii has attracted much attention from researchers.This study evaluated the aphid resistance of eighteen lilies in the field and greenhouse.The correlation between the number of aphids and lily leaf color parameters, stomatal density and size, thickness of leaf layers, leaf surface wax content, and leaf water content was assessed.The results show that 'Palazzo', 'Nymph' and L. lancifolium displayed higher resistance to A. gossypii.In contrast, 'Black Beauty' and 'Magnefique' were more susceptible to A. gossypii than the other lilies.There was a significant

Conclusions
As a limitation to the development of the lily industry, A. gossypii has attracted much attention from researchers.This study evaluated the aphid resistance of eighteen lilies in the field and greenhouse.The correlation between the number of aphids and lily leaf color parameters, stomatal density and size, thickness of leaf layers, leaf surface wax content, and

Figure 3 .
Figure 3. Scatterplot of correlation between greenhouse aphid population and field aphid population (Pearson correlation coefficient).

Figure 3 .
Figure 3. Scatterplot of correlation between greenhouse aphid population and field aphid population (Pearson correlation coefficient).

Insects 2023 ,
14, x FOR PEER REVIEW 12 of 16findings highlight the complexity of the role of surface waxes in the tertiary trophic relationship between plants, phytophagous insects, and natural enemy insects.

Table 2 .
Color parameters of the adaxial and abaxial leaves of eighteen lilies.

Table 3 .
Stomatal parameters of eighteen lilies.Values expressed as mean ± SE (n = 15); means in the same column followed by different lowercase letters are significantly different at the 0.05 level (Duncan's test).

Table 4 .
The thickness of each layer of eighteen lilies.SE (n = 3); means in the same column followed by different lowercase letters are significantly different at the 0.05 level (Duncan's test).