Through its natural or cultivated insular population distribution, Ilex aquifolium L. is a paramount species which is exceptionally suitable for studying phenotypic variability and plasticity through the assessment of morphological, physiological, biochemical and genomic features with respect to acclimation and/or adaptation efficiency. The current study is focused on four insular populations of Ilex aquifolium from Eastern Europe (i.e., in Romania, Hungary, Serbia and Bulgaria), and presents an initial evaluation of phenotypic variability in order to conclude our research on phylogenetic relationships and phytochemical profiles, including several descriptive and quantitative morphological traits. Taken together, the data from different methods in this paper indicate that the Bulgarian and Romanian populations can be distinguished from each other and from Serbian and Hungarian populations, while the latter show a higher level of resemblance with regards to their quantitative morphological traits. It is likely that these morphological traits are determined through some quantitative trait loci implicated in stress responses generated by light, temperature, soil water, soil fertility and salinity conditions that will need to be analysed in terms of their physiological, genomic and metabolomics traits in future studies.
Keywords: Ilex aquifolium L.; SE Europe; SEM; common holly; morphometry; phytochemicals; rbcL.