Supplement 1. Variance partitioning analysis. Breaking force (F brk) is expressed as the product of width (w), thickness (t), and tensile strength (σ brk) (eq.1): (eq.1) F brk = w • t • σ brk After logarithm transformation, eq.1 turns additive (eq.2): (eq.2) log 10 (F brk) = log 10 (w) + log 10 (t) + log 10 (σ brk) → F brk ' = w' + t' + σ brk ' The variance (Var) of the breaking force can be then expressed as the sum of the covariances (Cov) between each term and the mechanical resistance (eq.3): (eq.3) Var(F brk ') = Cov(F brk ',w') + Cov(F brk ',t') + Cov(F brk ',σ brk ') The contribution (Cont) of each component to F brk ' is expressed as its covariance divided by the variance of F brk ' (eqs.4): (eqs.4) Cont(w) = Cov(F brk ',w')/Var(F brk ') Cont(t) = Cov(F brk ',t')/Var(F brk ') Cont(σ brk) = Cov(F brk ',σ brk ') /Var(F brk ') Cont(w) + Cont(t) + Cont(σ brk) = 1 Therefore, the contribution of each component is translated as relative importance in mechanical resistance. 2 Supplement 2. Tables Table S1. Sampling sites by bioregion, with geographical coordinates, species collected, date(s) of collection (mm/yy), and habitat details (depth in m). Species ID: Amphibolis antarctica (Aa), A. griffithii (Ag), Cymodocea nodosa (Cn), C. rotundata (Cr), C. serrulata (Cs), Enhalus acoroides (Ea), Halodule pinifolia (Hp), H. uninervis (Hu), Halophila ovalis (Ho), H. ovalis ssp. bullosa (Hb), Posidonia australis (Pa), P. coriacea (Pc), P. oceanica (Po), P. sinuosa (Ps), Ruppia maritima (Rm), Syringodium isoetifolium (Si), Thalassia hemprichii (Th), Thalassodendron pachyrhizum (Tp), Zostera capricorni (Zc), Z. marina (Zm), Z. nigricaulis (Zg), and Z. noltei (Zn). Number of samples per species is included between brackets next to each species code at each site. Sites Geographical coordinates Species ID Date(s) of collection Habitat details TIP: Tropical Indo-Pacific region (8 species, 5 sites, 21 species-at-site sets)