Phylogenetic Diversity of Urban Floras in the Central Urals

Modern cities harbor a high diversity of plants, and urban floras are significantly different from non-urban floras especially when considering the proportion of alien species found in cities. However, it is not clear whether urban areas disproportionately select for species from relatively few evolutionary lineages or provide opportunities for species across the full spectrum of plant lineages. Here, we examined the taxonomic and phylogenetic diversity of the floras in four cities (Yekaterinburg, Kamensk-Uralsky, Krasnoufimsk, and Turinsk) in the understudied region of Central Urals (Russian Federation). We classified native species into indigenous and apophytic species, namely, those that are sensitive to anthropogenic disturbance and those that have expanded their range with human activity, respectively. Alien species were classified into archaeophytes and neophytes according to when they were introduced (i.e., before or after than 1800). Phylogenetic diversity was quantified using Faith’s index to reflect total evolutionary history in urban areas and mean phylogenetic distance (MPD) to reflect species dissimilarity. Phylogenetic diversity of native species was higher than that for alien species, and the standardized effect size (SES) of MPD for natives was positive, reflecting their general dissimilarity from one another, while it was very negative for aliens, showing that they were phylogenetically clustered. However, among natives, apophytes were significantly clustered, while indigenous species were overdispersed. For the aliens, MPD was higher for archaeophytes compared to neophytes, though both groups were significantly clustered. These results show that urbanization leads to a non-random selection of plants. Apophytes and alien plants were composed of closely related species, reflecting similar ecological traits and are likely to be pre-adapted to the environmentally altered and highly disturbed urban environment.


INTRODUCTION
Anthropogenic activity has weakened natural biogeographic barriers that limit the distributions of plants (Kueffer, 2017;Potgieter and Cadotte, 2020). These activities have intensified species migration and have contributed to the establishment of alien plant assemblages in the floras of many different geographical areas (Lonsdale, 1999;Aronson et al., 2016). The highest numbers of alien plants are concentrated in urban floras where purposeful introductions and environmental changes are most apparent (Aronson et al., 2014;Potgieter and Cadotte, 2020). Globally, it has been shown that cities generally harbor a high diversity of alien plant species; frequently between 25 and 50% of urban floras are classified as alien (Ricotta et al., 2009;Cadotte, 2020). The presence of alien plants influences the overall composition and structure of urban floras (Ricotta et al., 2009(Ricotta et al., , 2012Cadotte, 2020). The net effect of urban environmental conditions, human preferences, and intensified species immigration are hypothesized to be major drivers of biotic homogenization in cities (McKinney and Lockwood, 1999;Olden et al., 2004;McKinney, 2006;La Sorte et al., 2007).
However, taxonomic diversity does not directly reveal ecological patterns. It cannot reveal whether harsh urban environments select for similar species or whether species differences support coexistence (Cavender-Bares et al., 2009;Cadotte and Tucker, 2017). Biotic homogenization should be driven by two distinct mechanisms. First, species benefiting from urbanization should increase the similarity among urban floras. Second, if urbanization selects for specific traits and niches (e.g., Zhu et al., 2019), then relatively few, non-random groups should be driving this homogenization.
To account for such non-random species composition, different facets of biodiversity have been considered, including phylogenetic diversity, which takes into account the phylogenetic relationships between species based on the evolutionary relationships connecting all species together (Cadotte and Davies, 2016;Bitomský et al., 2020).
The use of phylogenies in ecology has expanded greatly over the past 20 years, and the reasons to use this information are twofold. First is that phylogenies provide direct information about how evolutionary history and specifically speciation events shape local diversity (Gerhold et al., 2015), thus revealing the ways in which human activities reshape the influence of this evolutionary history. Second is that it is often inferred that phylogenies are a good surrogate for functional, ecological, or niche diversity, especially when we conceive of functional diversity as including large numbers of traits in multivariate contexts (Tucker et al., 2018). While the relationship between phylogenetic and functional diversity can be complex and influenced by methodological decisions and ecological processes (Cadotte et al., 2019), this surrogacy has been shown to be the case in urban areas (Lososová et al., 2016), and so species expanding their ranges with urbanization likely result in reduced functional space relative to the native habitats that existed prior to human settlement and urban expansion. Furthermore, given the importance of ecosystem function and the delivery of ecosystem services in urban areas, reduced phylogenetic and functional diversity might result in lower functioning (Flynn et al., 2011;Cadotte, 2013Cadotte, , 2017. Phylogenetic analysis can reveal non-random patterns based on speciation and biogeography or from ecological processes that select for certain species. Evolutionarily, we might expect that the native flora is the product of relatively few successful clades and would therefore be relatively closely related, while alien species are a selection across many clades and regions and therefore represent higher phylogenetic diversity. Conversely, urban stressors and disturbances might select for species with the requisite adaptations allowing them to persist, resulting in both groups being phylogenetically underdispersed or clustered. A few studies have shown that alien species are phylogenetically less diverse than natives in cities (e.g., Čeplová et al., 2015;Knapp et al., 2017;Zhu et al., 2019), but there are surprisingly few such analyses.
However, native and alien floras themselves are heterogeneous assemblages and subdividing them in meaningful ways can further reveal how urbanization influences plant composition. The native flora can be separated into species that appear to be sensitive to human activity and mostly persist in natural ecosystems (i.e., what we refer to as indigenous) and into those that appear to have expanded their range and abundance by flourishing in human-dominated landscapes (i.e., apophytes) (Sukopp, 2006). We should expect that apophytes are more closely related than other species because they possess traits and strategies that are better suited to urban environments. The alien flora can be divided into species that assimilated into the regional flora long ago (i.e., archaeophytes) and those that were introduced more recently (i.e., neophytes) (Pysek et al., 2004). In an analysis of urban floras in Europe and the United States, Ricotta et al. (2009Ricotta et al. ( , 2012 showed that the phylogenetic diversity of alien species was significantly lower than that of native species, and archaeophytes were characterized by lower phylogenetic diversity compared to neophytes. A potential reason is that neophytes are more likely to be planted and supported by gardening in cities more recently (after 1800) because of their uniqueness, while archaeophytes are those species with preadaptations that allow them to persist or thrive after accidental introduction through early trade and agricultural expansion.
Despite the size of Russia, the diversity of cities, and its long history, analyses of urban biodiversity in Russian cities have not appeared in many international journals. Comparative analyses of the taxonomic structure of native and alien plants of urban floras in Russia have been performed and primarily published in Russian journals and monographs (Grigoryevskaya, 2000;Terekhina, 2000;Antipina, 2002;Panasenko, 2003;Sutkin, 2006;Tretyakova and Shurova, 2013;Senator et al., 2015;Antipova and Antipova, 2016;Golovanov and Abramova, 2017), and the patterns uncovered in this region can help inform our understanding of urbanization impacts on biodiversity globally. Across Russia, species richness of floras for cities with more than 1,000,000 people is about 1,200 species on average, about 800 for cities with 100,000-1,000,000, and about 500 species for towns with fewer than 100,000 people; the proportion of alien plants in Russian cities reaches 30-50% (Veselkin et al., 2017;Tretyakova et al., 2018). Urban floras are characterized by high taxonomic diversity, often exceeding diversity compared to the non-urban flora in the same regions (Tretyakova, 2016). There are obvious differences in the taxonomic composition of urban floras in Russian cities compared to non-urban floras of the same region. Species diversity of spore plants, gymnosperms, and monocotyledons is lower in urban floras, as well as the proportion of representatives of certain families of dicots that are typical in the boreal flora (Tolmachev, 1974;Khokhryakov, 2000). At the same time, families that tend to be affiliated with disturbed habitats are over-represented in urban floras (Terekhina, 2000;Antipina, 2002;Beresutsky and Panin, 2007;Antipova and Antipova, 2016;Tretyakova, 2016;Golovanov and Abramova, 2017).
In Russia, the phylogenetic diversity of urban flora has not yet been considered and so the goal of the present study was to characterize the phylogenetic diversity of native and alien plants in the urban floras of Sverdlovsk region in the Central Urals. This region was selected due to the high level of urbanization (proportion of urban population is about 88%) and its location on the transition zone between taiga and forest-steppe biomes. We expect (1) higher similarity of apophytes and archaeophytes among cities because these groups are composed of species adapted to disturbed habitats. We hypothesize that (2) alien species will exhibit lower phylogenetic diversity compared to native species since both alien and native species are nonrandom sets of species: aliens are able to excel under human preference and urban environmental filters and are thus likely drawn from fewer clades while natives are adapted to a wider range of conditions and are thus likely to be drawn from more disparate clades. We further hypothesize (3) apophytes to show similar patterns to aliens and exhibit low phylogenetic diversity compared to indigenous species, again because they are nonrandom subsets that can thrive under human influences. Finally, we expect (4) archaeophytes to have lower phylogenetic diversity than neophytes, since the latter are likely conscientiously drawn from more diverse species pools from around the world.

MATERIALS AND METHODS
Four cities were selected for our analysis of urban flora phylogenetic diversity: Krasnoufimsk, Yekaterinburg, Kamensk-Uralsky, and Turinsk (Figure 1). The selected cities differ from each other in terms of their geographic location, size, and level of economic development. Krasnoufimsk is located on the border of the western foothills of the Urals and the East European Plain, Yekaterinburg is on the border of the eastern foothills and the spine part of the Urals, Kamensk-Uralsky is on the border of the eastern foothills of the Urals and the West Siberian Plain, and Turinsk is on the West Siberian Plain. Together, they form a chain (from 56 • 37 N to 58 • 02 N and from 57 • 46 E to 63 • 42 E), which covers the Cis-Urals (geological features to the west of the Urals), mountain part of the Central Urals, Trans-Urals, and Western Siberia. In terms of population, the cities ranged from 17,000 to over 1.4 million residents ( Table 1).
The floras of these cities were studied from 2000 to 2016. The urban floras include all vascular plant entities (divisions Lycopodiophyta, Polypodiophyta, Pinophyta, and Magnoliophyta) that occur spontaneously within the administrative boundaries of the cities. The urban flora for the Central Urals is based on field surveys performed by the authors. Every type of habitat, including both natural/seminatural and humanmade (Figure 2), was thoroughly surveyed for the presence of species. Species lists are regularly updated based on new findings. These data were supplemented with the information from herbarium collections [Herbariums of the Institute of Plant and Animal Ecology of the Ural Branch of the Russian Academy of Sciences (SVER), Ural Federal University (UFU), and Kurgan State University] and published data on the flora of the Central Urals. The goal of these surveys and herbarium searches was to compile as complete a list as possible.
Native and alien species were distinguished in the urban flora using checklists of the flora of Sverdlovsk Region (Knyazev et al., 2016(Knyazev et al., , 2017(Knyazev et al., , 2018(Knyazev et al., , 2019a(Knyazev et al., ,b, 2020 in which species status was evaluated based on the appearance in the studied area due to human activity. We categorized native species as either indigenous or apophytic based on whether they exhibited expanded ranges because of their association with disturbed habitat. The indigenous plants included species confined to natural communities and that generally avoid anthropogenic habitats while apophyte plants included anthropogenic habitatassociated species from the native flora. The assignment of plants to the selected groups was based on the species distribution analysis in natural and anthropogenic habitats in the urban habitats (see Tretyakova, 2014;Baranova et al., 2018).
We considered alien species as those species unintentionally introduced into the city territory, as well as decorative or purposefully introduced species, for which the presence of seed or vegetative propagules were observed, or plants found outside their cultivation sites (Tretyakova and Shurova, 2013;Baranova et al., 2018). Alien species were divided into two groups based on their residence time: archaeophytes and neophytes (Pysek et al., 2004). Archaeophytes are the alien species that appeared in the Central Urals before the arrival of the Russian population prior to 1800. Neophytes appeared in the Central Urals after this date. The checklist of the flora of Sverdlovsk region, parts I-VI (Knyazev et al., 2016(Knyazev et al., , 2017(Knyazev et al., , 2018(Knyazev et al., , 2019a(Knyazev et al., ,b, 2020, served as a basic source for assignment to a groups of archaeophytes and neophytes. We included records from Flora Rossica by K. F. Ledebourg (von Ledebour et al., 1842-1853 that summarized information about the flora of the Urals, accumulated by the beginning of the 19th century.
Final lists of the studied urban floras, including grouping of species into indigenous plants, apophytes, neophytes, and archaeophytes, are available in the GBIF repository 1 .
When assessing the similarity of the species composition among urban floras, the Jaccard similarity coefficient was used, which was calculated as a/(a + b + c), where a is the number of species presented in both lists, and b and c are the numbers of unique species (Legendre and Legendre, 2012).
We used a combined dated phylogenetic tree published by Zanne et al. (2014), which includes 32,223 species of land plants to extract a phylogeny for our region. In the urban flora in the Central Urals, 1,035 plant species were recorded, and synonymy was aligned with The Plant List (The Plant List: URL 2 ). Two hundred fifty-six of our observed plant species were absent in the phylogenetic tree, and so we added them as polytomies to the corresponding genera (239), and in the absence of genera,  (17).
To characterize the phylogenetic diversity in these urban floras, we calculated Faith's PD, which sums the phylogenetic tree branch lengths that connect all species of a given urban flora or the corresponding grouping (e.g., native vs. alien). This index is a phylogenetic analog of species richness, and it is usually strongly correlated with it (Faith, 1992;Tucker et al., 2017). We also used a complementary phylogenetic diversity measure, the average pairwise phylogenetic distance between species [mean phylogenetic distance (MPD)], which is relatively insensitive to species richness and is well-suited to comparing phylogenetic and taxonomic diversity facets of different-sized species lists, as well as providing insights into the roles of different ecological mechanisms like environmental filtering and species interactions (Webb et al., 2002;Cadotte and Davies, 2016).
We compared the phylogenetic diversity in each of the species groups listed in Table 1 to null expectations generated from randomizations of species membership in each grouping for each city (Webb, 2000;Ricotta et al., 2009;Cadotte and Davies, 2016). The complete urban flora observed for each urban area was considered as a pool of species from which the plant classes were assembled. We estimated MPD for the random assemblages and repeated this 1,000 times. We then calculated a standardized effect size (SES) by subtracting the null model mean and dividing the difference by the null model standard deviation (Webb, 2000;Cadotte and Davies, 2016). Negative SES.MPD values correspond to phylogenetic clustering; in this case, MPD is lower compared to the null expectation, i.e., the group is composed of more related species. On the contrary, positive SES.MPD values indicate phylogenetic overdispersion. Null model analysis can also test significance of clustering or overdispersion effect. The corresponding p-value is calculated as a proportion of null distribution lesser than the empirical value of a metric. Clustering is considered to be significant when p < 0.025 and overdispersion is considered to be significant when p > 0.975 (two-tailed hypothesis; Swenson, 2014). All calculations were performed in R version 4.0.2 (R Core Team, 2020). We used package picante (Kembel et al., 2010) to calculate all phylogenetic metrics including PD, MPD, and SES.MPD.

RESULTS
The number of species recorded in these urban floras varied from 288 species in Turinsk to 960 species in Yekaterinburg. All major taxonomic groups of vascular plants are represented in these urban floras, but angiosperms (96.2-96.8%) and, specifically, dicotyledonous (74.3-81.9%) plants dominate these lists (see Table 1). Native plants constitute the majority of these urban floras. In the analyzed cities, with, on average, 67% of species being native, among them the indigenous class contained more species than apophytes. Alien plants account for about 30% of the urban species and were predominantly neophytes. Archaeophytes and neophytes account for 7% and 25% of all urban flora, respectively (see Table 1).
Flora similarity analysis revealed that species composition in the indigenous group was significantly more differentiated among cities than the apophytes (t = 5.24, p < 0.001). The Jaccard coefficients for indigenous plants varied from 0.14 to 0.53, with a mean of 0.32, and apophytes varied from 0.67 to 0.94, with a mean of 0.82 (see Table 2). Among alien plants, species similarity was marginally significantly higher in archaeophytes (archaeophyte mean = 0.69, neophyte mean = 0.48, t = 2.08, p = 0.071; see Table 2).
The results show that phylogenetic diversity, assessed by the Faith's index, was higher in the urban flora of Yekaterinburg and lower in Turinsk, as was predicted by differences in species richness. At the same time, all the considered urban floras had relatively similar MPD values (see Table 3) ranging from 297 to 312 Ma.
The MPD of native species in all our urban floras was significantly higher than that of alien species (native mean = 322.75, alien mean = 252.80, t = 10.67, p < 0.001). The MPD of the indigenous species was significantly higher than that of apophytic (indigenous mean = 347.90, apophytic mean = 281.65, t = 6.52, p = 0.006) and alien species (alien mean = 252.80, t = 8.33, p < 0.001) while the MPD of apophytes was also significantly higher than for alien species (t = 5.078, p = 0.01), albeit with a smaller magnitude of this difference. Within the alien species, neophytes had significantly higher MPD values than archaeophytes (neophyte mean = 257.12, archaeophyte mean = 235.42, t = 3.18, p = 0.023).
In general, the MPD of all species groups increased with the city size. It was highest in the urban flora of Yekaterinburg and smallest in the urban flora of Turinsk (see Table 3) with the exception of indigenous and apophytic species MPD being highest in the flora of Turinsk (the smallest city).
SES analysis showed (Figure 3) that in all four urban floras, the MPD of native species exceeded the random expectation for MPD (i.e., SES.MPD was positive), indicating phylogenetic overdispersion, while the SES.MPD of alien species, on the contrary, was negative, and exhibited substantial phylogenetic clustering. Interestingly, the two native plant groupings behaved in fundamentally different ways: the phylogenetic overdispersion of native plants was caused by indigenous species, while clustering was apparent for apophyte species (Figure 3). Both alien species classes exhibited phylogenetic clustering (Figure 3).

DISCUSSION
In the present study, we addressed the task of quantitative description of phylogenetic diversity and structure of urban floras of four cities in the understudied region in the Central Urals in Russia. We compared metrics of conventional taxonomic and phylogenetic diversity and also performed null-model analysis to reveal patterns of clustering and overdispersion in groups of native and alien species. Our results show that phylogenetic diversity of native species was higher than that for alien species, and the SES of MPD for natives was positive, reflecting their general dissimilarity from one another, while it was very negative for aliens, showing that they were phylogenetically clustered. However, among natives, apophytes were significantly clustered, while indigenous species were overdispersed. For the aliens, MPD was higher for archaeophytes compared to neophytes, though both groups were significantly clustered.
Our results also show that larger cities (both in terms of area and population size) support more species in total and in the various plant groupings, as has been shown previously (Pyšek, 1998;Veselkin et al., 2017;Tretyakova et al., 2018). The cities examined here harbor substantial numbers of alien taxa, more than 30% of the urban flora, which is remarkably similar to other urban plant analyses (Ricotta et al., 2009;Aronson et al., 2014;Cadotte, 2020).
Urban floras are undergoing fundamental changes, making it different from non-urban flora in diversity, composition, and structure. It is possible to identify these changes using a mix of approaches, and here we show these urban flora effects using measures of diversity, beta-diversity, and phylogenetic patterns. It is clear that species that expand their ranges in response to human activities are increasing similarity in the urban floras. Both apophytes and aliens had substantially higher Jaccard similarity values than for indigenous species, indicating that urbanization drives biotic homogenization (McKinney, 2006).
Our results show that the phylogenetic diversity of native species is higher than that of alien species, indicating that natives, and especially indigenous species, constitute a greater diversity of lineages. Among alien species, the MPD is higher in neophytes than in archaeophytes and similar observations have been made for the floras of European and American cities (Ricotta et al., 2009(Ricotta et al., , 2012. More importantly, though, indigenous species were significantly overdispersed while apophytes and aliens were significantly clustered (Figure 3). These results indicate that urbanization selects for species non-randomly. An analysis of the taxonomic, ecological, and geographical diversity of the species occurring in our region may help explain our findings of phylogenetic clustering among aliens and apophytes.
The decrease in the MPD of alien species is associated with a simplification of their taxonomic structure in comparison with  native species. Native species encompass all major clades of the flora of temperate latitudes (Lycopodiophyta, Polypodiophyta, Pinophyta, and Magnoliophyta). Among the alien plants, there are representatives of only two clades, Pinophyta and Magnoliophyta; spore plants are completely absent. Similarly, the decrease in MPD of apophytes in comparison with indigenous species can be explained similarly by a decrease in their taxonomic diversity. The apophytic group lacks representatives of the Lycopodiophyta, Polypodiophyta, and Pinophyta; the number of families is decreasing too; for example, representatives of the families Cyperaceae and Orchidaceae are absent. The strong phylogenetic clustering among alien species and apophytes implies that these closely related species possess a few suites of ecological traits that allow them to adapt to certain environmental conditions and coexist in urban environments (Weiher and Keddy, 1995;Knapp et al., 2008Knapp et al., , 2017Proche et al., 2008;Ricotta et al., 2009) and are the outcome of phylogenetic conservatism of niches (Losos, 2008). For example, there are higher ambient temperatures in the city compared to the surroundings, which allows thermophilic plants to exist. Urban soils are highly alkaline and highly saline, which provides plants adapted to saline or high pH soils (Sukopp, 2004;Godefroid et al., 2007;Thompson and McCarthy, 2008).
Conversely to the patterns for alien species and apophytes, indigenous species were statistically significantly phylogenetically overdispersed. Indigenous species have more functional differences, can use ecological space more efficiently, and exist in a wide range of conditions (MacArthur and Levins, 1967;Proche et al., 2008). The effect of phylogenetic overdispersion is usually interpreted as an evidence for the leading role of biotic interactions leading to the competitive exclusion of closely related species and thus increasing MPD (Webb, 2000;Mayfield and Levine, 2010). We believe that this is not the case with respect to the indigenous in our urban floras. These indigenous species are associated with a set of heterogeneous habitats and there is not likely to be direct interactions between species that occur in different habitats, so it is difficult to consider the effect of competitive exclusion. Furthermore, the scale we are examining is well above that for individual interactions, and this scale should better detect the influences of the environment and habitat diversity on phylogenetic patterns (Swenson, 2019), or the fact that the natives reflect the evolutionarily diverse lineages that evolved in this region or expanded into this region after the last glaciation event, meaning that underlying biogeographic histories of speciation and migration drive local ecophylogenetic patterns (Gerhold et al., 2015).
An analysis of the ecological (habitat types) and geographical (distribution area) diversity can help explain the effect of urbanization on plant composition. Simplification of ecological diversity through habitat loss and homogenization can explain the low phylogenetic diversity of apophytes in comparison with indigenous plants, as well as archaeophytes in comparison with neophytes. Lower habitat diversity results in selection for certain traits in the apophytes in comparison to indigenous plants and is evidenced by the absence of halophytic plants, as well as rocky, petrophytic, and steppe plants. Furthermore, the apophytes appear to have larger ranges; there are no Asian, North Asian, and Ural endemic species, as well as hypoarctic-boreal, arctoboreal, and steppe species.
Similarly, archaeophytes are mostly segetal and ruderal plants (Pyšek et al., 2002) mostly from South European (Mediterranean) and Central Asian steppe regions (Pyšek et al., 2002;Ricotta et al., 2009). Conversely, neophytes are a more diverse assemblage from a wider geographic origin, including American (North American, Central American, and South American), African, East Asian, and Siberian plants with multiple avenues of introduction, and among them, there are species of meadow, steppe, coastal, halophytic, and aquatic communities. Most prominently, these neophytes include many purposefully planted ornamental species, which should come from a more diverse sampling of evolutionary lineages, thus maintaining higher phylogenetic diversity , but in our findings, neophytes do not contribute to high phylogenetic diversity. Neophytes, like archaeophytes, are from disproportionately few lineages, perhaps because of not only the influences of urbanization, but also the fact that the location climate in the Ural region might limit successfully establishing species to those adapted to colder climates.
Conversely, the group of indigenous plants is highly variable both in the species number (from 71 in urban flora of Turinsk to 478 in urban flora of Yekaterinburg) and in species composition [the species similarity is 0.14-0.53 (see Tables 1, 3)]. The variability of the indigenous species group is primarily provided by the composition of natural habitats existing in the urban area. In the urban flora of Yekaterinburg, Krasnoufimsk, and Turinsk, forest plants are the richest group; the group of meadow plants is in second place for this indicator. The urban flora of Yekaterinburg contains a more diverse group of bog plants. In the urban flora of Kamensk-Uralsky, steppe species are in first place, the proportion of which increases due to groups of rocky, rocky petrophytic steppe, and petrophytic steppe plants, and forest and meadow plants are shifted to the second and third places. Urban flora of Kamensk-Uralsky is also distinguished by the presence of a large group of halophytic plants in its composition.
While we can attribute these phylogenetic patterns to the influences of urbanization, the reality is that urban environments can contain a complex mix of relictual and novel habitats (Aronson et al., 2017) that will undoubtedly select for different kinds of species. For example, land use history in urban areas can greatly impact phylogenetic diversity (Cheng et al., 2018), with relict forests containing higher phylogenetic diversity than recovered and secondary forests, which tend to be phylogenetically clustered (Borges et al., 2020). As a result, as long as urban areas have relatively few intact relict habitats, and more disturbed or other impacted sites, we should see a high prevalence of phylogenetically clustered communities. Furthermore, there are likely feedbacks between land use and plant invasion on phylogenetic patterns, especially for native communities. Gutiérrez-Cánovas et al. (2020) show that along the Spanish coast, urbanization had stronger negative impacts on native taxonomic and functional diversity when sites also contained invasive plants. We were not able to tease apart these interactions with our data, but future work should examine interactions among different drivers of biodiversity change in cities.
Finally, we wish to note the value in protecting not only indigenous species, but also the ecological diversity they represent. It is understood that the biotic homogenization caused by urban development and alien introductions serves to reduce resilience to future environmental change (McKinney, 2006). However, within cities, human populations rely on the benefits provided by healthy and diverse ecosystems (Gómez-Baggethun and Barton, 2013). Not only do alien species have the potential to directly impact ecosystem service delivery (Potgieter et al., 2017), but the loss of native diversity can also impact ecosystem services. The functioning of ecosystems has been shown to decline with both the loss of species richness (Tilman et al., 2014) and the loss of phylogenetic diversity (Cadotte, 2013). Thus, we encourage municipal policies that preserve and enhance diverse native habitats within urban areas.

CONCLUSION
The urban flora is characterized by high alien species richness, but relatively low phylogenetic diversity. This reflects the fact that natives represent high taxonomic, geographic, and ecological differentiation and that they persist in a broad set ecological conditions and provide diverse benefits. The alien species are not a random set of species; they include taxonomically closely related species with similar ecological properties that might provide adaptation to urban environmental conditions. Thus, the high species richness of the urban flora of the Ural cities, as well as European ones (Knapp et al., 2008), is provided mainly by closely related species that are functionally similar and adapted to an urban environment.

AUTHOR CONTRIBUTIONS
AT, BY, and MC conceived this work and designed the methodology. AT, BY, PK, and NG collected and analyzed the data. AT and BY created the original draft. MC supervised the analysis and writing process. All authors edited and revised subsequent drafts of the manuscript, approved the final version, and agreed to be held accountable for the work.

FUNDING
Funding for this collaboration was provided to MC by the