search for


Ecogeographical variations of the vegetative and floral traits of Lilium amabile Palibian
J Plant Biotechnol 2021;48:236-245
Published online December 31, 2021
© 2021 The Korean Society for Plant Biotechnology.

Viet Yen Nguyen ·Rameshwar Rai ·Jong-Hwa Kim ·Ji-Young Kim·Jong-Kuk Na

Department of Horticulture, Kangwon National University, Chuncheon 200-701, Korea
Oriental Bio-herb Research Institute, Kangwon National University, Chuncheon 200-701, Korea
Department of Controlled Agriculture, Kangwon National University, Chuncheon 200-701, Korea
Correspondence to: e-mail:,

The first two authors (VYN and RR) contributed equally to this work.
Received November 6, 2021; Revised December 13, 2021; Accepted December 13, 2021.
cc This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
This study aimed to characterize the morphological variations in the vegetative and floral traits of 73 wild Lilium amabile plants from six habitats in Korea. It was observed that L. amabile is distributed nationwide at any altitude from 300 m (Mt Mangdaeam) to 1550 m (Mt Halla). The majority of the natural habitats of L. amabile were found on mountain slopes, and some were found in rugged mountain regions. The down-facing flowers of this species not only had many blotches but also dense trichomes, and the flowering time was found to be from mid-June to mid-July. ANOVA revealed significant variations in vegetative and floral traits among the six habitats, indicating that the environment has substantial influences on the various growth parameters of L. amabile, such as plant height; number of leaves, bracts, papillae, and flowers; leaf angle; and lengths of the anther, longest blotch, and nectary of the petiole. In addition, the vegetative and floral traits were found closely correlated with each other under the direct impact of the environment. These findings will facilitate to find the appropriate environmental conditions for the conservation and development of L. amabile population as future lily-breeding materials.
Keywords : Population diversity, Morphological traits, Floral traits, Ecogeographical habitats, Lilium amabile

Lilies are one of the important geophytes and these floricultural crop plants are used for the production of cut flowers, for pot plant purposes at home, park, garden, landscape decoration, medicinal purposes, as food, cosmetic, and also used for various ethnic rituals (Dhiman et al. 2018; Grassotti and Gimelli 2011). The genus Lilium comprises nearly 110 species in 7 sections, and there are more than 10,000 commercially important cultivars (Bakhshaie et al. 2016; De Jong 1974; Mathews 2007, Mc Rae 1998). With recent addition of former genus Nomocharis to Lilium, now the genus Lilium includes about 125 species (Gao and Gao 2016), and they are widely distributed throughout the world, mainly in boreal and temperate zones in the Northern Hemisphere, including the Himalayas, Southeast Asia, Japan and China, the Caucasus, and North America (Gao and Gao, 2016; McRae 1998). The east coast of Asia, especially coastal regions of the southwest and north china, is considered as the natural home of almost 59 species of wild lily and is believed to be the center of wild lily diversity (De Jong 1974).

Korea is another major center of wild lily diversity, due to diverse climatic conditions that support Lilium prosperity (Lucidos et al. 2013). Wilson (1925) first reported the distribution of different lily species in Korea. However, the distribution of wild lily species in Korea had not been studied further until Lighty (1969) reported 11 lily species in South Korea. Among 11 wild lily species, seven species were classified into Sinomartagon section, one species into Daurolirion section, and three species into Martagon section (Kim 2008; Lee 2003). Seven species belonging to the Sinomartagon section include L. amabile, L. callosum, L. cernuum, L. concolor, L. lancifolium, and L. leichtlini var maximowizii. However, there are only few studies conducted for morphological characteristics, habitat ecology, and environmental impacts on wild lilies, including L. amabile, in Korea, and rather most previous studies were carried out to identify wild lily species and habitats.

L. amabile is a flowering bulbous perennial wild lily (Kim and Lee, 1990b) and distributed from the southernmost island Jeju (33°N 126°E) to the easternmost island Uleung, Korea (37°N 137°E) (Synge 1980; Wilson 1925). Also, it is distributed in Liaoning province (Haw 1986) and the Manchuria, China (Chung 1965). Although flower and leaf morphology of L. amabile were reported previously (Jeong 1991; Kim and Lee 1990a; Lee 1980; Song et al. 1989), these studies provided only general information (Wilson 1925) and lacked any information about phenotypic variations among different habitats. L. amabile plants can be found both sunny and shade areas with sufficient water supply or drought-prone environment (Chung et al. 2015; Jeong 1991). On the other hand, genetic studies of L. amabile were carried out recently, revealing karyotype and B chromosome variation, cytological variations, and long terminal repeat (LTR) retro transposon diversities among diploid and B-chromosome aneuploidy has been already analyzed in the previous studies (Lee et al. 2019; Nguyen et al. 2019).

The present study was undertaken (1) to determine the difference between the vegetative and flower traits under a given environment that directly impacts the appearance of L. amabile and to unravel (2) correlation between them—thereby understanding the phenotypic plasticity and the habitat’s geographical features of L. amabile—searching for the appropriate environmental type for the conservation and development of the L. amabile population.

Materials and Methods

Population survey and sample collection

The distribution sites of Lilium amabile Palibian in Korea were surveyed to investigate geo-ecological features and its diversities, and variation for different vegetative and floral traits. Six populations representing six different eco-geographical regions/mountains in South Korea were used for the evaluation of the traits. To examine environmental conditions, climatic components of natural habitats were collected using data loggers (Table 1). The populations were consisting of 73 individuals: 24 from mountain (Mt) Moonphilbong, 14 from Mt Gongjak, 6 from Mt Halla, 7 from Mt Maebong, 17 from Mt Mangdaeam, and 6 from Mt Kkotbong (Fig. 1, Table 2). Population habitats were denoted as MP for Mt Moonphilbong, GJ for Mt Gongjak, HL for Mt Halla, MB for Mt Maebong, MD for Mt Mangdaeam, and KB for Kkotbong, respectively (Table 1). The vegetative traits and the floral traits of the indicated number of individuals from each population were examined on the habitats and those individuals were collected for the further study.

Climate characteristics of the sites of L. amabile collection

ID Habitat locality Mean annual rainfall (mm) Mean temperature (°C) Mean relative humidity (%)
MP Mt. Moonphilbong, Danyang-gun, Gangwon-Do 1900 9.1 73.0
GJ Mt. Gongjak, Hongcheon-gun, Gangwon-Do 1405 10.3 69.5
HL Mt. Halla, JeJu-Do 4670 6.55 80.0
MB Mt Maebong, Hongcheon, Gangwon-Do 1000 10.5 71.0
MD* Mt Mangdaeam, Inje-gun, Gangwon-Do 1210 10.1 69.6
KB Mt Kkotbong, Yanggu-gun, Gangwon-Do 1000 10.5 71.0

MP, Moonphilbong; GJ, Gongjak; HL, Halla; MB, Maebong; MD, Mangdaem; KB, Kkotbong.

*The climate characteristics of MD were retrieved from Wikipedia (

Geographical and environmental features of the L. amabile habitats in Korea

ID Population size Number of the studied lines Longitude (°E) Latitude (°N) Altitude (m) Flowering time Habitat environment
MP 300 24 37.05 128.49 400 20-Jun Mountain roadside with herbs and woody plants (trees with the main stem throughout the plant, ~2–3 m high). Loam soil, 10% humidity, sunny weather with open sky.
GJ < 100 13 37.45 128.01 400–750 22-Jun Mountain roadside shaded by big trees. Dry sandy soil, 5% humidity.
HL < 20 6 33.2 126.32 1450–1550 11-Jul Very small population; shaded by azaleas, bamboos, and grasses. The soil of Mt.Halla is mostly composed of volcanic ash, volcanic sand, and lapilli (well-drained humus-enriched soil) and is not the best for farming; 35% humidity. Upper part receives full sun.
MB 100 7 37.81 127.81 400–600 16-Jun Mountain roadside with herbs and woody plants (trees with the main stem throughout the plant, ~2–3 m high). Sandy soil, 10% humidity. Sunny weather with open sky.
MD 200 17 38.05 128.5 300–400 22-Jun Flat hills co-inhabited with herbs and shrubs (any woody plant < 2.5 m high). Sandy loam. Sunny weather with open sky; receives full sun.
KB 100 6 38.12 127.51 339–450 18-Jun Mountain roadside with herbs and woody plants (trees with the main stem throughout the plant, ~3–4 m). Sandy soil; sunny with open sky.
Total 73 16-June–11-July

MP, Moonphilbong; GJ, Gongjak; HL, Halla; MB, Maebong; MD, Mangdaeam; KB, Kkotbong Mountain.

Fig. 1. Locations of the six L. amabile habitats surveyed in Korea (Blue triangles). The solid blue line represents the main mountain range “Baekdudaegan” in the Korean peninsula, and the dotted line represents the mountain range “Nakdongjeongmaek.”

Location of natural habitats

The geographical data of the six natural habitats, including latitude, longitude, and altitude, were recorded using the portable GPS (Garmin Colorado 300, Garmin, Kansas, TX, USA) navigator device. Their relative locations, directions, and angles of the mountain slopes were recorded by using Google maps.

Analysis of phenotypic characteristics

Twenty-five quantitative traits, 11 vegetative and 14 floral traits, were examined for 73 individuals from the six populations for three years, 2018-2020. Eleven vegetative growth characteristics include plant height (PHT), number of leaves (NOL), leaf length at middle (LLM), leaf width at middle (LWM), leaf length at top (LLT), leaf width at top (LWT), leaf length at bottom (LLB), leaf width at bottom (LWB), angle the leaf (AOL) from stem, stem diameter (STD), and number of bract (NOB) per plant. Fourteen floral traits were examined for number of flower (NOF) per stem, length of outer tepal (LOT), width of outer tepal (WOT), length of inner tepal (LIT), width of inner tepal (WIT), blotch number in inner tepal (BNIT), length of longest blotch (LOLB), length of filament (LOF), length of anther (LOA), length of style (LOS), length of ovary (LOO), length of nectary of petiole (NOP), number of papillae (NP), and diameter of flower (DOF).

Data analysis

The data of twenty five morphological traits were subjected to various statistical analyses; descriptive analysis, analysis of variance (ANOVA), Pearson correlations among and between the traits, and principal coordinate analysis (PCoA). To test differences in each trait among six L. amiable populations, we applied one‐way ANOVA followed by least significant difference (LSD) post hoc analyses using SPSS v.20. To analyze genotypic composition and environmental influence on the phenotype, the mean values were used to calculate genetic parameters such as the genotypic coefficient of variation (GCV), phenotypic coefficient of variation (PCV), and coefficient of variation (CV) based on Singh and Choudhary (1979). PCV and GCV were estimated based on the estimation of genotypic and phenotypic variance and expressed in percentage according to Burton and Devane (1953). Then, PCV and GCV were classified to 0~10% (Low), 10~20% (Moderate), and > 20% (High) as suggested by Shivasubramanian and Menon (1973).

Estimations of population subdivisions (QST) were calculated according to QST=VAP/ (VAP+VBP) under the report for a predominantly autogamous species (Badri et al. 2015; Morgan et al. 2001; Spitze 1993), where VAP is population variance component among them and VBP within populations variance component. For each population, broad-sense heritability (H2bs) of the traits was estimated as the ratio of the genetic variance on the sum of genetic variance and environmental variance estimated according to Badri (2007) and categorized as 0~10% (Low), 10~20 (Moderate), and >20% (High). The genetic advance as percent mean (GAM) was calculated as suggested by Johnson et al. (1955) and categorized as 0~10% (Low), 10~20% (Moderate), and >20% (High). The Pearson correlations between measured traits were estimated using SPSS software (SPSS Inc. Chicago, IL, USA) and used PCoA imbedded in GENALEX 6.51b2 software ( The clustering analysis of individual plants was performed based on dissimilarity matrix using Euclidean distances estimated on the mean lines' values with the ward's method in XLSTAT version 2015.1 (Addinsoft Inc., New York, USA).


Distribution, compositional features of the populations, and the natural growth environment of L. amabile in South Korea

To investigate distribution, features of population, and habitat environment of L. amabile in Korea, L. amabile populations were surveyed in six mountains, Mt Moonphilbong, Mt Gongjak, Mt Maebong, Mt Halla, and Mt Kkotbong (Fig. 1). Among the surveyed mountains, Mt Mangdaeam is the lowest habitat with altitude of 300 m and Mt Halla the highest with altitude of 1550 m. A large distribution of L. amabile was observed at an altitude of 300~600 m, in which population size was more than 100 individuals. In case of Mt Halla habitat at an altitude of 1450 m~1550 m, population size was diminished to less than 20 individuals, and habitats are mixed with azalea, bamboo, and grasses (Table 2). The soil of Mt Halla, mostly made out of volcanic ash, volcanic sand, and lapilli (good drain humus soil), is not suitable for farming due to low humidity of 35%. From 400 to 750 m altitude of Mt Gongjak, less than 100 individuals were scattered over the mountain completely covered by big trees and grasses on the dried sandy soil with 5% humidity. However, in an altitude of 300 to 450 m (Mt Moonphilbong, Mt Maebong, Mt Mangdaeam, and Mt Kkotbong), populations were large around 100~300 individuals. But the habitat of L. amablie population in the mountainous areas was very different. At rocky Mt Moonphilbong, habitats were full of herbaceous plant and woody plant (a tree with a straight main stem, about 2~3 m high) on loam soil with of 10% humidity. Habitat in Mt Maebong was very sunny and consisting of herbs and woody plants in sandy soil with 10% humidity. Mt Mangdaeam habitat is on flat sunny hills, and L. amabile grows with herbs and shrubs (woody plants less than 2.5 m high) on sandy loam soil (Table 2).

In this study, the peak of blooming of the L. Amabile in South Korea was between 16 June and 11 July (Table 2). The earliest blooming was observed at Mt Maebong and the latest at Mt Halla. Flowering time tended to be delayed with elevation of latitude. At the peak of blooming, the average temperature and humidity was between 18~23°C and 80~82%.

Diversities and variation in vegetative traits among the studied populations

Surveyed region was not large except for Mt Halla (Fig. 1), but significant variations in growth traits were detected from different geo-graphical regions (Table 3). The plant height (PHT) were ranged from 31.51 to 89.08 cm, the tallest one recorded from Kkotbong populations and the shortest one from Maebong populations. The highest number of leaves were found from Kkotbong populations (i.e 70), whereas the least number of leaves from Mangdaeam populations. Likewise, Kkotbong populations showed the thickest stem diameter with the average of 5.22 mm while the thinnest stem diameter was recorded from Mangdaeam populations with the average of 3.19 mm. Halla and Maebong populations showed high leaf density with 1.28 and 1.76 leaf per centimeter compared to the rest with 0.8~0.9. The leaf shape of L. amabile was a long and narrow, and leaves at the middle of plants were the largest in size, followed by those at bottom and top except for Halla population, in which leaves at the bottom was the largest leaf. Intriguingly, middle leaf size of Mt Gongjak population exhibited significantly bigger than those in upper and lower leaves (Table 3). The leaf angle (AOL) ranged from 38.46 to 85.83°, and the largest leaf angle was observed in Mt Halla and Mt Gongjak population, and the smallest in Mt Moonphillbong.

ANOVA for 11 vegetative traits among the six L. amabile populations

MP 49.33±8.37bc 40.25±10.49b 5.14±1.15cd 0.83±0.16a 2.64±1.10b 0.68±0.18bc 4.22±1.31b 0.74±0.17bc 38.46±7.36c 3.46±0.44b 2.46±1.53bc
GJ 52.10±8.12b 40.62±10.52b 8.31±1.46a 0.96±0.14a 4.85±1.66a 0.80±0.21b 5.58±1.74ab 0.80±0.18ab 81.23±9.24a 3.47±0.75b 1.08±0.28c
HL 41.00±6.60cd 52.50±11.59a 3.87±0.68cd 0.80±0.11a 3.77±0.52ab 1.11±0.08a 5.13±0.53ab 0.95±0.08ab 85.83±1.86a 3.55±0.97b 1.33±0.52c
MB 31.51±4.67d 55.57±17.9ab 4.74±0.62d 0.53±0.11b 2.44±0.0.75b 0.54±0.13c 4.00±1.35b 0.56±0.10c 49.71±7.67b 3.64±0.75b 3.57±1.13b
MD 42.16±7.37bcd 37.59±13.51b 6.12±1.02bc 0.76±0.17a 2.65±0.75b 0.69±0.14bc 5.01±1.04ab 0.77±0.20abc 49.06±7.95b 3.19±0.48b 4.18±1.29b
KB 89.08±10.83a 70.00±21.71a 7.10±1.18ab 0.97±0.19a 3.47±1.06ab 0.74±0.21b 6.22±0.72a 0.99±0.16a 57.17±6.49b 5.22±0.85a 6.50±1.87a

CV (%) 104.82 79.32 86.87 55.99 103.95 70.11 53.34 51.59 155.58 54.83 197.63
LSD5% 13.12 13.36 1.47 0.16 1.16 0.18 1.19 0.17 17.96 0.67 1.67
PCV 33.84 36.54 30.67 23.64 44.45 29.78 29.34 25.61 47.05 23.06 69.99
GCV 29.69 21.07 24.32 15.19 28.12 19.00 13.33 13.40 44.38 14.89 55.31
H2bs (%) 76.98 33.24 62.88 41.26 40.02 40.68 20.65 27.38 85.79 41.67 62.46
QST 0.98 0.88 0.96 0.91 0.90 0.90 0.81 0.85 0.99 0.91 0.96
GAM 53.66 25.02 39.73 20.10 36.65 24.96 12.48 14.45 67.70 19.79 90.04

PHT, plant height; NOL, number of leaves; LLM, leaf length at the middle; LWM, leaf width at the middle; LLT, leaf length at the top; LWT, leaf width at the top; LLB, leaf length at the bottom; LWB, leaf width at the bottom; AOL, angle of leaf; STD, stem diameter; and NOB, number of bracts

The values shown are means ± SD. Uppercase letters in trait data denote significant the differences based on the least significant difference post-hoc analysis (LSD; P < 0.05) after one-way ANOVA.

The phenotypic coefficient of variation (PCV) was greater than a genetic coefficient of variation (GCV) in all studied traits, indicating that habitat environment had an impact on all surveyed characteristics. However, the difference between PCV and GCV showed different levels of the influence of environmental factors on the traits. The biggest difference between PCV and GCV is observed in the LLT (44.45 vs 28.12), LLB (29.34 vs 13.33), NOL (36.54 vs 21.07), and NOB (69.99 vs 55.31), indicating that these traits are highly affected by geographical environment and that they can be improved under an appropriate environmental conditions. Heritability (H2bs %) of 11 growth traits were high as 20.65~85.79%, which can be improved by the selective method. Characteristics of PHT, LLM, LWM, LLT, LWT, AOL, STD, and NOB exhibited high heritability and genetic advance. Thus, these eight traits can be improved by growing conditions. On the other hand, NOL, LLB, and LWB may be regulated by an additive gene action.

Diversities and variations for floral traits among the studied populations

To investigate variation in floral traits, fourteen floral characteristics were examined for all individuals from six populations. The average number of flowers ranged 1~3, 3 in Mt Kkotbong population and 1 in Mt Gongjak population. The average length of outer tepal (LOT) was recorded as the longest with 5.11 mm in Mt Gongjak population and the shortest with 4.01 mm in Mt Kkotbong population. Average length of inner tepal (LIT) was recorded as the longest with 5.03 mm in Mt Gongjak population and the shortest with 4.05 mm in Mt Kkotbong. The blotch number in inner tepal (BNIT) ranged from 48.7 to 82.7, and most abundant on flowers in Mt Maebong and Mt Moonphilbong population. The length of the longest blotch (LOLB) ranged from 2.28 (Mt Halla) to 6.86 mm (Mt Maebong). The length of filament (LOF) ranged from 2.99 to 3.61 cm, and the variation of the length of anther (LOA) ranged from 1.19 to 3.10 cm. The length of style (LOS) varied from 2.41 to 2.97 cm. Mt Halla population exhibited the longest length of ovary (LOO) as 1.63 cm, whereas Mt Mangdaeam population showed the shortest LOO as 1.02. The length of nectary of petiole (NOP) ranged from 1.63 to 3.13 cm, with the longest recorded from Mt Halla population and the shortest from measured for Mt Mangdaeam. In case of the number of papillae (NP), the measured variation ranged from 3.00 to 19.43. The highest NP was recorded from Mt Maebong population and the lowest from Mt Halla. The flower diameter (DOF) varied from 4.82 to 6.08 cm.

As per estimation demonstrated, all floral characteristics showed significant variations among the habitats. These results showed that the flower characteristics in each region have their own specific characteristics. L. amabile in Mt Moonphilbong and Mt Maebong produced thick flowers with numerous black dots on the petals. Also, flowers of Mt Maebong population showed numerous papillae on the inner tepal. Mt Mangdaeam population showed the smallest indices of L. amabile in most floral traits. In contrast, most indices were bigger in Mt Kkotbong population. L. amabile in Mt Gongjak exhibited single cotton and little papillae on the tepal. One of the most intriguing findings among the surveyed populations was that L. amabile population at Mt Halla exhibited exceptionally long anther and the nectary of petiole (Table 4).

ANOVA for 14 floral traits among the six L. amabile populations

MP 1.42±0.58b 4.84±0.75ab 1.49±0.19a 4.63±0.66ab 1.89±0.30a 76.50±22.27a 3.90±2.24bc 3.38±0.43ab 1.29±0.28b 2.74±0.59ab 1.40±0.51ab 1.87±0.37b 13.46±4.39b 6.07±1.09a
GJ 1.15±0.38b 5.11±0.46a 1.34±0.20ab 5.03±0.47a 1.67±0.28a 62.69±20.24ab 2.59±1.24bc 3.50±0.42ab 1.19±0.23b 2.97±0.34ab 1.45±0.30ab 2.16±0.89b 5.49±2.02cd 5.72±0.91ab
HL 1.83±0.75b 4.62±0.25ab 1.45±0.08a 4.53±0.29ab 1.80±0.17a 58.00±12.00ab 2.28±0.55c 3.03±0.25b 3.10±0.18a 3.72±0.39a 1.63±0.08a 3.13±0.43a 3.00±1.79d 4.82±0.29b
MB 2.00±0.82b 4.61±0.62ab 1.36±0.16ab 4.44±0.56ab 1.81±0.19a 82.86±17.40a 6.86±1.35a 3.61±0.41a 1.38±0.23b 2.67±0.29b 1.26±0.21ab 1.87±0.28b 19.43±5.56a 5.87±0.63ab
MD 1.41±0.71b 4.13±0.61b 1.19±0.14b 4.14±0.51b 1.60±0.20a 48.71±16.67b 3.82±1.56bc 2.99±039b 1.21±0.27b 2.41±0.60b 1.02±0.10b 1.63±0.34b 7.18±3.64cd 4.97±0.82ab
KB 3.00±0.63a 4.01±0.63b 1.35±0.18ab 4.05±0.77b 1.79±0.26a 68.00±7.62ab 5.08±1.36ab 3.05±0.55ab 1.28±0.10b 2.43±0.45b 1.22±0.08ab 1.69±0.36b 18.33±4.32ab 6.08±0.72a

CV (%) 116.70 31.68 30.95 27.28 24.59 68.71 119.09 26.69 137.91 48.38 55.32 76.58 186.38 33.33
LSD5% 0.66 0.60 0.17 0.54 0.23 18.37 1.83 0.39 0.48 0.66 0.32 0.53 5.45 0.84
PCV 50.55 15.78 15.03 14.45 15.57 33.63 56.70 14.39 43.01 28.52 29.13 32.76 63.67 17.93
GCV 31.48 8.22 8.09 6.92 5.69 17.93 31.34 6.72 39.21 11.70 14.07 20.72 52.42 8.41
H2bs (%) 38.78 27.16 28.99 22.94 13.36 28.43 30.55 21.82 83.13 16.82 23.35 39.99 67.77 21.99
QST 0.90 0.85 0.85 0.82 0.74 0.85 0.86 0.81 0.98 0.78 0.82 0.90 0.96 0.82
GAM 40.83 8.83 8.98 6.83 4.29 19.69 35.69 6.47 73.65 9.88 14.01 26.99 88.89 8.12

NOF-Number of flowers, LOT-Length of outer tepal, WOT-Width of outer tepal, LIT-Length of inner tepal, WIT-Width of inner tepal, BNIT-Blotch number in inner tepal, LOLB-Length of longest blotch, LOF-Length of filament, LOA-Length of anther, LOS-Length of style, LOO-Length of ovary, NOP-Length of nectary of petiole, NP-Number of papillae, DOF-Diameter of flower.

Values shown are means ± SD. Uppercase letters in trait data denote significant differences based on least significant difference post hoc analysis (LSD; P < 0.05) after one‐way ANOVA

The phenotypic coefficient of variation (PCV) is greater than a genotypic coefficient of variation (GCV) in all studied floral traits. This showed the impact of environmental factors on the surveyed characteristics. However, the distance between PCV and GCV showed different levels of the influence of environmental factors on the traits. The biggest differences between PCV and GCV were observed in NOF (50.55; 31.48), BNIT (33.63; 17.93), LOLB (56.70; 31.34), LOS (28.52; 14.07), and LOO (29.13; 14.07), respectively. The heritability (H2bs %) estimation of all 14 floral traits ranged from 13.36 to 83.13%. These traits possessed high heritability except for WIT and LOS, which can be improved by the selective method. NOF, LOLB, LOA, NOP, and NP traits have possessed high heritability and high genetic advances. However, LOT, WOT, LIT, and DOF traits have demonstrated high heritability but low genetic advances. Thus, in the studied traits, NOF, LOLB, LOA, NOP, NP can be improved by growing condition. LOT, WOT, LIT, and DOF traits may be regulated by an additive gene action.

Correlation among the vegetative and floral traits

Analysis of the Pearson correlations was conducted for 300 possible correlations among the measured traits, resulting in 122 pairs of significant correlations among the traits. Among them, 90 pairs showed very significant correlations, in which 75 exhibited positive correlation while 15 pairs showed negative correlations (Table 5). Most correlations between the parameters of vegetative (aerial) growth were positive, and also not most but many pairs exhibited positive correlations among floral traits (Table 5). The flower numbers (NOF) were positively correlated with PHT, NOL, STD, NOB and NP. The LOT, WOT, LIT, and WIT were positively correlated with most of floral traits. The NP was negatively correlated with LLT, LWT, LLB, and AOL, but it was positively correlated with NOB, NOF, BNIT, LOLB, and DOF (Table 5).

Phenotypic correlation coefficient among 25 quantitative traits in L. amabile palibian

NOL 0.47** 1
LLM 0.47** 0.05ns 1
LWM 0.52** 0.17ns 0.45** 1
LLT 0.36** 0.09ns 0.63** 0.46** 1
LWT 0.35** 0.31** 0.18 ns 0.55** 0.58** 1
LLB 0.44** 0.16ns 0.62** 0.36** 0.55** 0.44** 1
LWB 0.51** 0.23* 0.26* 0.61** 0.43** 0.71** 0.52** 1
AOL 0.10ns 0.14ns 0.37** 0.26* 0.54** 0.45** 0.33** 0.27* 1
STD 0.61** 0.74** 0.16ns 0.29* 0.13ns 0.27* 0.14ns 0.25* 0.11ns 1
NOB 0.31** 0.46** 0.07ns -0.10ns -0.11ns 0.05ns 0.18ns 0.11ns -0.29* 0.41** 1
NOF 0.41** 0.74** -0.08ns 0.12ns -0.04ns 0.25* 0.10ns 0.26* 0.00ns 0.59** 0.52** 1
LOT 0.07ns -0.12ns 0.20ns 0.19ns 0.32** 0.18ns 0.07ns 0.06ns 0.17ns -0.06ns -0.44** -0.25* 1
WOT 0.13ns 0.04ns -0.17ns 0.09ns 0.06ns 0.20ns -0.12ns 0.08ns -0.05ns 0.14ns -0.32** 0.03ns 0.61** 1
LIT 0.01ns -0.15ns 0.20ns 0.13ns 0.33** 0.09ns 0.03ns 0.01ns 0.24* -0.04ns -0.47** -0.22ns 0.88** 0.51** 1
WIT 0.14ns 0.00ns -0.14ns 0.04ns -0.06ns 0.04ns -0.19ns 0.01ns -0.13ns 0.10ns -0.12ns 0.07ns 0.50** 0.76** 0.42** 1
BNIT 0.09ns 0.11ns -0.02ns 0.04ns 0.09ns 0.04ns -0.02ns 0.07ns -0.20ns 0.26* 0.02ns 0.14ns 0.39** 0.52** 0.36** 0.49** 1
LOLB 0.04ns 0.17ns -0.13ns -0.32** -0.29** -0.31** -0.18ns -0.23ns -0.24* 0.19ns 0.30** 0.25* 0.01ns 0.07ns 0.01ns 0.18ns 0.41** 1
LOF -0.02ns -0.04ns 0.10ns 0.04ns 0.06ns -0.10ns -0.05ns -0.04ns 0.20ns 0.07ns -0.33** -0.12ns 0.69** 0.49** 0.68** 0.44** 0.38** 0.10ns 1
LOA -0.08ns 0.22ns -0.31** -0.03ns 0.13ns 0.50** 0.04ns 0.28* 0.41** 0.05ns -0.23ns 0.12ns 0.19ns 0.30** 0.16ns 0.21ns 0.05ns -0.15ns -0.04ns 1
LOS 0.10ns 0.09ns -0.06ns 0.12ns 0.35** 0.43** 0.17ns 0.32** 0.33ns -0.03ns -0.29* -0.02ns 0.47** 0.42** 0.40** 0.37** 0.22ns -0.04ns 0.21ns 0.48** 1
LOO 0.00ns 0.00ns -0.10ns 0.12ns 0.20ns 0.14ns 0.07ns 0.03ns 0.24* -0.01ns -0.33** -0.08ns 0.32** 0.25* 0.30* -0.04ns 0.14ns -0.14ns 0.09ns 0.25* -0.10ns 1
NOP -0.10ns -0.02ns -0.16ns -0.10ns 0.09ns 0.25* -0.17ns 0.08ns -0.38** -0.15ns -0.42** -0.06ns 0.33** 0.20ns 0.30* 0.17ns -0.12ns -0.17ns 0.14ns 0.55** 0.45** 0.15ns 1
DOF 0.26* 0.17ns 0.08ns 0.16ns 0.18ns 0.09ns 0.07ns 0.15ns -0.17ns 0.27* 0.19ns 0.20ns 0.34** 0.25* 0.30* 0.23* 0.53** 0.18ns 0.21ns -0.13ns 0.04ns 0.19ns -0.12ns 1
NP 0.22ns 0.28* -0.19ns -0.20ns -0.34** -0.33** -0.30** -0.17ns -0.55** 0.30* 0.32** 0.32** -0.10ns 0.14ns -0.19ns 0.24* 0.44** 0.47** 0.06ns -0.23ns -0.22ns -0.13ns -0.24* 0.31** 1

PHT, plant height; NOL, number of leaves; LLM , leaf length at the middle; LWM, leaf width at the middle; LLT, leaf length at the top; LWT, leaf width at the top; LLB, leaf length at the bottom; LWB, leaf width at the bottom; AOL, angle of leaf; STD, stem diameter; NOB, number of bracts; NOF, number of flowers; LOT, length of the outer tepal; WOT, width of the outer tepal; LIT, length of the inner tepal; WIT, width of the inner tepal; BNIT, blotch number in the inner tepal; LOLB, length of the longest blotch; LOF, length of filament; LOA, length of anther; LOS, length of style; LOO, length of ovary; NOP, length of the nectary of petiole; NP, number of papillae; DOF, diameter of flower.

*P = 0.05, significant; **P = 0.01, highly significant; ns: not significant.

Principal coordinates analysis (PCoA) resulted in three principal coordinates, relevant to the morphological variations. These three principal coordinates accounted for 52.16% of the total variations. The first and second coordinates accounted for 21.21% and 19.66% of the total variations, respectively (Fig. 2). PCoA analysis revealed a significant separation between vegetative and floral traits except for NOF, number of flower, which exhibited close relationship with vegetative traits rather than floral traits (Fig. 2). There were positive correlations among NOF, PHT, NOL, STD, and NOB and among AOL, LLT, LWT, LWB, LLB, and LWM (Fig. 2). Similarly most floral traits except for NOF exhibited close relationships. Among floral traits, there were close relationships among LIT, LOT, LOS, LOA, LOO, NOP, LOF, WOT, and WIT, and among BNIT, DOF, LOLB, and NP (Fig. 2).

Fig. 2. PCoA analysis based on 25 morphological traits. Red and black rhombi depict vegetative and floral traits, respectively

Genetic diversity among the individuals within the sampled populations

Based on the characteristics of 25 morphological traits, a tree-shaped diagram was generated to show genetic relationships among L. amabile individuals from the six habitats (Fig. 1). The hierarchical cluster analysis was performed to examine the aggregation patterns of 73 lines of L. amabile, resulting in three groups (Fig. 3). The subclasses were classified based on morphology. The group 1 was consisting of 23 lines including each 3 from Moonphilbong and Gongjak, 16 from Mangdaeam, and 1 line from Kkotbong. This groups included those individuals with the smallest average of the traits. The group 2 was comprised of 41 lines, 21 from Moonphilbong, 4 from Gongjak, 6 from Halla, 7 from Maebong, 1 from Mangdaeam and 2 from Kkotbong. The second group is the group with the average of small vegetative traits but the largest floral traits.The group 3 consisted of 9 lines including 6 lines from Gongjak and 3 from Kkotbong. This group was characterized with a group with higher average in vegetative and floral traits. Taken together, the cluster analysis revealed that the individuals in the same group are closely related to each other in morphological characteristics.

Fig. 3. Dendrogram based on the morphological relationships among the L. amabile individuals from six different habitats in Korea

Lilium, as one of the major genus of geophytes, is grown worldwide for diverse ornamental purposes. This genus possesses a great diversity in many horticultural traits which is manifested in flower color, shape, size, fragrance, disease resistance, and many physiological characteristics (Dhiman et al. 2015). Lilium has comprised seven sections, among which three sections, Archelirion, Leucolirion, and Sinomartagon, are considered as the major economic groups likely due to the availability of the wide genetic diversities and various morphological characteristics. Wild lily species differ in crossability more or less, but the wild species (2n=2x=24) produces fertile hybrids within each sections (McRae 1998; Van Tuyl et al. 2011). L. amabile, one of the endemic wild lily species in Korea, belongs to the section of Sinomartagon, and it could provide not only a key for understanding evolutionary pass of Korean wild lilies but also an important material for the new novel breeding of the Lilium. However, many endemic species are getting threatened by global climate change as well as habitat destruction by human. Therefore, the understanding of endemic species, including the levels and forms of genetic variations, would be crucial for the successful protection, and utilization of rare and endangered species (Dhiman et al. 2018).

In this study, it was examined whether eco-geographical difference in the Korean Peninsula could result in correlated variations in morphological characteristics of L. amabile. The difference in the habitat environment caused significant morphological changes of L. amabile in both vegetative and floral traits, and some degree of distinctive regional characteristics was observed in some populations. L. amabile in the Mt Moonphilbong grows on steep hills with grasses and small trees, and it has small leaf angle (38.46°) and small leaf size that likely help to reduce direct exposure to sunlight and water loss in arid climates (Fig. 4A). L. amabile plants from Mt Halla habitat with grasses, azaleas (shrub) and bamboos showed a tendency to grow tilting towards sunlight. Their leaf blade spreads nearly perpendicular to the stem and the top leaf size is long and wide to receive maximum sunlight (Fig. 4B). In Mt Gongjak, L. amabile grows under the shade of big trees, resulting in a large leaf angle with tall stems and long and thin leaves likely to receive the maximum amount of sunlight (Fig. 4C). In the Kkotbong region, L. amabile populations showed a large variation in vegetative growth traits compared to other populations. Among floral traits of Mt Halla population, the length of anther, style, ovary, nectary of petiole is longer than those in other populations (Fig. 5A, Table 4). Population from Mt Gongjak exhibited a large leaf angle, long and thin leaves, and few bracts (Fig. 4C, table 3). Plants in Mt Moonphilbong was medium in height among six populations, but flower has large diameter, many blotches on the tepal, and vivid color (Fig. 5B, Table 3, 4). In Mt Maebong, L. amabile showed characteristics of arid area plants by small thick leaves and many big blotches in tepal (Fig. 4A Fig. 5B, Table 3~4).

Fig. 4. Leaf morphological variations of L. amabile among different geographical regions of Korea. A. Normal type (Mt Moonphilbong, Mt Mangdaeam, Mt Maebong, and Mt Kkotbong), B. Short-straight-leaf perpendicular to the stem (Mt Halla), C. Long curved leaf with a wide leaf angle (Mt Gongjak).

Fig. 5. Variations in L. amabile flower color and blotch among different habitats in Korea. The flowers showed distinctive regional characteristics. Flowers with A) orange-red color and small blotches from Mt Halla and Mt Mangdaeam; B) orange-red color and large blotches from Mt Moonphilbong, Mt Maebong, and Mt Kkotbong; and C) orange color with small blotches from Mt Gongjak.

Taken together, our results showed that habitat has strong effects on the morphological characteristics of L. amabile (Table 3~4), indicating the geographical characteristics of habitat contribute differently to ecological diversification most likely due to different environment. This may explains why PCV was higher than GCV in all surveyed traits (Table 3~4), which implies that L. amabile may be sensitive to environment changes and adapted to such environment. Similar to morphological changes observed in L. amabile, another Korean wild lily, L. tsingtauense, showed morphological changes in response to habitat environment. It developed verticillate leaves at the upper parts of plants, which could be another form of ecological adaptation to receive maximum sunlight because its habitats were under the shade of trees and shrubs (Lim 2013) as seen at some of L. amabile habitats. The geographical variation in morphology reveals phenotypic modifications to habitat environment and evolutionary pattern of plant species, and such morphological changes specify the local regional deviations (Du et al. 2014; Ellison et al. 2004). Also, the patterns of species distribution in different zones could be a consequence by various interacting factors, such as geographical features, competition and dynamic of regional species, regional soil conditions, and environmental variables (Criddle et al. 2003; Eriksson. 1996). Most of the lilies thriving in South Korea are colored, and those presenting varied flower colors and shapes receive the greatest ornamental and commercial value. The flowers of L. amabile have a turk’s-cap shape. The distinctive characteristic of L. amabile is the presence of white, stiff hairs covering the entire plant.


This study showed that morphological traits are useful tools for finding variation within and among populations of L. amabile. Our results suggest that the evolutionary trend of vegetative and floral traits could be related to the selective pressure to better adaptation to the given environmental conditions, especially extreme environment conditions. Intriguingly, it was observed that significant morphological variations in L. amabile occurred in relatively small region, suggesting that L. amabile may be a fast adapting species to local climates as seen many invasive plant species. Further population and genetic studies would give more insights of L. amabile morphological variations in a broad area with significant environmental changes.


This study was carried out with the support of the GSP Project No. 213007-05-5-SBM10, the Ministry of Agriculture, Food and Rural Affairs, Republic of Korea.

Conflict of Interests

The authors declare that there is no conflict of interests.

  1. Ellison AM, Buckley HL, Miller TE, Gotelli NJ (2004) Morphological variation in Sarracenia purpurea (Sarraceniaceae): geographic, environmental, and taxonomic correlates. Am J Bot. 91(11):1930-1935
    Pubmed CrossRef
  2. Badri M, Ilahi H, Huguet T, Aouani ME (2007). Quantitative and molecular genetic variation in sympatric populations of Medicago laciniata and M. truncatula (Fabaceae): relationships with ecogeographical geographical factors. Genet. Res. 89: 107-122
    Pubmed CrossRef
  3. Badri M, Bouhaouel I, Arraouadi S, Tammalli W, Huguet T, Aouani ME (2016) Variation in tolerance among Tunisian populations of Medicago truncatula. Plant Genetic Resources: Characterization and Utilization, 14:41- 49
  4. Burton GW and Devane EW (1953) Estimating heritability in tall fascue from replicated clonal material. Agron. J. 45:478-481
  5. Chung MY, López-Pujol J, Chung JM, Kim KJ, Park SJ, and Chung MG (2015). Polyploidy in Lilium lancifolium: evidence of autotriploidy and no niche divergence between diploid and triploid cytotypes in their native ranges. Flora 213:57-68. doi: 10.1016/j.flora.2015.04.002
  6. Chung T (1965) Illustrated encyclopedia of Fauna & Flora of Korea. Vol. 5. Tracheophyta. Genus Chrysosplenium. Ministry of Education. Seoul. pp. 497-483
  7. Criddle RS, Church JN, Smith BN, Hansen LD (2003) Fundamental Causes of the Global Patterns of Species Range and Richness1. Russian Journal of Plant Physiology 50, 192-199
  8. De Jong P (1974) Some notes on the evolution of lilies. Lily Yearbook, vol 27, pp. 23-28
  9. Dhiman MR, Moudgil S, Prakash C, Kumar R, and Kumar S (2015) Biodiversity in Lilium: A Review. International Journal of Horticulture. 8(8):83-97
  10. Du YP, He HB, Wang ZX, Wei C, Li S, Jia GX (2014) Investigation and evaluation of the genus Lilium resources native to China. Genet Resour Crop Evol 61:395-412
  11. Eriksson O (1996) Regional dynamics of plants: a review of evidence for remnant, source-sink and metapopulations. Oikos 77:248- 258
  12. Gao YD, Gao XF (2016) “Accommodating Nomocharis in Lilium (LIliaceae),” Phytotaxa. 277(2):205-210
  13. Grasotti A, Gimelli F (2011) Bulb and cut flower production in the genus Lilium: current status and the future. Acta Hortic 900:21-42
  14. Haw SG (1986) The lilies of China. B.T. Batsford Ltd., London
  15. Jeong JH, Kim KS, and Hong YP (1991). Distribution of Korean native lilies and environmental conditions of their native habitats. J. Kor. Soc. Hort. Sci. 32:270-277
  16. Johnson HW, Robinson HF, Comstock RE (1955) Estimates of genetic and environmental variability in soybeans. Agronomy Journal. 47(7):314-318
  17. Kim T (2008) Wild flowers and resources plants in Korea. 1-5, Seoul: SNU press
  18. Kim YS, Lee WB (1990a) A study of morphological characters on the genus Lilium L. in Korea. Korean Journal of Plant Taxonomy. 20(3):165-178
  19. Kim YS, Lee WB (1990b) A study of anatomical characters on the genus Lilium L. in Korea,” Korean Journal of Plant Taxonomy, vol. 20(3):179-189
  20. Lee C (2003) Coloured flora of Korea 1-2, Seoul: Hyangmunsa.
  21. Lee SI, Nguyen TX, Kim JH, Kim NS (2019) Cytological variations and long terminal repeat (LTR) retrotransposon diversities among diploids and B-chromosome aneuploids in Lilium amabile Palibin. Genes & genomics. 41(8):941-950
  22. Lee T (1980) Illustrated flora of Korea. Hyangmunsa, Seoul.
  23. Lighty RW (1969) The lilies of Korea. In: The lily yearbook, RHS, London. 31:31-39
  24. Lim KB (2013) Korean species lilies (Lilium tsingtauense). Lilies and related plants. RHS. 87-92
  25. Lucidos JGB, Kwang R, Younis A, Kim CK, Hwang YJ, Son BG and Lim KB (2013) Different day and night temperatures responses in Lilium hansonii in relation to growth and flower development. Hort Environ Biotech. 54:405-411
  26. Matthews V (2007) The International Lily Register and Checklist. RHS London. 4th Edition
  27. McRae E (1998) Lily species., In Lilies, Timber Press, Portland, Ore.USA, pp. 319-330
  28. Morgan KK, Hicks J, Spitze K, Latta L, Pfrender ME, Weaver CS et al (2001). Patterns of genetic architecture for life-history traits and molecular markers in a subdivided species. Evolution 55:1753-1761
    Pubmed CrossRef
  29. Nguyen TX, Kim BR, Park DR, Kim YK, Nguyen VY, Na JK, Kim NS, Kim JH (2019) Karyotype and B chromosome variation in Lilium amabile Palibin. Genes & genomics. 41(6):647-655
  30. Singh RK and Choudhary BD (1979), Biometrical methods in quantitative genetic analysis, New Delhi: Kalyani Publisher
  31. Sivasubramanian S, Menon M (1973) Heterosis and inbreeding depression in rice. Madras Agric. J. 60(7):1139-1140
  32. Song JT, Jeong HB, Kim BW, Jin HS (1989) The Sauras of Korean resourses plants, Seoul. Ilheung Co. pp. 512-513
  33. Spitze K (1993) Population-Structure in Daphnia obtusa - Quantitative genetic and allozymic variation. Genetics, 135(2):367-374
    Pubmed CrossRef
  34. Synge PM (1980) Lilies, London: B.T. Batsford
  35. Van Tuyl JM, Arens P (2011) Lilium breeding history of the mordern cultivar assortment. Acta Horticult 900:223-230
  36. Wilson EH (1925) The lilies of Eastern Asia. Dulau & Co. Ltd, London. pp 110

December 2021, 48 (4)
Full Text(PDF) Free

Social Network Service

Cited By Articles
  • CrossRef (0)

Funding Information
  • CrossMark
  • Crossref TDM