Skip to main content

Morphological variation in Cynodon dactylon (L.) Pers., and its relationship with the environment along a longitudinal gradient



Geographical variation in morphological traits may reflect evolutionary patterns of morphological adaptability along environmental gradients. Comprehensive information on longitudinal patterns of morphological trait variation is very meaningful to explore morphological diversity and evolutionary trends in widespread bermudagrass.


To explore the spatial patterns of morphological traits, we investigated 10 morphological traits of bermudagrass and 10 soil nutrient indexes and collected local climate data for 13 different regions from 119°E to 105°E along the latitude 34°N.


Considerable variations in morphological traits were observed at different longitudes, and the variations in most of the evaluated traits within populations were lower than those among populations. All of the 13 different longitudinal sites were divided into three groups based on morphological traits by cluster analysis. The major sources of diversity at the different longitudes were leaf length of the erect shoot, leaf width of the erect shoot, and the internode lengths of the erect shoot and stolon as determined by principal component analysis. Pearson correlation analysis also indicated that longitude was significantly and negatively correlated with these traits as well. Mean average rainfall was significantly correlated with leaf length of the erect shoot and the internode lengths of the erect shoot and stolon, while mean average temperature was only significantly correlated with internode length of the erect shoots. Available sulfur was significantly correlated with internode length of the erect shoot, plant height, and reproductive branch height, while the exchangeable Ca was significantly correlated with internode lengths of the erect shoot and stolon. Soil pH was significantly correlated with the internode length of the stolon. Longitude is an important factor that affects morphological trait variation in wild bermudagrass, and the leaves of the erect shoot and the internode length enlarged significantly with the collection sites moving from east to west.


Different combinations and interactions of environmental factors (soil and climate) along a longitudinal gradient may have strong effects on one or more morphological traits of bermudagrass.


Phenotypic plasticity is a widespread and important phenomenon that allows plants to adapt to the environment, and the adaptive mechanism may be related to the specific environments in which the plants live; that is, different environments may select for plants that have different phenotypic plasticity mechanisms [1,2,3,4]. It is generally accepted that fluctuating environmental conditions play a key role in the morphological variation in plants [5, 6]. Morphological variation in plants is generally a feedback for the changing climatic conditions and reflects the adaptive evolution [7,8,9,10,11,12], while the soil factor has apparently shown to be a driving force to determine the morphological traits of plants [13,14,15,16]. Morphological variation can evolve in response to environmental variability and may show clear patterns in geographical regions; for example, the distribution of Drosophila species is implicated in its innate thermal tolerance limits, and a majority of measured morphological traits of Sorghum bicolor Moench show significant regional differences, etc. [17,18,19,20]. However, few studies have addressed whether plant species would show appreciable changes or a significant pattern of morphology along a longitudinal gradient. The climatic variability hypothesis suggests that individuals of a species may exhibit larger morphological variation with increasing climatic variability, which could happen in areas of lower precipitation at decreased longitude [21, 22]. The leaf length and width of Pistacia atlantica Desf. subsp. Atlantica were found to be significantly related to precipitation, and precipitation is the major factor that determines plant functional traits in Stipa species [23, 24]. It is generally accepted that precipitation varies considerably with changes in longitude in many areas of the world. In the climate of China, precipitation also decreases linearly with increasing longitude from east to west in regions far from the sea [25]. At present, there is no information about the effects of climatic changes on morphological traits in widespread plant species, which makes it unclear how plants adapt to regular changes in climate, especially precipitation.

Bermudagrass, Cynodon dactylon (L.) Pers., is an important perennial grass species in the botanical family Poaceae [26]. Originally native to Africa, C. dactylon has become naturalized in many places around the world and is widely distributed in tropical, subtropical, and temperate coastal areas [27]. C. dactylon has the advantages of strong viability, rapid reproduction, and resistance to trampling, making it a suitable grass for lawns, livestock forage, soil conservation, and remediation of contaminated soils [12, 13]. As early as the 1960s, researchers at Oklahoma State University divided C. dactylon germplasm into six variants based on their growth habits, leaf color, plant type, and somatic chromosome count [14]. Bermudagrass was also divided into four biological types based on external morphological traits including leaf color, stem color, and shaped type in Mauritius [28]. Harlan et al. (1970) divided the C. dactylon collections from around the world into nine species and 10 variants, and Taliaferro reviewed and revised the taxonomy of C. dactylon [14, 29]. In the previous studies that explored the geographical clines in morphology for C. dactylon, most were focused on different habitats [30]. We can ask two questions; (1) how does the widespread C. dactylon adapt to the environmental changes along a longitudinal cline by morphological variation? And (2) what relationships are there between morphological traits and climate? At present, no studies have been conducted on morphological variation in C. dactylon that are strictly related to different longitudes to answer these questions. Thus, in our study, we measured 10 quantitative morphological traits in plants of C. dactylon collected from locations at 13 different longitudes and analyzed the relationships between these variations and the local soil and climatic conditions. The objectives of this research were to i) analyze the variability in 10 quantitative morphological traits of C. dactylon collected at sites representing 13 different longitudes, ii) explore the pattern of morphological variation among and within populations, and iii) evaluate the relationships between morphological traits, soil, and climate at the different longitudes.

Material and methods

Plant material

In August 2016, a total of 260 wild bermudagrass individuals were sampled as representatives of 13 populations growing in regions at different longitudes from 119°E to 105°E along latitude 34°N in southern China (Table 1). Meteorological data, including the mean average rainfall (MAR) and the mean average temperatures (MAT) over thirty years (1981–2010) were recorded for the 13 sampled regions (Table 1).

Table 1 Geographical and climatic factors of bermudagrass in 13 different longitude regions

Morphological trait data collection

The morphological characters, including leaf length of the erect shoot (LLES), leaf width of the erect shoot (LWES), internode length of the erect shoot (ILES), diameter of the erect shoot (DES), leaf length of the stolon (LLS), leaf width of the stolon (LWS), internode length of the stolon (ILS), diameter of the stolon (DS), plant height (PHT), and reproductive branch height (RBH) were measured for 20 individual plants per site. The sampling method was as described in our previous study [31]. The samples were measured with the aid of a caliper during the highest growth period in August. The fourth mature leaf from the top of the erect shoot to the base was selected for measurement of LLES and LWES, and the fourth section stems from the top of the erect shoot to the basal part were used for measurement of DES and ILES. The methods used to measure the stolon characters LLS, LWS, DS, and ILS were identical to those used to measure the erect shoot above. The natural height of bermudagrass was measured for PHT, and the vertical height from the surface to the natural top of the reproductive branch was measured for RBH. The values for 10 quantitative morphological traits were used for further data analysis. The climate data for the sampling sites was obtained from the website of the China meteorological data network.

Soil chemical characteristics

The soil samples (0–20 cm depth) were collected separately under the cover of the wild bermudagrass from 20 corresponding quadrats (10 m × 10 m) at each site, and the soil was collected using the principle of multi-point sampling. The soil samples were taken within different quadrats and mixed at each site. The samples were then sealed in polyethylene bags using the quarterly method to divide the soil into portions of approximately 1 kg. In order to facilitate the different analyses that require grinding and screening, the collected soil samples were transported to the laboratory and allowed to air-dry at 25 °C. The soil samples were analyzed for total nitrogen (TN), available nitrogen (AN), available sulfur (AS), available kalium (AK), available phosphate (AP), soil organic matter (SOM), exchangeable Ca, exchangeable Na, exchangeable Mg, and pH based on standard methodologies [32]. Briefly, the SOM was determined by the potassium dichromate wet combustion procedure. TN and AN were analyzed using the Kjeldahl method and the alkaline nitrogen diffusion method, respectively. The AS, AK, AP were determined using the barium chloride turbidimetric method, flame emission spectrometry and Olsen method, respectively. The exchangeable Ca, exchangeable Na, and exchangeable Mg were measured by continuum-source atomic-absorption spectrometry. Soil pH was tested with a glass electrode (soil to water ratio = 1:2.5).

Statistical analysis

The statistical parameters including maximum, minimum, mean, standard deviation (SD), and coefficient of variation (CV) for each trait at each longitude site were estimated by using EXCEL. Nested variation analysis was used to calculate the significance among populations and within populations for 10 morphological traits of wild bermudagrass, which were collected from 13 different sites. Within-population and among-population variance components were calculated to determine the percentage of variability owing to the longitude and the local micro-environmental changes (SPSS 22.0). We explored the relationships between morphological traits, longitude, soil, and climate by Pearson’s correlation analysis and presented resulted by HemI software. The 13 different longitude sites were clustered based on the 10 morphological traits using cluster analysis (CA). Principle component analysis (PCA) of the 10 morphological traits was used to evaluate the diversity of morphological traits in bermudagrass. Figures were drawn with Sigmaplot 10.0 (Systat Software Inc.) and the R package.


Descriptive statistics of morphological traits and soil nutrients along a longitudinal gradient.

The maximum, minimum, mean, SD, and CV of 10 morphological traits at each site are given in Table 2. The CVs for the 10 traits ranged from 13.78% (LWES) to 34.13% (LLES). The CV values for LLES, ILES, LLS, ILS, PHT, and RBH were > 20%, meaning that there were larger morphological variations in those traits along the longitudinal gradient. The CV values for LWES, DES, LWS, and DS were < 20% but > 10%, which means that there were smaller morphological variations in these traits along the longitudinal gradient, and the CV values for the length and height characters were greater than for the width characters. Observable morphological variation existed in bermudagrass along the longitudinal gradient, and the 10 quantitative traits showed moderate and higher morphological variation.

Table 2 Average values for morphological traits of bermudagrass germplasm collected at 13 different longitudinal sites along latitude 34oN

The mean values for the soil nutrient contents at the different longitude sites are shown in Table 3. There was no obvious distribution pattern for soil nutrients along the longitudinal gradient. The 0soil nutrient CVs were very high, and all indexes were > 27% along the longitudinal gradient. The CV values were in the following relative order: AP > exchangeable Mg > AS>exchangeable Na > AK > exchangeable Ca > SOM > TN > AN>pH. In addition, the variation in pH had almost no effect compared to the magnitude of the environmental changes along the longitudinal gradient.

Table 3 Soil nutrient contents and pH at 13 collection sites for wild bermudagrass at different longitudes along latitude 34oN

Variation in morphological traits among and within populations.

The term “population” means that all individuals of the same species occupy a certain space in a certain period of time. We considered the samples collected in each site as representative of 13 populations in this study. High and significant levels of variation were found among and within populations for most measured traits (P < 0.01), and only LWS and DS showed no significant differences within populations (Table 4). The values of the variance components for all morphological traits among populations were > 15% and the range of changes was from 15.66% (LLES) to 51.53% (PHT). The morphological trait variation was mainly controlled by longitude and environmental factors among the populations. LWS had the lowest variance component value (5.26%), while DS had the highest variance component (49.72%) within populations. The variance components of the other traits ranged from 10 to 23% within populations, which indicated that the traits with lower variance component values within populations showed less genetic variation. Results from nested variation analysis showed that the variance components of morphological traits among individuals within populations were generally lower than those among populations.

Table 4 ANOVA for 10 morphological traits among and within 13 populations of wild bermudagrass

Correlation analysis

There were significant correlations between the tested traits (Fig. 1). The results showed that leaf width traits were positively and significantly correlated with diameter traits in most cases. LLES and ILES were significantly correlated with LLS and ILS, separately. In addition, the height traits (PHT and RBH) showed positive and significant correlations with ILES and LLS. The results of correlation analyses between longitude and 10 morphological traits of bermudagrass showed that there were significant correlations for some morphological traits such as LLES (P < 0.05), LWES (P < 0.05), ILES (P < 0.01), and ILS (P < 0.01), and all tested traits were negatively correlated with longitude (Fig. 2). Thus, the larger morphological traits were associated with lower longitudes; the leaves on the erect shoots were wider, and internodes were longer in bermudagrass with decreasing longitude. LLES (P < 0.05), ILES (P < 0.01), and ILS (P < 0.01) were significantly and negatively correlated with MAR, and ILES (P < 0.05) was significantly and negatively correlated with MAT, which implied that bermudagrass has longer leaves and internodes in areas with colder temperatures with less rainfall.

Fig. 1

Heat map for the correlation coefficients between morphological traits of bermudagrass

Fig. 2

Heat map for the correlation of the soil nutrients, longitude, climate and morphological traits

ILES (P < 0.05) and ILS (P < 0.01) displayed positive and significant correlations with exchangeable Ca, and ILS (P < 0.05) was significantly and positively correlated with pH, which indicates that longer internodes are positively associated with higher contents of exchangeable Ca and higher pH. ILES (P < 0.01), PHT (P < 0.01), and RBH (P < 0.01) were significantly and negatively correlated with AS, which shows that plants of bermudagrass were taller with increasing AS content. In addition, the leaf width of the stolon and stolen diameter did not show significant correlations with any geo-climatic or soil parameters.

Cluster analysis and principal component analysis of morphological traits

Figure 3 shows a dendrogram based on the cluster analysis of the morphological distance matrix for the 13 different longitudes. These different regions were divided into three clusters: the morphological traits included in clusters A, B, and C were in accordance with the approximate longitudes. Also, the bermudagrass populations were divided into many sub-populations by the complex environmental and geographical factors, which were observed at both low and high longitudes. Sites 1, 2, 3, 4, 7, and 9 grouped together in cluster A. This cluster contained plants with the minimum values of the morphological traits; plants were short with and narrow leaves and stems, except for plant height, which belonged to the dwarf and thin type. Sites 8 and 11 were included in cluster B. The bermudagrass plants in cluster B differed from the others by their maximum leaf width and diameter characteristics, which had the best creeping characters and belonged to the type of middle and wide. Cluster C consists of populations 5, 6, 10, 12, and 13 that had best upright characters and belonged to the tall and stocky plant type, because the plants in Cluster C had the maximum values for LLES, LLES, ILES, ILS, PHT, and RBH.

Fig. 3

Cluster analysis of morphological traits

To estimate the relative significance of the evaluated traits and the longitudinal variation, we performed principal component analysis based on the 10 morphological traits. The PCA results showed that the first two PC axes (PC1 and PC2) explained 80.5% of the total multivariate variation (Table 5). The relatively large load in PC1 are LWES, ILES, LLES and ILS, and when the relatively large load in PC2 were PHT, DS, DES and RBH. The contribution of PC1 (60.03%) was much greater than that of PC2 (20.47%). Along the longitude gradient, morphological diversity in bermudagrass was mainly due to LWES, ILES, LLES, and ILS in PC1. In order to investigate the relationships between longitude and morphological traits and to determine the morphological variability, we calculated the Euclidean distances. The distribution of sites along the first two principle component axes partition similar groups, which confirmed the results of the cluster analysis (Fig. 4).

Table 5 Component matrix for principal component analysis (PCA) of the 10 morphological traits in bermudagrass
Fig. 4

The distribution of bermudagrass populations on the PC1 and PC2 by using principal component


The relationships between morphological traits and longitude

Previous studies have shown that morphological variations are significantly affected by geographical locations [33,34,35,36,37]. Morphological variation as a response to environmental changes at different geographic positions is common in C. dactylon [38, 39]. A study of 260 accessions of bermudagrass germplasm indicated that there is extensive morphological variation in C. dactylon populations along longitudinal gradients, which has been found in some other studies [31, 40, 41]. The correlation analysis showed that the traits of width and dimeter were seldom significantly correlated with the longitude, soil nuriture or climate. It seems these traits may be less susceptible to environmental changes and suggests that the mechanism for stress tolerance is found in bermudagrass, thereby reducing the change of irreversibly responding to varying environmental factors. In the contrast, most of other traits were significant correlated with each other and also significantly related with enviroment factors. These traits were significantly effected by environment factors including longitude variation so that they had high values in CV along the changes of longtitude. Leaf traits of the erect shoot and internode length were significantly affected by longitude, and these traits were also the main source of diversity in bermudagrass populations along a longitudinal gradient, which shows that longitude is an important factor that influences morphological variation in bermudagrass. In general, we observed a trend in which the morphological traits increased in size with decreasing longitude, suggesting that the larger sizes are an adaptation to drought, lower temperatures, and shorter growing seasons as well as longer and deeper periods of dormancy, enabling the plants to survive adverse conditions and to increase their competitive abilities. In many cases, plants are generally smaller and show reduced leaf trait sizes in order to reduce plant transpiration and maintain efficient utilization of the limited water resources under drought conditions [42,43,44]. Bermudagrass is a widespread species that is adaptable to adverse environments that shows a contrary strategy in morphological traits to adapt to dry and low temperature conditions. However, high temperatures did not cause increases in plant size along a latitudinal gradient in another study [31].

We have identified significant variation in morphological traits in bermudagrass among and within populations, which indicates that the different environments over a large scale had a remarkable influence on the morphological traits (Table 4). Variation in components of morphological traits within populations was lower than that among general populations, which means that longitude is the major factor influencing morphological variation in bermudagrass. Some studies have shown that genetic variation is the main determinant for adaptability of a species within and between populations [45,46,47,48]. Morphological variation is commonly influenced by genetic variation, environmental variation, or the interactions between them [4, 49]. Morphological variation in bermudagrass is mainly determined by the influence of longitude in this study, which can enhance the spreading to a large scale likewise.

Effect of climatic and soil nutrients on morphological traits in bermudagrass

Previous studies in species such as rice and Sarracenia purpurea (pitcher plant) have investigated the spatial patterns of morphological trait variation along geographical gradients, and showed that the geographical patterns are shaped by edaphic and climatic factors [50, 51]. Morphological variations are well documented in many plant species such as Afzelia africana (African mahogany) and black plum, which are affected by climatic conditions [52,53,54]. Fruits and leaflets of A. africana show significant morphological variation under different climatic conditions [52]. Bejiga et al. (1996) discovered that climatic factors are significantly correlated with seed weight [55]. In our study, longer internodes were found in those areas with lower MAR and MAT. Previous studies showed that internode length in bermudagrass is significantly and negatively related to rainfall and temperature [31]. In areas with low rainfall and low temperatures, one adaptation mechanism for drought could be an increase in the internode length that would allow the leaf to more effectively absorb resources.

Soil nutrient factors are the most important components of the soil ecosystem and are necessary elements for plant growth [56, 57]. Soil nutrients are affected by many factors, such as climate, topography, microorganisms, and artificial factors [58, 59]. Both exchangeable Ca and pH showed significant variation along the longitudinal gradient, but most soil nutrients did not show an obvious difference. Bermudagrass morphological traits such as height and internode length were significantly affected by AS, exchangeable Ca, and pH in our study.

The combinations of and interactions between environmental factors (soil and climate) and longitude may significantly influence morphological variation in bermudagrass [52, 60]. In this study, a large proportion of the internode lengths in regions 8 and 10 were longer than those in plants in the vicinity, possibly due to the fact that the MAR in these areas was much higher than in the neighboring regions. This suggests that morphological variation in internode length in bermudagrass would arise when the variation in precipitation is abrupt. The higher AS (available sulfur) content of region of 5, which was second only to region of 13, had maximum values for PHT and RBH. This indicates that mid-longitude sites may be accompanied by optimal temperature and precipitation, which could allow soil characteristics to limit and or alter plant productivity.

Impacts of longitude on possible applications, and adaptive evolution of C. dactylon populationsThe tested traits were separated into three groups by cluster analysis; the dwarf-thin type, the middle-wide type, and the tall-stocky type. The bermudagrass plants in cluster A belong to the dwarf-thin type and can be used for lawns. The bermudagrass plants in cluster B have the best creeping features and are ideal for use in soil conservation and slope protection. The plants in cluster C have the best upright characters and belong to the tall-stocky type, which can be used for livestock forage. The bermudagrass plants in cluster C are mainly concentrated in higher longitude sites with less rainfall, and using this type of bermudagrass in breeding can contribute to conservation and reasonable utilization of precious water resources in arid environments.

The level of morphological variation in bermudagrass was found to be high, which enhances its ability to adapt to different environments along a longitudinal gradient. Compared with the other regions, the morphological traits were obviously different in marginal regions and central areas. The morphological trait sizes were generally larger in the marginal and central longitudes than in the other sampling sites along the longitudinal gradient. This result is supported by the non-linear trend in morphological trait values with respect to longitude, and the larger sizes of the morphological traits will occur in areas with intermediate and edge longitudes. In Lianyungang, which was highest longitude site of all the selected locations, the sizes of the morphological traits were generally larger than in the neighboring areas. This phenomenon can be explained by the observation that a high level of environmental variability, especially for rainfall, is present in Lianyungang, and it is thought that the highest level of morphological plasticity would occur in peripheral regions rather than in central areas, because environmental variation is probably the highest there [61, 62]. Earlier studies in other species have shown that higher levels of variability would be evident in marginal populations of Acer platanooides and Picea abies rather than in central populations [62, 63]. We found that two characters, width and diameter, showed a high degree of stability and were rarely affected by longitude and environmental factors, and the maximum values for width and diameter in bermudagrass were found in the middle longitude site Sanmenxia. Higher levels of diversity exist in central populations compared to marginal populations [64]. We can explain the results of our study from the following two perspectives. First, two populations of a species may have expanded their territory during their long-term evolutionary history, and these two populations then form a contact zone at some point in time. If the two populations were genetically differentiated, the species would thus have a higher diversity in the contact area [65]. Second, it has been suggested that high levels of environmental stress in edge sites could reduce morphological plasticity compared to that in central areas [66].


The pronounced east-to-west longitudinal variations in morphological traits in bermudagrass show that longitude plays a major role in morphological variation. We observed that morphological trait sizes tended to increase with decreasing longitude in order for the plants to adapt to the changing environment, especially LLES, LWES, ILES, and ILS. However, the values for the morphological traits in plants from the edge and central sites were generally higher than in the other sampling locations along the longitudinal gradient because of the edge effects and the founder effects. Different combinations and interactions of environmental factors such as AS, exchangeable Ca, pH, MAT, and MAR at each site may obscure the general trends in trait changes in bermudagrass along the longitudinal gradient. Our research explored the evolutionary trends in morphological trait variation in wild bermudagrass populations along a longitudinal gradient, and provides abundant wild resources for breeding in bermudagrass. However, there are also shortcomings to this experiment; for example, we did not take into account that the content of soil nutrients could change with the changing of the seasons. Also, DNA molecular marker technology should be used to examine genetic variation in bermudagrass populations from China and other countries where C. dactylon grows wild.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Available kalium


Available nitrogen


Available phosphate


Available sulfur


Cluster analysis


Coefficient of variation


Diameter of erect shoot


Diameter of stolon


Internode length of erect shoot


Internode length of stolon


Leaf length of erect shoot


Leaf length of stolon


Leaf width of erect shoot


Leaf width of stolon


Mean annual rainfall


Mean annual temperature




Principal component analysis


Turf height


Reproductive branch height


Standard deviation


Soil organic matter


Total nitrogen


  1. 1.

    Bradshaw AD. Evolutionary significance of phenotypic plasticity in plants. Adv Genet. 1965;13(1):115–55.

    Article  Google Scholar 

  2. 2.

    Agrawal AA. Phenotypic plasticity in the interactions and evolution of species. Science. 2001;294:321–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Engel K, Tollrian R, Jeschke JM. Integrating biological invasions, climate change and phenotypic plasticity. Comm Integ Biol. 2011;4(3):247–50.

    Article  Google Scholar 

  4. 4.

    Rubio De Casas R, Vargas P, Pe’rez-Corona E, Manrique E, Quintana J, Garcı’a-Verdugo C, Balaguer L. Field patterns of leaf plasticity in adults of the long-lived evergreen Quercus coccifera. Ann Bot. 2007;100(2):325–34.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Milne JM, Murphy KJ, Thomaz SM. Morphological variation in Eichhornia azurea (Kunth) and Eichhornia crassipes (Mart.) Solms in relation to aquatic vegetation type and the environment in the floodplain of the Rio Paraná, Brazil. Hydrobiologia. 2006;570(1):19–25.

    Article  Google Scholar 

  6. 6.

    Cordell S, Goldstein G, Mueller-Dombois D, Webb D, Vitousek PM. Physiological and morphological variation in Metrosideros polymorpha, a dominant Hawaiian tree species, along an altitudinal gradient: the role of phenotypic plasticity. Oecologia. 1998;113(2):188–96.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Aitken SN, Yeaman S, Holliday JA, Wang T, Curtis-McLane S. Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol Appl. 2008;1(1):95–111.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Hoffmann AA, Sgrò CM. Climate change and evolutionary adaptation. Nature. 2011;470(7335):479–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Alberto FJ, Aitken SN, Alía R, Gonzáonzrto FJ, Z SC, Hänninen H, Kremer A, Lefèvre F, Lenormand T, Yeaman S, Whetten R, Savolainen O. Potential for evolutionary responses to climate change- evidence from tree populations. Glob Change Biol. 2013;19(6):1645–61.

    Article  Google Scholar 

  10. 10.

    Franks SJ, Weber JJ, Aitken SN. Evolutionary and plastic responses to climate change in terrestrial plant populations. Evol Appl. 2014;7(1):123–39.

    Article  Google Scholar 

  11. 11.

    Merilä J, Hendry AP. Climate change, adaptation, and phenotypic plasticity: the problem and the evidence. Evol Appl. 2014;7(1):1–14.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Beard JB. Turfgrass: science and culture. Englewood Cliffs: Prentice-Hall; 1973. p. 658

    Google Scholar 

  13. 13.

    Taliaferro CM. Bermudagrass (Cynodon (L.) Rich). In: Casler MD, Duncan RR, editors. Turfgrass biology, genetics and breeding. New Jersey: Wiley; 2003. p. 235–257.

  14. 14.

    Harlan JR. Sources of variation in Cynodon dactylon (L.)Pers. Crop Sci. 1969;9(6):774–8.

    Article  Google Scholar 

  15. 15.

    Holdaway RJ, Richardson SJ, Dickie IA, Peltzer DA, Coomes DA. Species- and community-level patterns in fine root traits along a 120000-year soil chronosequence in temperate rain forest. J Ecol. 2011;99(4):954–63.

    Article  Google Scholar 

  16. 16.

    Lambers H, Brundrett MC, Raven JA, Hopper SD. Plant mineral nutrition in ancient landscapes: high plant species diversity on infertile soils is linked to functional diversity for nutritional strategies. Plant Soil. 2010;334(1–2):11–31.

    CAS  Article  Google Scholar 

  17. 17.

    Overgaard J, Kristensen TN, Mitchell KA, Hoffmann AA. Thermal tolerance in widespread and tropical drosophila species: does phenotypic plasticity increase with latitude? Am Nat. 2011;178(1):80–96.

    Article  Google Scholar 

  18. 18.

    Ayana A, Bekele E. Geographical patterns of morphological variation in sorghum (Sorghum bicolor (L.) Moench) germplasm from Ethiopia and Eritrea: quantitative characters. Euphytica. 2000;115(2):91–104.

    Article  Google Scholar 

  19. 19.

    Li CY, Zhang XJ, Liu XL, Luukkanen O, Berninger F. Leaf morphological and physiological responses of quercus aquifolioides along an altitudinal gradient. Silva Fennica. 2006;40(1):5–13.

    Article  Google Scholar 

  20. 20.

    Chang SW, Hsiao JY. Altitudinal genetic differentiation and diversity of Taiwan Lily (Lilium longiflorum var. formosanum; Liliaceae) using RAPD markers and morphological characters. Int J Plant Sci. 2001;162(2):287–295.

  21. 21.

    Janzen DH. Why mountain passes are higher in the tropics. Am Nat. 1967;101(919):233–49.

    Article  Google Scholar 

  22. 22.

    Stevens GC. The latitudinal gradient in geographical range: how so many species coexist in the tropics. Am Nat. 1989;133(2):240–56.

    Article  Google Scholar 

  23. 23.

    El ZA, Benhassaini H. Morphological leaf variability in natural populations of Pistacia atlantica Desf. Subsp. atlantica along climatic gradient: new features to update Pistacia atlantica subsp. atlantica key. Int J Biometeorol. 2016;60(4):577–89.

    Article  Google Scholar 

  24. 24.

    Lü XM, Zhou G, Wang Y, Song XL. Effects of changing precipitation and warming on functional traits of zonal Stipaplants from inner Mongolian grassland. Acta Meteorologica Sinica. 2016;30(3):412–25.

    Google Scholar 

  25. 25.

    Shi XY. Regional climate meridional and zonal feature in China. China: Chinese academy of meteorological science; 2004.

    Google Scholar 

  26. 26.

    Gatschet MJ, Taliaferro CM, Anderson J, Porter DR, Anderson MP. Cold acclimation and alterations in protein synthesis in bermudagrass crowns. J Am Soc Hortic Sci. 1994;119(3):477–80.

    CAS  Article  Google Scholar 

  27. 27.

    Taliaferro CM. Diversity and vulnerability of Bermuda turf-grass species. Crop Sci. 1995;35(2):327–31.

    Article  Google Scholar 

  28. 28.

    Rochecouste E. Studies on the biotypes of Cynodon dactylon botanical investigation. Weed Res. 1962;2:1–23.

    Article  Google Scholar 

  29. 29.

    Harlan JR, JMJD W, Huffine WW, Deakin JR. A guide to the species of Cynodon (Gramineae). Bull Okla Agric Exp Stn Stillwater. 1970;B673:1–37.

  30. 30.

    Potenza G, Fascetti S, Castronuovo D, Lovelli S, Perniola M, Viggiani R, Rossi R, Marchione V, Candido V. Collection and preliminary characterisation of native turfgrass accessions of Cynodon dactylon L. in the Mediterranean area. J Food Agric Environ. 2014;1212(2):770–4.

    Google Scholar 

  31. 31.

    Zhang JX, Wang ML, Guo ZP, Guan YZ, Guo YX, Yan XB. Variations in morphological traits of bermudagrass and relationship with soil and climate along latitudinal gradients. Hereditas. 2018;155(1):31.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Lu RK. Chemical analysis method of agricultural soil. China: China agricultural science press; 2000. p. 106–7.

    Google Scholar 

  33. 33.

    Soolanayakanahally RY, Guy RD, Silim SN, Song M. Timing of photoperiodic competency causes phenological mismatch in balsam poplar (Populus balsamifera L.). Plant Cell Environ. 2013;36(1):116–27.

    PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Cooke JEK, Eriksson ME, Junttila O. The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms. Plant Cell and Environ. 2012;35(10):1707–28.

    CAS  Article  Google Scholar 

  35. 35.

    Anway LGHC. A morphological and chemical analysis of geographical variation in Tilia L. of eastern North America. Brittonia. 1972;24(1):2–8.

    Article  Google Scholar 

  36. 36.

    Naïma AA, Asbahani E, El Alem Y, El Finti A, Msanda F, El Mousadik A. Variation in oil content and fatty acid composition in preselected argan trees with morphological characters and geographical localization. Med J Nutrition Metab. 2013;6(3):217–25.

    Article  Google Scholar 

  37. 37.

    Antigoni K, Catarina P, Martínez-Freiría F. Where does diversity come from? Linking geographical patterns of morphological, genetic, and environmental variation in wall lizards. BMC Evol Biol. 2018;18(1):124.

    Article  Google Scholar 

  38. 38.

    Wofford DS. Heritability estimates for turfgrass characteristics in bermudagrass. Crop Sci. 1985;25(1):133–6.

    Article  Google Scholar 

  39. 39.

    Liu JX, Guo AG, Guo HL. Morphological variation and types of Cynodon dactylon. Acta Pratacultural Sci. 2003;12(6):99–104.

    CAS  Google Scholar 

  40. 40.

    Reasor EH, Brosnan JT, Staton ME, Lane T, Trigiano RN, Wadl PA, Conner JA, Schwartz BM. Genotypic and phenotypic evaluation of off-type grasses in hybrid Bermudagrass [Cynodon dactylon (L.) Pers.×C. transvaalensis Burtt-Davy] putting greens using genotyping-by-sequencing and morphological characterization. Hereditas. 2018;155(1):8.

    Article  Google Scholar 

  41. 41.

    Akbari M, Salehi H, Niazi A. Evaluation of diversity based on morphological variabilities and ISSR molecular markers in Iranian Cynodon dactylon (L.) Pers. accessions to select and introduce cold-tolerant genotypes. Mol Biotechnol. 2018;60(4):259–70.

    CAS  Article  Google Scholar 

  42. 42.

    Wang RZ, Gao Q. Photosynthesis, transpiration and water use efficiency in two divergent Leymus chinensis populations from North-Eastern China. Photosythetica. 2001;39(1):123–6.

    Article  Google Scholar 

  43. 43.

    Wang RZ, Ripley EA, Zu YG, Nie SQ. Demography of reproductive and biomass allocation of grassland and dune Leymus chinensis on the Songnen plain, North-Eastern China. J Arid Environ. 2001;49(2):289–99.

    Article  Google Scholar 

  44. 44.

    Wang R, Gao Q. Morphological responses of Leymus chinensis (Poaceae) to the large-scale climatic gradient along the north-East China transect (NECT). Divers Distrib. 2004;10(1):65–73.

    Article  Google Scholar 

  45. 45.

    Kongkiatngam P, Waterway MJ, Fortin MG, Coulman BE. Genetic variation within and between two cultivars of red clover (Trifolium pratense L.): comparison of morphological, isozyme, and RAPD markers. Euphytica. 1995;84(3):355–61.

    Article  Google Scholar 

  46. 46.

    Julier B, Huyghe C, Ecalle C. Within- and among cultivar genetic variation in alfalfa: forage quality, morphology, and yield. Crop Sci. 2000;40:365–9.

    Article  Google Scholar 

  47. 47.

    Bolanos-Aguilar ED, Huyghe CH, Julier B, Ecalle CH. Genetic variation for seed yield and its components in alfalfa (Medicago sativa L.) populations. Agronomie. 2002;20(3):333–46.

    Article  Google Scholar 

  48. 48.

    Rawat K, Bakshi M. Provenance variation in cone, seed and seedling characteristics in natural populations of Pinus wallichiana a.B. jacks (blue pine) in India. Ann For Res. 2011;54(1):39–55.

    Google Scholar 

  49. 49.

    Rehfeldt GE. A model of genetic variation for Pinus ponderosa in the inland northwest (USA): applications in gene resource management. Can J For Res. 1991;21(10):1491–500.

    Article  Google Scholar 

  50. 50.

    Ellison AM, Buckley HL, Miller TE, Gotelli NJ. Morphological variation in Sarracenia purpurea (Sarraceniaceae): geographic, environmental, and taxonomic correlates. Am J Bot. 2004;91(11):1930–5.

    Article  Google Scholar 

  51. 51.

    Kondo M, Pablico PP, Aragones DV, Agbisit R, Abe J, Morita S, Courtois B. Genotypic and environmental variations in root morphology in rice genotypes under upland field conditions. Plant Soil. 2003;255(1):189–200.

    CAS  Article  Google Scholar 

  52. 52.

    Houehanou TD, Prinz K, Hellwig F, Assogbadjo AE, Gebauer J, Glele Kakaï RL, Sinsin B. Morphological trait variation and relationships of Afzelia africana Sm. Caused by climatic conditions and anthropogenic disturbance in Benin (West Africa). Genet Resour Crop Evol. 2019;66(5):1091–105.

    CAS  Article  Google Scholar 

  53. 53.

    Hounkpevi A, Azihou AF, Kouassi EK, Porembski S, Glèlè KR. Climate-induced morphological variation of black plum (Vitex doniana Sw.) in Benin, West Africa. Genet Resour Crop Evol. 2016;63(3):1073–84.

    Article  Google Scholar 

  54. 54.

    Loidi J. Plant eco-morphological traits as adaptations to environmental conditions: some comparisons between different biomes across the world. Geographical changes in vegetation and plant functional types. Cham: Springer; 2018. p. 59–71.

    Google Scholar 

  55. 55.

    Bejiga G, Tsegaye S, Tullu A, Erskine W. Quantitave evaluation of Ethiopian landraces of lentil (Lens culinaris). Genet Resour Crop Evol. 1996;43(4):293–301.

    Article  Google Scholar 

  56. 56.

    Robertson GP, Vitousek PM. Nitrification potentials in primary and secondary succession. Ecology. 1981;62(2):376–86.

    Article  Google Scholar 

  57. 57.

    Vitousek PM, Matson PA, Cleve KV. Nitrogen availability and nitrification during succession: primary, secondary, and old-field seres. Plant Soil. 1989;115(2):229–39.

    Article  Google Scholar 

  58. 58.

    Jian N. Carbon storage in terrestrial ecosystems of China: estimates at different spatial resolutions and their responses to climate change. Clim Chang. 2001;49(3):339–58.

    Article  Google Scholar 

  59. 59.

    Post WM, Izaurralde RC, Mann LK, Bliss N. Monitoring and verifying changes of organic carbon in soil. Clim Chang. 2001;51(1):73–99.

    Article  Google Scholar 

  60. 60.

    Okello J, Okullo JBL, Eilu G, Nyeko P, Obua J. Morphological variations in Tamarindus indica LINN. Fruits and seed traits in the different agroecological zones of Uganda. Int J Ecol. 2018;99:9–21.

    Google Scholar 

  61. 61.

    Volis S, Mendlinger S, Olsvig-Whittaker L, Safriel UN, Orlovsky N. Phenotypic variation and stress resistance in core and peripheral populations of Hordeum spontaneum. Biodivers Conserv. 1998;7(6):799–813.

    Article  Google Scholar 

  62. 62.

    Rusanen M, Vakkari P, Blom A. Genetic structure of Acer platanoides and Betula pendula in northern Europe. Can J For Res. 2003;33(6):1110–5.

    Article  Google Scholar 

  63. 63.

    Muona O, Paule L, Szmidt AE, Kärkkäinen K. Mating system analysis in a central and northern European population of picea abies. Scand J Forest Res. 1990;5(1–4):97–102.

    Article  Google Scholar 

  64. 64.

    Barrett SCH, Husband BC. The genetics of plant migration and colonization. In: Brown AHD, Clegg MT, Kahler AL, Weir BS, editors. Plant population genetics, breeding and genetic resources. Sunderland: Sinauer Associates; 1989. p. 254–77.

    Google Scholar 

  65. 65.

    Walter R, Epperson BK. Geographic pattern of genetic diversity in Pinus resinosa: contact zone between descendants of glacial refugia. Am J Bot. 2005;92:92–100.

    Article  Google Scholar 

  66. 66.

    Mägi M, Semchenko M, Kalamees R, Zobel K. Limited phenotypic plasticity in range-edge populations: a comparison of co-occurring populations of two Agrimonia species with different geographical distributions. Plant Biol. 2011;13(1):177–84.

    Article  Google Scholar 

Download references


I would like to thank Ming-Zhang for their help in sample collection.


No funding was received.

Author information




ML analyzed and interpreted the patient data and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yuxia Guo or Xuebing Yan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, M., Zhang, J., Guo, Z. et al. Morphological variation in Cynodon dactylon (L.) Pers., and its relationship with the environment along a longitudinal gradient. Hereditas 157, 4 (2020).

Download citation


  • Cynodon dactylon
  • Wild beomudagrass
  • Longitude
  • Phenotypic plasticity
  • Soil nutrients
  • Climate