Schematic representation of the morphological traits measured in the samples of Erica andevalensis and E. mackayana. Leaf traits (a) (LF = Leaf length; LW = Leaf width; LP = Petiole length), variation on leaf morphology (b) (LR = Leaf rolling degree), and indumentum (c) (LSG = length of glandular hairs; LNR = presence of non-glandular hairs). Stamen traits (d) (AL = Anther length; AO = Anther pore; AAL = Anther appendix length; AAC = Anther appendix curvature; AK = Anther knob), flower morphology (e) (FL = Corolla length; FW = Corolla width; FWC = Corolla opening; FS = Style exertion; FHP = Pedicel hairiness, FLP = Pedicel length), and arrangement of a flowering branch (f) (FI = Flowers per inflorescence; LIA = Leaf insertion angle, L5 = 5nate leaves, LVD = Density of whorls). The ovary (g) (OL = Ovary length; OW = Ovary width; ON = Nectary size) and detail of ovary surface (h) (OH = Ovary hairiness) and stem (i) (BLH = Branch longest hair). Coding, measures, and additional information in Suppl. material 1: table S2.