- Open Access
Geographic range size and speciation in honeyeaters
BMC Ecology and Evolution volume 22, Article number: 86 (2022)
Darwin and others proposed that a species’ geographic range size positively influences speciation likelihood, with the relationship potentially dependent on the mode of speciation and other contributing factors, including geographic setting and species traits. Several alternative proposals for the influence of range size on speciation rate have also been made (e.g. negative or a unimodal relationship with speciation). To examine Darwin’s proposal, we use a range of phylogenetic comparative methods, focusing on a large Australasian bird clade, the honeyeaters (Aves: Meliphagidae).
We consider the influence of range size, shape, and position (latitudinal and longitudinal midpoints, island or continental species), and consider two traits known to influence range size: dispersal ability and body size. Applying several analytical approaches, including phylogenetic Bayesian path analysis, spatiophylogenetic models, and state-dependent speciation and extinction models, we find support for both the positive relationship between range size and speciation rate and the influence of mode of speciation.
Honeyeater speciation rate differs considerably between islands and the continental setting across the clade’s distribution, with range size contributing positively in the continental setting, while dispersal ability influences speciation regardless of setting. These outcomes support Darwin’s original proposal for a positive relationship between range size and speciation likelihood, while extending the evidence for the contribution of dispersal ability to speciation.
How speciation rates vary among geographic regions and clades remains one of evolutionary biology’s most pivotal questions [1,2,3]. The influences on speciation of climate [4, 5], ecological opportunity [6, 7], latitudinal position [8, 9], traits , and geographic range size [11, 12] have been widely investigated. Yet conflict remains between theoretical expectations and empirical data. Particularly, the role of geographic range has yet to be resolved. Theory suggests that geographic range size should have a positive influence on the likelihood of speciation [1, 13, 14], although the form of this relationship may be influenced by the mode of speciation [15, 16]. Species with larger range sizes are expected to have a higher chance of speciation because of the larger area in which geographic barriers may form [17, 18]. Larger ranges tend to be associated with larger population sizes and genetic variability, and are more robust to localised environmental disturbances, also lowering the chance of extinction [11, 19,20,21].
Empirical investigations have found varying evidence for both positive [11, 19] and negative [10, 12, 21, 22] relationships between range size and speciation rate. Some theoretical arguments suggest that the species with the largest ranges may have ranges large enough to engulf barriers [18, 23], reducing speciation. However, smaller ranges are associated with higher extinction risk, also reducing speciation probability [20, 24, 25]. Thus, a further proposal has been made that geographic range size and speciation likelihood should show a unimodal relationship, in which intermediate ranges have the highest rate of speciation [14, 18, 23, 26].
Although current empirical support for a unimodal relationship is limited  and the idea has been described as unsatisfying , it remains appealing from a theoretical perspective because it also considers traits that influence the relationship between range size and speciation, such as dispersal ability and body size. Dispersal ability can either limit or promote gene flow, therefore affecting the probability of population divergence [22, 29,30,31], though this can depend on the mode of speciation . Greater dispersal ability is associated with large range sizes [33,34,35], but is also thought to reduce speciation probability [19, 31, 36], complicating assessment of its effects on the relationship between geographic range size and speciation likelihood. Moreover, both range size and dispersal ability are also influenced by body size [33, 37], which in turn has an independent influence on speciation probability [38, 39] (Fig. 1).
Further complexity arises because the shape and geographic position of a species range can influence both range size and speciation likelihood. Shape can impact speciation because narrow or elongated ranges are more likely to give rise to dispersal events  and are thought to be prone to bisection by geographic barriers [13, 18], thereby increasing speciation probability, especially when geographic barriers are small . Geographic position may act as a constraint on range size and is related to broader geographic gradients in both range size and body size. Acknowledging variation among taxa, range size and body size tend to decrease toward the tropics [41,42,43] (Fig. 1). Furthermore, species found in large biogeographic regions with a large spatial extent can have larger ranges than those found on small areas such as some islands [44, 45]. Moreover, island species derived from a continental relative are unlikely to recolonise the mainland  which can limit opportunity for range expansion in island species.
Any study seeking to understand the relationship between geographic range and speciation likelihood or rate has to take the interrelationship between traits and biogeographical constraints into account. The overall form of the relationship between geographic range size and speciation rate often emerges as weak, and often weakly negative [e.g. 10]. This could be the outcome of interactions between interrelated drivers of speciation, variation among taxa in traits such as body size and dispersal ability, and variation in biogeographic constraints between different major regions of the world (Fig. 1). One useful approach involves focusing on a large monophyletic group in order to understand the underlying causes and drivers of speciation at a range of scales [47,48,49].
We applied such an approach, using honeyeaters (Aves: Meliphagidae), to investigate the influence of geographic range size on speciation rate, taking into account the influences of range position, dispersal ability and body size. The honeyeaters are a large (192 species), diverse, monophyletic clade of passerine birds that are currently restricted to the east of Wallace’s line (except for one species found in Bali) and have diversified throughout Australasia in the past 25 million years [50,51,52,53]. Ranging from arid Australia to tropical Indonesia, the honeyeaters have diverged to fill a wide variety of niches and have considerable diversity in body size, dispersal ability and range size [50, 51]. The most common mode of speciation in birds is allopatric speciation in the broadest sense [17, 54]. This is evident in honeyeaters [55,56,57], and dispersal events throughout Indonesia and the South Pacific are thought to have played a prominent role [53, 58], indicative of peripatric speciation. Honeyeaters are well studied and described; distributions are known, and trait data are prevalent throughout the literature. Moreover, honeyeaters are distributed across both a large continent and many islands, and are monophyletic: ideal for examining the geographic range size idea in the expectation also of an influence of predominant modes of speciation. Here, we examine the hypotheses that: (i) the relationship between geographic range size and speciation rate is positive, and (ii) that the form of the relationship is likely to be influenced by the mode of speciation, while accounting for other factors that are likely to influence the geographic range size-speciation rate relationship (Fig. 1a).
Phylogenetic reconstruction and range size
We included ultra-conserved elements, mitochondrial and nuclear loci, alongside topological constraints from previous phylogenomic studies to generate a honeyeater phylogeny of 192 species (Fig. 2). The topology of genera within the tree is concordant with Andersen et al.  and relationships within the genus Meliphaga follow McCullough et al. . In general, relationships between species within genera broadly follow previous reconstructions (; see Additional file 1: Fig. S4 for a comparison).
The species with the smallest range size (42.4 km2) is the Rotuma myzomela (Myzomela chermesina), from Rotuma Island, Fiji. By contrast, the species with the largest range, the singing honeyeater (Gavicalis virescens), has a geographic range that spans nearly all of Australia (7,306,701 km2; mean honeyeater range size = 650,335.2 km2, median honeyeater range size = 78,121.45 km2).
Phylogenetic generalised least squares regressions revealed a positive relationship between range size and speciation rate (p = 0.01; Fig. 3a; Additional file 2: Table S1). In this model, speciation rate also increases with dispersal ability (p = 0.002) and declines with increasing latitudinal midpoint (p = 0.02; Fig. 3a). Species with less elongated or disjunct ranges are associated with higher rates of speciation (p = 0.02; Additional file 2: Table S1).
Geographic range size decreases with body size (p = 0.03; Fig. 3b). Range size and dispersal ability are unrelated, but dispersal ability and range shape show a positive relationship (p = 0.02; Additional file 2: Table S1). Range position is strongly correlated with range size, shape and dispersal ability (Additional file 1: Figs. S5–S11; Additional file 2: Table S1).
The geographic setting term, islands/continent, had a significant interaction with range size (p = 0.019), but not with dispersal ability (p = 0.124) or body size (p = 0.145; Additional file 2: Table S1). In the continental species only, higher rates of speciation are associated with larger ranges (p = < 0.0001; Additional file 1: Figs. S12–S15; Additional file 2: Table S1)
Bayesian phylogenetic path analysis
Here, the only trait found to influence speciation rate is dispersal ability: species with greater dispersal abilities have higher rates of speciation. Range size is influenced by range position, with range size increasing with latitude and decreasing with longitude, and with continental range sizes larger than those on islands (Additional file 1: Fig. S16; Additional file 2: Tables S2, S3). Relationships between dispersal ability, range position, shape, and body size were in keeping with the phylogenetic regressions (Additional file 2: Table S1).
Due to the influence of islands on species range size, analysis was repeated separately on island and continental species. In island species speciation rate is related to high dispersal ability and small body size (Fig. 4a; Additional file 2: Tables S4, S5). By contrast, in the continental species, higher speciation rates are significantly associated with larger range sizes and greater dispersal ability (Fig. 4b; Additional file 2: Tables S4, S5). In both island and continental species, range size increases with latitude and decreases with longitude. Body size is also found to be influenced by position, decreasing with latitude and increasing with longitude (Fig. 4).
The spatiophylogenetic model revealed considerable spatial and phylogenetic structure to the variation in speciation rate (Fig. 5). Speciation rates are variable across both the continental and island settings, as the previous analyses indicated, and range size retains a positive relationship with speciation rate. So does dispersal ability, and the negative relationship with latitudinal midpoint, indicating higher rates in the tropics, is also retained (Fig. 5c; Additional file 2: Table S6).
State-dependent speciation and extinction analyses
The state-dependent speciation and extinction models verified the influence of range size on speciation rate, though at first the outcomes seem to contradict the previous model outcomes. The FiSSE model closest to significance for range size identified small ranges less than 100,000 km2 as having higher rates of speciation than range sizes larger than this value (λsmall range = 5.96, λlarge range = 4.61, p = 0.058; Additional file 2: Tables S7, S8). Four HiSSE models, with Akaike weights 0.66, 0.067, 0.056, and 0.051, respectively, occupy the bulk of the model weight (wi = 0.834), with the remainder of the weight associated with a CID-4 model (wi = 0.143; Additional file 2: Tables S9, S10). Importantly, however, the FiSSE models also confirmed that island species have higher speciation rates than their continental counterparts (λislands = 6.47, λcontinents = 4.68, p = 0.049; Additional file 2: Tables S7, S8). This was supported by the HiSSE modelling, where the best HiSSE model (wi = 0.952) outperformed the four CID models which were included as null models (All have wi ≤ 0.001; Additional file 2: Tables S9, S10). Differences in the range size distributions of the species on islands and the continent provides the explanation (Fig. 2c). When classified on a binary basis, ~ 96% (65 of 68) of the island-dwelling species have range sizes less than 100,000 km2 in extent, whereas ~ 72% (89 of 124) of the continental species have ranges larger than this threshold.
Such a finding is in keeping with the phylogenetic regressions, spatiophylogenetic modelling and the structural equation models demonstrating the importance of separating the continental from the island settings. If the differences between island and continental settings are significant, then a state-dependent speciation and extinction analysis differentiating large and small range-sized species from those with intermediate-sized ranges should reveal an effect with the former having the higher values. A FiSSE model using the quartile range to make such an adjudication (i.e. grouping the lower and upper quartiles) is indicative that this may be the case (FiSSE: λintermediate range = 4.92, λsmall + large range = 5.71, p = 0.081; Additional file 2: Tables S7, S8). Perhaps more noteworthy, the HiSSE model adopting the same classification outperforms (wi = 0.688) the CID models (best CID model wi = 0.107; Additional file 2: Tables S9, S10).
Theory suggests that larger ranges should have the highest speciation rates, with these relationships being potentially modified both by the mode of speciation and by other contributing factors and species traits [18, 34, 36, 37]. The latter, such as dispersal ability and range spatial position [12, 35, 40], may have independent effects on both range size and the likelihood of speciation [30, 36]. When this complexity is taken into account by applying several analytical approaches, the positive relationship between range size and speciation rate is supported, and a significant interaction found between range size and islands indicates that the patterns vary depending on geographic setting, which likely reflects mode of speciation. Specifically, honeyeater speciation rate differs considerably between islands and the continental setting across the clade’s distribution, with dispersal ability influencing speciation rate regardless of setting and range size contributing positively in the continental setting. The higher rates in the tropics are, in part, a reflection of the differential distribution of the continental and island settings. Therefore, in the Meliphagidae it is clear that dispersal associated speciation across islands and range size related allopatric speciation in the continental setting have both contributed to diversification. Moreover, we have been unable to reject the two hypotheses we sought to test, lending support both to theory and to previous findings for the honeyeaters.
In the case of the first hypothesis (a positive relationship between range size and speciation likelihood) [1, 17, 18], our results support the idea that, where vicariant speciation is likely to be the dominant mode of allopatric speciation, a positive relationship with range size is to be expected. Darwin  originally proposed that larger ranges would have a greater chance of allopatric speciation events. These ideas have persisted [13, 23] and empirical studies have found positive relationships between range size and speciation [11, 19]. By contrast, our results do not support the hypothesis of a largely negative or unimodal relationship between range size and speciation [10, 27].
Among Australian continental vertebrates, including birds, patterns of vicariant speciation associated with changing climates in particular, but also landform evolution, have played a major role in species diversification [61,62,63]. Diversification patterns and their outcomes are similar among mammals, amphibians, and birds  and our findings for continental honeyeaters are congruent with these patterns (Fig. 5). Examples include high speciation rates of honeyeaters on the east coast, with a hotspot around the Atherton Tableland in eastern Queensland (Fig. 5). Additionally, environmental heterogeneity and aridification of Australia has resulted in increased diversification of arid-zone clades in vertebrates . The increasing aridification of Australia has also resulted in shrinking mesic biomes, restricting species ranges, driving extinction and vicariant speciation. This environmental heterogeneity and the shifting of biomes has had widespread effects and facilitated allopatric speciation in mammals , reptiles , fish , and amphibians  across Australia. Although we have not investigated the specific historical environmental drivers of speciation, our results appear to be consistent with these patterns (Fig. 5) and previous investigations of barriers influencing honeyeaters [56, 57].
Nonetheless, specific identification of vicariant speciation as important for the honeyeaters is new. Honeyeaters are well established as good dispersers and have longer and more projected wings compared to their meliphagoid relatives [29, 70]. Vicariance events have been suggested to have acted in this group [53, 71], but to date, evidence for speciation in honeyeaters has predominantly supported the role of dispersal events throughout Indonesia and the South Pacific [53, 55, 72], associated with their nectarivorous lifestyle [35, 70]. Although we confirm the results of previous honeyeater studies which have suggested that dispersal has been key to honeyeater speciation especially in island species [53, 55], we find that vicariant patterns associated with a large range size mechanism in the continental setting is also important. Differences between the continental and island settings are supported by all the analytical approaches, including the state-dependent diversification models. In both settings, dispersal is an important contributor to diversification rate, but in the continental setting, range size contributes too.
The finding of smaller ranges associated with higher rates of speciation is consistent with previous studies that have investigated range size and speciation in birds [10, 12]. This result is, however, at least partially driven by island endemism. Higher rates of speciation are evident in species from the genera Philemon, Myzomela, and Lichmera, all of which include a greater proportion of island than continental species (Fig. 2). Such patterns of island speciation became apparent when plotting occurrence records coloured according to speciation rate (Fig. 5a), and small islands such as those throughout the Lesser Sunda Islands and the Banda Arc are still visually associated with higher rates of speciation in the predicted speciation rate map from the spatiophylogenetic analysis (Fig. 5b). This finding of islands as speciation hotspots is not surprising given that islands frequently play key roles in avian speciation [e.g. 12, 73, 74]. The island setting is clearly characterised by dispersal-associated speciation. High dispersal ability typically improves the ability of species to reach oceanic islands, increasing rates of diversification, which frequently play key roles in avian speciation , as found for the honeyeaters [53, 55, 72].
Latitudinal gradients in speciation rates have been well investigated and the subject of much theoretical discussion . Most recently, empirical studies have often uncovered patterns in which speciation rates decrease towards the tropics [8, 9]. Although we find evidence for a latitudinal gradient in honeyeater speciation rates, suggesting higher rates in the tropics, the outcome appears to have more to do with the island versus continental setting than a tropical versus extratropical setting, though obviously the two cannot entirely be distinguished. The outcome is wholly in keeping with the geographic complexity that may be expected from the differences between clades and environmental settings and their interactions . It also highlights the benefits of a focus on such settings and specific taxa for distinguishing the underlying likely influences on macroecological and macroevolutionary patterns .
Inferences of modes of speciation from current ranges should be cautious, as current ranges may not always be representative of past evolutionary history. However, simulations have demonstrated that even if species ranges have changed over time, speciation history is often detectable from current distributions .
Our results indicate that the relationship between geographic range size and speciation rate is positive and is influenced by mode of speciation. However, relationships between other factors are likely to influence the geographic range size-speciation rate relationship, including dispersal ability, and accounting for covariance due to shared biogeographic occupancy is also critical. Our ability to understand the proximal basis of speciation rate variation across the Tree of Life can be greatly aided by elucidating these complex multi-trait relationships.
We used recent phylogenomic analyses [59, 60] in conjunction with traditional nuclear and mitochondrial markers  to construct a comprehensive phylogeny of honeyeaters. We followed the IOC world bird list (version 10.2; ), which recognises 191 species of honeyeater and follows recent taxonomic revisions suggested by genomic work [60, 69], including the recently described Myzomela species [78,79,80]. We included Myzomela obscura rubrotincta as a distinct species, Myzomela rubrotincta, due to previous work suggesting that this represents a distinct species based on morphology and vocalisations  and molecular work finding this species has been more closely related to M. cruentata than to M. obscura .
All available DNA sequences for Meliphagidae species were downloaded from GenBank. This included sequences for eight nuclear genes (Fib-5, GAPDH, RAG-1, RAG-2, FIB-B17, c-mos, BDNF and TGBF2) and five mitochondrial genes (12 S, cyt-b, COI, ND2 and ND3). The coverage for each of these genes varied from 3 species to 187 species (Additional file 2: Table S12). We included the following nine taxa as outgroups; Gerygone chloronota, Gerygone chrysogaster, Acanthiza apicalis, Acanthiza chrysorrhoa, Dasyornis broadbenti, Pardalotus striatus, Pardalotus punctatus, Sericornis frontalis and Sericornis perspicillatus (Additional file 2: Table S12). Sequences were aligned for each loci using MUSCLE . From here, individual gene trees were generated for each of the nuclear and mitochondrial genes using IQ-TREE . This was implemented in the IQ-TREE web server with an automatic substitution model in which IQ-TREE determines the best-fit substitution model for the data [83, 84].
For the phylogenomic data we used information from two recent phylogenomic studies based on ultraconserved elements (UCE); one resolved genus-level relationships of honeyeaters  and the second resolved some species-level relationships within the genus Meliphaga . Unfortunately, sequence data for ultraconserved elements were not publicly available from the genus-level study . However to include the consideration of honeyeater genus-level relationships resolved from this study we extracted the topology of the ASTRAL tree presented as Fig. 3 by Andersen et al.  to use as a constraint tree. UCEs were downloaded for 12 species of honeyeater (; BioProject PRJNA509981; see Additional file 2: Table S12 for UCE accession numbers and the species included in the constraint tree). Ultraconserved elements were clustered based on locus, aligned using MAFFT , and trees for each locus were generated with IQ-TREE  using a mixed GTR model. A total of 4676 gene trees was generated for ultraconserved elements.
The 13 gene trees and 4676 UCE trees were used as input for ASTRAL (Accurate Species TRee ALgorithm; ), alongside the genus-level topology of honeyeaters from Andersen et al.  which was used as a constraint tree. ASTRAL provides a statistically consistent estimation of the true species tree from unrooted gene trees, under the multi-species coalescent model . The resulting ASTRAL topology tree (Additional file 1: Fig. S1) was then calibrated using a penalized likelihood approach implemented in the program treePL [87, 88] (see Additional file 3 for further details on an additional calibration method).
For the penalised likelihood approach, we first needed to estimate branch lengths of the tree in substitutions per site. To do so we used RAxML-NG . Here, we used the ASTRAL topology tree and the concatenated alignment of the 13 nuclear and mitochondrial loci, to optimize branch lengths and free model parameters on a fixed topology. The concatenated alignment was partitioned by loci and codon position, and we assigned a separate GTR + I model of rate heterogeneity to each locus. The resulting tree was then calibrated using penalized likelihood implemented in the program treePL [87, 88]. TreePL explicitly allows for rate variation across branches but penalises rate differences after cross-validating initial analyses. Unfortunately, fossil coverage of the honeyeaters is poor; only one honeyeater fossil is known . Since the primary hypotheses being examined in this study do not require accurate dating of the tree, this was not done, and the root of the phylogeny was fixed to one.
Finally, to have a fully sampled phylogenetic tree for all honeyeater species the two missing species, Melipotes carolae and Myzomela dammermani were added into the phylogeny using TACT . The placement of these species was based on taxonomic suggestions from Birds of the World available through the Cornell Lab of Ornithology .
Range size and phenotypic trait variables
We obtained range data for all honeyeater species from BirdLife International . Range size based on interpolated range polygons, such as the BirdLife International range maps, tends to over-estimate the true extent of occurrence, which may be a limitation to this study, but which is an unbiased one relative to the hypotheses being tested. Following previous studies which have also used these maps [12, 35], we only considered extant resident and breeding ranges to avoid introducing bias between migratory versus non migratory species and areas where species have been introduced. We used Quantum Geographic Information System (QGIS; ) to extract information relating to range size, shape and position. Geographic range size was extracted (km2). Range shape was quantified as the total range area divided by the total perimeter for each species. Using this index, species with smaller values (smaller area and a larger perimeter) have more elongated or disjunct ranges . Here, we also considered latitudinal and longitudinal extent as important aspects of shape, which are thought to be related to probability of range bisection [13, 18]. Extent was calculated by generating a bounding box around each species range and extracting the length and width of the range extent in kilometres. We characterised range position using two approaches; centroids were extracted from polygons for each species range, which were calculated in QGIS and are influenced by the shape of a species’ range, and we also extracted latitudinal and longitudinal midpoints from the range extent. Both approaches are commonly used and despite being slightly different, these two measures were highly correlated (Additional file 1: Fig. S5) and so midpoints were only considered for analysis.
Here, we also noted whether species were classified as island endemic species or continental species. Species whose ranges overlapped with Australia, Tasmania and New Guinea were considered continental species (n = 124) and all other species as island species (n = 68). This classification takes into consideration the historical biogeography of the region. Both Tasmania and New Guinea sit on continental crust and have been repeatedly connected to Australia via the Sahul Shelf throughout the Pleistocene , which has enabled considerable biotic interchange and gene flow in the Australo-Papuan region.
To examine factors that may contribute to or influence the geographic range size speciation relationship we included measures of body size and dispersal ability. Measures of body size, using mass in grams, were extracted from the literature . We used a similar approach to previous datasets , where for species without direct measurements of body mass we substituted missing values using genus averages, considering the same species list used for the phylogeny (; Additional file 2: Table S11). Dispersal ability was assessed by using the hand-wing index (HWI), an index of wing shape commonly used as an indicator of flight performance and dispersal capability [29, 30]. HWI was calculated as Kipp’s distance (the distance between the tip of the first secondary feather to the tip of the longest primary feather, measured on the folded wing) divided by wing length, then multiplied by 100. Measurements of wing length and Kipp’s distance were extracted from Marki et al. , these measurements were averaged for each species and hand-wing index was calculated. Additional values of HWI were then substituted from Sheard et al. . From here, HWI values were missing for five species of honeyeater (Melipotes carolae, Myzomela prawiradilagae, Myzomela irianawidodoae, Myzomela wahe, and Microptilotis imitatrix) and these were substituted in using phylogenetic imputation implemented in R  using the package ‘phytools’  (Additional file 2: Table S11).
Speciation rate estimates
To obtain a continuous estimate of speciation rate for use in phylogenetic regressions, phylogenetic path analysis and the spatiophylogenetic model, we used the diversification rate (DR) statistic. This was calculated in R  using the DR_statistic function from the ‘fisse’ R package . The DR statistic was selected estimates tip-specific rates of diversification without a formal parametric model, has been shown to provide a better estimate of speciation than net diversification, and performs as well as other commonly used statistics [100, 101]. We additionally used hidden-state dependent speciation and extinction models as a complementary approach, which estimate their own rates of speciation .
To test the relationship between range size and speciation, while accounting for other factors that are likely to influence this relationship, measures of range size, shape, and position for each species were treated as traits and used in phylogenetic generalised least-squares regressions (PGLS) alongside measures of body size and dispersal ability . Continuous variables were log transformed to improve normality and absolute values of latitude and longitude were used. PGLS regression tests were run using the ‘ape’ (v.5.3 ) and ‘nlme’  packages in R . An estimated lambda parameter was used to control for the amount of phylogenetic effect in the model residuals [106, 107]. To examine whether relationships differed between islands and continents, we tested an interaction term in our analysis. We reasoned that if the interaction term is significant, examination of our second hypothesis, that the form of the geographic range size-speciation relationship is influenced by mode of speciation, could be furthered by phylogenetic regression undertaken separately for the island and continental species. Previous studies of other bird clades have employed separate examinations of island and mainland species , and we adopted this rationale (additional split regression analyses, including those with nonsignificant interaction terms, are present in Additional file 2: Table S1).
Phylogenetic bayesian structural equation models
To further uncover the relationship between range size and speciation, we implemented a Bayesian version of phylogenetic path analysis . This analysis was undertaken in R  using ‘brms’ , for Bayesian multilevel models using the probabilistic programming language STAN [110, 111]. Unlike simple phylogenetic regressions, the use of structural equation models enables us to test the multiple interactions between all variables and traits in a single model, an ideal approach to determine the relationship between range size and speciation while accounting for other factors that are likely to influence this relationship (Fig. 1a). Specific model details can be found in Additional file 3. Path analyses were performed on all species, and separately on island and continental species following the rationale above, enabling us to test our second hypothesis that the form of the geographic range size-speciation relationship is influenced by mode of speciation.
Spatiophylogenetic modelling  was then used to estimate the spatial and phylogenetic structure in speciation rate, and to examine the contributions of range size, dispersal ability, body size and range position to variation in speciation rate (DR statistic). To model the spatial distributions of each species, we downloaded available occurrence records for Meliphagidae from eBird via GBIF [113, 114] (see Additional file 2: Table S13 and Additional file 3 for further details), our analysis included a total of 254,886 honeyeater occurrence records.
The spatiophylogenetic model was implemented in R using the package “INLA” which uses an integrated nested Laplace approximation to estimate the joint posterior distribution of model parameters [115, 116]. INLA is an approach for latent Gaussian Markov random field models, providing a significantly faster alternative to Markov chain Monte Carlo [116, 117]. We tested the influence of range size, position (latitudinal and longitudinal midpoints), dispersal ability and body size on speciation rate and included both a phylogenetic random effect and spatial random effect in the model. For the phylogenetic effect, INLA requires a phylogenetic precision matrix which is the inverse of a phylogenetic covariance matrix. Prior to inverting, the phylogenetic covariance matrix was standardised by dividing by its determinant raised to the power of 1/Nspecies . To model the spatial effect across the entire landscape, we used a spatial mesh, which was constructed over the point occurrence records and averaged across the spatial random fields of each occurrence record for each species. This approach integrates a spatial random field across each species’ distribution, so each species contributes a single datapoint to the likelihood, reducing any bias that could arise from biased data in occurrence records such as clustered data due to differential sampling effort, or when the number of occurrence records varies between common and rare species (see  for further model explanation). We used a mesh consisting of 877 vertices and used default priors for both the phylogenetic and spatial effects. The mesh selected provided a suitable coverage of honeyeater distributions and allowed reasonable computational times. Prior to analysis, all predictor variables were standardised by subtracting the mean and dividing by the standard deviation. The significance of effects in the models were based on whether the 95% credible intervals of the effect size overlapped with zero.
State-dependent speciation and extinction analyses
To verify which factors influence speciation rates in honeyeaters, we applied two methods of state-dependent speciation and extinction analyses: FiSSE (Fast, intuitive State-dependent Speciation-Extinction analysis ) and HiSSE (Hidden State-dependent Speciation and Extinction ). We first used the R package “fisse”  to explore trait-dependent speciation. We then used HiSSE  to verify the FiSSE outcome. HiSSE represents an improvement on previous models [118, 119] as it not only models the influence of a binary character on observed speciation rates, but also models an unobserved hidden character that could potentially impact species diversification rates.
The FiSSE analyses explored a range of thresholds to determine whether an association might exist between each trait and diversification rates. Thresholds were classified based on frequency distributions of the data, and the median, 25% and 75% quartiles (Additional file 2: Tables S8, S9). HiSSE models were then investigated where for any given trait, the FiSSE analyses revealed significant, or near significant (p < 0.07) outcomes (Additional file 2: Tables S8, S9). Using the R package “hisse” , twenty different models were fitted. Sixteen of these are HiSSE approaches that model the observed states (0 and 1A), but also assume a hidden state with each of the observed states (0B and 1B). These models estimate speciation and extinction for each of the four states and have transition rates between each of them, with each of the models representing various constraints on the transition rate matrix . The remaining four are null models, corresponding to various forms of character independent diversification (CID). Two CID-2 models are equivalent in complexity to a binary SSE model, and two CID-4 models are equivalent in complexity to the full HiSSE model . As our tree (n = 192) is below the threshold (n = 800) necessary to recover unbiased individual transition rates  we focus our attention on the AIC and relative model weights of HiSSE and CID models.
Availability of data and materials
Additional material is available from Figshare at https://doi.org/10.26180/19375091.
Darwin C. On the origin of species by natural selection. London: John Murray; 1859.
Coyne J, Orr H. Speciation. Sunderland: Sinauer; 2004.
Aguilée R, Gascuel F, Lambert A, Ferriere R. Clade diversification dynamics and the biotic and abiotic controls of speciation and extinction rates. Nat Commun. 2018;9:1–13.
Brennan IG, Keogh JS. Miocene biome turnover drove conservative body size evolution across Australian vertebrates. Proc R Soc B. 2018;285:20181474.
Morales-Barbero J, Gouveia SF, Martinez PA. Historical climatic instability predicts the inverse latitudinal pattern in speciation rate of modern mammalian biota. J Evol Biol. 2021;34:339–51.
Foster KL, Piller KR. Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae). BMC Evol Biol. 2018;18:1–12.
Testo WL, Sessa E, Barrington DS. The rise of the Andes promoted rapid diversification in Neotropical Phlegmariurus (Lycopodiaceae). New Phytol. 2019;222:604–13.
Rabosky DL, Chang J, Title PO, Cowman P, Sallan L, Friedman M, et al. An inverse latitudinal gradient in speciation rate for marine fishes. Nature. 2018;559:392–5.
Igea J, Tanentzap AJ. Angiosperm speciation cools down in the tropics. Ecol Lett. 2020;23:692–700.
Cally JG, Stuart-Fox D, Holman L, Dale J, Medina I. Male-biased sexual selection, but not sexual dichromatism, predicts speciation in birds. Evolution. 2021;75:931–44.
Cardillo M, Huxtable JS, Bromham L. Geographic range size, life history and rates of diversification in Australian mammals. J Evol Biol. 2003;16:282–8.
Conway M, Olsen BJ. Contrasting drivers of diversification rates on islands and continents across three passerine families. Proc R Soc B. 2019;286:20191757.
Rosenzweig ML. Geographical speciation: on range size and the probability of isolate formation. In: Proceedings of the Washington State University conference on biomathematics and biostatistics. Washington State University Press. 1978. p. 172–94.
Gaston KJ, Chown SL. Geographic range size and speciation. Oxford: Oxford University Press; 1999.
Losos JB, Glor RE. Phylogenetic comparative methods and the geography of speciation. Trends Ecol Evol. 2003;18:220–7.
Tobias JA, Ottenburghs J, Pigot AL. Avian diversity: speciation, macroevolution, and ecological function. Annu Rev Ecol Evol Syst. 2020;51:533–60.
Rosenzweig ML. Species diversity in space and time. Cambridge: Cambridge University Press; 1995.
Gaston KJ. The structure and dynamics of geographic ranges. Oxford: Oxford University Press; 2003.
Owens IP, Bennett PM, Harvey PH. Species richness among birds: body size, life history, sexual selection or ecology? Proc R Soc B. 1999;266:933–9.
Gaston KJ, Fuller RA. The sizes of species’ geographic ranges. J Appl Ecol. 2009;46:1–9.
Leão TC, Lughadha EN, Reich PB. Evolutionary patterns in the geographic range size of Atlantic Forest plants. Ecography. 2020;43:1–11.
Jablonski D, Roy K. Geographical range and speciation in fossil and living molluscs. Proc R Soc B. 2003;270:401–6.
Chown SL, Gaston KJ. Areas, cradles and museums: the latitudinal gradient in species richness. Trends Ecol Evol. 2000;15:311–5.
Jablonski D. Extinction and the spatial dynamics of biodiversity. Proc Natl Acad Sci USA. 2008;105:11528–35.
Castiglione S, Mondanaro A, Melchionna M, Serio C, Di Febbraro M, Carotenuto F, et al. Diversification rates and the evolution of species range size frequency distribution. Front Ecol Evol. 2017;5:1–10.
Chown SL. Speciation and rarity: separating cause from consequence. In: The biology of rarity. Dordrecht: Springer; 1997. p. 91–109.
Inostroza-Michael O, Hernandez CE, Rodriguez-Serrano E, Avaria-Llautureo J, Rivadeneira MM. Interspecific geographic range size-body size relationship and the diversification dynamics of neotropical furnariid birds. Evolution. 2018;72:1124–33.
Fine PV. Ecological and evolutionary drivers of geographic variation in species diversity. Annu Rev Ecol Evol Syst. 2018;46:369–92.
Claramunt S, Derryberry EP, Remsen JV, Brumfield RT. High dispersal ability inhibits speciation in a continental radiation of passerine birds. Proc R Soc B. 2012;279:1567–74.
Weeks BC, Claramunt S. Dispersal has inhibited avian diversification in Australasian archipelagoes. Proc R Soc B. 2014;281:20141257.
Ashby B, Shaw AK, Kokko H. An inordinate fondness for species with intermediate dispersal abilities. Oikos. 2020;129:311–9.
Gavrilets S, Li H, Vose MD. Patterns of parapatric speciation. Evolution. 2000;54:1126–34.
Brown JH, Stevens GC, Kaufman DM. The geographic range: size, shape, boundaries, and internal structure. Annu Rev Ecol Syst. 1996;27:597–623.
Lester SE, Ruttenberg BI, Gaines SD, Kinlan BP. The relationship between dispersal ability and geographic range size. Ecol Lett. 2007;10:745–58.
Sheard C, Neate-Clegg MH, Alioravainen N, Jones SE, Vincent C, MacGregor HE, et al. Ecological drivers of global gradients in avian dispersal inferred from wing morphology. Nat Commun. 2020;11:2463.
Birand A, Vose A, Gavrilets S. Patterns of species ranges, speciation, and extinction. Am Nat. 2012;179:1–21.
Agosta SJ, Bernardo J. New macroecological insights into functional constraints on mammalian geographical range size. Proc R Soc B Biol Sci. 2013;280:20130140.
Liow LH, Fortelius M, Bingham E, Lintulaakso K, Mannila H, Flynn L, et al. Higher origination and extinction rates in larger mammals. Proc Natl Acad Sci USA. 2008;105:6097–102.
Monroe MJ, Bokma F. Do speciation rates drive rates of body size evolution in mammals? Am Nat. 2009;174:912–8.
Pigot AL, Phillimore AB, Owens IP, Orme CD. The shape and temporal dynamics of phylogenetic trees arising from geographic speciation. Syst Biol. 2010;59:660–73.
Stevens GC. The latitudinal gradient in geographical range: how so many species coexist in the tropics. Am Nat. 1989;133:240–56.
Meiri S, Dayan T. On the validity of Bergmann’s rule. J Biogeogr. 2003;30:331–51.
Pintor AF, Schwarzkopf L, Krockenberger AK. 2015 Rapoport’s rule: do climatic variability gradients shape range extent? Ecol Monogr. 2003;85:643–59.
Gaston KJ. Species-range size distributions: products of speciation, extinction and transformation. Philos Trans R Soc Lond B Biol Sci. 1998;353:219–30.
Fortes RR, Absalão RS. The applicability of Rapoport’s rule to the marine molluscs of the Americas. J Biogeogr. 2004;31:1909–16.
Ricklefs RE, Bermingham E. Taxon cycles in the lesser Antillean avifauna. Ostrich. 1999;70:49–59.
Marquet PA, Fernández M, Navarette AA, Valdovinos C. Diversity emerging: towards a deconstruction of biodiversity patterns. In: Lomolino M, Heaney L, editors. Frontiers of biogeography: new directions in the geography of nature. Sunderland: Sinauer Associates; 2004. p. 191–209.
Ryan PG, Bloomer P, Moloney CL, Grant TJ, Delport W. Ecological speciation in South Atlantic Island finches. Science. 2007;315:1420–3.
Lamichhaney S, Han F, Webster MT, Andersson L, Grant BR, Grant PR. Rapid hybrid speciation in Darwin’s finches. Science. 2018;359:224–8.
Ford HA. Why does the distribution of the Honeyeaters (Meliphagidae) conform so well to Wallace’s Line? Lisse: Balkema Publishers; 2001.
Higgins PJ, Ford HA, Christidis L. Family meliphagidae. In: del Hoyo J, Elliott A, Christie DA, editors. Handbook of the birds of the world: penduline-tits to shrikes. Barcelona: Lynx Edicions; 2008. p. 498–691.
Joseph L, Toon A, Nyári ÁS, Longmore NW, Rowe KM, Haryoko T, et al. A new synthesis of the molecular systematics and biogeography of honeyeaters (Passeriformes: Meliphagidae) highlights biogeographical and ecological complexity of a spectacular avian radiation. Zool Scr. 2014;43:235–48.
Marki PZ, Jonsson KA, Irestedt M, Nguyen JM, Rahbek C, Fjeldså J. Supermatrix phylogeny and biogeography of the Australasian Meliphagides radiation (Aves: Passeriformes). Mol Phylogenet Evol. 2017;107:516–29.
Price TD. Speciation in birds. Greenwood Village: Roberts and Co; 2008.
Norman JA, Rheindt FE, Rowe DL, Christidis L. Speciation dynamics in the Australo-Papuan Meliphaga honeyeaters. Mol Phylogenet Evol. 2007;42:80–91.
Toon A, Hughes JM, Joseph L. Multilocus analysis of honeyeaters (Aves: Meliphagidae) highlights spatio-temporal heterogeneity in the influence of biogeographic barriers in the Australian monsoonal zone. Mol Ecol. 2010;19:2980–94.
Dolman G, Joseph L. Evolutionary history of birds across southern Australia: structure, history and taxonomic implications of mitochondrial DNA diversity in an ecologically diverse suite of species. Emu. 2015;115:35–48.
Andersen MJ, Naikatini A, Moyle RG. A molecular phylogeny of Pacific honeyeaters (Aves: Meliphagidae) reveals extensive paraphyly and an isolated Polynesian radiation. Mol Phylogenet Evol. 2014;71:308–15.
Andersen MJ, McCullough JM, Friedman NR, Peterson AT, Moyle RG, Joseph L, et al. Ultraconserved elements resolve genus-level relationships in a major Australasian bird radiation (Aves: Meliphagidae). Emu. 2019;119:218–32.
McCullough JM, Joseph L, Moyle RG, Andersen MJ. Ultraconserved elements put the final nail in the coffin of traditional use of the genus Meliphaga (Aves: Meliphagidae). Zool Scr. 2019;119:1–8.
Byrne M, Yeates DK, Joseph L, Kearney M, Bowler J, Williams MA, et al. Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biota. Mol Ecol. 2008;17:4398–417.
Brennan IG, Oliver PM. Mass turnover and recovery dynamics of a diverse Australian continental radiation. Evolution. 2017;71:1352–65.
Pepper M, Keogh JS. Life in the “dead heart” of Australia: the geohistory of the Australian deserts and its impact on genetic diversity of arid zone lizards. J Biogeogr. 2021;48:716–46.
Powney GD, Grenyer R, Orme CDL, Owens IPF, Meiri S. Hot, dry and different: Australian lizard richness is unlike that of mammals, amphibians and birds. Glob Ecol. 2010;19:386–96.
Rabosky DL, Hutchinson MN, Donnellan SC, Talaba AL, Lovette IJ. Phylogenetic disassembly of species boundaries in a widespread group of Australian skinks (Scincidae: Ctenotus). Mol Phylogenet Evol. 2014;77:71–82.
Westerman M, Blacket MJ, Hintz A, Armstrong K, Woolley PA, Krajewski C. A plethora of planigales: genetic variability and cryptic species in a genus of dasyurid marsupials from northern Australia. Aust J Zool. 2016;64:303–11.
Rabosky DL, Donnellan SC, Talaba AL, Lovette IJ. Exceptional among-lineage variation in diversification rates during the radiation of Australia’s most diverse vertebrate clade. Proc R Soc B. 2007;274:2915–23.
Unmack PJ, Bagley JC, Adams M, Hammer MP, Johnson JB. Molecular phylogeny and phylogeography of the Australian freshwater fish genus Galaxiella, with an emphasis on dwarf galaxias (G. pusilla). PLoS ONE. 2012;7:e38433.
Catullo RA, Keogh JS. Aridification drove repeated episodes of diversification between Australian biomes: evidence from a multi-locus phylogeny of Australian toadlets (Uperoleia: Myobatrachidae). Mol Phylogenet Evol. 2014;79:106–17.
Marki PZ, Kennedy JD, Cooney CR, Rahbek C, Fjeldså J. Adaptive radiation and the evolution of nectarivory in a large songbird clade. Evolution. 2019;73:1226–40.
Driskell AC, Christidis L. Phylogeny and evolution of the Australo-Papuan honeyeaters (Passeriformes, Meliphagidae). Mol Phylogenet Evol. 2004;31:943–60.
Mapel XM, Gyllenhaal EF, Modak TH, DeCicco LH, Naikatini A, Utzurrum RB, et al. Inter- and intra-archipelago dynamics of population structure and gene flow in a Polynesian bird. Mol Phylogenet Evol. 2021;156:107034.
Kennedy JD, Borregaard MK, Marki PZ, Machac A, Fjeldså J, Rahbek C. Expansion in geographical and morphological space drives continued lineage diversification in a global passerine radiation. Proc R Soc B. 2018;285:20182181.
Valente L, Phillimore AB, Melo M, Warren BH, Clegg SM, Havenstein K, Tiedelmann R, Illera JC, Thébaud C, Aschenbach T, Etienne RS. A simple dynamic model explains the diversity of island birds worldwide. Nature. 2020;579:92–6.
Vasconcelos T, O’Meara BC, Beaulieu JM. Retiring “cradles” and “museums” of biodiversity. Am Nat. 2022;199:194–205.
Skeels A, Cardillo M. Reconstructing the geography of speciation from contemporary biodiversity data. Am Nat. 2019;193:240–55.
Gill F, Donsker D, Rasmussen P. IOC world bird list (v10.2). 2020. https://doi.org/10.14344/IOC.ML.10.2
Prawiradilaga DM, Baveja P, Suparno, Ashar H, Ng NS, Gwee CY, et al. A colourful new species of Myzomela honeyeater from rote island in eastern Indonesia. Treubia. 2017;44:77–100.
Irham M, Ashari H, Suparno, Trainor CR, Verbelen P, Wu MY, et al. A new Myzomela honeyeater (Meliphagidae) from the highlands of Alor Island, Indonesia. J Ornithol. 2019;161:313–24.
Rheindt FE, Prawiradilaga DM, Ashari H, Suparno, Gwee CY, Lee GW, et al. A lost world in Wallacea: description of a montane archipelagic avifauna. Science. 2020;367:167–70.
Thibault M, Defos du Rau P, Pineau O, Pangimangen W. New and interesting records for the Obi archipelago (north Maluku, Indonesia), including field observations and first description of the vocalisation of Moluccan Woodcock Scolopax rochussenii. Bull Br Ornithol Club. 2013;133:83–112.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004;5:113.
Trifinopoulos J, Nguyen LT, von Haeseler A, Minh BQ. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44:232–5.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37:1530–4.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Mirarab S, Warnow T. ASTRAL-II: coalescent-based species tree estimation with many hundreds of taxa and thousands of genes. Bioinformatics. 2015;31:44–52.
Sanderson MJ. Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002;19:102–9.
Smith SA, O’Meara BC. treePL: divergence time estimation using penalized likelihood for large phylogenies. Bioinformatics. 2012;28:2689–90.
Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019;35:4453–5.
Boles WE. Fossil honeyeaters (Meliphagidae) from the Late Tertiary of Riversleigh, north-western Queensland. Emu. 2005;105:21–6.
Chang J, Rabosky DL, Alfaro ME. Estimating diversification rates on incompletely sampled phylogenies: theoretical concerns and practical solutions. Syst Biol. 2020;69:602–11.
Billerman SM, Keeney BK, Rodewald PG, Schulenberg TS. Birds of the World. Ithaca: Cornell Laboratory of Ornithology; 2020. https://birdsoftheworld.org/bow/home
BirdLife International and handbook of the birds of the world. Bird species distribution maps of the world. Version 2018.1. 2018. http://datazone.birdlife.org/species/requestdis
QGIS Development Team. QGIS geographic information system. Open Source Geospatial Foundation Project. 2016.
O’Connell JF, Allen J. Dating the colonization of Sahul (Pleistocene Australia-New Guinea): a review of recent research. J Archaeol Sci. 2004;31:835–53.
Dunning JB. CRC handbook of avian body masses. Boca Raton: Taylor & Francis Group, LLC; 2008.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2020. https://www.R-project.org/.
Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Rabosky DL, Goldberg EE. FiSSE: a simple nonparametric test for the effects of a binary character on lineage diversification rates. Evolution. 2017;71:1432–42.
Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491:444–8.
Title PO, Rabosky DL. Tip rates, phylogenies and diversification: what are we estimating, and how good are the estimates? Methods Ecol Evol. 2019;10:821–34.
Beaulieu JM, O’Meara BC. Detecting hidden diversification shifts in models of trait-dependent speciation and extinction. Syst Biol. 2016;65:583–601.
Grafen A. The phylogenetic regression. Philos Trans R Soc Lond B Biol Sci. 1989;326:119–57.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. 2013.
Pagel MD. Inferring the historical patterns of biological evolution. Nature. 1999;401:877–84.
Revell LJ. Phylogenetic signal and linear regression on species data. Methods Ecol Evol. 2010;1:319–29.
von Hardenberg A, Gonzalez-Voyer A. Disentangling evolutionary cause-effect relationships with phylogenetic conformatory path analysis. Evolution. 2012;67:378–87.
Bürkner PC. brms: an R package for Bayesian multilevel models using Stan. J Stat Softw. 2017;80:1–28.
Gelman A, Carlin JB, Stern HS, Dunson DB, Vehtari A, Rubin DB. Bayesian data analysis. Boca Raton: CRC Press; 2013.
Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: a probabilistic programming language. J Stat Softw. 2017;76:1–32.
Dinnage R, Skeels A, Cardillo M. Spatiophylogenetic modelling of extinction risk reveals evolutionary distinctiveness and brief flowering period as threats in a hotspot plant genus. Proc R Soc B. 2020;287:20192817.
GBIF.org (13 August 2020) GBIF occurrence download. https://doi.org/10.15468/dl.ee29jh.
eBird Basic Dataset. Version: EBD_relSep-2020. Ithaca: Cornell Lab of Ornithology; 2020.
Martins TG, Simpson D, Lindgren F, Rue H. Bayesian computing with INLA: new features. Comput Stat Data Anal. 2013;67:68–83.
Rue H, Martino S. Approximate Bayesian inference for hierarchical Gaussian Markov random fields models. J Stat Plan Inference. 2006;137:3177–92.
Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc. 2009;71:319–92.
Maddison WP, FitzJohn RG. The unsolved challenge to phylogenetic correlation tests for categorical characters. Syst Biol. 2015;64:127–36.
Rabosky DL, Goldberg EE. Model inadequacy and mistaken inferences of trait-dependent speciation. Syst Biol. 2015;64:340–55.
We thank Carla Sgrò and Craig White for comments on the manuscript. Melodie McGeoch and members of the ChownLab and McGeoch Lab provided comments on aspects of this work. We also thank two anonymous reviewers for providing helpful comments which improved the manuscript.
E.M.H. was supported by an Australian Government Research Training Program (RTP) Scholarship and Monash Graduate Excellence Scholarship (MGES). S.L.C. was funded by Australian Research Council Discovery Project DP170101046. M.D.M. was funded by Australian Research Council Discovery Early Career Research Award DE180101558.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Honeyeater topology tree estimated with ASTRAL, generated from 4676 UCE trees of 12 honeyeater species, 13 mitochondrial and nuclear gene trees (coverage varying from 3 species to 187 species), alongside a constraint tree extracted from Andersen et al.  that resolves genus-level relationships of honeyeaters. Figure S2. Dated phylogeny of honeyeaters, generated in MCMCtree using the topology tree from ASTRAL alongside the fully partitioned and concatenated dataset of all mitochondrial and nuclear genes, with the root of the phylogeny fixed to one. Node labels show the estimated age. Error bars represent 95% highest posterior densities for estimated node age. Figure S3. Dated phylogeny of honeyeaters, generated in TreePL using the topology tree from ASTRAL alongside the fully partitioned and concatenated dataset of all mitochondrial and nuclear genes, with the root of the phylogeny fixed to one. Node labels show the estimated age. Figure S4. Tanglegram comparing the topology of the phylogeny generated in this study (left) with that obtained from previous studies (right; ). Colours indicate similarities and the black lines indicate topological differences, with the lines matching the same species between the two trees. Note the Marki et al.  phylogeny contains 186 honeyeater species and whereas the tree from this studycontains 192 species, but here it has been trimmed to the same 186 species. Figure S5. Correlation plot indicating relationships between all traits and range size variables included in this study. Figure S6. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the TreePL calibration method. Figure S7. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the TreePL calibration method. Figure S8. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the TreePL calibration method. Figure S9. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the MCMCtree calibration method. Figure S10. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the MCMCtree calibration method. Figure S11. Biplots showing the relationship between traits and variables included in this study. Trendlines and p-values are generated from PGLS regressions, using the tree from the MCMCtree calibration method. Figure S12. Results from PGLS regressions showing the relationship between traits and variables included in this study, subset into island species (n = 68; top row—blue panels) and continental species (n = 124; bottom row—red panels). Trendlines and p-values are generated from PGLS regressions, using the tree from the TreePL calibration method. Figure S13. Results from PGLS regressions showing the relationship between traits and variables included in this study, subset into island species (n = 68; top row—blue panels) and continental species (n = 124; bottom row—red panels). Trendlines and p-values are generated from PGLS regressions, using the tree from the TreePL calibration method. Figure S14. Results from PGLS regressions showing the relationship between traits and variables included in this study, subsetinto island species (n = 68; top row—blue panels) and continental species (n = 124; bottom row—red panels). Trendlines and p-values are generated from PGLS regressions, using the tree from the MCMCtree calibration method. Figure S15. Results from PGLS regressions showing the relationship between traits and variables included in this study, subset into island species (n = 68; top row—blue panels) and continental species (n = 124; bottom row—red panels). Trendlines and p-values are generated from PGLS regressions, using the tree from the MCMCtree calibration method. Figure S16. Full phylogenetic Bayesian structural equation model testing the influence of range position (latitudinal and longitudinal midpoints) and islands, on body size and dispersal ability, as well as the interaction between body size and dispersal, and the influence of all these traits on range size, range shape, and ultimately speciation. Model results: Blue arrows indicate a significant negative relationship; Red arrows indicate a significant positive relationship. Grey arrows indicate no significant relationship.
PGLS results (both TreePLtree and MCMCtree tree). Table S2. SEM results TreePL tree. Table S3. SEM results MCMCtree tree. Table S4. SEM results island vs. continental species TreePL tree. Table S5. SEM results island vs continental species MCMCtree tree. Table S6. INLA model results (bothTreePL tree and MCMCtree tree). TableS7. FiSSE results (TreePL tree). TableS8. FiSSE results (MCMCtree tree). Table S9. HiSSE results (TreePL tree). Table S10. HiSSE results (MCMCtree tree). Table S11. Honeyeater trait dataset. Table S12. Accession table, contianing the GenBank accession numbers and list of species and outgroups included in this study. Table S13. Simulated occurrence data for missing honeyeater species.
Supplementary text file.
About this article
Cite this article
Hay, E.M., McGee, M.D. & Chown, S.L. Geographic range size and speciation in honeyeaters. BMC Ecol Evo 22, 86 (2022). https://doi.org/10.1186/s12862-022-02041-6
- Speciation mode
- Geographic range