Abstract Trait integration arises through both selection on functional coordination and shared developmental pathways. Different anatomical components must both work well and develop together to generate individuals with the appropriate physiology to survive and reproduce in their environment. In this study, we used a common garden experiment and Bayesian multilevel models to test whether stomatal anatomy coordinates leaf gas exchange, Rubisco kinetics, and leaf size across 10 closely related species of Limonium from the Balearic Islands. The results indicate that the anatomical determinants of maximum stomatal conductance, stomatal density and size, were functionally coordinated with Rubisco kinetics – species whose stomatal anatomy was correlated with low stomatal conductance have evolved Rubisco enzymes better adapted to low operational chloroplastic CO2 concentrations. Lower stomatal density was associated with greater leaf size, which can be explained by a greater proportion of pavement cells in large-leaved species. These results suggest that both selection for functional coordination (stomata and Rubisco kinetics) and shared development pathways (stomatal density and leaf area) likely shape patterns of trait integration between species.