Defining Biologically Meaningful Biomes Through Floristic, Functional, and Phylogenetic Data

While we have largely improved our understanding on what biomes are and their utility in global change ecology, conservation planning, and evolutionary biology is clear, there is no consensus on how biomes should be delimited or mapped. Existing methods emphasize different aspects of biomes, with different strengths and limitations. We introduce a novel approach to biome delimitation and mapping, based upon combining individual regionalizations derived from floristic, functional, and phylogenetic data linked to environmentally trained species distribution models. We define “core Biomes” as areas where independent regionalizations agree and “transition zones” as those whose biome identity is not corroborated by all analyses. We apply this approach to delimiting the neglected Caatinga seasonally dry tropical forest biome in northeast Brazil. We delimit the “core Caatinga” as a smaller and more climatically limited area than previous definitions, and argue it represents a floristically, functionally, and phylogenetically coherent unit within the driest parts of northeast Brazil. “Caatinga transition zones” represent a large and biologically important area, highlighting that ecological and evolutionary processes work across environmental gradients and that biomes are not categorical variables. We discuss the differences among individual regionalizations in an ecological and evolutionary context and the potential limitations and utility of individual and combined biome delimitations. Our integrated ecological and evolutionary definition of the Caatinga and associated transition zones are argued to best describe and map biologically meaningful biomes.

Calibrating phylogenies assuming bifurcation or budding alters inferred macroevolutionary dynamics in a densely sampled phylogeny of bivalve families

Analyses of evolutionary dynamics depend on how phylogenetic data are time-scaled. Most analyses of extant taxa assume a purely bifurcating model, where nodes are calibrated using the daughter lineage with the older first occurrence in the fossil record. This contrasts with budding, where nodes are calibrated using the younger first occurrence. Here, we use the extensive fossil record of bivalve molluscs for a large-scale evaluation of how branching models affect macroevolutionary analyses. We time-calibrated 91% of nodes, ranging in age from 2.59 to 485 Ma, in a phylogeny of 97 extant bivalve families. Allowing budding-based calibrations minimizes conflict between the tree and observed fossil record, and reduces the summed duration of inferred ‘ghost lineages’ from 6.76 billion years (Gyr; bifurcating model) to 1.00 Gyr (budding). Adding 31 extinct paraphyletic families raises ghost lineage totals to 7.86 Gyr (bifurcating) and 1.92 Gyr (budding), but incorporates more information to date divergences between lineages. Macroevolutionary analyses under a bifurcating model conflict with other palaeontological evidence on the magnitude of the end-Palaeozoic extinction, and strongly reduce Cenozoic diversification. Consideration of different branching models is essential when node-calibrating phylogenies, and for a major clade with a robust fossil record, a budding model appears more appropriate.

treedata.table: a wrapper for data.table that enables fast manipulation of large phylogenetic trees matched to data

The number of terminals in phylogenetic trees has significantly increased over the last decade. This trend reflects recent advances in next-generation sequencing, accessibility of public data repositories, and the increased use of phylogenies in many fields. Despite R being central to the analysis of phylogenetic data, manipulation of phylogenetic comparative datasets remains slow, complex, and poorly reproducible. Here, we describe the first R package extending the functionality and syntax of data.table to explicitly deal with phylogenetic comparative datasets. treedata.table significantly increases speed and reproducibility during the data manipulation steps involved in the phylogenetic comparative workflow in R. The latest release of treedata.table is currently available through CRAN (https://cran.r-project.org/web/packages/treedata.table/). Additional documentation can be accessed through rOpenSci (https://ropensci.github.io/treedata.table/).

Three new Russula species in sect. Ingratae (Russulales, Basidiomycota) from southern China

Three new species of Russulasection Ingratae, found in Guizhou and Jiangsu Provinces, southern China, are proposed: R. straminella, R. subpectinatoides and R. succinea. Photographs, line drawings and detailed morphological descriptions for these species are provided with comparisons against closely-related taxa. Phylogenetic analysis of the internal transcribed spacer (ITS) region supported the recognition of these specimens as new species. Additionally, R. indocatillus is reported for the first time from China and morphological and phylogenetic data are provided for the Chinese specimens.

A Natural Colonisation of Asia: Phylogenomic and Biogeographic History of Coin Spiders (Araneae: Nephilidae: Herennia)

Reconstructing biogeographic history is challenging when dispersal biology of studied species is poorly understood, and they have undergone a complex geological past. Here, we reconstruct the origin and subsequent dispersal of coin spiders (Nephilidae: Herennia Thorell), a clade of 14 species inhabiting tropical Asia and Australasia. Specifically, we test whether the all-Asian range of Herennia multipuncta is natural vs. anthropogenic. We combine Anchored Hybrid Enrichment phylogenomic and classical marker phylogenetic data to infer species and population phylogenies. Our biogeographical analyses follow two alternative dispersal models: ballooning vs. walking. Following these assumptions and considering measured distances between geographical areas through temporal intervals, these models infer ancestral areas based on varying dispersal probabilities through geological time. We recover a wide ancestral range of Herennia including Australia, mainland SE Asia and the Philippines. Both models agree that H. multipuncta internal splits are generally too old to be influenced by humans, thereby implying its natural colonisation of Asia, but suggest quite different colonisation routes of H. multipuncta populations. The results of the ballooning model are more parsimonious as they invoke fewer chance dispersals over large distances. We speculate that coin spiders’ ancestor may have lost the ability to balloon, but that H. multipuncta regained it, thereby colonising and maintaining larger areas.

Ancestral protein topologies draw the rooted bacterial tree of life

Many experimental and predicted observations remain unanswered in the current proposed trees of life (ToL). Also, the current trend in reporting phylogenetic data is based in mixing together the information of dozens of genomes or entire conserved proteins. In this work, we consider the modularity of protein evolution and, using only two domains with duplicated ancestral topologies from a single, universal primordial protein corresponding to the RNA binding regions of contemporary bacterial glycyl tRNA synthetase (bacGlyRS), archaeal CCA adding enzyme (arch-CCAadd) and eukaryotic rRNA processing enzyme (euk-rRNA), we propose a rooted bacterial ToL that agrees with several previous observations unaccounted by the available trees.

A new genus for the tiny hawk Accipiter superciliosus and semicollared hawk A. collaris (Aves: Accipitridae), with comments on the generic name for the crested goshawk A. trivirgatus and Sulawesi goshawk A. griseiceps

Multiple molecular phylogenetic studies have demonstrated that two Neotropical raptors, tiny hawk Accipiter superciliosus and its sister species semicollared hawk A. collaris, are not closely related to core Accipiter, and that A. superciliosus, at least, possesses osteological characters not replicated in the remainder of the genus. Based on these data, there is a need to recognise their distinctiveness at generic level. However, as recently noted in two global bird checklists, no name is available to accommodate them, so we provide a new nomen here. Furthermore, two Asian accipitrids, crested goshawk A. trivirgatus and its presumed closest relative Sulawesi goshawk A. griseiceps, are also phylogenetically distinctive; in this case the genus-group name Lophospiza is applicable. We also designate type species for two genus-group names (Hieraspiza and Eusparvius) currently in the synonymy of Accipiter, and, as an aid to future workers, we provide a synonymy of the genus Accipiter and a list of species currently included in Accipiter for which published molecular phylogenetic data are apparently lacking.

Hypothesis Testing With Rank Conditions in Phylogenetics

A phylogenetic model of sequence evolution for a set of n taxa is a collection of probability distributions on the 4n possible site patterns that may be observed in their aligned DNA sequences. For a four-taxon model, one can arrange the entries of these probability distributions into three flattening matrices that correspond to the three different unrooted leaf-labeled four-leaf trees, or quartet trees. The flattening matrix corresponding to the tree parameter of the model is known to satisfy certain rank conditions. Methods such as ErikSVD and SVDQuartets take advantage of this observation by applying singular value decomposition to flattening matrices consisting of empirical data. Each possible quartet is assigned an “SVD score” based on how close the flattening is to the set of matrices of the predicted rank. When choosing among possible quartets, the one with the lowest score is inferred to be the phylogeny of the four taxa under consideration. Since an n-leaf phylogenetic tree is determined by its quartets, this approach can be generalized to infer larger phylogenies. In this article, we explore using the SVD score as a test statistic to test whether phylogenetic data were generated by a particular quartet tree. To do so, we use several results to approximate the distribution of the SVD score and to give upper bounds on the p-value of the associated hypothesis tests. We also apply these hypothesis tests to simulated phylogenetic data and discuss the implications for interpreting SVD scores in rank-based inference methods.