Convergence assessment for Bayesian phylogenetic analysis using MCMC simulation

Robustness of divergence time estimation despite gene tree estimation error: a case study of fireflies (Coleoptera: Lampyridae)

Genomic data have become ubiquitous in phylogenomic studies, including divergence time estimation, but provide new challenges. These challenges include, among others, biological gene tree discordance, methodological gene tree estimation error, and computational limitations on performing full Bayesian inference under complex models. In this study, we use a recently published firefly (Coleoptera: Lampyridae) anchored hybrid enrichment data set (AHE; 436 loci for 88 Lampyridae species and 10 outgroup species) as a case study to explore gene tree estimation error and the robustness of divergence time estimation. First, we explored the amount of model violation using posterior predictive simulations because model violations are likely to bias phylogenetic inferences and produce gene tree estimation error. We specifically focused on missing data (either uniformly distributed or systematically) and the distribution of highly variable and conserved sites (either uniformly distributed or clustered). Our assessment of model adequacy showed that standard phylogenetic substitution models are not adequate for any of the 436 AHE loci. We tested if the model violations and alignment errors resulted indeed in gene tree estimation error by comparing the observed gene tree discordance to simulated gene tree discordance under the multispecies coalescent model. Thus, we show that the inferred gene tree discordance is not only due to biological mechanism but primarily due to inference errors. Lastly, we explored if divergence time estimation is robust despite the observed gene tree estimation error. We selected four subsets of the full AHE data set, concatenated each subset and performed a Bayesian relaxed clock divergence estimation in RevBayes. The estimated divergence times overlapped for all nodes that are shared between the topologies. Thus, divergence time estimation is robust using any well selected data subset as long as the topology inference is robust.

Evolutionary lability of a key innovation spurs rapid diversification

Rates of lineage diversification vary considerably across the tree of life, often as a result of evolutionary innovations1-5. Although the ability to produce new traits can vary between clades and may drive ecological transitions6-9, the impact of differences in the pace at which innovations evolve at macroevolutionary scales has been overlooked. Complex teeth are one innovation that contributed to the evolutionary success of major vertebrate lineages10-12. Here we show that evolutionary lability of tooth complexity, but not complexity itself, spurs rapid diversification across ray-finned fishes. Speciation rates are five times higher when transitions between simple and complex teeth occur rapidly. We find that African cichlids are unique among all fishes; they are dominated by lineages that transition between simple and complex teeth at unparalleled rates. This innovation interacted with the ecological versatility of complex teeth to spur rapid adaptive radiations in lakes Malawi, Victoria and Barombi Mbo. The marked effect on diversification stems from the tight association of tooth complexity with microhabitat and diet. Our results show that phylogenetic variation in how innovations evolve can have a stronger effect on patterns of diversification than the innovation itself. Investigating the impact of innovations from this new perspective will probably implicate more traits in causing heterogeneous diversification rates across the tree of life.

The effects of cryptic diversity on diversification dynamics analyses in Crocodylia

Incomplete taxon sampling due to underestimation of present-day biodiversity biases diversification analysis by favouring slowdowns in speciation rates towards the recent time. For instance, in diversification dynamics studies in Crocodylia, long-term low net-diversification rates and slowdowns in speciation rates have been suggested to characterize crocodylian evolution. However, crocodylian cryptic diversity has never been considered. Here, we explore the effects of incorporating cryptic diversity into a diversification dynamics analysis of extant crocodylians. We inferred a time-calibrated cryptic-species-level phylogeny using cytochrome b sequences of 45 lineages compared with the formally recognized 26 crocodylian species. Diversification rate estimates using the cryptic-species-level phylogeny show increasing speciation and net-diversification rates towards the present time, which contrasts with previous findings. Cryptic diversity should be considered in future macroevolutionary analyses; however, the representation of cryptic extinct taxa represents a major challenge. Additionally, further investigation of crocodylian diversification dynamics under different underlying genomic data is encouraged upon advances in population genetics. Our case study adds to the diversification dynamics knowledge of extant taxa and demonstrates that cryptic species and robust taxonomic assessment are essential to study recent biodiversity dynamics with broad implications for evolutionary biology and ecology.

How to Validate a Bayesian Evolutionary Model

Biology has become a highly mathematical discipline in which probabilistic models play a central role. As a result, research in the biological sciences is now dependent on computational tools capable of carrying out complex analyses. These tools must be validated before they can be used, but what is understood as validation varies widely among methodological contributions. This may be a consequence of the still embryonic stage of the literature on statistical software validation for computational biology. Our manuscript aims to advance this literature. Here, we describe, illustrate, and introduce new good practices for assessing the correctness of a model implementation with an emphasis on Bayesian methods. We also introduce a suite of functionalities for automating validation protocols. It is our hope that the guidelines presented here help sharpen the focus of discussions on (as well as elevate) expected standards of statistical software for biology.

Phylogenetic relationships of Neogene hamsters (Mammalia, Rodentia, Cricetinae) revealed under Bayesian inference and maximum parsimony

There is an ongoing debate about the internal systematics of today's group of hamsters (Cricetinae), following new insights that are gained based on molecular data. Regarding the closely related fossil cricetids, however, most studies deal with only a limited number of genera and statements about their possible relationships are rare. In this study, 41 fossil species from the Late Miocene to the Pliocene, belonging to seven extinct cricetine genera, Collimys, Rotundomys, Neocricetodon, Pseudocricetus, Cricetulodon, Apocricetus and Hattomys are analysed in a phylogenetic framework using traditional maximum parsimony and Bayesian inference approaches. Following thorough model testing, a relaxed-clock Bayesian inference analysis is performed under tip-dating to estimate divergence times simultaneously. Furthermore, so-called 'rogue' taxa are identified and excluded from the final trees to improve the informative value of the shown relationships. Based on these resulting trees, the fit of the topologies to the stratigraphy is assessed and the ancestral states of the characters are reconstructed under a parsimonious approach and stochastic character mapping. The overall topologies resulting from Bayesian and parsimonious approaches are largely congruent to each other and confirm the monophyly of most of the genera. Additionally, synapomorphies can be identified for each of these genera based on the ancestral state reconstructions. Only Cricetulodon turns out to be paraphyletic, while 'Cricetulodon' complicidens is a member of Neocricetodon. Lastly, this work makes a contribution to a debate that went on for decades, as the genus Kowalskia can be confirmed as junior synonym of Neocricetodon.

Phylogenetic analysis of the complete mitochondrial genome of the orange-winged sulphur butterfly Dercas nina Mell 1913 (Insecta: Lepidoptera: Pieridae: Coliadinae)

Dercas nina Mell 1913 (Pieridae) is a little-studied butterfly species endemic to China that flies primarily in the forest canopy. Genome skimming by Illumina sequencing allowed assembly of 146,702 reads for complete 1471.3-fold mean coverage of the circular 15,264 bp mitogenome from D. nina consisting of 82.1% AT nucleotides. A gene order typical of butterflies was recovered consisting of 13 protein-coding genes, 22 tRNAs, two rRNAs, and a predicted control region. The Dercas nina COX1 open reading frame begins with atypical start codon CGA. Six protein-coding genes (COX1, COX2, ND2, ND3, ND4, ND5) with single-nucleotide (T) stop codons, and two protein-coding genes (ATP6, ATP8) with two-nucleotide (TA) stop codons encoded in the DNA were inferred to be completed by adenine nucleotides from the Poly-A tail of the mRNA. Bayesian's phylogenetic reconstruction places the D. nina and D. lycorias mitogenomes as sister clades. Dercas mitogenomes were sister to those from genus Colias in the monophyletic subfamily Coliadinae. The mitogenome phylogeny is consistent with previous molecular phylogenetic hypotheses based on other markers, but differs somewhat from a morphology-based hypothesis that suggested that Dercas was more closely related to genus Gonepteryx. This may falsify the hypothesis or may instead reflect mitochondrial-nuclear phylogenetic discordance.

An Automated Convergence Diagnostic for Phylogenetic MCMC Analyses

Assessing convergence of Markov chain Monte Carlo (MCMC) based analyses is crucial but challenging, especially so in high dimensional and complex spaces such as the space of phylogenetic trees (treespace). In practice, it is assumed that the target distribution is the unique stationary distribution of the MCMC and convergence is achieved when samples appear to be stationary. Here we leverage recent advances in computational geometry of the treespace and introduce a method that combines classical statistical techniques and algorithms with geometric properties of the treespace to automatically evaluate and assess practical convergence of phylogenetic MCMC analyses. Our method monitors convergence across multiple MCMC chains and achieves high accuracy in detecting both practical convergence and convergence issues within treespace. Furthermore, our approach is developed to allow for real-time evaluation during the MCMC algorithm run, eliminating any of the chain post-processing steps that are currently required. Our tool therefore improves reliability and efficiency of MCMC based phylogenetic inference methods and makes analyses easier to reproduce and compare. We demonstrate the efficacy of our diagnostic via a well-calibrated simulation study and provide examples of its performance on real data sets. Although our method performs well in practice, a significant part of the underlying treespace probability theory is still missing, which creates an excellent opportunity for future mathematical research in this area.

Who cares? Elucidating parental care evolution in extant birds

Extant birds stand out among vertebrates in the diversity of parental care types they present, spanning absence of care to uniparental care by either sex, biparental care, or even cooperative care. Despite years of research, key questions remain regarding parental care evolution in birds. Firstly, the parental care type in the most recent ancestor of extant birds is a matter of controversy, with proposed ancestral states including no care, uniparental male or female care, and biparental care. Another unsolved question is the direction, order, and frequency of transitions between parental care types. We address these key questions using a database of 5,438 bird species (~50% of extant diversity) and modern phylogenetic comparative methods controlling simultaneously for model and phylogenetic uncertainty as well as potential confounding effects of state-dependent diversification. Our results indicate that the most likely ancestral state for extant birds is male-only care, with a posterior probability of 0.8. Transition rates across parental care types were generally low and heterogeneous; loss of parental care virtually never occurs and transitions away from female-only or cooperative care most often lead to biparental care. Given the low transition rates, future research should analyze the factors favoring the maintenance of care types.

MCMC Methods for Parameter Estimation in ODE Systems for CAR-T Cell Cancer Therapy

Chimeric antigen receptor (CAR)-T cell therapy represents a breakthrough in treating resistant hematologic cancers. It is based on genetically modifying T cells transferred from the patient or a donor. Although its implementation has increased over the last few years, CAR-T has many challenges to be addressed, for instance, the associated severe toxicities, such as cytokine release syndrome. To model CAR-T cell dynamics, focusing on their proliferation and cytotoxic activity, we developed a mathematical framework using ordinary differential equations (ODEs) with Bayesian parameter estimation. Bayesian statistics were used to estimate model parameters through Monte Carlo integration, Bayesian inference, and Markov chain Monte Carlo (MCMC) methods. This paper explores MCMC methods, including the Metropolis-Hastings algorithm and DEMetropolis and DEMetropolisZ algorithms, which integrate differential evolution to enhance convergence rates. The theoretical findings and algorithms were validated using Python and Jupyter Notebooks. A real medical dataset of CAR-T cell therapy was analyzed, employing optimization algorithms to fit the mathematical model to the data, with the PyMC library facilitating Bayesian analysis. The results demonstrated that our model accurately captured the key dynamics of CAR-T cell therapy. This conclusion underscores the potential of parameter estimation to improve the understanding and effectiveness of CAR-T cell therapy in clinical settings.

Practical guidelines for Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC)

Phylogenetic estimation is, and has always been, a complex endeavor. Estimating a phylogenetic tree involves evaluating many possible solutions and possible evolutionary histories that could explain a set of observed data, typically by using a model of evolution. Values for all model parameters need to be evaluated as well. Modern statistical methods involve not just the estimation of a tree, but also solutions to more complex models involving fossil record information and other data sources. Markov chain Monte Carlo (MCMC) is a leading method for approximating the posterior distribution of parameters in a mathematical model. It is deployed in all Bayesian phylogenetic tree estimation software. While many researchers use MCMC in phylogenetic analyses, interpreting results and diagnosing problems with MCMC remain vexing issues to many biologists. In this manuscript, we will offer an overview of how MCMC is used in Bayesian phylogenetic inference, with a particular emphasis on complex hierarchical models, such as the fossilized birth-death (FBD) model. We will discuss strategies to diagnose common MCMC problems and troubleshoot difficult analyses, in particular convergence issues. We will show how the study design, the choice of models and priors, but also technical features of the inference tools themselves can all be adjusted to obtain the best results. Finally, we will also discuss the unique challenges created by the incorporation of fossil information in phylogenetic inference, and present tips to address them.

A macroevolutionary analysis of European Late Upper Palaeolithic stone tool shape using a Bayesian phylodynamic framework

Phylogenetic models are commonly used in palaeobiology to study the patterns and processes of organismal evolution. In the human sciences, phylogenetic methods have been deployed for reconstructing ancestor-descendant relationships using linguistic and material culture data. Within evolutionary archaeology specifically, phylogenetic analyses based on maximum parsimony and discrete traits dominate, which sets limitations for the downstream role cultural phylogenies, once derived, can play in more elaborate analytical pipelines. Recent methodological advances in Bayesian phylogenetics, however, now allow us to infer evolutionary dynamics using continuous characters. Capitalizing on these developments, we here present an exploratory analysis of cultural macroevolution of projectile point shape evolution in the European Final Palaeolithic and earliest Mesolithic (approx. 15 000-11 000 BP) using a Bayesian phylodynamic approach and the fossilized birth-death process model. This model-based approach leaps far beyond the application of parsimony, in that it not only produces a tree, but also divergence times, and diversification rates while incorporating uncertainties. This allows us to compare rates to the pronounced climatic changes that occurred during our time frame. While common in cultural evolutionary analyses of language, the extension of Bayesian phylodynamic models to archaeology arguably represents a major methodological breakthrough.

Comparison of different agents and doses of anti-vascular endothelial growth factors (aflibercept, bevacizumab, conbercept, ranibizumab) versus laser for retinopathy of prematurity: A network meta-analysis

Laser photocoagulation (LPC) and/or intravitreal anti-vascular endothelial growth factor (anti-VEGF) injections constitute the current standard treatment for retinopathy of prematurity (ROP). This network meta-analysis focus on whether a ranking of interventions may be established for different dose levels of intravitreal injection of anti-VEGF agents (aflibercept, bevacizumab, conbercept, ranibizumab) as primary treatments for ROP versus laser in terms of retreatment rate as primary outcome, and time to retreatment and refractive error as secondary endpoints, since best anti-VEGF dosage remains under debate. Sixty-eight studies (15 randomized control trials and 53 nonrandomized studies) of 12,356 eyes of 6445 infants were retrieved from databases (2005 Jan. - 2023 June). Studies were evaluated for model fit, risk of bias and confidence of evidence in Network Meta-Analysis (CINeMA). Bayesian NMA showed that anti-VEGF drugs were not inferior to laser in terms of retreatment rate. For intravitreal bevacizumab (IVB), doses half of the conventional infant dose showed a low risk of retreatment rate (risk ratio (RR) of 1.43; 95% credible interval (CrI): 0.508, 4.03). On probability ranking as surface under the cumulative ranking curve (SUCRA) plot, half dose of bevacizumab had a better position than conventional and augmented (1.2-2 times the regular dose) doses. A similar probability trend was observed for half vs. conventional doses of aflibercept and ranibizumab. Conventional infant dose of conbercept showed the lowest risk for retreatment (RR 0.846; 95% CrI: 0.245, 2.91). For secondary endpoints, lower doses of anti-VEGF agents were associated with shorter times to retreatment. The largest changes were noted for the augmented doses of bevacizumab and ranibizumab (0.3 mg) with means of 14.1 weeks (95% CrI: 6.65, 21.6) and 12.8 weeks (95% CrI: 3.19, 20.9), respectively. Finally, NMA demonstrated better refractive profile for anti-VEGF than laser therapy, especially for the conventional infant doses of bevacizumab and ranibizumab which exhibited a significantly better refractive profile than LPC, with mean differences of 1.67 (spherical equivalent - diopters) (95% CrI: 0.705, 2.67) and 2.19 (95% CrI: 0.782, 3.59), respectively. In the SUCRA plots, LPC had a markedly different position with a higher probability for myopia. Further clinical trials comparing different intravitreal doses of anti-VEGF agents are needed, but our findings suggest that low doses of these drugs retain efficacy and may reduce ocular and systemic undesired events.

How Trustworthy Is Your Tree? Bayesian Phylogenetic Effective Sample Size Through the Lens of Monte Carlo Error

Bayesian inference is a popular and widely-used approach to infer phylogenies (evolutionary trees). However, despite decades of widespread application, it remains difficult to judge how well a given Bayesian Markov chain Monte Carlo (MCMC) run explores the space of phylogenetic trees. In this paper, we investigate the Monte Carlo error of phylogenies, focusing on high-dimensional summaries of the posterior distribution, including variability in estimated edge/branch (known in phylogenetics as "split") probabilities and tree probabilities, and variability in the estimated summary tree. Specifically, we ask if there is any measure of effective sample size (ESS) applicable to phylogenetic trees which is capable of capturing the Monte Carlo error of these three summary measures. We find that there are some ESS measures capable of capturing the error inherent in using MCMC samples to approximate the posterior distributions on phylogenies. We term these tree ESS measures, and identify a set of three which are useful in practice for assessing the Monte Carlo error. Lastly, we present visualization tools that can improve comparisons between multiple independent MCMC runs by accounting for the Monte Carlo error present in each chain. Our results indicate that common post-MCMC workflows are insufficient to capture the inherent Monte Carlo error of the tree, and highlight the need for both within-chain mixing and between-chain convergence assessments.

A new abelisaurid dinosaur from the end Cretaceous of Patagonia and evolutionary rates among the Ceratosauria

Gondwanan dinosaur faunae during the 20 Myr preceding the Cretaceous-Palaeogene (K/Pg) extinction included several lineages that were absent or poorly represented in Laurasian landmasses. Among these, the South American fossil record contains diverse abelisaurids, arguably the most successful groups of carnivorous dinosaurs from Gondwana in the Cretaceous, reaching their highest diversity towards the end of this period. Here we describe Koleken inakayali gen. et sp. n., a new abelisaurid from the La Colonia Formation (Maastrichtian, Upper Cretaceous) of Patagonia. Koleken inakayali is known from several skull bones, an almost complete dorsal series, complete sacrum, several caudal vertebrae, pelvic girdle and almost complete hind limbs. The new abelisaurid shows a unique set of features in the skull and several anatomical differences from Carnotaurus sastrei (the only other abelisaurid known from the La Colonia Formation). Koleken inakayali is retrieved as a brachyrostran abelisaurid, clustered with other South American abelisaurids from the latest Cretaceous (Campanian-Maastrichtian), such as Aucasaurus, Niebla and Carnotaurus. Leveraging our phylogeny estimates, we explore rates of morphological evolution across ceratosaurian lineages, finding them to be particularly high for elaphrosaurine noasaurids and around the base of Abelisauridae, before the Early Cretaceous radiation of the latter clade. The Noasauridae and their sister clade show contrasting patterns of morphological evolution, with noasaurids undergoing an early phase of accelerated evolution of the axial and hind limb skeleton in the Jurassic, and the abelisaurids exhibiting sustained high rates of cranial evolution during the Early Cretaceous. These results provide much needed context for the evolutionary dynamics of ceratosaurian theropods, contributing to broader understanding of macroevolutionary patterns across dinosaurs.

Comparison of Bayesian Coalescent Skyline Plot Models for Inferring Demographic Histories

Bayesian coalescent skyline plot models are widely used to infer demographic histories. The first (non-Bayesian) coalescent skyline plot model assumed a known genealogy as data, while subsequent models and implementations jointly inferred the genealogy and demographic history from sequence data, including heterochronous samples. Overall, there exist multiple different Bayesian coalescent skyline plot models which mainly differ in two key aspects: (i) how changes in population size are modeled through independent or autocorrelated prior distributions, and (ii) how many change-points in the demographic history are used, where they occur and if the number is pre-specified or inferred. The specific impact of each of these choices on the inferred demographic history is not known because of two reasons: first, not all models are implemented in the same software, and second, each model implementation makes specific choices that the biologist cannot influence. To facilitate a detailed evaluation of Bayesian coalescent skyline plot models, we implemented all currently described models in a flexible design into the software RevBayes. Furthermore, we evaluated models and choices on an empirical dataset of horses supplemented by a small simulation study. We find that estimated demographic histories can be grouped broadly into two groups depending on how change-points in the demographic history are specified (either independent of or at coalescent events). Our simulations suggest that models using change-points at coalescent events produce spurious variation near the present, while most models using independent change-points tend to over-smooth the inferred demographic history.

Nucleotide Substitution Model Selection Is Not Necessary for Bayesian Inference of Phylogeny With Well-Behaved Priors

Model selection aims to choose the most adequate model for the statistical analysis at hand. The model must be complex enough to capture the complexity of the data but should be simple enough not to overfit. In phylogenetics, the most common model selection scenario concerns selecting an adequate substitution and partition model for sequence evolution to infer a phylogenetic tree. Previously, several studies showed that substitution model under-parameterization can bias phylogenetic studies. Here, we explored the impact of substitution model over-parameterization in a Bayesian statistical framework. We performed simulations under the simplest substitution model, the Jukes-Cantor model, and compare posterior estimates of phylogenetic tree topologies and tree length under the true model to the most complex model, the $\text{GTR}+\Gamma+\text{I}$ substitution model, including over-splitting the data into additional subsets (i.e., applying partitioned models). We explored 4 choices of prior distributions: the default substitution model priors of MrBayes, BEAST2, and RevBayes and a newly devised prior choice (Tame). Our results show that Bayesian inference of phylogeny is robust to substitution model over-parameterization and over-partitioning but only under our new prior settings. All 3 current default priors introduced biases for the estimated tree length. We conclude that substitution and partition model selection are superfluous steps in Bayesian phylogenetic inference pipelines if well-behaved prior distributions are applied and more effort should focus on more complex and biologically realistic substitution models.

Utility of cytochrome c oxidase I for the deciphering of unstable phylogeny and taxonomy of gorals, genus Nemorhaedus Hamilton Smith, 1827 (Bovidae, Ovibovina)

Gorals represent ungulate mammals of the Palearctic and Indo-Malayan realms that face habitat destruction and intense hunting pressure. Their classification has been the subject of various (mainly genetic) assessments in the last decade, but some results are conflicting, hampering some conservation-based decisions. Genetic sampling of gorals has increased considerably in recent years, at least for mitochondrial (mt) DNA. Results based on two mt genes (cytochrome b and the D-loop) are currently available. Still, the utility of cytochrome oxidase subunit I remains unanalysed, even though it belongs among the gene markers that enable a correct species identification in mammals. This study examines phylogenetic relationships and species delimitation in gorals using all currently available cytochrome oxidase subunit I sequences, including the not yet analysed goral population from Pakistan. Our results of various phylogenetic approaches, such as maximum parsimony, likelihood and Bayesian inference, and exploration of species boundaries via species delimitation support the validity of six species of goral, namely N.baileyi, N.caudatus, N.cranbrooki, N.evansi, N.goral, and N.griseus. This result accords well with results based on other mt genes, especially the cytochrome b from the highly exhaustive data sampling. Our study also summarises common sources of errors in the assessment of goral phylogeny and taxonomy and highlights future priorities in understanding goral diversification.

Fidelity of hyperbolic space for Bayesian phylogenetic inference

Bayesian inference for phylogenetics is a gold standard for computing distributions of phylogenies. However, Bayesian phylogenetics faces the challenging computational problem of moving throughout the high-dimensional space of trees. Fortunately, hyperbolic space offers a low dimensional representation of tree-like data. In this paper, we embed genomic sequences as points in hyperbolic space and perform hyperbolic Markov Chain Monte Carlo for Bayesian inference in this space. The posterior probability of an embedding is computed by decoding a neighbour-joining tree from the embedding locations of the sequences. We empirically demonstrate the fidelity of this method on eight data sets. We systematically investigated the effect of embedding dimension and hyperbolic curvature on the performance in these data sets. The sampled posterior distribution recovers the splits and branch lengths to a high degree over a range of curvatures and dimensions. We systematically investigated the effects of the embedding space's curvature and dimension on the Markov Chain's performance, demonstrating the suitability of hyperbolic space for phylogenetic inference.

Exploring genome gene content and morphological analysis to test recalcitrant nodes in the animal phylogeny

An accurate phylogeny of animals is needed to clarify their evolution, ecology, and impact on shaping the biosphere. Although datasets of several hundred thousand amino acids are nowadays routinely used to test phylogenetic hypotheses, key deep nodes in the metazoan tree remain unresolved: the root of animals, the root of Bilateria, and the monophyly of Deuterostomia. Instead of using the standard approach of amino acid datasets, we performed analyses of newly assembled genome gene content and morphological datasets to investigate these recalcitrant nodes in the phylogeny of animals. We explored extensively the choices for assembling the genome gene content dataset and model choices of morphological analyses. Our results are robust to these choices and provide additional insights into the early evolution of animals, they are consistent with sponges as the sister group of all the other animals, the worm-like bilaterian lineage Xenacoelomorpha as the sister group of the other Bilateria, and tentatively support monophyletic Deuterostomia.

Gilbertaria, a first crustose genus in the Sphaerophoraceae (Lecanoromycetes, Ascomycota) for Catillaria contristans, Toninia squalescens and related species

Lecideoid lichen-forming fungi are a large, heterogeneous group that includes many species described during the nineteenth century that are of unclear taxonomic status. We revise such a group, the species of which have previously been treated under the much-misunderstood names Catillaria contristans or Toninia squalescens, and use a seven-locus phylogeny to determine its phylogenetic position. We found strong support for a previously unrecognized monophyletic lineage within the Sphaerophoraceae, comprising five phylogenetic species, and describe the new genus Gilbertaria to accommodate them. The new genus is characterized by a crustose growth form, 1-septate ascospores, thick ((1.5–)2–3(–4) μm wide) paraphyses and asci of the Biatora-type. We revise the nomenclature and give new delimitations and descriptions of the Northern Hemisphere species Gilbertaria contristans comb. nov., G. holomeloides comb. nov., G. squalescens comb. nov. and describe the new species G. astrapeana from the Falkland Islands.

Inferring the demographic history of the North American firefly Photinus pyralis

The firefly Photinus pyralis inhabits a wide range of latitudinal and ecological niches, with populations living from temperate to tropical habitats. Despite its broad distribution, its demographic history is unknown. In this study, we modelled and inferred different demographic scenarios for North American populations of P. pyralis, which were collected from Texas to New Jersey. We used a combination of ABC techniques (for multi-population/colonization analyses) and likelihood inference (dadi, StairwayPlot2, PoMo) for single-population demographic inference, which proved useful with our RAD data. We uncovered that the most ancestral North American population lays in Texas, which further colonized the Central region of the US and more recently the North Eastern coast. Our study confidently rejects a demographic scenario where the North Eastern populations colonized more southern populations until reaching Texas. To estimate the age of divergence between of P. pyralis, which provides deeper insights into the history of the entire species, we assembled a multi-locus phylogenetic data covering the genus Photinus. We uncovered that the phylogenetic node leading to P. pyralis lies at the end of the Miocene. Importantly, modelling the demographic history of North American P. pyralis serves as a null model of nucleotide diversity patterns in a widespread native insect species, which will serve in future studies for the detection of adaptation events in this firefly species, as well as a comparison for future studies of other North American insect taxa.

Variation in eye abundance among scallops reveals ontogenetic and evolutionary convergence associated with life habits

Eyes are remarkable systems to investigate the complex interaction between ecological drivers and phenotypic outcomes. Some animals, such as scallops, have many eyes for visual perception, but to date, the evolution of multiple-eye systems remains obscure. For instance, it is unclear whether eye number changes over a lifetime or varies among species. Scallops are a suitable model group to investigate these questions considering the interspecific variation of adult size and ecological diversity. We tested whether eye abundance scales with body size among individuals and species and whether it varies with life habits. We performed comparative analyses, including a phylogenetic ANCOVA and evolutionary model comparisons, based on eye count and shell height (as a proxy of body size) across 31 scallop species. Our analyses reveal that patterns of increasing relationship with body size are not concordant among taxa and suggest ontogenetic convergence caused by similar ecologies. Accordingly, selective optima in eye numbers are associated with shifts in life habits. For instance, species with increased mobility have significantly more eyes than less mobile species. The convergent evolution of greater eye abundance in more mobile scallops likely indicates a visual improvement based on increased levels of oversampling of the surrounding environment.