Leaf Phenotypic Variation of Chuanqian Crape Myrtle in Natural Populations
Click on author name to view affiliation information
Measurement of morphological factors. LA = leaf area (cm2), LBA = leaf base angle, LC = leaf circumference (cm), LL = leaf length (cm), LW = leaf width (cm), LVA = leaf vein angle.
Cluster analysis of leaf morphological characteristics. Using the coefficient of variation of four leaf traits—leaf vein angle (LVA), leaf area (LA, cm2), leaf base angle (LBA), and leaf length (LL, cm)—the 10 populations were initially grouped into two main categories.
Click on author name to view affiliation information
Leaves are the primary tissue for photosynthesis and transpiration, playing an irreplaceable role in the ecosystem. Understanding variations in leaf morphology across different populations is crucial for better insight into plant classification, evolution, and cultivation practices. From 10 populations of Lagerstroemia excelsa and 80 individual specimens, we measured key leaf traits, including leaf length (LL), leaf width (LW), leaf area (LA), leaf circumference (LC), leaf basal angle (LBA), and leaf vein angle (LVA), to determine the geographic variation trends in leaf morphology across natural populations. Our results revealed several important findings: 1) significant differences in leaf morphology were observed among the populations; 2) temperature appears to be a primary factor influencing the potential geographic distribution of L. excelsa; 3) the 10 populations were grouped into two main categories, indicating that similar germplasm tends to be distributed contiguously. These insights are valuable for future efforts in the introduction, breeding, selection, and cultivation of L. excelsa.
Leaves play a crucial role in photosynthesis and transpiration (Tian et al. 2016), making them indispensable components of the ecosystem. Studying leaf morphology offers valuable insights into how plants absorb water and nutrients, release photosynthetic products, and facilitate gas exchange and light absorption (Nobel 2009). Leaf characteristics are influenced by the interaction between plant genetics and environmental conditions (Falconer and Mackay 1960). Variations in leaf characteristics among plant populations (Wang et al. 2021) reflect their adaptive mechanisms to environmental conditions and serve as indicators of genetic diversity (Cao et al. 2014). Therefore, understanding changes in leaf morphological traits across various plant populations in diverse ecological environments is essential for insights into plant classification, evolution, and the selection of new varieties.
Chuanquan crape myrtle [Lagerstroemia excelsa (Dode) Chun ex S.K. Lee & L.F. Lau], a tall tree species in Lythraceae, is native to China and widely distributed across Sichuan, Guizhou, Hunan, Hubei, and Chongqing (Xiang et al. 2017). This species is popular as an ornamental plant and is notable for its strong disease and insect resistance, long flowering period, and adaptability to stone and pebble soils, making it an essential species for controlling rocky desertification (Zhou et al. 2018). Additionally, its graceful trunk has made it a favored choice for bonsai production. Despite its broad ecological range and potential applications, research on the geographic variation in the leaf morphology of L. excelsa remains limited. Investigating the differences in leaf characteristics among various populations and exploring the response of leaf phenotypic traits to environmental changes is imperative.
Methods such as dispersion degree analysis, principal component analysis (PCA), and cluster analysis of leaf circumference are commonly employed to identify leaf morphological characteristics (Fawad et al. 2023). These methods help researchers understand the variations and differences in leaf traits among plant populations, significantly contributing to leaf morphology classification (Pegiou et al. 2023). However, it is important to note that these analyses primarily focus on comparing plant populations and may not adequately capture the correlation between leaf characteristics and environmental factors. Therefore, incorporating environmental attributes into morphological studies of diverse plant populations is both necessary and valuable.
This study aims to investigate the geographic variation in leaf morphology among ten populations of L. excelsa through statistical analysis and correlation analysis between leaf morphology and environmental factors. Specifically, the research seeks to answer the following questions: (1) Are there significant differences in leaf morphology among different populations of L. excelsa? (2) If such differences exist, are they significantly related to environmental characteristics?
We first gathered distribution data for L. excelsa populations using resources such as the Global Biodiversity Information Facility (http://www.gbif.org), the Chinese Virtual Herbarium (http://www.cvh.org.cn), Plant Plus of China (http://www.iplant.cn), and the China National Knowledge Infrastructure (https://www.cnki.net). Subsequently, 10 populations from five provinces, each characterized by distinct climatic conditions and representing the primary distribution areas, were selected as the research subjects (Fig. 1, Table 1).
Citation: HortScience 60, 3; 10.21273/HORTSCI18375-24
Environmental data were collected from global climate and weather sources (WorldClim, http://www.worldclim.org), which included 19 environmental variables (bio1–bio19) (Table 2) and one terrain variable (Elevation, Ele).
From each population, eight healthy plants with similar growth stages were selected, ensuring more than 30 m between each plant. Twelve leaves were collected from different orientations in the middle section of the branches of each tree. To ensure data consistency, only fully mature leaves were collected during the flowering period, specifically excluding leaves proximal to buds. This strict sampling protocol guaranteed that all analyzed leaves had reached their complete developmental stage. The leaves were then flattened and scanned to measure the leaf vein angle (LVA), leaf area (LA, cm2), leaf base angle (LBA), leaf length (LL, cm), leaf circumference (LC, cm), and leaf width (LW, cm) using AutoCAD (Fig. 1). The aspect ratio (AR) is defined as the ratio of LL to LW. The morphometric factor (MF) is calculated using the formula 4π × LA/LC2, Both AR and MF were subsequently calculated using IBM SPSS 26 (IBM Corp., Armonk, NY, USA).
Four types of analyses were conducted to identify the geographic variation trends in L. excelsa leaf morphology: 1) analysis of intrapopulation leaf trait variation, 2) correlation analysis of leaf traits, 3) correlation analysis between leaf traits and environmental factors, 4) PCA of leaf trait variation among different populations, and 5) cluster analysis and genetic classification of leaf morphological characteristics. All data analyzes were performed using IBM SPSS 26 (IBM Corp.).
The statistical analysis results of phenotypic traits for 10 L. excelsa populations are presented in Table 3. Significant variations were observed among the eight analyzed leaf phenotypic traits, with coefficient of variation (CV) values ranging from 6.94% to 35.37%. The CV values for each trait, ranked from highest to lowest, are as follows: LVA, 35.37%; LA, 24.32%; LBA, 21.79%; LL, 14.62%; LC, 13.57%; LW, 11.00%; AR, 10.02%; and MF, 6.94%.
These results indicate that the MF exhibits the lowest variability, with a CV of only 6.94%, suggesting that the leaf morphology of L. excelsa is relatively stable. In contrast, traits related to leaf area (LVA, LA, LBA) showed significant variations, all exceeding 15%. This suggested that as latitude, longitude, and environmental conditions change, L. excelsa adjusted its leaf size to adapt to its environment. These traits demonstrated high sensitivity to environmental changes, resulting in substantial variations in their characteristics.
The results of the correlation analysis of leaf traits are presented in Table 4. LL showed significant positive correlations with LW, LA, LC, and AR. LW was significantly positively correlated with LL, LA, and LC. LA exhibited significant positive correlations with LL, LW, and LC. LC had significant positive correlations with LL, LW, LA, and AR. AR was significantly positively correlated with LL and LC but showed a significant negative correlation with MF. Notably, the correlation between LL and LC had the highest correlation coefficient of 0.990.
LL is significantly positively correlated with bio5 (maximum temperature of the warmest month), with a correlation coefficient of 0.688, suggesting that LL increases as the maximum temperature of the warmest month rises. Additionally, LL is significantly negatively correlated with elevation, with a coefficient of −0.659, indicating that LL decreases as elevation increases. Similarly, the leaf AR is significantly negatively correlated with elevation, with a coefficient of −0.701, implying that the length-to-width ratio of the leaf also decreases with increasing elevation (Table 5).
The PCA analysis identified two principal components, with contribution rates of 56.3% and 23.8%, respectively. In the analysis of the first principal component, LL and LC had high eigenvalues of 0.987, suggesting that this component can be considered the “leaf size factor.” In the analysis of the second principal component, LW and MF had significant eigenvalues of 0.657 and 0.667, respectively, indicating that the second component can be termed the “leaf shape factor” (Table 6). The cumulative contribution rate of the first two principal components is 80%, and the rapid increase in the cumulative contribution rate of each shape indicator suggests that the genetic diversity of leaf morphology in L. excelsa is relatively limited across multiple dimensions.
The Euclidean distance hierarchical clustering method was used to analyze the leaf characteristics of different Chuanqian crape myrtle populations, with the CV of four leaf traits, LVA, LA, LBA, and LL, as parameters. When the Euclidean distance was set to 25, the 10 populations were initially grouped into two main categories. TR and BL formed one category, whereas WL, TZ, ES, SZ, CB, MT, CN, and GZ formed the other (see Table 1 for population definitions). Notably, ES and SZ exhibited a close relationship, as did CB, MT, CN, and GZ (Fig. 2).
Citation: HortScience 60, 3; 10.21273/HORTSCI18375-24
Studying the phenotypic traits of plant leaves can enhance our understanding of the variation and diversity of L. excelsa, which is beneficial for its introduction, breeding, selection, and cultivation (Henry and Nevo 2014). Through an analysis of phenotypic diversity in 10 populations of L. excelsa, we found that the traits MF and AR changed in response to variations in environmental and geographic factors. Leaf size–related traits such as LVA, LA, and LBA exhibited significant changes, whereas morphological traits such as MF and AR showed minor variations. This suggests that the leaf size of L. excelsa was greatly influenced by environmental conditions, whereas its overall leaf morphology remained relatively stable. The findings also indicate that L. excelsa possessed strong adaptability to the environment, highlighting its breeding potential.
Correlation analysis revealed that LL and AR were significantly negatively correlated with elevation, indicating that the leaf length and aspect ratio of L. excelsa decreased as elevation increases. Previous research had shown that plant leaf size often varied with elevation (Liu et al. 2020; Milla and Reich 2011). As elevation increases, leaf size typically decreases due to significant changes in environmental conditions, such as lower temperatures, reduced air pressure, and diminished light intensity. These changes affect plant growth and survival, prompting adaptive responses. Leaf size is a crucial indicator of a plant’s adaptation to its environment (Wang et al. 2019). At high altitudes with lower light intensity, plants may reduce leaf area to minimize water loss through transpiration and photosynthetic inefficiency, enhancing leaf functional efficiency. This adaptive strategy helped plants survive and reproduce successfully in varied elevation conditions (Gonzalo-Turpin and Hazard 2009), optimizing leaf function under ecological pressures and limited resources (Niklas et al. 2023).
The results also showed that the leaf length of L. excelsa increased with the rise in the maximum temperature of the warmest month. This finding aligned with the distribution patterns of tree species in forests. In warmer, more humid lower-latitude regions, broadleaf forests with larger leaves dominated, whereas colder, drier higher-latitude regions were characterized by coniferous forests with smaller leaves (Kira 1991). The significant increase in leaf length with higher maximum temperatures indicates that L. excelsa develops longer leaves in warmer environments. In lower-latitude regions with higher temperatures, abundant rainfall, and high humidity, the vapor pressure gradient on leaf surfaces was low, which hindered transpiration. A narrower leaf shape, characterized by increased leaf length and minimal width changes, facilitates transpiration and reduces leaf temperature, reflecting an adaptation to high-temperature conditions.
PCA of L. excelsa leaf traits identified the leaf size factor and leaf shape factor as two major principal components influencing leaf morphology, with a cumulative contribution rate of 80%. PCA is a dimensionality reduction technique that transforms multiple correlated variables into a smaller set of uncorrelated variables, known as principal components (Salem and Hussein 2019). These components are arranged by the variance they explain, and the contribution rate shows their explanatory power for the overall data variability (Dunteman 1989). The high contribution rates of the first two principal components suggested that the leaf size and shape factors effectively captured leaf morphological variation, indicating that most variability was concentrated along these dimensions.
Cluster analysis grouped the 10 L. excelsa populations into two main clusters: TR and BL formed one category, and WL, TZ, ES, SZ, CB, MT, CN, and GZ comprised the other. This demonstrated a degree of phenotypic diversity among L. excelsa populations. On the basis of their geographic distribution, most similar germplasm sources are located in contiguous regions. This suggests that over long-term evolution, different populations have undergone directional differentiation under various natural selection pressures, resulting in distinct variations in leaf morphological traits. This study classifies L. excelsa populations based on the morphological traits of leaves from different geographic regions, shedding light on the geographic variation trends in leaf morphology.
This study examined the leaf morphological diversity of L. excelsa from 10 populations across 5 provinces in China, analyzing 960 leaf samples and 8 phenotypic traits using statistical methods. The results revealed significant variation and rich diversity in leaf phenotypic traits, with various degrees of correlation among the 8 traits. Leaf size–related traits of L. excelsa were notably influenced by environmental changes, whereas overall leaf morphology remained relatively stable. The rapid increase in cumulative contribution rates of all traits suggested relatively low genetic diversity in multiple directions of L. excelsa leaf morphology.
Cluster analysis grouped the 10 populations into two main categories, with similar germplasm found to be distributed contiguously. This highlights the distinct geographic characteristics of leaf morphology in L. excelsa and provides valuable insights for the geographic classification of its germplasm. These findings are essential for future efforts in the introduction, breeding, selection, and cultivation of L. excelsa.
Measurement of morphological factors. LA = leaf area (cm2), LBA = leaf base angle, LC = leaf circumference (cm), LL = leaf length (cm), LW = leaf width (cm), LVA = leaf vein angle.
Cluster analysis of leaf morphological characteristics. Using the coefficient of variation of four leaf traits—leaf vein angle (LVA), leaf area (LA, cm2), leaf base angle (LBA), and leaf length (LL, cm)—the 10 populations were initially grouped into two main categories.
Contributor Notes
Measurement of morphological factors. LA = leaf area (cm2), LBA = leaf base angle, LC = leaf circumference (cm), LL = leaf length (cm), LW = leaf width (cm), LVA = leaf vein angle.
Cluster analysis of leaf morphological characteristics. Using the coefficient of variation of four leaf traits—leaf vein angle (LVA), leaf area (LA, cm2), leaf base angle (LBA), and leaf length (LL, cm)—the 10 populations were initially grouped into two main categories.