The squirrel is in the detail: Anatomy and morphometry of the tail in Sciuromorpha (Rodentia, Mammalia).

In mammals, the caudal vertebrae are certainly among the least studied elements of their skeleton. However, the tail plays an important role in locomotion (e.g., balance, prehensility) and behavior (e.g., signaling). Previous studies largely focused on prehensile tails in Primates and Carnivora, in which certain osteological features were selected and used to define tail regions (proximal, transitional, distal). Interestingly, the distribution pattern of these anatomical characters and the relative proportions of the tail regions were similar in both orders. In order to test if such tail regionalization can be applied to Rodentia, we investigated the caudal vertebrae of 20 Sciuridae and six Gliridae species. Furthermore, we examined relationships between tail anatomy/morphometry and locomotion. The position of selected characters along the tail was recorded and their distribution was compared statistically using Spearman rank correlation. Vertebral body length (VBL) was measured to calculate the proportions of each tail region and to perform procrustes analysis on the shape of relative vertebral body length (rVBL) progressions. Our results show that tail regionalization, as defined for Primates and Carnivora, can be applied to almost all investigated squirrels, regardless of their locomotor category. Moreover, major locomotor categories can be distinguished by rVBL progression and tail region proportions. In particular, the small flying squirrels Glaucomys volans and Hylopetes sagitta show an extremely short transitional region. Likewise, several semifossorial taxa can be distinguished by their short distal region. Moreover, among flying squirrels, Petaurista petaurista shows differences with the small flying squirrels, mirroring previous observations on locomotory adaptations based on their inner ear morphometry. Our results show furthermore that the tail region proportions of P. petaurista, phylogenetically more basal than the small flying squirrels, are similar to those of bauplan-conservative arboreal squirrels.

Reevaluating claims of ecological speciation in Halichoeres bivittatus.

Allopatry has traditionally been viewed as the primary driver of speciation in marine taxa, but the geography of the marine environment and the larval dispersal capabilities of many marine organisms render this view somewhat questionable. In marine fishes, one of the earliest and most highly cited empirical examples of ecological speciation with gene flow is the slippery dick wrasse, Halichoeres bivittatus. Evidence for this cryptic or incipient speciation event was primarily in the form of a deep divergence in a single mitochondrial locus between the northern and southern Gulf of Mexico, combined with a finding that these two haplotypes were associated with different habitat types ("tropical" vs. "subtropical") in the Florida Keys and Bermuda, where they overlap. Here, we examine habitat assortment in the Florida Keys using a broader sampling of populations and habitat types than were available for the original study. We find no evidence to support the claim that haplotype frequencies differ between habitat types, and little evidence to support any differences between populations in the Keys. These results undermine claims of ecological speciation with gene flow in Halichoeres bivittatus. Future claims of this type should be supported by multiple lines of evidence that illuminate potential mechanisms and allow researchers to rule out alternative explanations for spatial patterns of genetic differences.

Climate change may induce connectivity loss and mountaintop extinction in Central American forests.

The tropical forests of Central America serve a pivotal role as biodiversity hotspots and provide ecosystem services securing human livelihood. However, climate change is expected to affect the species composition of forest ecosystems, lead to forest type transitions and trigger irrecoverable losses of habitat and biodiversity. Here, we investigate potential impacts of climate change on the environmental suitability of main plant functional types (PFTs) across Central America. Using a large database of occurrence records and physiological data, we classify tree species into trait-based groups and project their suitability under three representative concentration pathways (RCPs 2.6, 4.5 and 8.5) with an ensemble of state-of-the-art correlative modelling methods. Our results forecast transitions from wet towards generalist or dry forest PFTs for large parts of the study region. Moreover, suitable area for wet-adapted PFTs is projected to latitudinally diverge and lose connectivity, while expected upslope shifts of montane species point to high risks of mountaintop extinction. These findings underline the urgent need to safeguard the connectivity of habitats through biological corridors and extend protected areas in the identified transition hotspots.

Could climate change benefit invasive snakes? Modelling the potential distribution of the California Kingsnake in the Canary Islands.

The interaction between climate change and biological invasions is a global conservation challenge with major consequences for invasive species management. However, our understanding of this interaction has substantial knowledge gaps; this is particularly relevant for invasive snakes on islands because they can be a serious threat to island ecosystems. Here we evaluated the potential influence of climate change on the distribution of invasive snakes on islands, using the invasion of the California kingsnake (Lampropeltis californiae) in Gran Canaria. We analysed the potential distribution of L. californiae under current and future climatic conditions in the Canary Islands, with the underlying hypothesis that the archipelago might be suitable for the species under these climate scenarios. Our results indicate that the Canary Islands are currently highly suitable for the invasive snake, with increased suitability under the climate change scenarios tested here. This study supports the idea that invasive reptiles represent a substantial threat to near-tropical regions, and builds on previous studies suggesting that the menace of invasive reptiles may persist or even be exacerbated by climate change. We suggest future research should continue to fill the knowledge gap regarding invasive reptiles, in particular snakes, to clarify their potential future impacts on global biodiversity.

Comparing the dietary niche overlap and ecomorphological differences between invasive Hemidactylus mabouia geckos and a native gecko competitor.

Hemidactylus mabouia is one of the most successful, widespread invasive reptile species and has become ubiquitous across tropical urban settings in the Western Hemisphere. Its ability to thrive in close proximity to humans has been linked to the rapid disappearance of native geckos. However, aspects of Hemidactylus mabouia natural history and ecomorphology, often assumed to be linked with this effect on native populations, remain understudied or untested. Here, we combine data from ∂15N and ∂13C stable isotopes, stomach contents, and morphometric analyses of traits associated with feeding and locomotion to test alternate hypotheses of displacement between H. mabouia and a native gecko, Phyllodactylus martini, on the island of Curaçao. We demonstrate substantial overlap of invertebrate prey resources between the species, with H. mabouia stomachs containing larger arthropod prey as well as vertebrate prey. We additionally show that H. mabouia possesses several morphological advantages, including larger sizes in feeding-associated traits and limb proportions that could offer a propulsive locomotor advantage on vertical surfaces. Together, these findings provide the first support for the hypotheses that invasive H. mabouia and native P. martini overlap in prey resources and that H. mabouia possess ecomorphological advantages over P. martini. This work provides critical context for follow-up studies of H. mabouia and P. martini natural history and direct behavioral experiments that may ultimately illuminate the mechanisms underlying displacement on this island and act as a potential model for other systems with Hemidactylus mabouia invasions.

Eyes Wide Shut: the impact of dim-light vision on neural investment in marine teleosts.

Understanding how organismal design evolves in response to environmental challenges is a central goal of evolutionary biology. In particular, assessing the extent to which environmental requirements drive general design features among distantly related groups is a major research question. The visual system is a critical sensory apparatus that evolves in response to changing light regimes. In vertebrates, the optic tectum is the primary visual processing centre of the brain and yet it is unclear how or whether this structure evolves while lineages adapt to changes in photic environment. On one hand, dim-light adaptation is associated with larger eyes and enhanced light-gathering power that could require larger information processing capacity. On the other hand, dim-light vision may evolve to maximize light sensitivity at the cost of acuity and colour sensitivity, which could require less processing power. Here, we use X-ray microtomography and phylogenetic comparative methods to examine the relationships between diel activity pattern, optic morphology, trophic guild and investment in the optic tectum across the largest radiation of vertebrates-teleost fishes. We find that despite driving the evolution of larger eyes, enhancement of the capacity for dim-light vision generally is accompanied by a decrease in investment in the optic tectum. These findings underscore the importance of considering diel activity patterns in comparative studies and demonstrate how vision plays a role in brain evolution, illuminating common design principles of the vertebrate visual system.

Testing ontogenetic patterns of sexual size dimorphism against expectations of the expensive tissue hypothesis, an intraspecific example using oyster toadfish (Opsanus tau).

Trade-offs associated with sexual size dimorphism (SSD) are well documented across the Tree of Life. However, studies of SSD often do not consider potential investment trade-offs between metabolically expensive structures under sexual selection and other morphological modules. Based on the expectations of the expensive tissue hypothesis, investment in one metabolically expensive structure should come at the direct cost of investment in another. Here, we examine allometric trends in the ontogeny of oyster toadfish (Opsanus tau) to test whether investment in structures known to have been influenced by strong sexual selection conform to these expectations. Despite recovering clear changes in the ontogeny of a sexually selected trait between males and females, we find no evidence for predicted ontogenetic trade-offs with metabolically expensive organs. Our results are part of a growing body of work demonstrating that increased investment in one structure does not necessarily drive a wholesale loss of mass in one or more organs.

Identifying early modern human ecological niche expansions and associated cultural dynamics in the South African Middle Stone Age.

The archaeological record shows that typically human cultural traits emerged at different times, in different parts of the world, and among different hominin taxa. This pattern suggests that their emergence is the outcome of complex and nonlinear evolutionary trajectories, influenced by environmental, demographic, and social factors, that need to be understood and traced at regional scales. The application of predictive algorithms using archaeological and paleoenvironmental data allows one to estimate the ecological niches occupied by past human populations and identify niche changes through time, thus providing the possibility of investigating relationships between cultural innovations and possible niche shifts. By using such methods to examine two key southern Africa archaeological cultures, the Still Bay [76-71 thousand years before present (ka)] and the Howiesons Poort (HP; 66-59 ka), we identify a niche shift characterized by a significant expansion in the breadth of the HP ecological niche. This expansion is coincident with aridification occurring across Marine Isotope Stage 4 (ca. 72-60 ka) and especially pronounced at 60 ka. We argue that this niche shift was made possible by the development of a flexible technological system, reliant on composite tools and cultural transmission strategies based more on "product copying" rather than "process copying." These results counter the one niche/one human taxon equation. They indicate that what makes our cultures, and probably the cultures of other members of our lineage, unique is their flexibility and ability to produce innovations that allow a population to shift its ecological niche.

RWTY (R We There Yet): An R Package for Examining Convergence of Bayesian Phylogenetic Analyses.

Bayesian inference using Markov chain Monte Carlo (MCMC) has become one of the primary methods used to infer phylogenies from sequence data. Assessing convergence is a crucial component of these analyses, as it establishes the reliability of the posterior distribution estimates of the tree topology and model parameters sampled from the MCMC. Numerous tests and visualizations have been developed for this purpose, but many of the most popular methods are implemented in ways that make them inconvenient to use for large data sets. RWTY is an R package that implements established and new methods for diagnosing phylogenetic MCMC convergence in a single convenient interface.

Estimating the Effective Sample Size of Tree Topologies from Bayesian Phylogenetic Analyses.

Bayesian phylogenetic analyses estimate posterior distributions of phylogenetic tree topologies and other parameters using Markov chain Monte Carlo (MCMC) methods. Before making inferences from these distributions, it is important to assess their adequacy. To this end, the effective sample size (ESS) estimates how many truly independent samples of a given parameter the output of the MCMC represents. The ESS of a parameter is frequently much lower than the number of samples taken from the MCMC because sequential samples from the chain can be non-independent due to autocorrelation. Typically, phylogeneticists use a rule of thumb that the ESS of all parameters should be greater than 200. However, we have no method to calculate an ESS of tree topology samples, despite the fact that the tree topology is often the parameter of primary interest and is almost always central to the estimation of other parameters. That is, we lack a method to determine whether we have adequately sampled one of the most important parameters in our analyses. In this study, we address this problem by developing methods to estimate the ESS for tree topologies. We combine these methods with two new diagnostic plots for assessing posterior samples of tree topologies, and compare their performance on simulated and empirical data sets. Combined, the methods we present provide new ways to assess the mixing and convergence of phylogenetic tree topologies in Bayesian MCMC analyses.

Projecting future expansion of invasive species: comparing and improving methodologies for species distribution modeling.

Modeling the distributions of species, especially of invasive species in non-native ranges, involves multiple challenges. Here, we developed some novel approaches to species distribution modeling aimed at reducing the influences of such challenges and improving the realism of projections. We estimated species-environment relationships for Parthenium hysterophorus L. (Asteraceae) with four modeling methods run with multiple scenarios of (i) sources of occurrences and geographically isolated background ranges for absences, (ii) approaches to drawing background (absence) points, and (iii) alternate sets of predictor variables. We further tested various quantitative metrics of model evaluation against biological insight. Model projections were very sensitive to the choice of training dataset. Model accuracy was much improved using a global dataset for model training, rather than restricting data input to the species' native range. AUC score was a poor metric for model evaluation and, if used alone, was not a useful criterion for assessing model performance. Projections away from the sampled space (i.e., into areas of potential future invasion) were very different depending on the modeling methods used, raising questions about the reliability of ensemble projections. Generalized linear models gave very unrealistic projections far away from the training region. Models that efficiently fit the dominant pattern, but exclude highly local patterns in the dataset and capture interactions as they appear in data (e.g., boosted regression trees), improved generalization of the models. Biological knowledge of the species and its distribution was important in refining choices about the best set of projections. A post hoc test conducted on a new Parthenium dataset from Nepal validated excellent predictive performance of our 'best' model. We showed that vast stretches of currently uninvaded geographic areas on multiple continents harbor highly suitable habitats for parthenium. However, discrepancies between model predictions and parthenium invasion in Australia indicate successful management for this globally significant weed.

Mistaking geography for biology: inferring processes from species distributions.

Over the past few decades, there has been a rapid proliferation of statistical methods that infer evolutionary and ecological processes from data on species distributions. These methods have led to considerable new insights, but they often fail to account for the effects of historical biogeography on present-day species distributions. Because the geography of speciation can lead to patterns of spatial and temporal autocorrelation in the distributions of species within a clade, this can result in misleading inferences about the importance of deterministic processes in generating spatial patterns of biodiversity. In this opinion article, we discuss ways in which patterns of species distributions driven by historical biogeography are often interpreted as evidence of particular evolutionary or ecological processes. We focus on three areas that are especially prone to such misinterpretations: community phylogenetics, environmental niche modelling, and analyses of beta diversity (compositional turnover of biodiversity).

Molecular phylogenetics of squirrelfishes and soldierfishes (Teleostei: Beryciformes: Holocentridae): reconciling more than 100 years of taxonomic confusion.

Squirrelfishes and soldierfishes (Holocentridae) are among the most conspicuous species in the nocturnal reef fish community. However, there is no clear consensus regarding their evolutionary relationships, which is reflected in a complicated taxonomic history. We collected DNA sequence data from multiple single copy nuclear genes and one mitochondrial gene sampled from over fifty percent of the recognized holocentrid species and infer the first species-level phylogeny of the Holocentridae. Our results strongly support the monophyly of the clades Myripristinae (soldierfishes) and Holocentrinae (squirrelfishes). The molecular phylogenies differ with regard to previous hypotheses of relationships within the Myriprisitinae, resolving a clade of cryptic reef associated and deep water non-reef dwelling lineages (Corniger+Plectrypops+Ostichthys) that is the sister lineage to a monophyletic Myripristis. Within Holocentrinae, Neoniphon and Sargocentron are strongly supported as paraphyletic, while Holocentrus is nested within Sargocentron. Using Bayesian ancestral state reconstruction methods, we demonstrate the taxonomically diagnostic characters for Neoniphon and Sargocentron likely represent character states with a complex evolutionary history that is not reflective of shared common ancestry. We propose a new classification for Holocentrinae, recognizing four lineages that are treated as genera: Sargocentron Fowler, 1904, Holocentrus Scopoli, 1777, Flameo Jordan and Evermann, 1898, and Neoniphon Castelnau, 1875.

In defense of 'niche modeling'.

There is a growing awareness of problems with the estimation of the ecological tolerances of species through correlative modeling approaches. These problems have led some investigators to argue for abandoning terms such as 'ecological niche model' and 'environmental niche model' in favor of the ostensibly more value-neutral 'species distribution model', as the models are thought to frequently be poor estimators of the niche. Here, I argue that most applications to which these models are put require the assumption that they do estimate the niche, however imperfectly, and that obscuring this inescapable and potentially flawed assumption in the terminology may only serve to hinder the development of the field.

Uterine gene expression in the live-bearing lizard, Chalcides ocellatus, reveals convergence of squamate reptile and mammalian pregnancy mechanisms.

Although the morphological and physiological changes involved in pregnancy in live-bearing reptiles are well studied, the genetic mechanisms that underlie these changes are not known. We used the viviparous African Ocellated Skink, Chalcides ocellatus, as a model to identify a near complete gene expression profile associated with pregnancy using RNA-Seq analyses of uterine transcriptomes. Pregnancy in C. ocellatus is associated with upregulation of uterine genes involved with metabolism, cell proliferation and death, and cellular transport. Moreover, there are clear parallels between the genetic processes associated with pregnancy in mammals and Chalcides in expression of genes related to tissue remodeling, angiogenesis, immune system regulation, and nutrient provisioning to the embryo. In particular, the pregnant uterine transcriptome is dominated by expression of proteolytic enzymes that we speculate are involved both with remodeling the chorioallantoic placenta and histotrophy in the omphaloplacenta. Elements of the maternal innate immune system are downregulated in the pregnant uterus, indicating a potential mechanism to avoid rejection of the embryo. We found a downregulation of major histocompatability complex loci and estrogen and progesterone receptors in the pregnant uterus. This pattern is similar to mammals but cannot be explained by the mammalian model. The latter finding provides evidence that pregnancy is controlled by different endocrinological mechanisms in mammals and reptiles. Finally, 88% of the identified genes are expressed in both the pregnant and the nonpregnant uterus, and thus, morphological and physiological changes associated with C. ocellatus pregnancy are likely a result of regulation of genes continually expressed in the uterus rather than the initiation of expression of unique genes.

Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria.

Maxent, one of the most commonly used methods for inferring species distributions and environmental tolerances from occurrence data, allows users to fit models of arbitrary complexity. Model complexity is typically constrained via a process known as L1 regularization, but at present little guidance is available for setting the appropriate level of regularization, and the effects of inappropriately complex or simple models are largely unknown. In this study, we demonstrate the use of information criterion approaches to setting regularization in Maxent, and we compare models selected using information criteria to models selected using other criteria that are common in the literature. We evaluate model performance using occurrence data generated from a known "true" initial Maxent model, using several different metrics for model quality and transferability. We demonstrate that models that are inappropriately complex or inappropriately simple show reduced ability to infer habitat quality, reduced ability to infer the relative importance of variables in constraining species' distributions, and reduced transferability to other time periods. We also demonstrate that information criteria may offer significant advantages over the methods commonly used in the literature.

Evolution of cold-tolerant fungal symbionts permits winter fungiculture by leafcutter ants at the northern frontier of a tropical ant-fungus symbiosis.

The obligate mutualism between leafcutter ants and their Attamyces fungi originated 8 to 12 million years ago in the tropics, but extends today also into temperate regions in South and North America. The northernmost leafcutter ant Atta texana sustains fungiculture during winter temperatures that would harm the cold-sensitive Attamyces cultivars of tropical leafcutter ants. Cold-tolerance of Attamyces cultivars increases with winter harshness along a south-to-north temperature gradient across the range of A. texana, indicating selection for cold-tolerant Attamyces variants along the temperature cline. Ecological niche modeling corroborates winter temperature as a key range-limiting factor impeding northward expansion of A. texana. The northernmost A. texana populations are able to sustain fungiculture throughout winter because of their cold-adapted fungi and because of seasonal, vertical garden relocation (maintaining gardens deep in the ground in winter to protect them from extreme cold, then moving gardens to warmer, shallow depths in spring). Although the origin of leafcutter fungiculture was an evolutionary breakthrough that revolutionized the food niche of tropical fungus-growing ants, the original adaptations of this host-microbe symbiosis to tropical temperatures and the dependence on cold-sensitive fungal symbionts eventually constrained expansion into temperate habitats. Evolution of cold-tolerant fungi within the symbiosis relaxed constraints on winter fungiculture at the northern frontier of the leafcutter ant distribution, thereby expanding the ecological niche of an obligate host-microbe symbiosis.

Ecological and geographic modes of species divergence in wild tomatoes.

Understanding the role of geography and ecology in species divergence is central to the study of evolutionary diversification. We used climatic, geographic, and biological data from nine wild Andean tomato species to describe each species' ecological niche and to evaluate the likely ecological and geographical modes of speciation in this clade. Using data from >1000 wild accessions and publicly available data derived from geographic information systems for various environmental variables, we found most species pairs were significantly differentiated for one or more environmental variables. By comparing species' predicted niches generated by species distribution modeling (SDM), we found significant niche differentiation among three of four sister-species pairs, suggesting ecological divergence is consistently associated with recent divergence. In comparison, based on age-range correlation (ARC) analysis, there was no evidence for a predominant geographical (allopatric vs. sympatric) context for speciation in this group. Overall, our results suggest an important role for environmentally mediated differentiation, rather than simply geographical isolation, in species divergence.

Conclusions about niche expansion in introduced Impatiens walleriana populations depend on method of analysis.

Determining the degree to which climate niches are conserved across plant species' native and introduced ranges is valuable to developing successful strategies to limit the introduction and spread of invasive plants, and also has important ecological and evolutionary implications. Here, we test whether climate niches differ between native and introduced populations of Impatiens walleriana, globally one of the most popular horticultural species. We use approaches based on both raw climate data associated with occurrence points and ecological niche models (ENMs) developed with Maxent. We include comparisons of climate niche breadth in both geographic and environmental spaces, taking into account differences in available habitats between the distributional areas. We find significant differences in climate envelopes between native and introduced populations when comparing raw climate variables, with introduced populations appearing to expand into wetter and cooler climates. However, analyses controlling for differences in available habitat in each region do not indicate expansion of climate niches. We therefore cannot reject the hypothesis that observed differences in climate envelopes reflect only the limited environments available within the species' native range in East Africa. Our results suggest that models built from only native range occurrence data will not provide an accurate prediction of the potential for invasiveness if applied to areas containing a greater range of environmental combinations, and that tests of niche expansion may overestimate shifts in climate niches if they do not control carefully for environmental differences between distributional areas.