Phylogenetic Models Research Articles

Comparing Partitioned Models to Mixture Models: Do Information Criteria Apply?

The use of information criteria to distinguish between phylogenetic models has become ubiquitous within the field. However, the variety and complexity of available models are much greater now than when these practices were established. The literature shows an increasing trajectory of healthy skepticism with regard to the use of information theory-based model selection within phylogenetics. We add to this by analyzing the specific case of comparison between partition and mixture models. We argue from a theoretical basis that information criteria are inherently more likely to favor partition models over mixture models, and we then demonstrate this through simulation. Based on our findings, we suggest that partition and mixture models are not suitable for information-theory based model comparison. [AIC, BIC; information criteria; maximum likelihood; mixture models; partitioned model; phylogenetics.]

An Improved Codon Modeling Approach for Accurate Estimation of the Mutation Bias.

Phylogenetic codon models are routinely used to characterize selective regimes in coding sequences. Their parametric design, however, is still a matter of debate, in particular concerning the question of how to account for differing nucleotide frequencies and substitution rates. This problem relates to the fact that nucleotide composition in protein-coding sequences is the result of the interactions between mutation and selection. In particular, because of the structure of the genetic code, the nucleotide composition differs between the three coding positions, with the third position showing a more extreme composition. Yet, phylogenetic codon models do not correctly capture this phenomenon and instead predict that the nucleotide composition should be the same for all three positions. Alternatively, some models allow for different nucleotide rates at the three positions, an approach conflating the effects of mutation and selection on nucleotide composition. In practice, it results in inaccurate estimation of the strength of selection. Conceptually, the problem comes from the fact that phylogenetic codon models do not correctly capture the fixation bias acting against the mutational pressure at the mutation–selection equilibrium. To address this problem and to more accurately identify mutation rates and selection strength, we present an improved codon modeling approach where the fixation rate is not seen as a scalar, but as a tensor. This approach gives an accurate representation of how mutation and selection oppose each other at equilibrium and yields a reliable estimate of the mutational process, while disentangling the mean fixation probabilities prevailing in different mutational directions.

Pleistocene mitogenomes reconstructed from the environmental DNA of permafrost sediments

Traditionally, paleontologists have relied on the morphological features of bones and teeth to reconstruct the evolutionary relationships of extinct animals.1 In recent decades, the analysis of ancient DNA recovered from macrofossils has provided a powerful means to evaluate these hypotheses and develop novel phylogenetic models.2 Although a great deal of life history data can be extracted from bones, their scarcity and associated biases limit their information potential. The paleontological record of Beringia3-the unglaciated areas and former land bridge between northeast Eurasia and northwest North America-is relatively robust thanks to its perennially frozen ground favoring fossil preservation.4,5 However, even here, the macrofossil record is significantly lacking in small-bodied fauna (e.g., rodents and birds), whereas questions related to migration and extirpation, even among well-studied taxa, remain crudely resolved. The growing sophistication of ancient environmental DNA (eDNA) methods have allowed for the identification of species within terrestrial/aquatic ecosystems,6-12 in paleodietary reconstructions,13-19 and facilitated genomic reconstructions from cave contexts.8,20-22 Murchie etal.6,23 used a capture enrichment approach to sequence a diverse range of faunal and floral DNA from permafrost silts deposited during the Pleistocene-Holocene transition.24 Here, we expand on their work with the mitogenomic assembly and phylogenetic placement of Equus caballus (caballine horse), Bison priscus (steppe bison), Mammuthus primigenius (woolly mammoth), and Lagopus lagopus (willow ptarmigan) eDNA from multiple permafrost cores spanning the last 40,000 years. We identify a diverse metagenomic spectra of Pleistocene fauna and identify the eDNA co-occurrence of distinct Eurasian and American mitogenomic lineages.

Nucleotide Usage Biases Distort Inferences of the Species Tree.

Despite the importance of natural selection in species’ evolutionary history, phylogenetic methods that take into account population-level processes typically ignore selection. The assumption of neutrality is often based on the idea that selection occurs at a minority of loci in the genome and is unlikely to compromise phylogenetic inferences significantly. However, genome-wide processes like GC-bias and some variation segregating at the coding regions are known to evolve in the nearly neutral range. As we are now using genome-wide data to estimate species trees, it is natural to ask whether weak but pervasive selection is likely to blur species tree inferences. We developed a polymorphism-aware phylogenetic model tailored for measuring signatures of nucleotide usage biases to test the impact of selection in the species tree. Our analyses indicate that although the inferred relationships among species are not significantly compromised, the genetic distances are systematically underestimated in a node-height-dependent manner: that is, the deeper nodes tend to be more underestimated than the shallow ones. Such biases have implications for molecular dating. We dated the evolutionary history of 30 worldwide fruit fly populations, and we found signatures of GC-bias considerably affecting the estimated divergence times (up to 23%) in the neutral model. Our findings call for the need to account for selection when quantifying divergence or dating species evolution.

Phylogenetic analysis of combined mitochondrial genome and 32 nuclear genes provides key insights into molecular systematics and historical biogeography of Asian warty newts of the genus Paramesotriton (Caudata: Salamandridae).

The Paramesotriton Chang, 1935 genus of Asian warty newts is the second most diverse genus in the family Salamandridae, currently containing 14 recognized species from northern Vietnam to southwest-central and southern China. Although species of this genus have been included in previous phylogenetic studies, the origin and interspecific relationships of the genus are still not fully resolved, especially at key nodes in the phylogeny. In this study, we sequenced mitochondrial genomes and 32 nuclear genes from 27 samples belonging to 14 species to reconstruct the interspecific phylogenetic relationships within Paramesotriton and explore its historical biogeography in southern China. Both Bayesian inference and maximum-likelihood analyses highly supported the monophyly of Paramesotriton and its two recognized species groups (P. caudopunctatus and P. chinensis groups) and further identified five hypothetical phylogenetic cryptic species. Biogeographic analyses indicated that Paramesotriton originated in southwestern China (Yunnan-Guizhou Plateau/South China) during the late Oligocene. The time of origin of Paramesotriton corresponded to the second uplift of the Himalayan/Qinghai-Xizang (Tibetan) Plateau (QTP), rapid lateral extrusion of Indochina, and formation of karst landscapes in southwestern China. Principal component analysis (PCA), independent sample t-tests, and niche differentiation using bioclimatic variables based on locations of occurrence suggested that Paramesotriton habitat conditions in the three current regions (West, South, and East) differ significantly, with different levels of climatic niche differentiation. Species distribution model (SDM) predictions indicated that the most suitable distribution areas for the P. caudopunctatus and P. chinensis species groups are western and southern/eastern areas of southern China. This study increases our knowledge of the taxonomy, biodiversity, origin, and suitable distribution areas of the genus Paramesotriton based on phylogenetic, biogeographic, and species distribution models.

RAxML Grove: an empirical phylogenetic tree database.

SummaryThe assessment of novel phylogenetic models and inference methods is routinely being conducted via experiments on simulated as well as empirical data. When generating synthetic data it is often unclear how to set simulation parameters for the models and generate trees that appropriately reflect empirical model parameter distributions and tree shapes. As a solution, we present and make available a new database called ‘RAxML Grove’ currently comprising more than 60 000 inferred trees and respective model parameter estimates from fully anonymized empirical datasets that were analyzed using RAxML and RAxML-NG on two web servers. We also describe and make available two simple applications of RAxML Grove to exemplify its usage and highlight its utility for designing realistic simulation studies and analyzing empirical model parameter and tree shape distributions.Availability and implementationRAxML Grove is freely available at https://github.com/angtft/RAxMLGrove.Supplementary information Supplementary data are available at Bioinformatics online.

Both source‐ and recipient‐range phylogenetic community structure can predict the outcome of avian introductions

Competing phylogenetic models have been proposed to explain the success of species introduced to other communities. Here, we present a study predicting the establishment success of birds introduced to Florida, Hawaii and New Zealand using several alternative models, considering species' phylogenetic relatedness to source‐ and recipient‐range taxa, propagule pressure and traits. We find consistent support for the predictive ability of source‐region phylogenetic structure. However, we find that the effects of recipient‐region phylogenetic structure vary in sign and magnitude depending on inclusion of source‐region phylogenetic structure, delineation of the recipient species pool and the use of phylogenetic correction in the models. We argue that tests of alternative phylogenetic hypotheses including both source and recipient community phylogenetic structure, as well as important covariates such as propagule pressure, are likely to be critical for identifying general phylogenetic patterns in introduction success, predicting future invasions and for stimulating further exploration of the underlying mechanisms of invasibility.

Climate change aggravates non-target effects of pesticides on dragonflies at macroecological scales.

Critical gaps in understanding how species respond to environmental change limit our capacity to address conservation risks in a timely way. Here, we examine the direct and interactive effects of key global change drivers, including climate change, land use change, and pesticide use, on persistence of 104 odonate species between two time periods (1980-2002 and 2008-2018) within 100 × 100 km quadrats across the USA using phylogenetic mixed models. Non-target effects of pesticides interacted with higher maximum temperatures to contribute to odonate declines. Closely related species responded similarly to global change drivers, indicating a potential role of inherited traits in species' persistence or decline. Species shifting their range to higher latitudes were more robust to negative impacts of global change drivers generally. Inherited traits related to dispersal abilities and establishment in new places may govern both species' acclimation to global change and their abilities to expand their range limits, respectively. This work is among the first to assess effects of climate change, land use change, and land use intensification together on Odonata, a significant step that improves understanding of multispecies effects of global change on invertebrates, and further identifies conditions contributing to global insect loss.

You are more than what you eat: potentially adaptive enrichment of microbiome functions across bat dietary niches.

BackgroundAnimals evolved in a microbial world, and their gut microbial symbionts have played a role in their ecological diversification. While many recent studies report patterns of phylosymbiosis between hosts and their gut bacteria, fewer studies examine the potentially adaptive functional contributions of these microbes to the dietary habits of their hosts. In this study, we examined predicted metabolic pathways in the gut bacteria of more than 500 individual bats belonging to 60 species and compare the enrichment of these functions across hosts with distinct dietary ecologies.ResultsWe found that predicted microbiome functions were differentially enriched across hosts with different diets. Using a machine-learning approach, we also found that inferred microbiome functions could be used to predict specialized host diets with reasonable accuracy. We detected a relationship between both host phylogeny and diet with respect to microbiome functional repertoires. Because many predicted functions could potentially fill nutritional gaps for bats with specialized diets, we considered pathways discriminating dietary niches as traits of the host and fit them to comparative phylogenetic models of evolution. Our results suggest that some, but not all, predicted microbiome functions may evolve toward adaptive optima and thus be visible to the forces of natural selection operating on hosts over evolutionary time.ConclusionsOur results suggest that bats with specialized diets may partially rely on their gut microbes to fulfill or augment critical nutritional pathways, including essential amino acid synthesis, fatty acid biosynthesis, and the generation of cofactors and vitamins essential for proper nutrition. Our work adds to a growing body of literature suggesting that animal microbiomes are structured by a combination of ecological and evolutionary processes and sets the stage for future metagenomic and metabolic characterization of the bat microbiome to explore links between bacterial metabolism and host nutrition.

Field surveys in Croatia and North Macedonia reveal two novel phleboviruses circulating in sandflies.

Sandfly-borne phleboviruses are distributed widely throughout the Mediterranean Basin, presenting a threat to public health in areas where they circulate. However, the true diversity and distribution of pathogenic and apathogenic sandfly-borne phleboviruses remains a key issue to be studied. In the Balkans, most published data rely on serology-based studies although virus isolation has occasionally been reported. Here, we report the discovery of two novel sandfly-borne phleboviruses, provisionally named Zaba virus (ZABAV) and Bregalaka virus (BREV), which were isolated in Croatia and North Macedonia, respectively. This constitutes the first isolation of phleboviruses in both countries. Genetic analysis based on complete coding sequences indicated that ZABAV and BREV are distinct from each other and belong to the genus Phlebovirus, family Phenuiviridae. Phylogenetic and amino acid modelling of viral polymerase shows that ZABAV and BREV are new members of the Salehabad phlebovirus species and the Adana phlebovirus species, respectively. Moreover, sequence-based vector identification suggests that ZABAV is mainly transmitted by Phlebotomus neglectus and BREV is mainly transmitted by Phlebotomus perfiliewi. BREV neutralizing antibodies were detected in 3.3% of human sera with rates up to 16.7% in certain districts, demonstrating that BREV frequently infects humans in North Macedonia. In vitro viral growth kinetics experiments demonstrated viral replication of both viruses in mammalian and mosquito cells. In vivo experimental studies in mice suggest that ZABAV and BREV exhibit characteristics making them possible human pathogens.

Tube-dwelling in early animals exemplified by Cambrian scalidophoran worms

BackgroundThe radiation of ecdysozoans (moulting animals) during the Cambrian gave rise to panarthropods and various groups of worms including scalidophorans, which played an important role in the elaboration of early marine ecosystems. Although most scalidophorans were infaunal burrowers travelling through soft sediment at the bottom of the sea, Selkirkia lived inside a tube.ResultsWe explore the palaeobiology of these tubicolous worms, and more generally the origin and evolutionary significance of tube-dwelling in early animals, based on exceptionally preserved fossils from the early Cambrian Chengjiang Lagerstätte (Stage 3, China) including a new species, Selkirkia transita sp. nov. We find that the best phylogenetic model resolves Selkirkia as a stem-group priapulid. Selkirkia secreted a protective cuticular thickening, the tube, inside which it was able to move during at least part of its life. Partly based on measured growth patterns, we construe that this tube was separated from the trunk during a moulting process that has no direct equivalent in other scalidophorans. Although the ontogeny of Selkirkia is currently unknown, we hypothesize that its conical tube might have had the same ecological function and possibly even deep development origin as the lorica, a protective cuticular thickening found in larval priapulids and adult loriciferans. Selkirkia is seen as a semi-sedentary animal capable of very shallow incursions below the water/sediment interface, possibly for feeding or during the tube-secreting phase. Brachiopod epibionts previously reported from the Xiaoshiba Lagerstätte (ca. 514 Ma) also presumably occur in Selkirkia sinica from Chengjiang (ca. 518 Ma).ConclusionsOur critical and model-based approach provides a new phylogenetic framework for Scalidophora, upon which to improve in order to study the evolution of morphological characters in this group. Tube-dwelling is likely to have offered Selkirkia better protection and anchoring to sediment and has developed simultaneously in other Cambrian animals such as hemichordates, annelids or panarthropods. Often lost in modern representatives in favour of active infaunal lifestyles, tube-dwelling can be regarded as an early evolutionary response of various metazoans to increasing environmental and biological pressure in Cambrian marine ecosystems.

Multi-model approaches to phylogenetics: Implications for idealization

Phylogenetic models traditionally represent the history of life as having a strictly-branching tree structure. However, it is becoming increasingly clear that the history of life is often not strictly-branching; lateral gene transfer, endosymbiosis, and hybridization, for example, can all produce lateral branching events. There is thus motivation to allow phylogenetic models to have a reticulate structure. One proposal involves the reconciliation of genealogical discordance. Briefly, this method uses patterns of disagreement – discordance – between trees of different genes to add lateral branching events to phylogenetic trees of taxa, and to estimate the most likely cause of these events. I use this practice to argue for: (1) a need for expanded accounts of multiple-models idealization, (2) a distinction between automatic and manual de-idealization, and (3) recognition that idealization may serve the meso-level aims of science in a different way than hitherto acknowledged.

Evolvability in the fossil record

AbstractThe concept of evolvability—the capacity of a population to produce and maintain evolutionarily relevant variation—has become increasingly prominent in evolutionary biology. Paleontology has a long history of investigating questions of evolvability, but paleontological thinking has tended to neglect recent discussions, because many tools used in the current evolvability literature are challenging to apply to the fossil record. The fundamental difficulty is how to disentangle whether the causes of evolutionary patterns arise from variational properties of traits or lineages rather than being due to selection and ecological success. Despite these obstacles, the fossil record offers unique and growing sources of data that capture evolutionary patterns of sustained duration and significance otherwise inaccessible to evolutionary biologists. Additionally, there exist a variety of strategic possibilities for combining prominent neontological approaches to evolvability with those from paleontology. We illustrate three of these possibilities with quantitative genetics, evolutionary developmental biology, and phylogenetic models of macroevolution. In conclusion, we provide a methodological schema that focuses on the conceptualization, measurement, and testing of hypotheses to motivate and provide guidance for future empirical and theoretical studies of evolvability in the fossil record.

Environmental Water and Sediment Microbial Communities Shape Intestine Microbiota for Host Health: The Central Dogma in an Anthropogenic Aquaculture Ecosystem.

From increasing evidence has emerged a tight link among the environment, intestine microbiota, and host health status; moreover, the microbial interaction in different habitats is crucial for ecosystems. However, how the environmental microbial community assembly governs the intestinal microbiota and microbial communities of multiple habitats contribute to the metacommunity remain elusive. Here, we designed two delicate experiments from temporal and spatial scales in a shrimp culture pond ecosystem (SCPE). Of the SCPE metacommunity, the microbial diversity was mainly contributed to by the diversity of–βIntraHabitats and βInterHabitats, and water and sediment communities had a large contribution to the shrimp intestine community as shown by SourceTracker and Sloan neutral community model analyses. Also, phylogenetic bin-based null model results show that microbial assembly of three habitats in the SCPE appeared to be largely driven by stochastic processes. These results enrich our understanding of the environment–intestinal microbiota–host health closely linked relationship, making it possible to be the central dogma for an anthropogenic aquaculture ecosystem. Our findings enhance the mechanistic understanding of microbial assembly in the SCPE for further analyzing metacommunities, which has important implications for microbial ecology and animal health.

Genome-Wide Identification and Transcriptional Analysis of Arabidopsis DUF506 Gene Family.

The Domain of unknown function 506 (DUF506) family, which belongs to the PD-(D/E)XK nuclease superfamily, has not been functionally characterized. In this study, 266 DUF506 domain-containing genes were identified from algae, mosses, and land plants showing their wide occurrence in photosynthetic organisms. Bioinformatics analysis identified 211 high-confidence DUF506 genes across 17 representative land plant species. Phylogenetic modeling classified three groups of plant DUF506 genes that suggested functional preservation among the groups based on conserved gene structure and motifs. Gene duplication and Ka/Ks evolutionary rates revealed that DUF506 genes are under purifying positive selection pressure. Subcellular protein localization analysis revealed that DUF506 proteins were present in different organelles. Transcript analyses showed that 13 of the Arabidopsis DUF506 genes are ubiquitously expressed in various tissues and respond to different abiotic stresses and ABA treatment. Protein-protein interaction network analysis using the STRING-DB, AtPIN (Arabidopsis thaliana Protein Interaction Network), and AI-1 (Arabidopsis Interactome-1) tools indicated that AtDUF506s potentially interact with iron-deficiency response proteins, salt-inducible transcription factors, or calcium sensors (calmodulins), implying that DUF506 genes have distinct biological functions including responses to environmental stimuli, nutrient-deficiencies, and participate in Ca(2+) signaling. Current results provide insightful information regarding the molecular features of the DUF506 family in plants, to support further functional characterizations.

Influence of historical changes in tropical reef habitat on the diversification of coral reef fishes

Past environmental changes are expected to have profoundly impacted diversity dynamics through time. While some previous studies showed an association between past climate changes or tectonic events and important shifts in lineage diversification, it is only recently that past environmental changes have been explicitly integrated in diversification models to test their influence on diversification rates. Here, we used a global reconstruction of tropical reef habitat dynamics during the Cenozoic and phylogenetic diversification models to test the influence of (i) major geological events, (ii) reef habitat fragmentation and (iii) reef area on the diversification of 9 major clades of tropical reef fish (Acanthuridae, Balistoidea, Carangoidea, Chaetodontidae, Haemulinae, Holocentridae, Labridae, Pomacentridae and Sparidae). The diversification models revealed a weak association between paleo-habitat changes and diversification dynamics. Specifically, the fragmentation of tropical reef habitats over the Cenozoic was found to be a driver of tropical reef fish diversification for 2 clades. However, overall, our approach did not allow the identification of striking associations between diversification dynamics and paleo-habitat fragmentation in contrast with theoretical model’s predictions.

Human density is associated with the increased prevalence of a generalist zoonotic parasite in mammalian wildlife.

Macroecological approaches can provide valuable insight into the epidemiology of globally distributed, multi-host pathogens. Toxoplasma gondii is a zoonotic protozoan that infects any warm-blooded animal, including humans, in almost every ecosystem worldwide. There is substantial geographical variation in T. gondii prevalence in wildlife populations and the mechanisms driving this variation are poorly understood. We implemented Bayesian phylogenetic mixed models to determine the association between species’ ecology, phylogeny and climatic and anthropogenic factors on T. gondii prevalence. Toxoplasma gondii prevalence data were compiled for free-ranging wild mammal species from 202 published studies, encompassing 45 079 individuals from 54 taxonomic families and 238 species. We found that T. gondii prevalence was positively associated with human population density and warmer temperatures at the sampling location. Terrestrial species had a lower overall prevalence, but there were no consistent patterns between trophic level and prevalence. The relationship between human density and T. gondii prevalence is probably mediated by higher domestic cat abundance and landscape degradation leading to increased environmental oocyst contamination. Landscape restoration and limiting free-roaming in domestic cats could synergistically increase the resiliency of wildlife populations and reduce wildlife and human infection risks from one of the world's most common parasitic infections.

Spatially resolved assembly, connectivity and structure of particle-associated and free-living bacterial communities in a high Arctic fjord.

ABSTRACTThe assembly processes that underlie the composition and connectivity of free-living (FL) and particle-associated (PA) bacterial communities from surface to deep waters remain little understood. Here, using phylogenetic null modeling, we quantify the relative influence of selective and stochastic mechanisms that assemble FL and PA bacterial communities throughout the water column in a high Arctic fjord. We demonstrate that assembly processes acting on FL and PA bacterial communities are similar in surface waters, but become increasingly distinct in deep waters. As depth increases, the relative influence of homogeneous selection increases for FL but decreases for PA communities. In addition, dispersal limitation and variable selection increase with depth for PA, but not for FL communities, indicating increased residence time of taxa on particles and less frequent decolonization. As a consequence, beta diversity of PA communities is greater in bottom than in surface waters. The limited connectivity between these communities with increasing depth leads to highly distinct FL and PA bacterial communities in bottom waters. Finally, depth-related trends for FL and PA beta diversity and connectivity in this study are consistent with previous observations in the open ocean, suggesting that assembly processes for FL and PA bacterial communities may also be distinct in other aquatic environments.

Convergent evolution of skin surface microarchitecture and increased skin hydrophobicity in semi-aquatic anole lizards.

ABSTRACTAnimals that habitually cross the boundary between water and land face specific challenges with respect to locomotion, respiration, insulation, fouling and waterproofing. Many semi-aquatic invertebrates and plants have developed complex surface microstructures with water-repellent properties to overcome these problems, but equivalent adaptations of the skin have not been reported for vertebrates that encounter similar environmental challenges. Here, we document the first evidence of evolutionary convergence of hydrophobic structured skin in a group of semi-aquatic tetrapods. We show that the skin surface of semi-aquatic species of Anolis lizards is characterized by a more elaborate microstructural architecture (i.e. longer spines and spinules) and a lower wettability relative to closely related terrestrial species. In addition, phylogenetic comparative models reveal repeated independent evolution of enhanced skin hydrophobicity associated with the transition to a semi-aquatic lifestyle, providing evidence of adaptation. Our findings invite a new and exciting line of inquiry into the ecological significance, evolutionary origin and developmental basis of hydrophobic skin surfaces in semi-aquatic lizards, which is essential for understanding why and how the observed skin adaptations evolved in some and not other semi-aquatic tetrapod lineages.

Characterizing Amino Acid Substitution with Complete Linkage of Sites on a Lineage.

Amino acid substitution models are commonly used for phylogenetic inference, for ancestral sequence reconstruction, and for the inference of positive selection. All commonly used models explicitly assume that each site evolves independently, an assumption that is violated by both linkage and protein structural and functional constraints. We introduce two new models for amino acid substitution which incorporate linkage between sites, each based on the (population-genetic) Moran model. The first model is a generalized population process tracking arbitrarily many sites which undergo mutation, with individuals replaced according to their fitnesses. This model provides a reasonably complete framework for simulations but is numerically and analytically intractable. We also introduce a second model which includes several simplifying assumptions but for which some theoretical results can be derived. We analyze the simplified model to determine conditions where linkage is likely to have meaningful effects on sitewise substitution probabilities, as well as conditions under which the effects are likely to be negligible. These findings are an important step in the generation of tractable phylogenetic models that parameterize selective coefficients for amino acid substitution while accounting for linkage of sites leading to both hitchhiking and background selection.