Endangered Plants in Novel Urban Ecosystems Are Filtered by Strategy Type and Dispersal Syndrome, Not by Spatial Dependence on Natural Remnants

Understanding the contribution of cities to nature conservation is gaining increasing importance with a globally accelerating urbanization and requires insights into the mechanisms that underlie urban distribution patterns. While a considerable number of endangered plant species have been reported for cities, the spatial dependence of populations of these species on natural remnants versus anthropogenic ecosystems is critically understudied due to deficiencies in population distribution data. To which extent endangered species in anthropogenic ecosystems spatially rely on natural remnants is thus an open question. We used a unique dataset of 1,742 precisely mapped populations of 213 endangered plant species in the city of Berlin and related these point data to habitat patches that had been assigned to natural remnants, hybrid ecosystems and novel ecosystems according to the novel ecosystem approach. By applying point pattern analyses (Ripley’s K function, cross K function, cross pair correlation function) we unraveled the spatial dependence of the populations toward the different ecosystem types. Moreover, we tested how plant traits related to plant strategy and dispersal filter for species occurrence across ecosystems. Differentiating populations on anthropogenic sites revealed that populations in hybrid ecosystems spatially depended on natural remnants, but populations in novel ecosystems (i.e. more than a third of all populations) surprisingly didn’t. A conditional inference tree showed that endangered plant species in novel ecosystems are filtered for ruderal strategy type and wind dispersal syndrome, while competitive and stress-tolerant species were mainly confined to natural remnants. Our results highlight the importance of conserving natural remnants as habitats and seed sources of endangered plants. Yet novel urban ecosystems can support many populations of endangered plant species beyond the adjacency to natural remnants, with hybrid ecosystems likely acting as stepping stones. This indicates a specific contribution of urban ecosystems to biodiversity conservation. Since different filters modulate the species pools of different ecosystem types, novel urban ecosystems are not supposed to substitute fully the habitat functions of natural remnants. Our study thus highlights promising opportunities for involving the total range of urban ecosystem types into urban conservation approaches.

Understanding the contribution of cities to nature conservation is gaining increasing importance with a globally accelerating urbanization and requires insights into the mechanisms that underlie urban distribution patterns. While a considerable number of endangered plant species have been reported for cities, the spatial dependence of populations of these species on natural remnants versus anthropogenic ecosystems is critically understudied due to deficiencies in population distribution data. To which extent endangered species in anthropogenic ecosystems spatially rely on natural remnants is thus an open question. We used a unique dataset of 1,742 precisely mapped populations of 213 endangered plant species in the city of Berlin and related these point data to habitat patches that had been assigned to natural remnants, hybrid ecosystems and novel ecosystems according to the novel ecosystem approach. By applying point pattern analyses (Ripley's K function, cross K function, cross pair correlation function) we unraveled the spatial dependence of the populations toward the different ecosystem types. Moreover, we tested how plant traits related to plant strategy and dispersal filter for species occurrence across ecosystems. Differentiating populations on anthropogenic sites revealed that populations in hybrid ecosystems spatially depended on natural remnants, but populations in novel ecosystems (i.e., more than a third of all populations) surprisingly didn't. A conditional inference tree showed that endangered plant species in novel ecosystems are filtered for ruderal strategy type and wind dispersal syndrome, while competitive and stress-tolerant species were mainly confined to natural remnants. Our results highlight the importance of conserving natural remnants as habitats and seed sources of endangered plants. Yet novel urban ecosystems can support many populations of endangered plant species beyond the adjacency to natural remnants, with hybrid ecosystems likely acting as stepping stones. This indicates a specific contribution of urban ecosystems to biodiversity conservation.

INTRODUCTION
Cities can harbor a high biological richness (Kühn et al., 2004;Aronson et al., 2014), also including a considerable richness of endangered plant species (Schwartz et al., 2002;Lawson et al., 2008;Lenzen et al., 2009;Shwartz et al., 2014;Ives et al., 2016;Kowarik and von der Lippe, 2018;Planchuelo et al., 2019). Urban areas in Australia, for example, host 30% of the country's endangered plant species (Ives et al., 2016), with some of them occurring today exclusively in urban environments (Soanes and Lentini, 2019). While these studies indicate a high potential of cities to contribute to biodiversity conservation at a global scale (McKinney, 2008;Kowarik, 2011), survival in urban environments is a challenge for many species as demonstrated by establishment failures (Kowarik and von der Lippe, 2018), extinction debts (Hahs and McDonnell, 2014) and completed extinctions (Knapp et al., 2010;Duncan et al., 2011). Developing effective nature conservation strategies for cities thus requires insights into the mechanisms that underlie urban biodiversity patterns in general, and those of endangered species in particular.
Plants are not randomly distributed in cities, and biodiversity patterns were early identified in relation to urban structure, the level of human pressure, urbanization, or specific land use (Kowarik, 1990;McKinney, 2002McKinney, , 2008Shochat et al., 2006;Swan et al., 2011;Schmidt et al., 2014), thus resulting from an interaction of ecological and socio-economic parameters (Swan et al., 2011). Some major filters have been hypothesized to predict the composition of urban species assemblages by filtering for species with a specific suite of biological traits that support the colonization and persistence of species in cities Aronson et al., 2016;Kowarik and von der Lippe, 2018). Yet, despite increasing insights into the mechanisms of urban species assembly (Knapp et al., 2009;Beninde et al., 2015;Williams et al., 2015;Aronson et al., 2016), the responses of endangered plant species to urban filters are critically understudied (Shwartz et al., 2014). This is an essential gap of knowledge because responses of large species sets to urbanization cannot be simply translated to endangered species, as indicated by contrasting species' traits that have gone extinct or, alternatively, persist in cities (Knapp et al., 2010;Duncan et al., 2011;Williams et al., 2015).
Cities were often founded in landscapes of high geomorphological and biological richness (Kühn et al., 2004), leading to the inclusion of natural remnants into the urban matrix (e.g., Godefroid and Koedam, 2003;Rebelo et al., 2011;Schwartz et al., 2013). The importance of such remnants for endangered species is well established. Yet, an important part of the specific urban contribution to biodiversity conservation -beyond conserving natural remnants -depends on the role that anthropogenic urban ecosystems can play as habitats of endangered species. Such ecosystems are omnipresent in cities, resulting from a stepwise transformation of natural to novel urban ecosystems, and being modulated by different landuse types and associated legacies (Kowarik, 2011;Ramalho et al., 2014). Endangered species are filtered by their ability to colonize these usually highly fragmented anthropogenic sites through dispersal, and to negotiate with novel ecological features of urban sites.
Few recent studies have reported a considerable amount of endangered plant species also from habitats of specific urban origin (Kowarik and von der Lippe, 2018;Planchuelo et al., 2019;Soanes and Lentini, 2019), including two studies from Berlin, Germany. The first revealed that a similar share of endangered plant species was able to establish populations in natural remnants and anthropogenic ecosystems (Kowarik and von der Lippe, 2018). The second study showed that populations of a large set of plant species of highest conservation priority stretch over a broad range of natural and anthropogenic ecosystems; surprisingly with highest numbers of populations in novel ecosystems, but highest density of populations in natural remnants .
These studies highlight conservation opportunities for novel anthropogenic ecosystems. Yet the role of the latter as habitats of endangered species might be compromised by a spatial dependence to remnant ecosystems. Species typical of ancient forests, for example, are critically dispersal-limited (Jacquemyn et al., 2001). Urban filters that determine the species' assemblages of fragmented habitat patches Aronson et al., 2016) can thus be expected to determine as well the colonization success of endangered species in novel urban habitats, with natural remnants as source habitats. Thus, one critical question for understanding the urban contribution to biodiversity conservation is to which extent endangered species are spatially confined to natural remnants and their vicinity (source-sink dynamic ;Pulliam, 1988;Schreiber, 2010).
As dispersal is an ecological process that causes spatial autocorrelation in distribution patterns (Epperson, 2005;Bahn et al., 2008), spatial statistics can be used to assess whether or not populations located in anthropogenic urban ecosystems are still spatially dependent to populations in natural remnants. Point pattern analysis has already been widely applied in nonurban areas to assess the dependence of endangered plant species on specific habitat patches of other species (Wehenkel et al., 2015). Further examples involve work on the effects of animalmediated dispersal of an endangered shrub (Rodríguez-Pérez et al., 2012), or on the minimum areas for the conservation of an endangered plant (Lim et al., 2008). Even if the importance of spatial patterning for biodiversity conservation in cities is increasingly acknowledged (Rastandeh et al., 2017), insights into related spatial patterns and processes are scarce because exact spatial information on the locations of endangered species in urban landscapes is often missing.
Many of the plants that pass dispersal filters in cities tend to be ruderal species in terms of the CSR scheme (Grime, 1977), as they can rapidly colonize new urban habitats while tolerating a wide range of urban disturbances (Williams et al., 2015). Other species might be pre-adapted or adapting to urban environments, e.g., in terms of advantageous dispersal traits (McDonnell and Hahs, 2015). For instance, some seed traits related to dispersal increase the probability of human-mediated long-distance dispersal in cities (von der Lippe and Kowarik, 2008). Correspondingly, species with human vectors as an additional dispersal pathway increased in abundance over time in densely urbanized neighborhoods (Johnson et al., 2017) and were less affected by urban dispersal limitations (Chytrý et al., 2012).
Here, we tested for the spatial dependence between populations of endangered plant species in anthropogenic habitats and natural remnants in the urban landscape of a European metropolis. We used a unique dataset with the exact geographical position of 1,742 populations of 213 endangered plant species in Berlin. These species had been identified as species of highest conservation priority due to their endangerment in Berlin or at larger spatial scales, and their populations had been subsequently precisely mapped by experts (Meißner and Seitz, 2010). While we have previously analyzed  the affiliation of these populations to different types of biotopes and three types of ecosystems according to the novel ecosystem approach (i.e., natural remnants, hybrid ecosystems, novel ecosystems; Hobbs et al., 2009), we here focus on the spatial patterns among populations and the related biological features of the endangered species.
We hypothesized that populations of endangered plant species are not randomly distributed but spatially clustered in Berlin. We expected that populations in novel ecosystems occur independently from source populations in natural remnants since novel ecosystems filter for species with enhanced capacities of long-distance dispersal and for those that belong to a ruderal strategy type (Grime, 1977). Since hybrid ecosystems result from a less profound transformation of natural ecosystems compared to novel ecosystems, we assume that populations located in hybrid ecosystems still spatially relate to seed sources in natural remnants, with species pools filtered by different features compared to species that were able to colonize novel urban ecosystems.
In detail, we addressed the following research questions: (1) Are populations (n = 1,742) of endangered plant species (n = 213) located within natural remnants, hybrid ecosystems and novel ecosystems in Berlin spatially clustered? (2) Do populations located in the two types of anthropogenic ecosystems (i.e., hybrid ecosystems and novel ecosystems) spatially depend on natural remnants? If populations in novel ecosystems are not dependent on natural remnants, are they instead dependent on hybrid ecosystems? (3) Are endangered plant species with populations in natural remnants, hybrid ecosystems and novel ecosystems filtered by (a) dispersal-related species traits and (b) different strategy types according to the CSR scheme (Grime, 1977)?

Study Area and Population Data
This study was carried out in the city of Berlin, the capital and most populous city of Germany with 3.6 million inhabitants as of 2017. Berlin has a total area of 891km 2 , from which 59% is covered by built-up areas and streets, while green and blue spaces occupy 41% of the area -including forests (17.7%), lakes and rivers (6.1%), parks (5.6%), allotment gardens (5.3%), fields (5%), and grasslands (1.3%) (SenStadtUm, 2016).
Since the 19th century, the flora of Berlin has been thoroughly studied (Sukopp, 1987). Information on populations dynamics has been collected and updated numerous times in Red Lists of endangered species (Seitz et al., 2018). Until today, 17% of Berlin's flora has gone extinct since the mid-19th century and 29% is currently being endangered. In this study, we used a detailed dataset with the precise geographical location of 1,742 populations of 213 highly endangered plant species (Figure 1 and Supplementary Appendix S1) from Berlin's Flora Protection Program (Berliner Florenschutzkonzept). These species have been selected for the Flora Protection Program according to their: (1) red list status in Berlin and/or at higher spatial scales, (2) long term population trends (50-150 years), (3) short term population trends (1-15 years), and (4) biogeographical responsibility, i.e., when the range of a species in Berlin constitutes an important part of the global range of this species (Seitz, 2007). The populations of these species are being monitored since 2009 by several experts (Meißner and Seitz, 2010).

Biotope Types and Ecological Novelty
A methodological approach was developed in the 1980s to assign different biotope types to the whole area of Berlin by incorporating various combinations of land uses and environmental conditions (Sukopp and Weiler, 1988). The mapping was done as part of a combined work of 64 different projects between 2001 and 2014, where areas with high conservation value or low detectability in aerial photos such as forests were covered by field biotope mapping while built up areas and transportation corridors were primarily covered by data of the urban planning authorities (SenStadtUm, 2014). Today, the Berlin Environmental Atlas includes a very detailed area wide biotope mapping 1 that is updated regularly (SenStadtUm, 2014). Currently, the city of Berlin is divided into 79,268 patches, with each of these patches being assigned to one of 12 biotope classes. These are forests, grasslands, ruderal sites, standing bodies of waters, built-up areas, bogs and marshes, hedges, green spaces, fields, flowing moving bodies of waters, heaths, and other types. These biotope classes are in turn further subdivided in 7,483 biotope types at several hierarchical levels.
To relate the geographical position of each of the 1,742 populations of endangered plant species to remnant, hybrid and novel ecosystems and assess their spatial correlation toward natural remnants, we used the ecological novelty map of Berlin . Each of the 79,268 patches that had been discerned in the biotope map were assigned to either natural remnants, hybrid or novel ecosystems -see Planchuelo et al. (2019) for details: • Natural remnant ecosystems are relicts of natural ecosystems that are often slightly affected by urban impacts but still remain within their historical range of modifications. These include biotopes from mires and wetlands to near-natural forests or moderately used dry or wet grasslands (Figure 1).
• Hybrid ecosystems are human mediated ecosystems that have been modified from their historical counterparts but still have the potential to return to historical conditions. Examples range from managed grasslands to urban green spaces, young tree plantations in forests and parks or intermediate succession stages following abandonment. • Novel ecosystems are human shaped ecosystems that have potentially been irreversibly changed by large modifications to their abiotic conditions or biotic composition. They include built up areas, rooftops, vacant lots, abandoned industrial areas or intensively managed agricultural fields.

Spatial Analyses
To enable spatial point pattern analyses to assess the correlation between populations of endangered plants and patches of different novelty categories, the latter were represented by a point positioned at the centroid of each patch (Baddeley, 2008). Using the "polygon centroid" function in QGIS, we converted all 79,268 ecosystem patches that had been differentiated in Berlin's biotope map into point data.
To respond to the first research question, we first tested whether populations within natural remnants, hybrid ecosystems and novel ecosystems were spatially clustered. This was done by means of Ripley's L function (Dixon, 2002), which allows to verify if, and at which spatial scales, populations are clustered, regularly or randomly distributed. This test (function Linhom of the spatstat package in R) compares the spatial distribution of a point pattern against a generated set of points with a random distribution (which is named Poisson distribution, see Stoyan and Stoyan, 1994). If the curve of the L function of a point pattern lies above the Poisson line, the populations are clustered; if it lies below, the populations are regularly distributed; if it overlaps it, the populations are randomly distributed (Dixon, 2002). L plots of Ripley's K functions were used (Dixon, 2002), as they stabilize the variance and are more powerful than simple plots of the K function (Baddeley et al., 2015).
To respond to the second research question, we analyzed if populations in the two types of anthropogenic ecosystems (i.e., hybrid ecosystems with 632 populations and novel ecosystems with 687 populations) are spatially dependent on remnant ecosystem patches. Because one of our hypotheses was that populations in novel ecosystems are not related to natural remnants, we also tested if populations in novel ecosystems might instead be spatially dependent on hybrid ecosystem patches. Statistically significant spatial correlations between populations of endangered plant species and the patches of different novelty categories were assessed with global envelopes of cross L functions (function Lcross.inhom of the spatstat package in R). This spatial analysis compares the observed correlation between the two point patterns of the empirical data against a null model, which is generated by keeping the same location of the ecosystem patches while creating simulations with random locations for the populations (Baddeley et al., 2015). If the curve of the cross L function for the empirical point pattern transgresses outside of the global envelope of the null model at any moment, there is a statistically significant spatial correlation between the populations and the novelty category (Baddeley et al., 2014(Baddeley et al., , 2015. Alternatively, if the curve does not transgress outside of the global envelope, we consider that there is a missing spatial correlation between the populations and the novelty category. Global envelopes were chosen because they avoid the problem of data snooping of pointwise envelopes and they allow to achieve an exact significance level: α = 1/(1 + n • simulations), see Baddeley et al., 2014. Thus, in order to achieve a 95% confidence interval, we generated global envelopes from 19 simulations of the null model. L plots of cross K functions were used (Dixon, 2002), as they stabilize the variance and are more powerful than the K plots (Baddeley et al., 2015).
In many cases, spatial correlations can be simply fruit of the spatial configuration of the urban matrix. Populations in hybrid ecosystems, for example, might be spatially dependent on remnant ecosystems simply because patches of hybrid ecosystems tend to be located in the proximity of patches of remnant ecosystems. We thus assessed with cross L functions if the patches of the different novelty categories were spatially correlated or not (see Supplementary Appendix S2). This was done only for patches where positive spatial relationships with populations of endangered plant species were previously shown.
The cross L function can reveal a spatial correlation between populations and novelty categories, but not if this correlation is positive or negative. We thus used a cross pair correlation function in a successive step (Penttinen et al., 1992;Baddeley, 2008; function pcfcross.inhom of the spatstat package in R) to assess the type of correlation at any given distance between the populations and the different novelty categories. There is a spatial dependence between the populations and a novelty category when the values of the cross pair correlation function lay over 1, with values below 1 indicating spatial repulsion (Penttinen et al., 1992;Baddeley et al., 2015). The higher or lower the values, the stronger the dependence or the repulsion.
All spatial statistics were performed in the R package spatstat (Baddeley et al., 2015). Ripley's isotropic edge corrections (Ripley, 1991) were applied to all calculations. A prior Monte Carlo test for spatial homogeneity proved our point data to be inhomogeneous (the number of points per area is not constant throughout the whole space), so we used statistical functions in "spatstat" that were adapted for inhomogeneous data.

Dispersal Traits and CSR Strategy
To respond to the third research question on the importance of plant dispersal traits and CSR strategy types in explaining the spatial distribution of populations of endangered plant species across the different types of ecosystems, we first gathered information on plant dispersal traits and CSR strategy (Grime, 1977) from the BIOFLOR database (Klotz et al., 2002; Table 1). We used this database because it focuses on the German flora and it achieves a high level of detail by ranking the various dispersal syndromes of each species, allowing us to select the most important one in relation to our study area.
In a second step, we performed a conditional inference tree analysis with the predictors shown in Table 1. A conditional inference tree is a non-parametric type of decision tree where the dataset is recursively split into dichotomous subsets which  are discriminated by the most significant predictor (Hothorn et al., 2006). To analyze if populations attributed to different categories of ecological novelty in each terminal node contribute disproportionally to the different types of ecosystems, we performed a log-linear model using the R-function loglin and computed the standardized residuals to reveal any positive or negative deviations from the expected frequencies of population distribution across the different types of ecological novelty. These deviations were then included as symbols in the resulting nodes in the conditional inference tree, displaying over-or under-representation of populations compared to the expected frequencies. All analyses were performed with the statistical and programming software R version 3.5.2 (R Core Team, 2018).

Spatial Clustering of Populations of Endangered Plant Species
The analysis of the distribution patterns of the endangered plant populations by means of Ripley's L function revealed a spatial clustering of populations in all three novelty categories over the entire range of distances between the populations, as all curves of the Ripley's L function exceed the dashed line that represents the null model of complete spatial randomness (Figure 2). Furthermore, populations that are located in natural remnant ecosystems have a higher level of clustering than those located in hybrid or novel ecosystems (Figure 2).

Spatial Relation Between Populations of Endangered Plant Species and Different Novelty Categories
Results from the cross L function show that from the populations located in anthropogenic ecosystems, only those in hybrid ecosystems are spatially correlated to natural remnant ecosystems ( Figure 3A, green line), as the observed curve of the point process crosses outside of the boundaries of the global envelope of the null model (Baddeley et al., 2014). An additional cross pair correlation function indicates that this correlation is positive (dependence), especially at distances from 0 to 700 m ( Figure 3B, green line), where values lie far above 1 (Penttinen et al., 1992;Baddeley et al., 2015). In contrast, populations located in novel ecosystems are spatially independent from remnant ecosystems (Figure 3A, blue line). Yet, they are spatially correlated to hybrid ecosystems ( Figure 3A, red line). In the latter case, the cross pair correlation function indicates that this correlation is positive (dependence), with a lower intensity, especially at shorter distances from 0 to 400 m ( Figure 3B, red line).

Dispersal Traits and CSR Strategy
From all the predictor variables, the conditional inference tree retained CSR strategy (p < 0.001) and dispersal syndrome (p = 0.001) as significant predictors of population occurrence across the different types of ecological novelty (Figure 4). Endangered species assigned to a ruderal strategy type (r, rc, rs, rcs) and with anemochory as main dispersal syndrome were exceedingly represented in novel ecosystems (Node 3, Figure 4). In contrast, competitive and stress-tolerant species (c, s, cs) were much more common in remnant ecosystems than elsewhere (Node 5, Figure 4). Finally, populations of species assigned to a ruderal strategy type in combination with other dispersal syndromes other than anemochory were represented equally amongst natural remnants, hybrid and novel ecosystems (Node 4, Figure 4).

DISCUSSION
It is well established that plants in cities are not randomly distributed in space, but follow specific spatial patterns, driven by feedbacks between ecological and socio-economic factors (Swan et al., 2011). While many studies have highlighted the importance of natural remnants in hosting species of conservation concern in urban regions (e.g., Godefroid and Koedam, 2003;Celesti-Grapow et al., 2013), recent studies also revealed a high relative importance of anthropogenic habitats for endangered plant species in cities (Kowarik and von der Lippe, 2018;Planchuelo et al., 2019;). Yet to which extent populations of endangered plant species in anthropogenic habitats rely on the adjacency to natural remnants as seed sources is a largely open question, since point data on the exact location of endangered plant species are mostly missing for cities. This study makes a step forward by linking the location of each of 1,742 populations of 213 plant species of highest conservation concern in Berlin to one of three ecosystem types that represent different transformation stages of natural to urban landscapes, i.e., natural remnants, hybrid ecosystems and novel ecosystems, following the novel ecosystem approach (Hobbs et al., 2009). Our methodological approach advances the understanding of spatial patterns of endangered species at the city scale, because the scarce studies on this topic are often based on species inventories of grid cells (e.g., Kühn et al., 2004;Schwartz et al., 2013;Schmidt et al., 2014). While spatial relationships between anthropogenic habitats and natural remnants can be masked by the environmental heterogeneity of grid cells, our study relies on point data for populations and on spatially explicit data from about 80,000 biotope patches covering the entire area of Berlin, which in a previous study were classified as natural remnants, hybrid or novel ecosystems . This enabled us to test for the habitat function of different transformation stages of natural ecosystems in cities for endangered plant species.

Spatial Relationship Between Natural Remnants and Populations of Endangered Plant Species
A first insight of this study is that the total populations of endangered plant species are not randomly distributed throughout Berlin but are overall clustered inside their respective habitats, as commonly happens in highly dynamic landscapes (Johst et al., 2002). Spatial clustering of endangered plant populations was most pronounced in natural remnants. It is likely that clustering of populations in this habitat type merely reflects the rarity and scattered spatial distribution of natural remnants in Berlin. Natural remnants only occupy 7% of the surface area within the city borders but host 24% of all populations of endangered plant species .
As a major and unexpected result of our study we found contrasting spatial dependences of endangered plant populations among the two anthropogenic types of ecosystems and natural remnants ( Figure 3A). Populations in hybrid ecosystems spatially depended on natural remnants -while populations in novel ecosystems did not. As dispersal is the ecological phenomenon that most often causes spatial dependence (Epperson, 2005;Bahn et al., 2008), our results suggest an exchange of species between natural remnants and hybrid ecosystems. This could happen through a source-sink dynamic (Pulliam, 1988;Schreiber, 2010), in which remnant ecosystems are acting as seed sources for adjacent hybrid ecosystems. The importance of natural remnants as relict habitats of endangered plant species is well known (Rebelo et al., 2011;Kowarik and von der Lippe, 2018). Correspondingly, Jarošík et al. (2011) demonstrated that the number of endangered species decreased with distance to natural areas. In Christchurch, New Zealand, a forest remnant within the city functioned as seed source for the colonization of adjacent gardens by native bush species, yet mostly within only 250 m of the forest margin (Doody et al., 2010). Our results corroborate these findings, as the strength of the spatial interaction between natural remnants and populations in hybrid ecosystems was highest in a spatial range between 0 and approximately 400 m (Figure 3B), suggesting that dispersal capabilities prevent many species from directly dispersing further than that (Jacquemyn et al., 2003).
As an alternative but not mutually exclusive explanation, some populations in hybrid ecosystems might be relicts of former remnants that have persisted despite land use change or increasing levels of urbanization. Grasslands, for example, are a largely anthropogenic type of vegetation in Berlin (Sukopp, 1990), but some endangered grassland species stem from habitats in open forests, clearings or forest edges (e.g., Festuca psammophila or Botrychium lunaria) and have survived the transformation of forests to grassland as hybrid ecosystems. To untangle the relative role of dispersal-related versus persistencerelated mechanisms underlying the presence of endangered plants in hybrid ecosystems is a future research direction which requires spatial information about land-use legacies.
Since novel urban ecosystems have experienced much deeper changes to the abiotic conditions and species assemblages (Hobbs et al., 2009;Kowarik, 2011), relict populations or seed banks of endangered plant species likely do not persist in these systems. However, novel ecosystems can offer habitat analogs for a range of native species (Lundholm and Richardson, 2010), i.e., new environmental conditions to which some species may already have been pre-adapted in advance (McDonnell and Hahs, 2015). In Berlin, a considerable amount of 687 populations (40% of the total) encompassing 142 highly endangered species occur in novel ecosystems. One important insight of this study is that in Berlin these populations are spatially independent from natural remnant ecosystems, likely due to dispersal barriers associated with higher levels of urbanization in novel ecosystems. In this context, hybrid ecosystems could partly act as stepping stones between natural remnants and novel ecosystems, as populations in novel ecosystems were slightly dependent on hybrid ecosystems -while at the same time populations in hybrid ecosystems were spatially dependent on remnant ecosystems. However, species exchange between hybrid ecosystems and novel ecosystems alone cannot explain the vast number of populations of endangered plant species in the latter, as populations in Berlin's novel ecosystems are only weakly correlated to hybrid ecosystems, and only so in distances over 400 m. Thus, the high number of populations here could be explained by additional natural or human mediated long-distance dispersal from far away sources and a subsequent establishment of meta-communities within ecologically interconnected novel ecosystems (Loreau et al., 2003). In our dataset, for instance, some species occur across the urban matrix in novel ecosystems that are very distant from each other (e.g., Cystopteris fragils or Dianthus carthusianorum). Conducting genetic analyses to assess whether these populations originate from the same source would be an interesting path for future research.

Strategy Type and Dispersal Syndrome Filter Species Pools of Different Habitat Types
Shifts of biological traits in urban floras compared to rural floras are well established (Knapp et al., 2008a(Knapp et al., , 2010Williams et al., 2015). Our study revealed that biological features (i.e., strategy type, dispersal syndrome) filtered endangered plant species differently for natural remnants, hybrid and novel ecosystems.
As expected, the conditional inference tree revealed that species assigned to a ruderal strategy (component) according to Grime (1977) were more present in novel urban ecosystems than anticipated by chance, while species related to competitive or stress strategies were more confined to remnant ecosystems (Figure 4). This goes along with the review by Williams et al. (2015) who found that many studies report increased numbers of ruderal species in urban floras -despite diverging results existing as well. Our study helps understand ambiguous previous results since different strategy types matter for different types of ecosystems. Thus, the transformation from natural ecosystems to novel urban ecosystems filters for species adhering to the ruderal strategy types, while no specifications were detected for hybrid ecosystems. Natural remnants, in contrast, are characterized by a higher presence of competitive and stress-tolerant species, most probably because of their limited dispersal as a consequence of a trade-off between competition and dispersal (Pellissier, 2015) and the low reproductive effort of stress-tolerant species (Chapin et al., 1993).
As spatial dependence in distribution patterns of plants is closely related to dispersal processes (Epperson, 2005;Bahn et al., 2008), we expected that dispersal-related plant traits filtered the suite of endangered plant species across the different types of ecosystems. A major insight of this study is that dispersal syndrome filtered for species composition of natural remnants versus novel ecosystems. Anemochory (i.e., wind dispersal) was significantly related to the latter, but only for species with a ruderal strategy component (Figure 4). The revealed relevance of wind dispersal for endangered species in novel ecosystems seems to contradict results of Williams et al. (2015) who found that most of the few studies on species extinction in cities showed that wind-dispersed species were more likely to disappear than species dispersed by other vectors (e.g., Williams et al., 2005;Knapp et al., 2010); yet studies on the total urban flora did not find consistent patterns in dispersal modes (Williams et al., 2015). These ambiguous results can be explained by the context dependence of relevant species traits in relation to different urban ecosystems as illustrated by Figure 4, and contribute to the understanding of the spatial independence of populations in novel ecosystems from natural remnants, since wind dispersal often translates to long-distance dispersal (Cain et al., 2000).
While our study revealed the importance of long-distance dispersal by wind for filtering species pools of endangered species throughout different urban ecosystem types, the role of humanmediated long-distance dispersal might be underestimated in our analyses. Indeed, dispersal by humans is considered one of the most important mechanisms of long-distance dispersal (Nathan, 2006), and species with human vectors as an additional dispersal pathway are increasingly found over time in highly urbanized neighborhoods (Johnson et al., 2017). Correspondingly, urban roads function as dispersal highways for a broad range of native species (von der Lippe and Kowarik, 2008), and likely more for anthropogenic than for remnant ecosystems as the former are usually more exposed to transportation corridors. Given that vehicles can move native species over long-distances along urban-rural gradients (von der Lippe and Kowarik, 2008), human-mediated dispersal might thus link novel with other ecosystem types. Furthermore, high seed production, which is a feature of many ruderal species (Grime, 1977), increases the chance of dispersal by vehicles, especially for seeds with a low terminal velocity like those of wind dispersed species (von der Lippe and Kowarik, 2012). In consequence, seeds of wind dispersed species might be introduced from distant sources (i.e., hybrid ecosystems) into novel ecosystems through the road network, while the same mechanism could support meta-communities of species already established in novel ecosystems. Work on pollen dispersal correspondingly indicates a high level of exchange between different types of urban ecosystems (Van Rossum, 2010;Van Rossum and Triest, 2012).

Future Perspectives
The present study suggests that meta-communities of plants might exist in novel ecosystems within cities. A future direction is to examine the spatial structure and dynamics of populations in novel urban environments (and probably related sources) also with genetic analyses.
On a different note, our research shows that differences in interspecific plant traits are related to the distribution of endangered plants across different kinds of urban ecosystems. Previous research has shown rapid changes and evolutionary adaptations in some urban plants (McDonnell and Hahs, 2015;Dubois and Cheptou, 2017), and that trait frequencies can shift as a consequence of different levels of urbanization (Knapp et al., 2008b) or urban stressors (Williams et al., 2015). Thus, it would be most interesting to assess patterns of intraspecific trait adaptations throughout populations of the same species located in ecosystems with varying levels of ecological novelty.
Finally, though Berlin shares similar land use legacies (Ramalho and Hobbs, 2012) and configuration of urban elements (Louf and Barthelemy, 2014) with other historical European cities, results of this study cannot be generalized for all urban contexts. The classification of the urban ecosystems according to their ecological novelty is a first step in providing a general scheme to compare species distributions between cities of different history, and comparable studies in other cities could be used to verify to which extent biodiversity patterns vary amongst them.

CONCLUSION
The urban contribution to biodiversity conservation is widely debated (Shwartz et al., 2014;, with one important question being the role of urban ecosystems of anthropogenic origin for endangered plant species. While previous studies have revealed an important relative contribution of human-shaped hybrid and novel ecosystems as habitats of endangered plant species (Kowarik and von der Lippe, 2018;Planchuelo et al., 2019) this study adds further insights into the underlying mechanisms, based on fine-grained georeferenced data for a high number of populations of endangered plant species and urban biotope patches.
The spatial analyses allow three major conclusions to be drawn. First, the spatial dependence of populations in hybrid ecosystems on natural remnants shows the latter are important habitats and seed sources for endangered plant species and highlights the importance of conserving these remains of natural landscapes in urban regions. Second, the missing spatial dependence of populations in novel ecosystems demonstrates that novel urban ecosystems can support populations of a range of endangered species beyond the adjacency to natural remnants. Third, our results suggest that hybrid ecosystems may act as essential stepping stones by mediating the spread of endangered plant species into novel urban ecosystems.
The existence of natural remnants in cities can partly be explained by the fact that urban agglomerations have preferably been founded in naturally rich landscapes (Kühn et al., 2004). Yet at least in Berlin, these legacies do not ultimately determine the spatial patterns of endangered plant species in profoundly changed novel urban ecosystems. This is a strong argument for including the total range of natural and anthropogenic ecosystems in urban conservation strategies.
Discerning different biological features (i.e., strategy type, dispersal syndrome) related to the presence of endangered plants in natural remnants, hybrid ecosystems and novel ecosystems demonstrates that different filters modulate the presence of endangered plants across the trajectory of urban habitat transformation. Identifying traits of species that successfully colonize novel urban ecosystems demonstrates environmental barriers for other endangered species. This has important implications for conservation efforts, e.g., for conserving and managing those species with highest priority in natural remnants that are likely unable to colonize novel urban ecosystems; or by strictly selecting dispersal-limited species for introduction to novel ecosystems in restoration measures.
Understanding the spatial relationships between populations of endangered plant species in different urban ecosystems thus supports conservation approaches that cover the total range of urban ecosystems and associated opportunities for biodiversity conservation in cities.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
All persons entitled to authorship have been so named. All authors conceived the ideas, designed methodology, and contributed to writing. GP processed the data. GP and ML analyzed the data.