An evaluation of new parsimony‐based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony‐based approaches. Here, we compare a parametric method, dispersal–extinction–cladogenesis (DEC), against a parsimony‐based method, dispersal–vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes‐DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study. Location World‐wide. Methods Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes‐DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years. Results Despite differences in the underlying biogeographical model, Bayes‐DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events – which in Bayes‐DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node. Main conclusions By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes‐DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide‐spanning spatial and temporal model proposed here could prove useful for testing large‐scale biogeographical patterns in plants.     

A new african fossil caprin and a combined molecular and morphological bayesian phylogenetic analysis of caprini (Mammalia: Bovidae)

Given that most species that have ever existed on Earth are extinct, no evolutionary history can ever be complete without the inclusion of fossil taxa. Bovids (antelopes and relatives) are one of the most diverse clades of large mammals alive today, with over a hundred living species and hundreds of documented fossil species. With the advent of molecular phylogenetics, major advances have been made in the phylogeny of this clade; however, there has been little attempt to integrate the fossil record into the developing phylogenetic picture. We here describe a new large fossil caprin species from ca. 1.9-Ma deposits from the Middle Awash, Ethiopia. To place the new species phylogenetically, we perform a Bayesian analysis of a combined molecular (cytochrome b) and morphological (osteological) character supermatrix. We include all living species of Caprini, the new fossil species, a fossil takin from the Pliocene of Ethiopia (Budorcas churcheri), and the insular subfossil Myotragus balearicus. The combined analysis demonstrates successful incorporation of both living and fossil species within a single phylogeny based on both molecular and morphological evidence. Analysis of the combined supermatrix produces superior resolution than with either the molecular or morphological data sets considered alone. Parsimony and Bayesian analyses of the data set are also compared and shown to produce similar results. The combined phylogenetic analysis indicates that the new fossil species is nested within Capra, making it one of the earliest representatives of this clade, with implications for molecular clock calibration. Geographical optimization indicates no less than four independent dispersals into Africa by caprins since the Pliocene.     

Fish species introductions provide novel insights into the patterns and drivers of phylogenetic structure in freshwaters

Despite long-standing interest of terrestrial ecologists, freshwater ecosystems are a fertile, yet unappreciated, testing ground for applying community phylogenetics to uncover mechanisms of species assembly. We quantify phylogenetic clustering and overdispersion of native and non-native fishes of a large river basin in the American Southwest to test for the mechanisms (environmental filtering versus competitive exclusion) and spatial scales influencing community structure. Contrary to expectations, non-native species were phylogenetically clustered and related to natural environmental conditions, whereas native species were not phylogenetically structured, likely reflecting human-related changes to the basin. The species that are most invasive (in terms of ecological impacts) tended to be the most phylogenetically divergent from natives across watersheds, but not within watersheds, supporting the hypothesis that Darwin''s naturalization conundrum is driven by the spatial scale. Phylogenetic distinctiveness may facilitate non-native establishment at regional scales, but environmental filtering restricts local membership to closely related species with physiological tolerances for current environments. By contrast, native species may have been phylogenetically clustered in historical times, but species loss from contemporary populations by anthropogenic activities has likely shaped the phylogenetic signal. Our study implies that fundamental mechanisms of community assembly have changed, with fundamental consequences for the biogeography of both native and non-native species.     

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae)

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.     

Probabilistic methods outperform parsimony in the phylogenetic analysis of data simulated without a probabilistic model

To understand patterns and processes of the diversification of life, we require an accurate understanding of taxon interrelationships. Recent studies have suggested that analyses of morphological character data using the Bayesian and maximum likelihood Mk model provide phylogenies of higher accuracy compared to parsimony methods. This has proved controversial, particularly studies simulating morphology‐data under Markov models that assume shared branch lengths for characters, as it is claimed this leads to bias favouring the Bayesian or maximum likelihood Mk model over parsimony models which do not explicitly make this assumption. We avoid these potential issues by employing a simulation protocol in which character states are randomly assigned to tips, but datasets are constrained to an empirically realistic distribution of homoplasy as measured by the consistency index. Datasets were analysed with equal weights and implied weights parsimony, and the maximum likelihood and Bayesian Mk model. We find that consistent (low homoplasy) datasets render method choice largely irrelevant, as all methods perform well with high consistency (low homoplasy) datasets, but the largest discrepancies in accuracy occur with low consistency datasets (high homoplasy). In such cases, the Bayesian Mk model is significantly more accurate than alternative models and implied weights parsimony never significantly outperforms the Bayesian Mk model. When poorly supported branches are collapsed, the Bayesian Mk model recovers trees with higher resolution compared to other methods. As it is not possible to assess homoplasy independently of a tree estimate, the Bayesian Mk model emerges as the most reliable approach for categorical morphological analyses.     

Colonization of the Southern Hemisphere by fur seals and sea lions (Carnivora: Otariidae) revealed by combined evidence phylogenetic and Bayesian biogeographical analysis

Fur seals and sea lions (Carnivora: Otariidae) evolved in the North Pacific and later dispersed throughout the Southern Hemisphere. However, the timing and number of dispersals into the Southern Hemisphere still remain poorly understood. To determine the biogeographical patterns of dispersal within fur seals and sea lions, we conducted cladistic analyses using combined evidence incorporating morphological and molecular data. The phylogeny produced in this study was then incorporated into Bayesian biogeographical analyses to reconstruct ancestral points of origin and dispersal patterns for otariid clades. Combined evidence analyses supported Callorhinus as the earliest diverging extant otariid, and a strongly supported northern sea lion clade (Zalophus, Eumetopias, and Proterozetes) as the sister group to a southern clade comprising the remainder of Otariidae. Fossil data constrained the timing and location of this dispersal as occurring between 6 and 7 Mya during a period of unusually cool sea surface temperatures and high productivity in the eastern equatorial Pacific, far older than suggested by prior studies. Our study indicates that the distribution of fur seals and sea lions is tightly linked to sea surface temperature and productivity, and suggests that otariids may be vulnerable to future anthropogenic climate change. © 2014 The Linnean Society of London     

The phylogenetic systematics of blue‐tailed skinks (Plestiodon) and the family Scincidae

Blue‐tailed skinks (genus Plestiodon) are a common component of the terrestrial herpetofauna throughout their range in eastern Eurasia and North and Middle America. Plestiodon species are also frequent subjects of ecological and evolutionary research, yet a comprehensive, well‐supported phylogenetic framework does not yet exist for this genus. We construct a comprehensive molecular phylogeny of Plestiodon using Bayesian phylogenetic analyses of a nine‐locus data set comprising 8308 base pairs of DNA, sampled from 38 of the 43 species in the genus. We evaluate potential gene tree/species tree discordance by conducting phylogenetic analyses of the concatenated and individual locus data sets, as well as employing coalescent‐based methods. Specifically, we address the placement of Plestiodon within the evolutionary tree of Scincidae, as well as the phylogenetic relationships between Plestiodon species, and their taxonomy. Given our sampling of major Scincidae lineages, we also re‐evaluate ‘deep’ relationships within the family, with the goal of resolving relationships that have been ambiguous in recent molecular phylogenetic analyses. We infer strong support for several scincid relationships, including a major clade of ‘scincines’ and the inter‐relationships of major Mediterranean and southern African genera. Although we could not estimate the precise phylogenetic affinities of Plestiodon with statistically significant support, we nonetheless infer significant support for its inclusion in a large ‘scincine’ clade exclusive of Acontinae, Lygosominae, Brachymeles, and Ophiomorus. Plestiodon comprises three major geographically cohesive clades. One of these clades is composed of mostly large‐bodied species inhabiting northern Indochina, south‐eastern China (including Taiwan), and the southern Ryukyu Islands of Japan. The second clade comprises species inhabiting central China (including Taiwan) and the entire Japanese archipelago. The third clade exclusively inhabits North and Middle America and the island of Bermuda. A vast majority of interspecific relationships are strongly supported in the concatenated data analysis, but there is nonetheless significant conflict amongst the individual gene trees. Coalescent‐based gene tree/species tree analyses indicate that incongruence amongst the nuclear loci may severely obscure the phylogenetic inter‐relationships of the primarily small‐bodied Plestiodon species that inhabit the central Mexican highlands. These same analyses do support the sister relationship between Plestiodon marginatus Hallowell, 1861 and Plestiodon stimpsonii (Thompson, 1912), and differ with the mitochondrial DNA analysis that supports Plestiodon elegans (Boulenger, 1887) + P. stimpsonii. Finally, because the existing Plestiodon taxonomy is a poor representation of evolutionary relationships, we replace the existing supraspecific taxonomy with one congruent with our phylogenetic results. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 165 , 163–189.     

Morphological phylogenetic analysis of Ophrys (Orchidaceae): insights from morpho‐anatomical floral traits into the interspecific relationships in an unresolved clade

Reconstructing the phylogeny of the sexually deceptive orchid genus Ophrys is crucial to our understanding of the evolution of its complex floral morphology. Molecular phylogenetic analyses showed that section Pseudophrys forms a well supported clade with Ophrys bombyliflora, O. tenthredinifera and O. speculum, but were unable to elucidate the relationships between these four groups of taxa. Here we conduct a morphological phylogenetic analysis of this unresolved clade of Ophrys based on a data matrix of 45 macro‐ and micromorphological and anatomical floral characters, using maximum parsimony and Bayesian inference. Our cladistic analysis yielded a single most parsimonious tree and a Bayesian 50% majority‐rule consensus tree which differed in their overall topology but agreed that O. tenthredinifera and O. bombyliflora are not sister groups. The phylogenetic placement of O. tenthredinifera was ambiguous since it shares six valid synapomorphies each with the cluster of O. speculum–O. bombyliflora and with section Pseudophrys. In contrast, O. bombyliflora is most likely the sister group to O. speculum, a finding that rejects an earlier morphological phylogenetic hypothesis and favours the existing molecular trees based on nuclear ITS rather than plastid data. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 454–476.     

Biogeography of the grasses (Poaceae): a phylogenetic approach to reveal evolutionary history in geographical space and geological time

The grasses (Poaceae) are the fifth most diverse family of angiosperms, including 800 genera and more than 10 000 species. Few phylogenetic studies have tried to investigate palaeo‐biogeographical and palaeo‐ecological scenarios that may have led to present‐day distribution and diversity of grasses at the family level. We produced a dated phylogenetic tree based on combined plastid DNA sequences and a comprehensive sample of Poaceae. Furthermore, we produced an additional tree using a supermatrix of morphological and molecular data that included all 800 grass genera so that ancestral biogeography and ecological habitats could be inferred. We used a likelihood‐based method, which allows the estimation of ancestral polymorphism in both biogeographical and ecological analyses for large data sets. The origin of Poaceae was retrieved as African and shade adapted. The crown node of the BEP + PACCMAD clade was dated at 57 Mya, in the early Eocene. Grasses dispersed to all continents by approximately 60 million years after their Gondwanan origin in the late Cretaceous. PACCMAD taxa adapted to open habitats as early as the late Eocene, a date consistent with recent phytolith fossil data for North America. C₄ photosynthesis first originated in Africa, at least for Chloridoideae in the Eocene at c. 30 Mya. The BEP clade members adapted to open habitats later than PACCMAD members; this was inferred to occur in Eurasia in the Oligocene. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 162 , 543–557.     

Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.     

Morphology‐based phylogenetic analysis and classification of the family Rhinocryptidae (Aves: Passeriformes)

The family Rhinocryptidae comprises an assemblage of 12 genera and 55 species confined to the Neotropical region. Here we present the first morphology‐based phylogenetic study of the Rhinocryptidae, using 90 anatomical characters (62 osteological, 28 syringeal) scored for all genera of the family and representatives of all families of the infraorder Furnariides. Parsimony analysis of this dataset recovered 7428 equally most‐parsimonious trees. The strict consensus of those trees was completely resolved at the genus level, with the topology (Liosceles (Psilorhamphus ((Eleoscytalopus + Merulaxis) (Acropternis ((Teledromas + Rhinocrypta) ((Pteroptochos + Scelorchilus) (Eugralla (Myornis + Scytalopus)))))))). The monophyly of the Rhinocryptidae as presently understood was recovered with strong support [eight synapomorphies and Bremer support (BS) = 6). Strongly supported internal arrangements included the basal position of the Amazonian genus Liosceles relative to the rest of the family (four synapomorphies, BS = 4), a clade containing Acropternis through Scytalopus (six synapomorphies, BS = 4), and other less inclusive nodes. The main points of congruence between the present morphological phylogeny and previous molecular phylogenetic work on the family were clades supported by six or more synapomorphies and Bremer values of 6–7: Eleoscytalopus + Merulaxis (eight synapomorphies, BS = 6), Scelorchilus + Pteroptochos (seven synapomorphies, BS = 7), Rhinocrypta + Teledromas (seven synapomorphies, BS = 7), and Eugralla + Myornis + Scytalopus (six synapomorphies, BS = 6). A classification derived from the morphological phylogeny is proposed, with new suprageneric taxa being named and diagnosed. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 377–432.     

Insights into plant biodiversity conservation in large river valleys in China: A spatial analysis of species and phylogenetic diversity

Large river valleys (LRVs) are heterogeneous in habitat and rich in biodiversity, but they are largely overlooked in policies that prioritize conservation. Here, we aimed to identify plant diversity hotspots along LRVs based on species richness and spatial phylogenetics, evaluate current conservation effectiveness, determine gaps in the conservation networks, and offer suggestions for prioritizing conservation. We divided the study region into 50 km × 50 km grid cells and determined the distribution patterns of seed plants by studying 124,927 occurrence points belonging to 14,481 species, using different algorithms. We generated phylogenies for the plants using the “V. PhyloMaker” R package, determined spatial phylogenetics, and conducted correlation analyses between different distribution patterns and spatial phylogenetics. We evaluated the effectiveness of current conservation practices and discovered gaps of hotspots within the conservation networks. In the process, we identified 36 grid cells as hotspots (covering 10% of the total area) that contained 83.4% of the species. Fifty‐eight percent of the hotspot area falls under the protection of national nature reserves (NNRs) and 83% falls under national and provincial nature reserves (NRs), with 42% of the area identified as conservation gaps of NNRs and 17% of the area as gaps of NRs. The hotspots contained high proportions of endemic and threatened species, as did conservation gaps. Therefore, it is necessary to optimize the layout of current conservation networks, establish micro‐nature reserves, conduct targeted conservation priority planning focused on specific plant groups, and promote conservation awareness. Our results show that the conservation of three hotspots in Southwest China, in particular, is likely to positively affect the protection of biodiversity in the LRVs, especially with the participation of the neighboring countries, India, Myanmar, and Laos.     

Seasonally different response of photosynthetic activity to daytime and night‐time warming in the Northern Hemisphere

Over the last century the Northern Hemisphere has experienced rapid climate warming, but this warming has not been evenly distributed seasonally, as well as diurnally. The implications of such seasonal and diurnal heterogeneous warming on regional and global vegetation photosynthetic activity, however, are still poorly understood. Here, we investigated for different seasons how photosynthetic activity of vegetation correlates with changes in seasonal daytime and night‐time temperature across the Northern Hemisphere (>30°N), using Normalized Difference Vegetation Index (NDVI) data from 1982 to 2011 obtained from the Advanced Very High Resolution Radiometer (AVHRR). Our analysis revealed some striking seasonal differences in the response of NDVI to changes in day‐ vs. night‐time temperatures. For instance, while higher daytime temperature (T_max) is generally associated with higher NDVI values across the boreal zone, the area exhibiting a statistically significant positive correlation between T_max and NDVI is much larger in spring (41% of area in boreal zone – total area 12.6 × 10⁶ km²) than in summer and autumn (14% and 9%, respectively). In contrast to the predominantly positive response of boreal ecosystems to changes in T_max, increases in T_max tended to negatively influence vegetation growth in temperate dry regions, particularly during summer. Changes in night‐time temperature (T_min) correlated negatively with autumnal NDVI in most of the Northern Hemisphere, but had a positive effect on spring and summer NDVI in most temperate regions (e.g., Central North America and Central Asia). Such divergent covariance between the photosynthetic activity of Northern Hemispheric vegetation and day‐ and night‐time temperature changes among different seasons and climate zones suggests a changing dominance of ecophysiological processes across time and space. Understanding the seasonally different responses of vegetation photosynthetic activity to diurnal temperature changes, which have not been captured by current land surface models, is important for improving the performance of next generation regional and global coupled vegetation‐climate models.