Biologic interactions determining geographic range size: a one species response to phylogenetic community structure

Range size variation in closely related species suggests different responses to biotic and abiotic heterogeneity across large geographic regions. Species turnover generates a wide spectrum of species assemblages, resulting in different competition intensities among taxa, creating restrictions as important as environmental constraints. We chose to adopt the widely used phylogenetic relatedness (NRI) measurement to define a metric that depicts competition strength (via phylogenetic similarity), which one focal species confronts in its environment. This new approach (NRIfocal) measures the potential of the community structure effect over performance of a single species. We chose two ecologically similar Peucaea sparrows, which co-occur and have highly dissimilar range size to test whether the population response to competition intensity is different between species. We analyzed the correlation between both Peucaea species population sizes and NRIfocal using data from point counts. Results indicated that the widespread species population size was not associated with NRIfocal, whereas the population of restricted-sized species exhibited a negative relationship with competition intensity. Consequently, a species' sensitivity to competition might be a limiting factor to range expansion, which provides new insights into geographic range analysis and community ecology.


Introduction
Biologic interactions are fundamental to the understanding of ecological patterns and processes. Webb et al. (2002) set the framework of phylogenetic community ecology to elucidate the effect of ecological interactions depicted via phylogenetic structure, that is, the degree of relatedness among species forming an assemblage. This approach is strongly supported by the tested observation of the inheritance of niche-related traits from ancestors, termed niche conservatism (Wiens et al. 2010; for birds see Lovette and Hochachka 2006), which reinforces the concept of a link between phylogenetic distance and ecological similarity reflected in behavioral (Houle 1997) and life-history traits (Burns and Strauss 2011). Consequently, demonstrating statistical support for phylogenetic and community assemblage patterns (i.e., clustering or evenness) suggests a potential process exists structuring species assemblages (Webb et al. 2002). Questions addressed under this method have resolved several issues related to ecology (ecosystem stability Cadotte et al. 2012; species lost and climate change, Willis et al. 2008), providing useful insights into the internal structure of phylogenetic and ecological relationships (Lovette andHochachka 2006, G omez et al. 2010). For example, Graham et al. (2009) demonstrated that biologic interactions among hummingbird species represented the leading factor in tropical lowland community assembly processes, even more important than environmental variation (i.e., habitat filtering). However, a bond must exist between community phylogenetic structure and the performance of each species, a fact that is largely overlooked. The fitness of particular species might consequently be explained by the composition of the community in which the species is found. Furthermore, Ricklefs (2004) reported that community level processes generated population level changes, driving current ecological patterns. Yang et al. (2013) evaluated the phylogenetic diversity (PD) of assemblages surrounding target individuals, showing that most individuals had a neutral tendency regarding the PD of adjacent plots. However, by exploring the phylogenetic distance of each individual in a focal species, the direct effect of the surrounding community structure may be measured. For instance, Jiang et al. (2010) designed an experiment to assess the success of an invader species in bacterial communities, based on phylogenetic distance between invader and native species. The authors found a positive relationship between phylogenetic distance and the probability of a species becoming established.
The central role of biotic interactions is considered contingent on a species geographic range. For example, Brown et al. (1996) reported that biotic interactions tend to limit the distribution and abundance of species at lower latitudes. Differences in a species geographic range size do not just suggest variation in response to environmental variables or niche breadth (Gaston and Spicer 2001), but can reflect a species response to biologic interactions, which were illustrated in classic experiments by Connell (1983), and more recent studies by Bullock et al. (2000). However, previous studies reported the influence of biologic interactions in two-species systems distributions, without evaluating the effects of the entire community. Recently, Villalobos et al. (2013) introduced a novel and interesting approach, in which the phylogenetic structure of species co-occurrence of a focal species is used to study broad coexistence patterns.
We hypothesize that the sensitivity of species confronting negative interactions is reflected in the species population attributes: species inhabiting different assemblages, consequently experience different levels of competition throughout the species geographic range, which results in different rates of change in a species overall fitness. Gaston (2009) indicated that population size is the outcome of several population structure attributes (e.g., levels of births, deaths, and migration). For instance, population density has been linked to species richness, producing higher density in areas where richness is low (i.e., density compensation; MacArthur et al. 1972), which is a pattern that was first described for island systems compared with mainland systems. Under these conditions, it is expected that populations of a species that occupy large geographic ranges (widespread) are not as influenced by co-occurring species with which they compete; alternatively, populations of species exhibiting restricted geographic ranges are more affected by increased potential competition. Support for range size heritability (Waldron 2007; but also see Webb and Gaston 2005) facilitates the expectation that related species would have similar range size; however, differences in the range of closely related species might serve as a viable system to test whether this dissimilarity is provided by a differential response to competition. The present study included two components: (1) we modified a widely used metric of phylogenetic structure (NRI; Webb et al. 2002) to center the attention toward a focal species and (2) two sympatric species in the genus Peucaea (Emberizidae), which exhibit very dissimilar range sizes, were used to evaluate whether a population size response to potential competition (through the modified metric) differs between the two sparrow species.

Peucaea sparrows and fieldwork
We conducted this study in southeast Mexico, in a region called the Tehuantepec Isthmus (Huidobro et al. 2006). This region is located in the narrowest stretch of land between the Gulf of Mexico and the Pacific Ocean, which is represented by the municipality of Juchit an, state of Oaxaca. Two sparrow species co-occur in the area: Peucaea ruficauda and P. sumichrasti. The former is a more widespread species (2.6 9 10 5 km 2 ; see next section) compared with P. sumichrasti, which is one of the most range-restricted avian species in Mexico (9.7 9 10 3 km 2 ; Wolf 1977). This species is endemic to the Tehuantepec Isthmus, and its entire range overlaps with P. ruficauda. Both sparrow species are common birds in the region and are similar in shape, behavior, and ecology (Wolf 1977). We selected 17 monitoring sites across the region, which had comparable vegetation type, human perturbation, and size, covering almost the entire geographic range of P. sumichrasti. The species composition (land birds) of each assemblage and Peucaea sparrow abundance were described by sampling 24 fixed-radius point counts separated by 200 m to avoid double counting individuals at each monitoring site. Each point count was sampled eight times by the same team of observers for one year, with each monitoring site being visited every 6 weeks on average.

Phylogeny and geographic range size
We conducted a phylogenetic reconstruction, which included all land bird species we observed during our fieldwork. We queried the GenBank database (NCBI, September-October 2011) for the mitochondrial COI gene (Hebert et al. 2003;Alif et al. 2011; Appendix A1), representing each species we identified from our 17 monitoring sites. Sequences of the species not deposited in GenBank were replaced with an available sequence from the closest relative, according to previous studies reported in the literature. Sequence alignments were made in Clustal W (Larkin et al. 2007), a BioEdit 7.0.9.0 (Hall 1999) accessory tool. The reconstruction was performed by Bayesian method using BEAST v1.7.5 (Drummond and Rambaut 2007). Nucleotide substitution model employed was HKY+G+I being identified as the more adequate by JModelTest 0.1.1. The Yule speciation process was set to model the tree prior. In order to calibrate the root node of the tree, we used the date from Jetz et al. (2012) for the divergence of Anatidae (lognormal mean 4.27, SD 1, zero offset 0). One independent 5,000,000 generation run was performed sampling at every 1000 generations. The outcome was analyzed in TreeAnnotator v1.7.5 discarding the 10% of trees and visualized in Mesquite v2.74 (Maddison and Maddison 2011)  We looked for range size phylogenetic signal as an Emberizidae family trait by conducting a randomization test (999 randomizations; Blomberg et al. 2003) using the comprehensive phylogeny of Emberizidae reported by DaCosta et al. (2009), and the range size of 49 sparrow species. The species distributions area calculations were generated in ArcGIS 10 and shapefiles from the Nature-Serve compilation (Ridgely et al. 2003) using the appropriate geographic coordinates regarding species distributions in North or Central America. Phylogenetic signal is detected when random distributions exhibit significant differences from observed values. Although similar geographic range sizes have been observed among close relatives in birds (Waldron 2007), we decided to test this hypothesis in the Emberizidae family because the outcome was scale dependent.

Data analysis
The traditional metric of the phylogenetic community (Net Relatedness Index -NRI-) is a standardized measure of the mean pairwise phylogenetic distance (MPPD, Webb et al. 2002), which is the phylogenetic distance among all possible pairs of species within a community. We modified the MPPD, so that our metric (MPPD focal ) did not reflect the distance among all pairs, rather it measured the distance from the focal species and each species included in the assemblage, and averaged those values (Fig. 1). The calculation was performed using the cophenetic distance between focal species and all others members of the community, from which the mean was calculated. The species abundances are included in the weighed version of this metric by calculating the weighed mean instead of the arithmetic mean; in this way, the metric reflects the real composition of the community. Like the NRI, we constructed a standardization named NRI focal , in which the observed MPPD focal values were compared with null distributions that were generated by creating communities of identical size by random draws from species pool (Kraft et al. 2007). In a single value, the NRI focal describes focal species relatedness and the set of species that co-occur with the focal species. In this way, NRI focal may be defined as the phylogenetic (ecological) similarity of Peucaea with all other observed bird species and serves as a measure of potential competition (Fig. 1).
We calculated NRI focal (both weighed and nonweighed) for P. ruficauda and P. sumichrasti for each assemblage and performed a Pearson′s product-moment correlation between NRI focal values and sparrow abundance. Peucaea abundance was normalized by root square transformation (Sokal and Rohlf 1987). The modified metric was based on the "mpd" function implemented in the PICANTE package (Kembel et al. 2010 Figure 1. The NRI focal calculation using a phylogenetic tree and its phylogenetic distance matrix. Notice that the difference with traditional NRI appears in the MPPD calculation, where the traditional is estimated by averaging the distances among all the possible pairs, while MPPD focal is estimated by averaging the distances among the focal species (species C in this example) and each species included in the assemblage. communities, measuring both the phylogenetic metrics of each assemblage and calculating the correlation between them.

Results and Discussion
The geographic range among Emberizidae sparrows represents a trait with significant phylogenetic signal (P < 0.001), indicating that similar species have a similar range size (Fig. 2). However, the difference in the geographic range size of sympatric and closely related Peucaea sparrows is noticeably large (2.5 9 10 5 km 2 ). This range size difference is of interest because of the ecological similarity between the sparrow species and the apparent lack of any form of geographic barrier; hence, this scenario may be explained by biologic interactions. We found that the widespread species (P. ruficauda) was more abundant compared with P. sumichrasti in 80% of the species assemblages. The relationship between range size and local abundance has been well established in several taxa and was demonstrated in passerine birds (Bock and Ricklefs 1983). Several mechanisms have been proposed to determine the positive correlation between range size and local abundance (revised and discussed by Gaston et al. 1997); for example, Holt et al. (1997) conferred special importance to among-species differences based on species differential responses to density-independent factors influencing population attributes. Here, results emphasized that biologic interactions were an important element in this process. For the resident species assemblages, the number of P. sumichrasti individuals showed a negative correlation with NRI focal (r = À0.592, P < 0.05; Fig. 3) when considering all species abundances (weighed NRI focal ). This suggests that competition strength (implied by increasing phylogenetic similarity) limits abundance in P. sumichrasti. The competition strength approached by the phylogenetic similarity leads the population volume in this restricted-sized species. Consequently, the community structure (composition of the species assemblage) is an important factor affecting one species population size. This result is analogous with that reported by Jiang et al. (2010), who found a positive relationship between the phylogenetic distance of invaders to bacterial communities and invader abundance. While the context of this preceding study differs to ours, it is worth noting the similar outcomes, even for very different organisms and differently sized communities. However, we did not find evidence of a relationship between widespread P. ruficauda abundance and NRI focal (Table 1). Dissimilarity in range size among ecologically similar species might be related to differences in sensitivity to competition among species. This difference might be linked to the categories defined by Yang et al. (2013) in a lightly similar context, who grouped phylogenetic diversity into neutral, repeller, and accumulators species. Sensitivity to competition as a trait has the potential to limit species distributions across large geographic regions, because a large area exhibits high species turnover; consequently, several communities may form levels of competition so great that a sensitive species could not compete successfully. Consequently, regions at low latitudes with high species turnover can be defiant to susceptible species producing small-sized range species and contributing (besides other factors; Arita et al. 2005) to the observed and proved pattern of Rapoport′s rule. Yang et al. (2013) identified only a few phylogenetic diversity repellers, which might indicate that few sensitive species produce a reduced number of species with small-sized ranges and elevated numbers of species with large-sized ranges. However, the opposite pattern is obtained in reality; whereby, many species have small-to moderate-sized ranges, while only a few have very large ranges (Brown et al. 1996). Thus, the complex interaction between biologic interactions and environmental variation is reaffirmed. The internal distributional range structure (sensu Brown et al. 1996), and the variability in species abundance throughout a species distribution (Brown et al. 1995), might provide a way of elucidating the response of species to potential competition across different assemblages. In even a relatively small area (Tehuantepec Isthmus), the observed sites showed enough species turnover (see Appendix A2) to depict dissimilar phylogenetically structured assemblages, which suggested a gradient of potential competition (Fig. 4). Moreover, the geographic boundaries of sensitive species might have been established by communities in which ecological (i.e., phylogenetic) similarity exceeds the level of similarity (i.e., competition) that a focal species could withstand. In other words, certain assemblages function as "stakes," limiting species geographic ranges. These stakes acted as biologic barriers setting a threshold that could not be trespassed by sensitive species. For instance, the P. sumichrasti population size reached low levels when the community phylogenetic similarity was high, preventing the dispersal of individuals to new areas. Although some individuals might cross these stakes, their numbers might be reduced, due to their being unable to establish stable populations. Areas where phylogenetic similarity is low would facilitate the free transit of sensitive species, until another type of barrier (i.e., environmental or biologic) is reached. Even though our monitoring sites were distributed throughout the geographic range of the restricted species, it would have been more accurate to evaluate our hypothesis by identifying the exact boundary of the range to assess the presence of overwhelming competition to which sensitive species are subjected. These boundaries are not distinct for birds; however, it might be viable to apply such studies that exhibit lower levels of movement to improve the experimental design.
In our study, the specific assemblage composition alone was not enough to explain the variability in population size among sites; the positive relationship between P. sumichrasti abundance and the nonweighted NRI focal was not statistically supported (P > 0.05). Communities might appear similar to each other when species' abundances are neglected. In the Appendix (A2), we show  that the dissimilarity among communities was low when only taking species composition into account, but become highly dissimilar when abundance was taken into account.
Thus, our findings are based on species turnover, which is only evident when it is analyzed through abundance. Individuals of one species interacted with individuals of  several species, and the level of competition differed in each interaction. We should think that individuals from one species are interacting with individuals of other species, rather than assuming that individuals interact with species; thus, abundances should be included in the analysis of phylogenetic structure whenever such data are available (Norden et al. 2012).Therefore, competition intensity among organisms must be weighed by phylogenetic distance. Pragmatically, when the scale and data are appropriate, patterns and processes are, respectively, revealed and inferred by weighing phylogenetic distance by the number of individuals. The analyses were conducted using Peucaea population size as the species response; nonetheless, competition strength represented by the entire community could affect other population parameters, including birth, recruitment, and fecundity rates, among others. Gaston (2009) reviewed the population structure attributes necessary as a framework to continue research on the effects of community strength on population range dynamics, and methods similar to the present study show promise. The phylogenetic community ecology approach as a means of studying a single species clarifies patterns and processes that otherwise might be masked. Although NRI focal and NRI are related methodologies (r = 0.626, P < 0.05), the metrics told different stories about the same assemblages. The NRI focal we applied here provided a new method to evaluate the potential community effect over the focal species. Both metrics are superficially similar and derived from the same theoretical framework; nonetheless, we caution that metric application and results differ and should be applied under different contexts.