We used linear mixed models (R library lme4, function lmer) to examine the effects of interactions and to estimate the variance of subtype richness associated with region and selected the best model based on Bayesian Information Criterion (Supplementary Table 3).