Local phenotypic adaptation
Common garden experiments were conducted to investigate the phenotypic divergence evidence for genome differentiation and adaptation of theA. viridiflora complex. The mean FST for the four lineages was 0.1519 (95% HPD: 0.14949-0.15074), and the overall QST of the six traits was higher than the mean FST, including the corolla diameter, petal length, spur length, pistil length, inflorescence number and leaf area, indicating that local adaptation was driven by these traits (Table 1). The highest QST was procured by measuring the corolla diameter, followed by the spur length, petal length, pistil length, inflorescence number and leaf area, which also showed a higher QST of floral characteristics than nutritional traits. Among the four lineages, CN showed the smallest corolla diameter and the shortest petal length, spur length and pistil length, followed by EL, while NW and NE showed the largest values of these phenotypes (the difference between NW and NE was not significant) (Figure S10A - D). Interestingly, the number of inflorescences showed the reverse order of the above traits (Figure S10E). EL showed the largest leaf area, followed by CN, while the difference in leaf area between NW and NE was not significant (Figure S10F).