Contribution of non‐native plants to the phylogenetic homogenization of U.S. yard floras

Cultivation and spread of non‐native plant species may result in either phylogenetic homogenization (increasing similarity) or differentiation (decreasing similarity) of urban floras. However, it is unknown how non‐native species influence homogenization of cultivated versus spontaneously occurring species in cities, and which traits are associated with species that promote homogenization versus differentiation. In this study, we compared homogenization effects of cultivated and spontaneous non‐native species in yard floras across and within seven widely distributed U.S. cities. Additionally, we explored which traits explained their particular contribution to homogenization. We recorded plant presence/absence in 178 private yards distributed among seven metropolitan statistical areas in the United States. We compared phylogenetic homogenization effects of non‐native species within both the cultivated and spontaneous species pools using phylogenetic dissimilarities and the homogenization index. Then, we expressed contributions of non‐native species to the homogenization of each pool as a function of two different sets of plant functional traits using phylogenetic generalized least square (PGLS) models across and within cities. Across cities, spontaneous non‐native species homogenized, and cultivated non‐native species differentiated, yard floras. Within the spontaneous pool, short, small‐seeded non‐native plants and non‐native grasses significantly homogenized yard floras. Within the cultivated pool, species contribution to homogenization was best predicted by plant height, presence of showy flowers, and growth form, with non‐native grasses significantly homogenizing cultivated yard floras. Within cities, non‐native species—whether they were cultivated or spontaneous—consistently homogenized yard floras of the three northern cities and differentiated yard floras of three of the four southern cities, suggesting that homogenization processes are context‐ and scale‐dependent. Likewise, traits explaining homogenization differed substantially among cities. The inconsistent patterns among cities in the plant traits that promoted homogenization of both cultivated and spontaneous species suggest that local environmental and anthropogenic conditions of individual cities imposed strong constraints on trait selection. Linking plant functional traits that promote homogenization with residents’ preferences for vegetation may further enhance understanding of how yard plant communities assemble at regional and local scales.


INTRODUCTION
In the United States, 75 million houses have a yard (US Census Bureau 2018). Together, these residential landscapes encompass a wide range of environmental and structural conditions. Yard floras typically contain both native and nonnative species, and spontaneously occurring and cultivated species (Smith et al. 2006, Marco et al. 2008, Bigirimana et al. 2012. Historical immigration from Europe and cultural similarities among U.S. regions have resulted in similar landscapes across widely dispersed cities (McKinney 2006, Ignatieva and Stewart 2009, Groffman 2014. However, the extent to which U.S. yard floras are homogeneous (i.e., taxonomically, phylogenetically, or functionally similar among regions) and the contributions of spontaneous and cultivated non-native species to the process of homogenization (i.e., promoting similarity) at different spatial scales are not clear.
Urban areas are frequently described as homogenized when compared to natural areas (i.e., urban areas in different regions are more similar to each other than the natural areas they replace; K€ uhn and Klotz 2006, McKinney 2006, Grimm et al. 2008, Lososov a et al. 2012, Groffman 2014. For example, previous research in the United States has demonstrated urban homogenization by showing that yard floras are compositionally and structurally more similar than the corresponding floras in surrounding natural areas (Wheeler 2017, Pearse 2018. However, biotic homogenization might also be assessed within and across cities, by comparing the relative similarities of different components of vegetation among yards or other urban habitats (Lososov a et al. 2012). In practice, biotic homogenization has not been tested using the latter approach across urban yards from a wide range of socio-ecological conditions. Species growing in yards can be considered as arising from two separate pools (Knapp et al. 2012, Groffman 2017): human-cultivated and a spontaneously occurring (establishing without assistance) species pool. Each of these pools is subject to distinct ecological assembly processes, such as habitat filtering and species interactions, as well as anthropogenic factors, that operate across different spatial scales (Williams 2009, Aronson et al. 2016, Groffman 2017. Across cities, the occurrence of spontaneous native species is mainly driven by biogeographical processes and habitat transformations that filter out species from the regional pool, while spontaneous non-native and cultivated species are mainly driven by climate, global trade, and policy regulations (Williams 2009, Kendal et al. 2012b. Within cities, both native and non-native spontaneous species are mainly driven by landuse identity and diversity, local conditions, socioeconomic and cultural factors, and species interactions, while cultivated species are mainly driven by socioeconomic and cultural factors alone (Cook et al. 2012, Aronson et al. 2016. As a result of all these interactions, a large proportion of yard floras (about 68-80%) are non-native (Marco et al. 2008, Bigirimana et al. 2012. Non-native species are one of the main drivers of ecosystem transformation (Vitousek et al. 1996), and the consequences of species introduction are currently a matter of study throughout the world (Gurevitch et al. 2011). Increasing trends in globalization are associated with the expansion of non-native species (Kuebbing et al. 2013). As a consequence, many ecosystems now harbor multiple introduced species, increasing the complexity of interspecific interactions. In urban environments, although it has been hypothesized that non-native species tend to homogenize urban plant communities (McKinney and Lockwood 1999), some studies suggest that non-natives can drive within-city floral differentiation (Olden andPoff 2003, La Sorte andMcKinney 2006). For example, Lososov a et al. (2012) found that biotic homogenization across 32 cities in central Europe depended on nonnative plants' residence time (the time a species had been present in the novel environment), with more recently introduced non-native species that had not achieved their potential range, increasing differentiation, and those that had had sufficient time to disperse into most suitable habitats increasing homogenization. Increasing numbers of non-native species in urban communities may lead to a decrease in phylogenetic diversity (Ceplov a et al. 2015) that might in turn reduce overall community functioning, stability, and resistance to environmental change by narrowing the available range of species-specific responses (Cadotte et al. 2012). On the other hand, non-native species also have the potential to increase phylogenetic diversity and species richness. Cultivated and spontaneous non-native species in yards are affected by different filtering and sorting processes (Groffman 2017, Pearse 2018) and may thus have contrasting effects on the composition of yard floras. A critical gap remains in our understanding of the influence of non-native species on community assembly in urban plant communities and their role in shaping composition and diversity.
Species traits influence how species interact with various filters in the process of community assembly (Williams et al. 2015, Aronson et al. 2016. Literature on biotic homogenization suggests that the common features of cities and human communities create conditions that promote species with similar traits and life histories (Groffman 2014, Williams et al. 2015. For spontaneous species, common functional traits of plants growing in human settlements include fast growth with a high demand for nutrients; wind pollination; small, readily dispersed seeds; and a tendency toward shorter lifespan than found in natural areas (Lososov a et al. 2006, Thompson and McCarthy 2008, Knapp et al. 2012, Williams et al. 2015. For cultivated species, predominant plant traits are heavily influenced by human choices based on the benefits they provide, ranging from aesthetic (flowers, color, texture) to functional (ground cover, shade provision, food provision; Kendal et al. 2012a, Larson 2016. In this regard, Pataki et al. (2013) refer to traits that are related to human attitudes and preferences as ecosystem service-based traits, to distinguish them from conventional categories of functional traits developed for natural ecosystems, and which are subject to natural selection and are linked to ecological processes. Nevertheless, few clear patterns have emerged from existing studies of plant traits in individual cities or urban habitats because of the complexity of plant-environment interactions in urban systems and limitations of existing data and methods (Williams et al. 2015). Therefore, determining which traits of non-native species make them successful in yards, and how their persistence and spread translates into biotic homogenization, may enhance our understanding of how yard plant communities assemble at regional and local scales.
In this study, we build on previous research on biotic homogenization of residential yards in the United States (Groffman 2014, Polsky 2014, Wheeler 2017, Pearse 2018 to quantify floristic compositional similarities across and within cities and determine what plant traits promote homogenization. Unlike previous studies, we specifically explored the extent to which nonnative species contribute to the homogenization/ differentiation of phylogenetic composition (using phylogenetic dissimilarities), thus incorporating phylogenetic information into traditional methods of studying communities like ordination that treat all species as independent and equally dissimilar. Moreover, we expand on previous studies that assessed homogenization of urban plant communities by reporting the particular contributions of cultivated versus spontaneously occurring non-native species to homogenization, and by relating species effects on homogenization to plant traits. In doing so, we specifically addressed the following questions: (1) Are cultivated and spontaneous non-native species similarly homogenizing urban yards across and within cities? and (2) What traits of these species explain their particular contribution to homogenization at different scales (across versus within cities)? Unlike previous studies that assessed biotic homogenization in urban areas using data from a variety of habitats and sampling schemes, we collected original data with the same sampling scheme for 178 private yards distributed among seven major metropolitan areas in the United States.
We hypothesized that non-native spontaneous (weedy) species would contribute to biotic homogenization because they have been in residence in urban areas for a longer time and are adapted to a broader range of biogeographical conditions than cultivated species (La Sorte et al. 2007, Lososov a et al. 2012, Ricotta et al. 2014. In contrast, we hypothesized that non-native cultivated species, which are more likely to be filtered by human preferences from commercial sources (Williams 2009), would have the opposite effect, increasing differentiation in floristic composition among yards. Across all yards and cities, we expected traits of non-native species to reflect the nature and magnitude of environmental and anthropogenic filters acting on each pool of species. Within cities, we did not expect a clear association between patterns of homogenization and non-native species or traits, given that local biophysical conditions, homeowners' management choices, and idiosyncratic aspects of cities (e.g., development history, demographics, cultural attributes, types of infrastructure, and connectivity) may have strong impacts on the dispersal and persistence of non-native species.

Yard selection
We selected seven major U.S. metropolitan statistical areas (hereafter "cities"): Boston, Massachusetts (BOS); Baltimore, Maryland (BAL); Los Angeles, California (LA); Miami, Florida (MIA); Minneapolis-St. Paul, Minnesota (MSP); Salt Lake City, Utah (SLC); and Phoenix, Arizona (PHX), to represent seven different ecological biomes and major climatic regions across the United States (Trammell et al. 2016). We stratified the cities by population density and socioeconomic status using the PRIZM marketing classification scheme (CLARITAS 2013, Polsky 2014, which classifies each Census Block Group in the United States into a single group based on analysis of the areal unit's population density, affluence, and life stage. Using this primary selection, we first contacted >100,000 households and identified >13,500 where the respondent was over 18 yr of age and their home had either a front or back yard (Polsky 2014). From this group, we conducted a telephone survey of 9480 residents, equally distributed among six of the seven cities (all but SLC), to identify individual residential properties as potential sampling sites (Polsky 2014). From the 9480 households surveyed, we randomly selected 21-30 residential properties per city for field sampling. In SLC, we randomly sent 50 letters to homeowners in each of six targeted PRIZM categories asking permission to visit their home for sampling. Within each PRIZM category, we visited five homes randomly selected from the positive respondents, for a total of 30 yards.

Vegetation data
We recorded plant species presence/absence in 178 yards across the seven cities (Data S1). We sampled all cities during the season of peak diversity (spring for LA and PHX; summer for BAL, BOS, MIA, MSP, and SLC). We sampled locations in BAL, BOS, MIA, MSP, and PHX in 2012, and those in LA and SLC in 2013. We surveyed the entire area of each yard. In areas of the yard with unmanaged vegetation or a woodland/woodlot, we sampled plant species presence/absence within a 2 m wide transect across the full yard or 100 m, whichever was shorter. Yard plants are often subspecies or cultivars, but we did not attempt to classify plants below the species level. For those plants that could not be identified at the species level (~15%), we only recorded the genus. We designated species as spontaneous or cultivated based on homeowner interviews, observations of placement, land use, and land-use history. For example, species in woodlots and unmanaged vegetation components were generally considered spontaneous. Any given species could be documented as both spontaneous and cultivated if different individuals of that species fell into different categories. Species that were not designated as either cultivated or spontaneous in the field were later classified in the laboratory based on records for the same species in other sampled yards. We matched species names to The Plant List (http:// www.theplantlist.org) version 1.1, using R package Taxonstand (Cayuela et al. 2017).

Native status and trait data
We classified species according to their status as either native or non-native for each U.S. state where the cities were located following the USDA PLANTS (http://www.plants.usda.gov) and EOL (http://www.eol.org) databases (see Appendix S1 in the Supporting Information for more details on native status classification). We assigned traits to the non-native species (1634 species) to examine whether their functional attributes could help explain why non-natives contribute to homogenization or differentiation. Our choice of traits to include depended on data availability, and different criteria were used for the cultivated and spontaneous pools. For example, for the cultivated pool we selected ecosystem service-based traits that have previously been shown to be associated with homeowners' needs and preferences: growth form, plant height, longevity, presence of showy flowers, edibility, and water use (Kendal et al. 2012a, Goodness et al. 2016, Avolio et al. 2018; Table 1). It has been suggested that traits related to overall ❖ www.esajournals.org plant size, such as plant height and larger growth forms (usually perennial), can influence aesthetic preferences. However, there is no wide consensus on which growth forms may be universally preferred, nor whether there is a preference for a single growth form. Edible plants grown in private yards are hypothesized to be similar across locations and thus to contribute to homogenization. The presence of showy (large, colorful, and conspicuous) flowers has been correlated with preferences for attractive plants over non-flowering vegetation in a variety of settings. Plant maintenance requirements largely depend on the local conditions where the plant is grown, and thus affect homeowners' preferences. In this study, we used water use to refer to the amount of water needed to achieve an acceptable level of plant health and aesthetics. We used the Water Use Classification of Landscape Species (WUCOLS) classification method and database developed by Costello and Jones (2014) to provide a trait for individual plant water needs. The WUCOLS method and database, to our knowledge, has not been previously used to understand community assembly in urban environments. Costello and Jones (2014) tabulated water requirement values (ks parameter) for more than 3000 ornamental species in six different areas of California. Species were classified as presenting very low requirements (VL; ks < 0.10), low requirements (L; 0.10 < ks < 0.30), moderate requirements (M; 0.40 < ks < 0.60), and high requirements (H; 0.70 < ks < 0.90). Average ks values across the six different areas of California were used in this study as a proxy for plant water requirements extended to cities included in our study.
For the spontaneous pool, we selected functional traits hypothesized to be selected for or against in urban environments, including traits that facilitate species' establishment and persistence, including specific leaf area (SLA), plant height, photosynthetic pathway, longevity, growth form, and seed mass (Weiher et al. 1999, Thompson and McCarthy 2008, Williams et al. 2015; Table 1). Specific leaf area correlates with species relative growth rates, palatability, and response to nutrient and moisture gradients, plant height relates to competitive ability (particularly for light) and resource allocation, and seed mass Notes: For continuous traits, the unit of measurement is provided in parentheses following the trait name. % Missing data refers to the percentage of species in the study lacking available trait information from accessible databases. Specific references for the TRY database can be found in Appendix S3.
† Reference category used in analysis is indicated with Ã . ‡ Cultivated/spontaneous. § Annual and biennial plants were combined in the same category. ¶ C 4 and CAM plants were combined in the same category.
relates to dispersal ability. Photosynthetic pathway reflects habitat affinity (with C 4 species potentially profiting from warmer temperatures due to urban heat island effects and being more frequent in yards in arid regions), whereas growth form and longevity correlate with extinction risk (with woody perennial vegetation less likely to be entirely extirpated in yards). We collected all functional traits from available databases: USDA PLANTS Database; BIEN database (Maitner 2017); TRY database (try-db.org; see Appendix S3 for specific references); the WUCOLS list (Costello and Jones 2014); Missouri Botanical Garden database (missouribotanicalgarden.org); Engemann et al. (2016); specialized literature; and expert knowledge (Table 1). We tested for correlations among traits using a Pearson correlation coefficient (r). Because no traits were correlated with r > 0.55, all traits were retained for analyses (Appendix S1: Table S2).
Data for continuous plant traits and water use were not available for some species (Table 1). Deleting species with missing data from the analysis would reduce the number of data points substantially (~45%), and probably bias the results because of the selective removal of species that were less well known. These cases were accommodated by estimating the missing values using phylogenetic information from species with available data. Some characteristics inherent to functional traits such as a strong phylogenetic trait signal and structural trade-offs between traits support statistical gap-filling of sparse trait matrices using this approach (Swenson 2014). For this, we first compared available trait data across four alternative evolutionary models (Early-Burst, Brownian motion, Ornstein-Uhlenbeck, and multivariate Ornstein-Uhlenbeck) using the R package Rphylopars (Goolsby et al. 2017). The best-fitting model on the basis of the lowest AIC value (i.e., Ornstein-Uhlenbeck) was then used to impute trait data (see Appendix S2 for more details on missing trait data imputation).
Although online databases that collate sampled trait data have provided a powerful tool for addressing mechanistic questions about multispecies assemblages McCarthy 2008, Duncan et al. 2011), they may not represent urban populations exposed to unique selection pressures (Alberti et al. 2017). To mitigate this concern, the majority of the traits we focused on for this study were categorical and, thus, unlikely to shift substantially across populations. Thus, we used the mean value in analysis for continuous traits where multiple values occurred for any given species.

Data analysis
A final set of 171 yards was included in the analysis after removing seven households with no non-native species in their yards. For the cultivated pool, only 167 yards were included in the analysis because four yards did not contain nonnative cultivated species.
We used an updated version of the Zanne (2013) phylogeny produced by Qian and Jin (2016) for all phylogenetic metrics, and we added species missing from this phylogeny at the genus level using the congeneric.merge function in the R package pez (Pearse et al. 2015). Hybrids were reduced to the genus level, and we excluded species for which there were no phylogenetic data (~1.5%) from the analysis.
To determine if a pair of yards had been homogenized or differentiated because of the introduction of non-native species, we calculated the homogenization (H) index (Qian and Ricklefs 2006, Qian et al. 2008, Lososov a et al. 2012; Box 1, Fig. 1). H index ranges from À1 to 1. Positive H index values indicate that non-native species contribute to phylogenetic differentiation. Negative H index values indicate that non-native species contribute to the phylogenetic homogenization of species assemblages. To calculate the H index, we measured pairwise phylogenetic dissimilarities in total and native species composition among yards using the evodiss_family function in the R package adiv (Pavoine 2017) and the formulas provided in Box 1. We further explored contributions of particular non-native species to the phylogenetic homogenization of yard floras (DJ T in Box 1) by measuring the change in the H index when each non-native species was removed from each pairwise comparison. We repeated the same procedure independently for the cultivated and spontaneous species pools across all yards and within the seven cities. To evaluate associations between DJ T and species relative frequency across all yards, we calculated Pearson correlation coefficients.
We used Cohen's d in R package lsr (Navarro 2015) as a measure of effect size when comparing the mean H index value of all yard pairwise comparisons to zero (no homogenization effect) across and within cities (Cohen's d around 0.20 = small effect, d around 0.50 = medium effect, and d around 0.80 = large effect; Cohen 1998). Unlike the t-test statistic, the effect size aims to estimate a population parameter and is not affected by the sample size. Also, we compared the same mean value of the H index for all yard pairwise comparisons to a null model to investigate whether observed homogenization could have arisen by chance. We interpreted departures from the null model as a signal that biotic homogenization was not a consequence of random stochastic processes, but rather of taxonomically selective anthropogenic factors. The null model was generated by a two-step process. First, 999 permutations of the species community matrix were created using the R package picante (Kembel et al. 2010) and the independent swap algorithm (Gotelli 2000), which maintains species occurrence frequency and sample species richness. Second, 999 random H index values were obtained from the randomized community matrices and their means used for the null model. We repeated this part of the analysis for each species pool across and within cities. We calculated P-values based on a two-tailed test.
Spatial autocorrelation in ecological data is a common phenomenon as observations at closer locations are usually more similar than would be expected by chance (Legendre 1993). To account for spatial autocorrelation between the H index of pairwise comparisons between yards across cities and geographical distance, we performed Mantel correlations. Geographical distances in kilometers between each location were calculated Box 1.

Homogenization index and species contribution to homogenization
Sites can be compared in a pairwise fashion according to their composition in evolutionary units (i.e., units of the branch lengths of a phylogenetic tree) exactly the same way they can be compared based on their composition in species (Pavoine 2016). Hence, traditional species dissimilarity indices may be adapted to account for phylogenetic dissimilarities. For example, the phylogenetic equivalent of the Jaccard index (Jaccard 1901) is evoD Jaccard (Pavoine 2016(Pavoine , 2017. evoD Jaccard ranges from 0 to 1, where 0 indicates that the pair of sites are identical in terms of evolutionary units and 1 indicates that sites have no evolutionary units in common (see Pavoine (2016) for details on the formula).
Consider two sites numbered 1 and 2 (Fig. 1A), containing five native (N i , in black) and four non-native (A i , in red) species between the two. evoD Jaccard (hereafter only J) can be calculated considering the whole pool of species (J total ) and only natives (J native ). These two dissimilarity indices can be further used to calculate the homogenization (H) index (Qian and Ricklefs 2006, Qian et al. 2008, Lososov a et al. 2012: H = J total -J native . H index ranges from À1 to 1. Positive H index values indicate that dissimilarity of native species is lower than dissimilarity of all species; therefore, non-native species contribute to phylogenetic differentiation. Negative H index values indicate that the introduction of non-native species contributes to the phylogenetic homogenization of species assemblages. H index thus reflects the combined effect of all non-native species on the biotic homogenization of the two sites. However, we can calculate contributions of particular non-native species to the phylogenetic homogenization of both sites by measuring the change in H index when each non-native species is removed from the pairwise comparison. To this end, we create additional pairwise site comparisons (Fig. 1B-E) where each non-native species is removed at a time (species A i in blue) and then recalculate H index for each pair of sites following the H index equation. The change in H index is therefore calculated as: DH Ai ¼ H T À H Ai , where DH Ai is the change in H index, H T is the H index of the original pairwise comparison (Fig. 1A), and H Ai is the H index of the newly created pairwise comparison without species A i (Fig. 1B-E). Given that J native is the same across all pairwise comparisons, then DH Ai ¼ DJ Ai , where DJ Ai ¼ J total À J Ai , and J Ai is the evoD Jaccard index considering all species except species A i .
Finally, DJ Ai can be summed up across all (n) pairwise site comparisons to obtain species A i total contribution to phylogenetic homogenization: DJ TAi ¼ P n j¼1 ðDJ Ai Þ j . In general, we refer to this term as DJ T .
using the great-circle method in the R package fields (Nychka et al. 2017). The significances of matrix correlations were tested using a Monte Carlo permutation test with 999 permutations, using the R package ade4 (Dray and Dufour 2007). Species cannot be assumed to be independent data points in regression models because of their phylogenetic relationship (Felsenstein 1985, Garland and Ives 2000, Ives and Zhu 2006. Therefore, we used the phylogenetic generalized least square (PGLS) method to test associations between DJ T and plant traits. The analysis was performed using the gls function in R package nlme (Pinheiro et al. 2018). PGLS explicitly incorporates the expected covariance among species into a statistical model fit by generalized least squares and allows for flexibility in the underlying evolutionary assumptions. We modeled DJ T for the cultivated and spontaneous non-native pools both across and within cities using a set of eight candidate models, including (1) a null model (only the intercept term), (2) six models testing each trait separately, and (3) the full model (all traits; Appendix S4: Tables S5, S6). We excluded growth form from the within-cities analysis to guarantee an appropriate ratio between the number of observations and predictor variables. In all models, we normalized continuous predictors between 0 and 1 before being entered into the model so that the coefficient values were comparable (Zuur et al. 2007). Plant height, SLA, and seed mass were log-transformed to reduce skewness and improve the Fig. 1. Diagram exemplifying calculations of species particular contribution to phylogenetic homogenization. First, the original pairwise site comparison containing five native (N i ) and four non-native (A i ) species is transformed into four new pairwise site comparisons (B, C, D, and E) where each non-native species is removed at a time (species A i in blue). For each pairwise site comparison, the evoD Jaccard (J) and H index are calculated. Finally, the particular contribution of species A i to phylogenetic homogenization is measured as DJ Ai . In the example, two non-native species are shared between the two sites (homogenizing), and two non-native species are only present in one of them (differentiating). Also, non-native species are separated from natives in the phylogeny. As a result, non-native species contribute to a different extent to phylogenetic homogenization. normality of the residuals. Following Anderson (2008), we calculated the AICc (the sample sizecorrected Akaike information criterion) values for each candidate model set and ranked them by their Akaike weights using the R package AICcmodavg (Mazerolle 2017). We used the same package to compute the standardized model-averaged coefficients.
We performed all statistical procedures in R version 3.4.1 (R Core Team 2017) and established significance at a < 0.05.

RESULTS
There were 2417 total cultivated and spontaneous species in all the yards sampled across the seven cities, of which 1634 (67.6%) were classified as non-native in the United States (Table 2, Fig. 2). In total, 1648 were classified as cultivated, of which 1310 (79.5%) were non-native. Among the spontaneous pool (1287 species), 653 (50.7%) were classified as non-native. Within cities, the number of pooled species varied widely, from 297 species in PHX to 964 in BOS ( Table 2  Notes: For each city, the number of sampled yards, the total number of sampled species, and the proportion of non-natives are shown. Note that total number of species in the study does not equal the sum of the number of species within each city because the same species are frequently found in different cities. Additional information on the biome, climatic region, mean annual temperature, and rainfall of each city can be found in Appendix S4: Table S1.  (Cohen's d = 0.94). These H indices were not different from their respective null model (Fig. 3). None of the Mantel tests comparing H index matrices and the geographical distance matrices were significant.
Among cultivated non-native species, those that contributed the most to the homogenization of yards across cities were Cynodon dactylon, Rosa hybrid, Hosta hybrid, Festuca arundinacea, and Hemerocallis hybrid, and those that contributed the most to differentiating yards were Buxus sempervirens, Berberis thunbergii, Vitis spp., Lamprocapnos spectabilis, and Cycas revoluta (Appendix S4: Tables S3, S4). Among spontaneous non-native species, those that contributed the most to the homogenization of yards across cities were Taraxacum campyloides, Acer platanoides, Cynodon dactylon, Poa Notes: Estimates whose 95% unconditional confidence intervals do not overlap zero are shown in bold. Positive estimates promote differentiation, while negative estimates promote homogenization. See Table 2 for corresponding reference categories in categorical variables. ALL, across all cities; BAL, Baltimore; BOS, Boston; LA, Los Angeles; MIA, Miami; MSP, Minneapolis-Saint Paul; PHX, Phoenix; and SLC, Salt Lake City. Table 4. Model-averaged estimates and unconditional SE (in parentheses) for each trait predicting non-native species contribution to homogenization (DJ T ) for the spontaneous pool within the seven cities. Notes: Estimates whose 95% unconditional confidence intervals do not overlap zero are shown in bold. Positive estimates promote differentiation, while negative estimates promote homogenization. We excluded growth form from the within-cities analysis to guarantee an appropriate ratio between the number of observations and predictor variables. See Table 2 for corresponding reference categories in categorical variables. ALL, across all cities; BAL, Baltimore; BOS, Boston; LA, Los Angeles; MIA, Miami; MSP, Minneapolis-Saint Paul; PHX, Phoenix; SLC, Salt Lake City; and SLA, specific leaf area.
pratensis, and Trifolium repens, and those that contributed the most to the differentiation of yards were Celastrus orbiculatus, Chelidonium majus, Commelina communis, Berberis thunbergii, and Rhamnus cathartica. Non-parametric correlations between species contributions to phylogenetic homogenization (DJ T ) and species relative frequency in yards were significant for both the cultivated (Pearson's r = À0.46; P < 0.001) and the spontaneous pool (Pearson's r = À0.60; P < 0.001).
Models within two AIC c points of the top model predicting DJ T for the cultivated nonnative species across cities included plant height, showy flowers, growth form, and the null model (Appendix S4: Tables S5, S6). For the spontaneous pool, only the model including plant height was ranked within two AIC c points. Nonnative grasses significantly contributed to the phylogenetic homogenization of the cultivated and spontaneous pools (model-averaged confidence intervals did not overlap zero; Table 3). Shorter and smaller-seeded spontaneous nonnative plants also significantly contributed to yard homogenization across cities (Table 4).
Within cities, traits of non-native species that significantly promoted homogenization among the cultivated pool included grass growth form (in LA and PHX) and showy flowers (LA) ( Table 3). For the spontaneous pool, traits of nonnatives that significantly differentiated yard floras within cities included seed mass (BAL and PHX), plant height (MIA and SLC), and perenniality (LA and SLC; Table 4). In contrast, traits of non-natives that promoted homogenization for the same pool included C 4 photosynthetic pathway (LA and PHX) and SLA (PHX).

Phylogenetic homogenization of U.S. yard floras
Understanding the extent to which non-native species contribute to the homogenization of biodiversity in urban areas is of great interest because globalization has contributed to the intentional and unintentional movement of species and because urbanization has increased globally in the last century and is projected to continue to increase. Here, we compared homogenization effects of cultivated and spontaneous non-native plant species on U.S. yard floras. As predicted, cultivated non-native species promoted phylogenetic differentiation of the cultivated yard floras across cities, while spontaneous non-native species promoted the phylogenetic homogenization of the spontaneous yard flora. Within cities, cultivated and spontaneous non-natives consistently promoted homogenization in the three northern cooler cities: BAL, BOS, and MSP.
Our findings further show that despite the significant effect of non-natives on the phylogenetic homogenization of yard floras at the continental scale, they are far from being completely uniform because of the cultivation and spread of nonnative plants (Cohen's d < 0.25). Similarly, biotic homogenization within cities was also limited; in particular, cultivated and spontaneous nonnative species strongly differentiated yard floras in PHX and LA, respectively, when compared to the other cities. Although other studies have reported mixed effects of non-native species on urban biotic homogenization at regional scales (McKinney 2004, 2006, K€ uhn and Klotz 2006, Qian et al. 2008, Lososov a et al. 2012, Ricotta et al. 2014, none examined their effect in highly managed urban areas such as private yards and at the local (within-metropolitan region) scale.
Our result that spontaneous non-native species contributed to homogenization of flora across the United States supports the idea that, at local scales, cultivated species are more yard-specific, whereas spontaneous species are more cosmopolitan. Most of the non-native spontaneous species recorded in our study correspond to weeds that were typically introduced to the United States decades ago (Ridley 1930) and that have had sufficient time to occupy most of the suitable habitats and to spread across larger geographical areas. Likewise, in Europe, archaeophytes (species introduced before AD 1500) have been found to homogenize urban floras (La Sorte et al. 2007, Lososov a et al. 2012, Ricotta et al. 2014. These cosmopolitan weeds tend to tolerate a wide range of environmental conditions and are thus able to invade available niches within broad climatic limits. The contribution of widespread weeds to the homogenization of yard floras was also supported by the significant correlation between our H index and the relative frequency of spontaneous non-native species in yards (Appendix S4: Table S3). The convergence between the observed mean H index and that of the null model within each city also highlights the relevance of random colonization processes in determining the homogenization of the spontaneous pool at local scales. However, across cities the significant (P = 0.034) deviation of the observed mean H index from that of the null model indicates that homogenization as a result of the dispersal and establishment of weeds in yards across the country has not been completely random.
For the cultivated pool, the observed mean H index also deviated from the null model both across and within some of the cities (BOS, LA, MIA, and PHX), which indicates that variation in homogenization patterns, and thus species composition, cannot be explained solely by habitat qualities, species natural dispersion, and macroclimatic variables. Rather, the geographic context of the studied yards must be also taken into account. For example, non-native cultivated species are strongly subject to local temperature and homeowners' preferences and management (Williams 2009, Kendal et al. 2012b, Aronson et al. 2016. Moreover, because they often have a limited potential range, they tend to increase floristic differentiation at the continental scale. In the three warmer southern cities (LA, MIA, and PHX), there was a clear tendency for cultivated non-native species to differentiate yard floras. In contrast, in the cooler northern cities (MSP, BOS, and BAL) cultivated non-native species tended to homogenize yard floras (note that in SLC, nonnative species promoted differentiation although with the lowest effect size across all cities; Cohen's d < 0.15). The proportion of non-native species was also consistently higher in the southern than in the northern cities (excluding again SLC; see Table 2, Fig. 2). These findings reinforce the idea that latitudinal and climatic factors influence regional urban floras (Williams 2009, Aronson et al. 2016, Jenerette 2016, Padull es Cubino 2018, and suggest that warmer urban areas have a larger available pool of horticultural non-native species adapted to local conditions than the available pool of natives. This pattern could reflect modern consequences of anthropogenic activities operating at an increasingly global extent that potentially promote the introduction of non-native species from all over the world (La Sorte and McKinney 2006). Although previous studies have shown that biotic dissimilarities increase with geographic distance (McKinney 2004, Qian et al. 2008, Lososov a et al. 2012, our H index showed no sign of spatial autocorrelation. Furthermore, our results contradicted the idea that biotic homogenization increases with increasing non-native species richness (McKinney 2004), at least in urban yards at the city scale.

Plant traits promoting phylogenetic homogenization
The consequences of environmental and anthropogenic influences on biotic homogenization at different scales were also reflected in plant functional traits, especially among the spontaneous pool. For example, among non-native spontaneous species, shorter, smaller-seeded plants and grasses (both widespread across the seven cities; see Appendix S4: Table S4) were significantly associated with homogenization of yard floras across cities. Frequent disturbance, commonly associated with urban areas (e.g., yard weeding), may favor shorter species with small seeds that reproduce and disperse more efficiently (Gilbert 1989). For example, spontaneous tree saplings are often removed by yard owners because these species are perceived as undesirable, untidy, or potentially dangerous (Kirkpatrick et al. 2007). Furthermore, seed mass of non-native spontaneous species has been found to be negatively related to their frequency in urban environments, such that smaller-seeded species occur more frequently (Thompson and McCarthy 2008), likely reflecting the trade-off between seed size, seed production, and dispersal distance. In contrast, species with greater seed mass have been proposed to perform better in stressful environments that are subject to drought, competition for light or water, deep shade, or burial (Westoby et al. 2002).
Literature on biotic homogenization suggests that the common features of urban areas create conditions that select for species with similar traits and life histories. Yet, our study also supports the idea that particular characteristics of individual cities-such as development history, human demographics, and cultural attributesimpose strong constraints on plant traits (McKinney 2006, Groffman 2014) that cause divergent patterns of association between plant traits and homogenization among cities. In fact, for the spontaneous pool, no plant trait consistently promoted homogenization in more than two cities. For example, annual plants only promoted homogenization in LA and SLC, and C 4 /CAM species only did so in LA and PHX.
Grasses were found to homogenize cultivated and spontaneous yard floras across cities, indicating that the similarity of lawn species across yards is driven, at least in part, by the influence of non-natives (Wheeler 2017). This pattern supports the idea that lawn landscaping styles across the residential United States have caused ecological homogenization (McKinney 2006, Groffman 2014, Wheeler 2017) that maintains similar ecosystems across a broad range of ecoregions. In this regard, the non-native turfgrass that contributed the most to the homogenization of yard floras and that was present in all seven cities was Bermuda grass (Cynodon dactylon), a medium-to fine-textured warm-season grass that spreads by rhizomes and stolons and has excellent heat, drought, and salt tolerance. Previous research assessing lawn and overall plant composition in the same seven cities concluded that species in the same four cooler northern cities (BAL, BOS, MSP, and SLC) were generally similar to each other, and species in the same three warmer cities (LA, MIA, and PHX) were distinct from the cool cities but less similar to each other (Wheeler 2017, Pearse 2018. Here, we have expanded on these findings by reporting that non-native grasses, despite being taxonomically different across regions, may have similar phylogenetic homogenizing effects on yard floras, as shown, for example, in the two warmest cities (LA and PHX).
Plant height was the best model predicting homogenization of the cultivated flora, with taller non-native cultivated plants causing differentiation across cities, further indicating that these species are more city-specific than shorter species. In this regard, a majority of trees in the southwestern U.S. urban areas were shown to be non-native (Nowak 2012, Avolio et al. 2018) because the original vegetation prior to urbanization was predominantly shrubland. Pincetl et al. (2013) quantified the tree species available for purchase in Los Angeles nurseries throughout the 20th century and found that there has been an increasing range of commercially available tree species sold by local nurseries, which may help explain this differentiation potential.
Showy flowers of non-native plants were also among the top-ranked traits in models for species contribution to homogenization. Accordingly, flowering is a conspicuous characteristic that is highly aesthetically regarded and preferred over non-flowering vegetation in private landscapes (Kendal et al. 2012a, Goodness et al. 2016. Although preferences for flower size, shape, and color may vary along gardener's personal landscape styles, our results suggest that flower showiness of non-native plants promoted homogenization of yard floras not only at the continental scale but also at the regional scale (particularly in LA), revealing the importance of aesthetics in determining ecological processes in urban landscapes.
Our water use index derived from Costello and Jones' (2014) WUCOLS list showed no effect on the homogenization of yard floras. This indicates that non-native plants with different levels of water requirement were evenly distributed across all yards, providing further evidence that irrigation and management practices might be compensating for variation in precipitation and thereby reducing climatic stress (Kendal et al. 2012b, Padull es Cubino 2018. Plant cold hardiness, although not included in this study because of unavailability of consistent data, is also a major factor that explains species distributions in urban environments at the continental scale (Jenerette 2016) and may influence the role that cultivated non-natives play in differentiating the flora across cities.
Although our study provides evidence that plant functional traits help explain phylogenetic homogenization of yard floras, the cities included in our analysis contained different numbers of sampled yards (between 21 and 30), which might have favored representation of certain locally abundant species and ultimately their contribution to biotic homogenization/differentiation. In contrast, some nonnative species might have not been detected in our vegetation surveys, as shown by nonsaturated accumulation curves (Appendix S4: Fig. S3). We believe that such limitations, although not affecting our conclusions, might be addressed by a more exhaustive and balanced sampling of yard floras at the continental scale.

CONCLUSIONS
Our study compared the degree of phylogenetic homogenization among yards at different spatial scales, discriminated between cultivated and spontaneous plant species pools, and specifically examined the role that non-native species play. Our results demonstrated that homogenization of yard floras by non-native species was prevalent for the spontaneous (weedy) flora but was context-and scale-dependent. Although non-native spontaneous species generally increased homogenization among yards across cities, non-native cultivated species tended to differentiate yard floras at the same scale. Importantly, such patterns had a limited effect on the mean value of the homogenization index, thus suggesting that urban yard floras at the continental scale still maintain their local idiosyncrasy and are far from being completely phylogenetically uniform; yet they have enough in common floristically to keep them from being completely heterogeneous. Furthermore, contrasting patterns of biotic homogenization by non-natives among and within cities suggest that local environmental and anthropogenic factors played a significant role in determining species pools. Hence, homogenization must be assessed for each particular location and species pool and further research should elucidate how homogenization varies across multiple environmental and anthropogenic gradients.
The inclusion of locally measured plant functional traits in the study of urban biotic homogenization may help overcome the lack of available trait information from accessible databases-as in this study-and enhance our understanding of the processes affecting ecosystem function. An assessment of the best way of computing and interpreting homogenization measures is also required to establish comparable results. Lists of species promoting homogenization and/or differentiation at different scales can inform urban planners and land managers interested in influencing the configuration of floras in urban landscapes in a context of global change. Further research is needed to elucidate the linkages between residents' preferences for vegetation and landscaping styles, and the plant functional traits that promote biotic homogenization.