Parsimony analysis of endemicity and analysis of endemicity: A fair comparison

Parsimony Analysis of Endemicity (PAE) has been compared with other methods regarding its performance to identify areas of endemism. It is frequently compared with the Analysis of Endemicity (AE), which seems to perform better than PAE to identify these areas. Here I compare PAE and AE considering the sympatric taxa diagnosed as endemic, being as strictly close as possible to sympatry, and using previously published data of Sciobius (Coleoptera: Curculionidae). AE identified more candidate areas of endemism than PAE, but the number of highly restricted endemic taxa to these areas was insufficient to support them as areas of endemism. Considering strictly sympatry (homopatry), PAE performed better than AE; however, both methods may identify areas with some grade of sympatry, but the recognition of which areas constitute real areas of endemism in the strict sense depends on the interpretation of the researcher.     

Comparison between Parsimony Analysis of Endemicity (PAE), Endemicity Analysis (EA), and an alternative coding of Three-Distribution Statements based on hypothetical distributions

Areas of endemism are identified by a variety of methods, none of which is universally accepted. Performance of each method depends upon the variables chosen. Here, we compare Parsimony Analysis of Endemicity (PAE), Endemicity Analysis (EA), and a new coding method that we propose, Three-Distribution Statements (3DS). We rate performance based on the ability to identify hypothetical predefined patterns that represent non-conflicting, nested, and overlapping areas of endemism. Additionally, we also compared properties commonly used in analyses, such as shape and size of the area and the number of taxa involved. We found that 3DS has the best performance in retrieving predefined areas. EA is the only method that resolved a completely overlapping pattern, but it also found spurious patterns. Resolution with PAE always had intermediate precision and efficiency and so is not the best option for analysis of endemism. We recommend the use of 3DS together with EA as the best available option for hypothesizing areas of endemism.     

Application of parsimony analysis of endemicity in Amazonian biogeography: an example with primates

The distributions of 51 non-human primate species are used for Parsimony Analysis of Endemicity (PAE) to determine the relationships among 14 interfluvial regions in the Amazon basin, South America. Two most parsimonious cladograms were found. The strict consensus tree of these cladograms suggests an early separation between Lower Amazonia (eastern) and Upper Amazonia (western). The major clusters of interfluvial regions identified in the PAE cladogram are congruent with the areas of endemism delimited for birds. When interfluvial regions are converted into avian areas of endemism, the PAE cladogram is congruent with one of the two general areas cladograms suggested for Amazonia based on phylogenies of several clades of forest birds. Our analysis suggests that PAE can be used as a tool to objectively identify areas of endemism at an intra-continental scale as well as to make historical inferences. However, the value of a PAE cladogram in this latter application should be always evaluated by congruence with area cladograms built upon cladistic biogeography procedures.     

Sympatry inference and network analysis in biogeography

A new approach for biogeography to find patterns of sympatry, based on network analysis, is proposed. Biogeographic analysis focuses basically on sympatry patterns of species. Sympatry is a network (= relational) datum, but it has never been analyzed before using relational tools such as Network Analysis. Our approach to biogeographic analysis consists of two parts: first the sympatry inference and second the network analysis method (NAM). The sympatry inference method was designed to propose sympatry hypothesis, constructing a basal sympatry network based on punctual data, independent of a priori distributional area determination. In this way, two or more species are considered sympatric when there is interpenetration and relative proximity among their records of occurrence. In nature, groups of species presenting within-group sympatry and between-group allopatry constitute natural units (units of co-occurrence). These allopatric units are usually connected by intermediary species. The network analysis method (NAM) that we propose here is based on the identification and removal of intermediary species to segregate units of co-occurrence, using the betweenness measure and the clustering coefficient. The species ranges of the units of co-occurrence obtained are transferred to a map, being considered as candidates to areas of endemism. The new approach was implemented on three different real complex data sets (one of them a classic example previously used in biogeography) resulting in (1) independence of predefined spatial units; (2) definition of co-occurrence patterns from the sympatry network structure, not from species range similarities; (3) higher stability in results despite scale changes; (4) identification of candidates to areas of endemism supported by strictly endemic species; (5) identification of intermediary species with particular biological attributes.     

Distribution patterns and biogeographic analysis of Austral Polychaeta (Annelida)

Aim We investigate the biogeography of Austral Polychaeta (Annelida) using members of the families Eunicidae, Lumbrineridae, Oenonidae, Onuphidae, Serpulidae and Spionidae and Parsimony Analysis of Endemicity (PAE). We determine whether observed polychaete distribution patterns correspond to traditional shallow-water marine areas of endemism, estimate patterns of endemism and relationships between areas of endemism, and infer the biological processes that have caused these patterns. Location The study is concerned with extant polychaete taxa occupying shallow-water areas derived from the breakup of the Gondwana landmass (i.e. Austral areas). Methods Similarity was assessed using a significance test with Jaccard's indices. Areas not significantly different at 0.99 were combined prior to the PAE. Widespread species and genera (155 taxa) were scored for presence/absence for each area of endemism. PAE was used to derive hypotheses of area relationships. Hierarchical patterns in the PAE trees were identified by testing for congruence with patterns derived from cladistic biogeographic studies of other Gondwanan taxa and with geological evidence. Results The polychaete faunas of four area-pairs were not significantly different and the areas amalgamated: South-west Africa and South Africa, New Zealand South Island and Chatham Islands, Macquarie Island and Antipodean Islands, and West Antarctica and South Georgia. Areas with the highest levels of species endemism were southern Australia (67.0%), South-east South America (53.2%) and South Africa (40.4%). About 60% of species and 7.5% of genera occupied a single area of endemism. The remainder were informative in the PAE. Under a no long-distance dispersal assumption a single minimal-length PAE tree resulted (l=367; ci=0.42); under dispersal allowed, three minimal-length trees resulted (l=278; ci=0.56). In relation to the sister grouping of the New Zealand areas and Australia we find congruence between our minimal-length trees and those derived from a biogeographic study of land plants, and with area relationships predicted by the Expanding Earth Model. Main conclusions The polychaete distribution patterns in this study differ slightly from the classical areas of endemism, most notably in being broader, thereby bringing into question the value of using single provincial system for marine biogeographic studies. The Greater New Zealand region is found to be ‘monophyletic’ with respect to polychaetes, that is comprising a genuine biogeographical entity, and most closely related to the polychaete fauna of southern Australia. This finding is consistent with studies of land plants and with the Expanding Earth model, but disagrees with conventional geology and biogeographic hypothesis involving a ‘polyphyletic’ New Zealand. Both vicariance and concerted range expansion (=biotic dispersion) appear to have played important roles in shaping present-day distribution patterns of Austral polychaetes. Shallow-water ridge systems between the Australian and Greater New Zealand continental landmasses during the Tertiary are thought to have facilitated biotic dispersion.     

Biogeography of the Lizard Genus Tropidurus Wied-Neuwied, 1825 (Squamata: Tropiduridae): Distribution,Endemism, and Area Relationships in South America

Based on comprehensive distributional records of the 23 species currently assigned to the lizard genus Tropidurus, we investigated patterns of endemism and area relationships in South America. Two biogeographic methods were applied, Parsimony Analysis of Endemicity (PAE) and Brooks Parsimony Analysis (BPA). Two areas of endemism were detected by PAE: the first within the domains of the semiarid Brazilian Caatinga, which includes seven endemic species, and the second in the region of the Serranía de Huanchaca, eastern Bolivia, in which three endemic species are present. The area cladograms recovered a close relationship between the Atlantic Forest and areas of the South American open corridor. The results revealed a close relationship among the provinces Caatinga (Cerrado, Parana Forest (Pantanal+Chaco)). The uplift of the Brazilian Central Plateau in the Late Pliocene-Early Pleistocene (4-2 Myr BP) has been interpreted as a major event responsible for isolation and differentiation of biotas along these areas. However, we emphasize that without the establishment of a temporal framework concerning the diversification history of Tropidurus it is premature to correlate cladogenetic events with specific time periods or putative vicariant scenarios. The limiting factors hampering the understanding of the biogeographic history of this genus include (1) the absence of temporal references in relation to the diversification of distinct clades within Tropidurus; (2) the lack of an appropriate taxonomic resolution of the species complexes currently represented by widely distributed forms; and (3) the need for a comprehensive phylogenetic hypothesis. We suggest that these three important aspects should be prioritized in future investigations.     

Phytogeographic analysis of taxa endemic to the Yucatán Peninsula using geographic information systems, the domain heuristic method and parsimony analysis of endemicity

Abstract. The distribution of plant taxa endemic to the Yucatán Peninsula was studied using Parsimony Analysis of Endemicity (PAE). The known distribution of 162 endemic plant taxa was plotted and the DOMAIN method together with environmental data were used to model the potential distribution for each taxon. The Peninsula was divided into a grid of quarter‐degree cells for the purpose of identifying distribution patterns. A total of 294 cells were analysed using known collection records and potential distribution of endemic taxa data. Two data matrices were constructed, a matrix of known distribution and a matrix of both the known and potential distribution. The two matrices were included in the PAE to identify areas of endemism. The areas determined with the known distribution were restricted and almost half of them remained unresolved, whereas with the potential distribution, approximately 90% of the cells were assigned to any one of the endemicity areas. Four endemism areas were identified: the Yucatán dry zone, Yucatán, El Petén and Belize. The areas of Yucatán and El Petén could be explained by current and Pleistocene climatic conditions and their congruence with other biological groups. Analysis of the potential distribution identified areas with patterns that share current environmental characteristics and a palaeoclimate history. Potential distribution modelling can eliminate uncertainties in biogeographical analysis caused by lack of data distribution and sample variation and produce information about the relationships between areas and taxa as well as the environmental affinities of taxa.     

Biogeographic analysis of endemic cacti of the Sierra Madre Oriental, Mexico

The distribution of cacti species that inhabit the Sierra Madre Oriental (SMO) was analysed. Grid-cells were analysed using parsimony analysis of endemicity (PAE) and endemism indices. Areas characterized by their diagnostic species were determined, aiming to propose areas for the conservation of threatened cacti. Distributional data were obtained from 1936 herbarium specimens, electronic information, and from field collections. Eight areas of endemism and three main clades were obtained from the grid-cell analysis. Areas obtained from the endemism indices are very similar to those obtained with the PAE, but differ in the association of grid-cells. PAE showed endemism patterns indicating that southern and central sections of the SMO province are the areas richest in geographically-restricted species. The results obtained with different endemism indices detected more or less the same areas, although the importance level is different. The corrected weighted endemism index can be considered as a reliable measure of endemism because it is unrelated to species richness. A regionalization of the SMO in three subprovinces is suggested, supported by characteristic cacti taxa and the existence of natural barriers.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 373–389.     

Areas of endemism of Mexican mammals: reanalysis applying the optimality criterion

In order to test Mexican areas of endemism of mammals identified by previous parsimony analyses of endemicity (PAEs), we applied the optimality criterion to three data matrices (based on point records, potential distributional models and the fill option in software NDM). We modelled the ecological niches of 429 terrestrial mammal species using the genetic algorithm for rule-set prediction (GARP) and models were projected as potential distributional areas. We overlapped the point occurrence data and the individual maps of potential distributions to a grid of 1° latitude–longitude. Three matrices of 247 grid cells (areas) and 429 species were built: (1) a binary matrix with '0' for absence and '1' for presence of at least one record of the species inside the grid-cell; (2) a three-state matrix similar to (1) but assigning the state '2' to the assumed presence in the model of potential distribution; and (3) a three-state matrix similar to (2), but applying the fill option of software NDM instead of using a model. The optimality criterion was performed in NDM version 2.7 and results were examined with VNDM version 2.7. The first and second matrices showed 13 areas of endemism and the third identified 16 areas of endemism. NDM provided a better resolution than PAE, allowing us to identify several new areas of endemism, previously undetected. Ecological niche models, projected as potential distributional areas, and the optimality criterion are very useful to identify areas of endemism, although they should be used with caution because they may overpredict potential distributional areas. PAE seems to underestimate the areas of endemism identified.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 468–478.     

Biogeographical patterns of the avifaunas of the Caribbean Basin Islands: a parsimony perspective

We analyzed the avifaunas of the Caribbean islands and nearby continental areas and their relationships using Parsimony Analysis of Endemicity (PAE), in order to assess biogeographical patterns and their concordance with geological and phylogenetic evidence. Using distributional information of birds obtained from published literature, a presence/absence matrix for 695 genera and 2026 species of land and freshwater birds was constructed and analyzed. Three different analyses were performed: for species, for genera, and for species and genera combined. In the combined analysis, the Lesser Antilles appear paraphyletic at the base of the cladogram. Then, two major clades are identified: South America (Andes, Venezuelan lowlands, Dutch West Indies and Trinidad and Tobago) and North America, including the Greater Antilles in a clade that is the sister area to Yucatan and the Central American countries nested from north to south. PAE results support Caribbean vicariant models and cladistic biogeographical hypotheses on area relationships, and show relative congruence with available phylogenetic data. Bird biogeography on the Caribbean islands appears to have been caused by both vicariance and dispersal processes. © The Willi Hennig Society 2007.     

Distributional patterns of the Neotropical humid montane forest avifaunas

The study of biogeographical patterns of the Neotropical humid montane forest avifauna has been prevented due to the lack of phylogenetic hypotheses for most taxa, and to the paucity of detailed studies about the geographical distribution for most of the species. Distributional patterns of this avifauna were explored by integrating predictive distributional models and Parsimony Analysis of Endemicity (PAE). Distributional maps were generated using the Genetic Algorithms for Rule Set Prediction for 442 species; this information was transformed into a data matrix for analysis with PAE. Hierarchical information corresponding to the genus level was incorporated to help resolve the relationships between areas. A strict consensus cladogram showed a clear separation between the Mesoamerican and South American avifaunas. Within the Mesoamerican clade, three geographically structured groups were obtained, whereas the South American clade showed a polytomy of three groups, with two of them, the South Andean Yungas and the Tepuis, located outside the main clade. The third group is a well resolved clade, including areas from eastern Panama and northern Venezuela to central Bolivia. Area relationships suggest a mixed history of dispersal and vicariant events, with the latter being the most important for explaining the biogeographical patterns found.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 175–194.     

Distributional modelling and parsimony analysis of endemicity of Senecio in the Mediterranean-type climate area of Central Chile

Aim Floristic blocks and areas of endemism resulting from a parsimony analysis of endemicity (PAE) using raw floristic data versus data generated from distributional modelling for 130 species in the genus Senecio Tourn. ex L. distributed in the Mediterranean‐type climate area of Central Chile were compared, and the results were used to identify conservation priorities for the flora of the region. Location Central Chile, between 30° and 38° S. Methods Using herbarium records, a species × area matrix consisting of presence/absence data was constructed from a 0.5° × 0.5° grid. Distributional modelling techniques incorporating vegetation formations, elevation and the contagion index were used to interpolate floristic composition of poorly known areas. Parsimony analysis of endemicity was used to identify floristic blocks and areas of endemism. Results Using the number of most parsimonious trees as an index, distributional modelling greatly optimized the results of the PAE analysis. Three floristic blocks and four areas of endemism were suggested based on the PAE results using potential distribution data not incorporating the contagion index, while four blocks and two areas of endemism were suggested from the PAE results using potential distribution data incorporating the contagion index. Floristic blocks for the northern coast, southern Andes, and northern/central Andes were found, with some blocks showing divisions within them representing distinct geographic subunits. Major breaks between and within floristic blocks were identified at 32.5°–33° S and 34.5°–35° S. Main conclusions The floristic blocks identified with the distributional modelling and PAE correspond well to results from some previous studies and support hypothesized biogeographic divisions within Central Chile. The results were similar to those obtained from parallel analysis of the entire tree flora of Central Chile. The vegetative formation‐based distributional modelling produced robust and reproducible results when used along with PAE, especially when the contagion index was incorporated, and is a useful technique for area classification. The results demonstrate the utility of Senecio as an indicator genus for biogeography and conservation in southern South America.     

一种改进的PAE方法及其在锦鸡儿组(锦鸡儿属)特有性分析中的应用

针对特有性简约性分析（PAE）不足之处,提出一个新的改进分析方法。主要区别是,新方法对分布区内分类群区分了原始和演化,相应地编码为0／1;为了获取分类群原始和演化的特性,分析以分类群分支图为基础。新方法是系统发育与地理分布相统一原理的一个具体的定量化探讨。用新方法分析了豆科锦鸡儿属锦鸡儿组15种,结果表明新方法优于以前的PAE方法。用新的改进方法分析得到的可能的祖先分布区是原始类群树锦鸡儿和一大类属内原始类群的分布区。与目前属的起源问题的一般观点相一致。     

Areas of endemism and patterns of diversity for aphids of the Qinghai-Tibetan Plateau and the Himalayas

Aim  The study aimed to identify areas of endemism for aphids in the Qinghai-Tibetan Plateau and the Himalayas (QTPH), and to test congruence between patterns of endemism and patterns of overall species richness identified in a previous study.
Location  The QTPH.
Methods  A distribution data base of 326 endemic aphids in the QTPH was compiled. The study area was divided into a grid of 2°× 2° operative geographical units. Parsimony analysis of endemicity (PAE) was used to identify areas of endemism, and the diversity patterns of endemic species were then mapped using GIS.
Results  We identified 326 endemic species belonging to 138 genera within Adelgidae and 14 subfamilies of Aphididae. Five areas of endemism were found using PAE analysis: the eastern Himalayas, the western Himalayas, north-western Yunnan, southern Tibet and the eastern QTPH. Maps of patterns of endemism identified four major centres for endemic aphids, namely the western Himalayas, the eastern Himalayas (or Sikkim-Assam Himalayas), north-western Hengduan Mountains and the mountains of southern Gansu Province, and three minor centres, southern Tibet, south-eastern Tibet and the eastern Qinghai Province in the north-eastern QTPH.
Main conclusions  Our study identifies major centres of aphid endemism. Furthermore, there is a noticeable congruence between patterns of endemism and patterns of species richness. The patterns of endemism were most likely influenced by the recent uplift of the QTPH.     

Application of parsimony analysis of endemicity to Mexican gymnosperm distributions: grid-cells, biogeographical provinces and track analysis

Parsimony analysis of endemicity (PAE) was used to analyse the distributional patterns of 124 species of Mexican gymnosperms, using two different sample units: grid-cells and biogeographical provinces. PAE analyses were based on distributional data from herbarium specimens and specialized literature. Two data matrices were constructed for 60 grid-cells of 2° and 14 biogeographical provinces. The analysis of the 2° grid-cell matrix led to 7084 cladograms. The strict consensus cladogram showed several clades equivalent to the results obtained with the biogeographical provinces. Three clades agree with some principal regions of distribution of Mexican pines, previously identified by several authors, located at the northern portion of the Baja California peninsula, the Sierra Madre Occidental, and the Sierra Madre Oriental. These areas represent important centres of species diversity and endemism for Mexican gymnosperms. The analysis of the province matrix led to two most parsimonious cladograms, which only differed in the position of the Sierra Madre Occidental province. The iterative procedure PAE with progressive character elimination was applied to identify generalized tracks, where clades of provinces were considered equivalent to generalized tracks, and each time a cladogram was obtained, species defining its clades were deleted and a new run was undertaken. We found five generalized tracks, mainly located in montane provinces. The distribution patterns of gymnosperms agree with the existence of several Mexican biogeographical provinces, and a different historical biogeography of the Mexican peninsulas from the rest of the country is evident.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 405–417.     

Cladistic analysis of distributions and endemism (CADE): using raw distributions of birds to unravel the biogeography of the South American aridlands

Aim To use published distributional and taxonomic information on Neotropical birds in a cladistic analysis of distributions and endemism (CADE) to generate a testable hypothesis of area‐relationships for the arid areas of endemism, particularly those of Central South America (the ‘arid diagonal’), and to clarify the different methodologies commonly associated with parsimony analysis of endemicity (PAE). Location South America. Methods Cladistic analysis of distributions and endemism. Results We obtain a tree where the relationships of most areas are resolved, and we find support for an exclusive clade of Central South American areas, with the Caatinga as sister to both the Chaco and Cerrado. Main conclusions There is a substantial amount of historical signal in avian distributions, when large numbers of taxa and multiple taxonomic levels are considered. Ecological noise and historical information are more easily distinguished in CADE analyses than they would be in PAE analyses. Based on our results we predict that among aridland birds, the Cerrado and Chaco species will be more closely related to each other than to Caatinga species.     

Parsimony Analysis of Endemicity (PAE) of Mexican hydrological basins based on helminth parasites of freshwater fishes

Aim The distributional patterns of helminthological fauna of freshwater fishes were analysed to postulate a general hypothesis on the relationships of some Mexican hydrological systems. Location Eight hydrological systems of central and eastern Mexico were studied and compared with records from Nicaragua. Methods A Parsimony Analysis of Endemicity (PAE) was applied to the presence/absence of ninety‐two helminth parasite taxa (Monogeneans, Digeneans, Cestodes, Acanthocephalans and Nematodes) of freshwater fishes, from eight Mexican hydrological systems, using the Hennig86 program. Results The results represent the first attempt for a biogeographical analysis through application of the PAE method to the distributional patterns of helminth parasites of freshwater fishes in Mexico. A single most parsimonious cladogram was obtained, which grouped all the Neotropical systems in accordance with previous proposals based on other plant and animal taxa. Main conclusions The most basal systems were Santiago and Lerma basins, which exhibited Nearctic affinities. The remaining areas of the cladogram showed Neotropical affinities. All the southeastern systems were grouped in a clade with the Nicaragua system, providing support for a ‘Mesoamerican province’ based on helminth parasites of cichlids. The cladogram also suggests that the treatment of the Lerma‐Santiago basin as a single biogeographical unit is inaccurate and that they should be treated as separate systems.     

Both Recent and Pleistocene geography determine animal distributional patterns in the Tuscan Archipelago

Cluster analyses, Parsimony Analysis of Endemicity, species–area relationships, Mantel tests, partial Mantel tests and levels of endemism were used to investigate the influence of current geography and paleogeography on the composition of the fauna of the Tuscan Islands, a small archipelago in the Mediterranean. The composition of terrestrial mollusc and tenebrionid beetle faunas was prevalently influenced by Pleistocene conditions. The composition of butterfly (and possibly reptile) island fauna was influenced by both present (Holocene) and Pleistocene factors. Finally, chrysids were influenced only by present factors. On the whole, both present and historical factors concurred to determine present distributional patterns, but their relative influence varied according to the mobility of the group considered. Results provided here give some insights into the treatment of endemics to elucidate historical relationships among areas and add further momentum to a general model of biotic dynamics in landbridge archipelagos.     

Temporal area approach for distributional data in biogeography

A structural approach to temporality in distributional data for use in palaeobiogeography is described herein. Pre-established areas in the distributional data matrix are split temporally, allowing a single geographical space to have multiple iterations [e.g. Area A (Lower Devonian), Area A (Middle Devonian)]. The resulting temporal matrix will allow the representation and capture of any differing relationships through time. Designed primarily for Parsimony Analysis of Endemicity (PAE) and biotic similarity analyses, this approach simply structures distributional data within a temporal partition, meaning that numerical methods can be used to assess relationships between areas to find a branching diagram. Created through the application of the temporal matrix to a given analysis, Temporal Area Approach (TAAp) is a structural approach that facilitates exploration of the data rather than being a hypothesis-driven model following analysis. Understanding the behaviour of non-phylogenetic palaeobiogeographical data and reducing the prevalence of temporal artefacts will lead to more robust area classifications.     

Biogeographic regions of North American mammals based on endemism

Since the 19th Century, two regions have been recognized for North American mammals, which overlap in Mexico. The Nearctic region corresponds to the northern areas and the Neotropical region corresponds to the southern ones. There are no recent regionalizations for these regions under the criterion of endemism. In the present study, we integrate two methods to regionalize North America, using species distribution models of mammals: endemicity analysis (EA) and parsimony analysis of endemicity (PAE). EA was used to obtain areas of endemism and PAE was used to hierarchize them. We found 76 consensus areas from 329 sets classified in 146 cladograms, and the strict consensus cladogram shows a basal polytomy with 14 areas and 16 clades. The final regionalization recognizes two regions (Nearctic and Neotropical) and a transition zone (Mexican Transition Zone), six subregions (Canadian, Alleghanian, Californian‐Rocky Mountain, Pacific Central America, Mexican Gulf‐Central America, and Central America), two dominions (Californian and Rocky Mountain), and 23 provinces. Our analysis show that North America is probably more complex than previously assumed. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 485–499.