Phylogenetic Methods Research Articles

Phylogenetic tree instability after taxon addition: empirical frequency, predictability, and consequences for online inference.

Online phylogenetic inference methods add sequentially arriving sequences to an inferred phylogeny without the need to recompute the entire tree from scratch. Some online method implementations exist already, but there remains concern that additional sequences may change the topological relationship among the original set of taxa. We call such a change in tree topology a lack of stability for the inferred tree. In this paper, we analyze the stability of single taxon addition in a Maximum Likelihood framework across 1, 000 empirical datasets. We find that instability occurs in almost 90% of our examples, although observed topological differences do not always reach significance under the AU-test. Changes in tree topology after addition of a taxon rarely occur close to its attachment location, and are more frequently observed in more distant tree locations carrying low bootstrap support. To investigate whether instability is predictable, we hypothesize sources of instability and design summary statistics addressing these hypotheses. Using these summary statistics as input features for machine learning under random forests, we are able to predict instability and can identify the most influential features. In summary, it does not appear that a strict insertion-only online inference method will deliver globally optimal trees, although relaxing insertion strictness by allowing for a small number of final tree rearrangements or accepting slightly suboptimal solutions appears feasible.

Arrangements of intramembranous muscles of wings are influenced by body mass across the radiation of phyllostomid bats.

Extensive research into bat flight mechanisms has highlighted the complex functional and evolutionary dynamics of their wing structures, yet the anatomical details of certain wing muscles remain elusive. In particular, the intramembranous plagiopatagiales proprii muscles, located within the plagiopatagium-an area of the wing lacking direct joint connections-exhibit remarkable variation across bat families. These muscles, which extend anteroposteriorly in macroscopic bundles, play a crucial role in wing stiffening, modulating membrane tension, and reducing wing curvature during flight. Since larger bats tend to have higher wing loading (WL; the ratio of body mass [BMa] to wing area) and may therefore experience increased patagial curvature and resultant drag, we hypothesized that body size significantly influences the evolutionary development of the plagiopatagiales proprii muscles. This study investigates the relationship between BMa and the morphology of the plagiopatagiales proprii in New World leaf-nosed bats (Phyllostomidae), employing bivariate allometry, multivariate analysis, and comparative phylogenetic methods across 24 species from eight phyllostomid subfamilies. Our findings reveal a significant phylogenetic signal in muscle architecture, along with positive evolutionary allometry in muscle area. This suggests an adaptive increase in muscle size in larger species, likely to counterbalance the increased WL, reduce wing curvature, and minimize drag. This research enhances our understanding of the functional and adaptive morphological evolution of intramembranous wing muscles in phyllostomid bats, underscoring their evolutionary significance.

The female and emended male description of the USA Interior Highlands endemic Allocapnia oribata Poulton & Stewart, 987 (Plecoptera: Capniidae)

The female of Allocapnia oribata Poulton & Stewart, 1987 is newly associated by tree- and genetic distance-based phylogenetic methods using the mitochondrial cytochrome c oxidase subunit I gene barcode fragment sequence data. The first description of the female and an emended description of the male of A. oribata are presented. Supporting data includes color images, scanning electron micrographs, and genetic analysis of DNA barcodes. Corrections to the holotype and paratype localities are presented with the help of B.C. Poulton. Additional notes on the phylogeny and endemism of Allocapnia Claassen, 1928 are provided.

Whole Genome Analysis Reveals Evolutionary History and Introgression Events in Bale Monkeys

Background/Objective: The Bale monkey (Chlorocebus djamdjamensis) is a threatened primate species endemic to Ethiopia and, in contrast to other members of the genus Chlorocebus, lives at high altitudes and feeds mainly on bamboo. Two populations of the species are present, one in continuous bamboo forest (CF) in the eastern part of the species’ range, and the other in fragmented forest (FF) in the western part. Based on mitochondrial DNA and phenotypic characteristics, previous studies have suggested introgression by parapatric congeners into the FF population but not into the CF population. The objective of this study was to gain insights into the evolutionary history of Bale monkeys and their potential genetic adaptations to high altitudes and for bamboo consumption. Methods: We sequenced the whole genomes of individuals from both populations and compared their genomes with those of the other five Chlorocebus species. We applied phylogenetic methods and conducted population demographic simulations to elucidate their evolutionary history. A genome-wide analysis was conducted to assess gene flow and identify mutations potentially associated with adaptations to high altitudes and for bamboo metabolism. Results: Our analyses revealed Bale monkeys as the sister clade to Chlorocebus aethiops and showed that gene flow occurred between C. aethiops and FF but not between C. aethiops and CF. In addition, we detected non-synonymous mutations in genes potentially associated with the adaptation to high altitudes (EPAS1) in both populations and with the adaptation for bamboo metabolism (TAS2R16, MPST, and TST) mainly in the CF population. Conclusions: Our study provides insights into the evolutionary history of a threatened primate species and reveals the genetic basis for its adaptions to unique environments and for diet specialization.

The meaning and measure of concordance factors in phylogenomics.

As phylogenomic datasets have grown in size, researchers have developed new ways to measure biological variation and to assess statistical support for specific branches. Larger datasets have more sites and loci, and therefore less sampling variance. While we can more accurately measure the mean signal in these datasets, lower sampling variance is often reflected in uniformly high measures of branch support-such as the bootstrap and posterior probability-limiting their utility. Larger datasets have also revealed substantial biological variation in the topologies found across individual loci, such that the single species tree inferred by most phylogenetic methods represents a limited summary of the data for many purposes. In contrast to measures of statistical support, the degree of underlying topological variation among loci should be approximately constant regardless of the size of the dataset. "Concordance factors" and similar statistics have therefore become increasingly important tools in phylogenetics. In this review, we explain why concordance factors should be thought of as descriptors of topological variation rather than as measures of statistical support, and argue that they provide important information about the predictive power of the species tree not contained in measures of support. We review a growing suite of statistics for measuring concordance, compare them in a common framework that reveals their interrelationships, and demonstrate how to calculate them using an example from birds. We also discuss how measures of topological variation might change in the future as we move beyond estimating a single "tree of life" towards estimating the myriad evolutionary histories underlying genomic variation.

Biodiversity of symbiotic gastropods (Mollusca, Eulimidae) in Vietnamese waters barcoded

ABSTRACT Eulimidae is a prominent group within the class Gastropoda, comprising 1087 species that have developed specialised symbiotic relationships with all extant echinoderm classes. At least 47 species of eulimids can be found in the waters of Vietnam. Representatives of the genus Pulicicochlea, the species Vitreobalcis tripneusticola, and four new undescribed species were observed for the first time. In this study, we update the existing checklists of Vietnamese eulimids, emphasising the use of molecular barcoding to determine the phylogenetic relationships among the collected species. The classical taxonomy of molluscs relies heavily on the morphology of the shell and radula. The family Eulimidae presents challenges as in many species these features are poorly developed, and molecular phylogenetic methods have proven to be more informative for the taxonomy and systematics of this family. To understand the evolutionary pathways in Eulimidae, we constructed a COI-based phylogenetic tree, which included original sequences primarily from Vietnam and all available sequences for Eulimidae from GenBank. The newly obtained genetic information is crucial for enhancing our understanding of this family of gastropod mollusсs, as well as for any future revisions that may be necessary.

Color polymorphic carnivores have faster speciation rates

Variation in coat color is a prominent feature in carnivores, thought to be shaped by environmental factors. As new traits could allow populations to occupy novel niches and habitats, color polymorphism may be maintained by balancing selection. Consequently, color polymorphic species may speciate more rapidly and can give rise to monomorphic daughter species. We thus predicted that, within the Carnivora, (i) speciation rate is higher in polymorphic lineages, (ii) divergence between color polymorphic lineages is more recent, and (iii) within closely related groups, polymorphic lineages are ancestral and monomorphic lineages derived. We also tested whether accelerated speciation rates relate to niche breadth, measured by the number of occupied habitats and range size. We collected data of 48 polymorphic and 192 monomorphic carnivore species, and assessed speciation rates using phylogenetic comparative methods. We found that polymorphic carnivores had higher speciation rates (λ1 = 0.29, SD = 0.13) than monomorphic species (λ0 = 0.053, SD = 0.044). Hidden and quantitative state speciation and extinction models inferred that color polymorphism was the main contributing factor, and that niche breadth was not of influence. Therefore, other selective forces than spatial niche segregation, such as predator-prey coevolution, may contribute to color polymorphism in wild carnivores.

Comparative Population Genomic Diversity and Differentiation in Trapdoor Spiders and Relatives (Araneae, Mygalomorphae).

Although patterns of population genomic variation are well-studied in animals, there remains room for studies that focus on non-model taxa with unique biologies. Here we characterise and attempt to explain such patterns in mygalomorph spiders, which are generally sedentary, often occur as spatially clustered demes and show remarkable longevity. Genome-wide single nucleotide polymorphism (SNP) data were collected for 500 individuals across a phylogenetically representative sample of taxa. We inferred genetic populations within focal taxa using a phylogenetically informed clustering approach, and characterised patterns of diversity and differentiation within- and among these genetic populations, respectively. Using phylogenetic comparative methods we asked whether geographical range sizes and ecomorphological variables (behavioural niche and body size) significantly explain patterns of diversity and differentiation. Specifically, we predicted higher genetic diversity in genetic populations with larger geographical ranges, and in small-bodied taxa. We also predicted greater genetic differentiation in small-bodied taxa, and in burrowing taxa. We recovered several significant predictors of genetic diversity, but not genetic differentiation. However, we found generally high differentiation across genetic populations for all focal taxa, and a consistent signal for isolation-by-distance irrespective of behavioural niche or body size. We hypothesise that high population genetic structuring, likely reflecting combined dispersal limitation and microhabitat specificity, is a shared trait for all mygalomorphs. Few studies have found ubiquitous genetic structuring for an entire ancient and species-rich animal clade.

Temporal lobe evolution in Hominidae and the origin of human lobe proportions.

Objectives Evolutionary changes in hominin social complexity have been associated with increases in absolute brain size. The temporal lobes are nestled in the middle cranial fossae (MCF) of the skull, the dimensions of which allow estimation of temporal lobe volume (TLV) in extant and fossil taxa. Materials and Methods The main aim of this study is to determine where along the hominid phylogeny, major temporal lobe size transitions occurred. We used computed tomography (CT) scans of crania, 3D photogrammetry data, and laser surface scans of endocranial casts to measure seven MCF metrics in 11 extant anthropoid taxa using multiple regressions to estimate TLV in 5 extant hominids and 10 fossil hominins. Phylogenetic comparative methods mapped temporal lobe size, brain size, and temporal lobe proportions onto phylogenetic trees broadly for Hominidae and specifically for Hominini. Results Extant Homo sapiens were not an outlier in relative brain size, temporal lobe size, or proportions of the temporal lobes, but some proportions within the lobe were uniquely altered. The most notable changes in relative temporal lobe size and proportions saw a decrease in relative temporal lobe size and proportions in the genus Pan compared to other extant great apes and fossil hominins while there was a relative increase in the temporal lobe width and length in Australopithecus-Paranthropus clade compared to the genus Homo and other extant great apes including modern humans. Discussion We do not find support for the social brain, environmental or functional craniology hypotheses alone but think it prudent to consider the implications of cerebral reorganization between the temporal lobes and other regions of the brain within the context of these hypotheses and with future investigation is warranted.

Construction of core genome multi-locus sequence typing schemes for population structure analyses of Nocardia species

Nocardia, a member of the Actinobacteria phylum, populates diverse habitats globally, with certain species being the cause of various clinical infections in humans. There is paucity of data regarding the population structure of this genus and of established genomic-based phylogenetic methods. We examined the whole genome sequences of 193 isolates spanning five major pathogenic Nocardia species sourced from public databases, encompassing diverse geographic regions. Using the chewBBACA pipeline, a species-specific core genome multilocus sequence typing (cgMLST) schema was created for N. cyriacigeorgica, N. farcinica, N. brasiliensis, N. wallacei, and N. abscessus. Additional genomic features that were examined included virulence factor (VF) profile, total length and open-reading frame count, the core genome length and core gene count, and GC content. Our findings indicated that: (i) N. brasiliensis diverges significantly from the other four species, underscoring its distinct evolutionary trajectory; (ii) the population structures of all species were polyclonal, with phylogenetic clustering occurring in the minority of isolates; (iii) clonal complexes were largely restricted to specific geographical locations, rather than demonstrating a global distribution, and (iv) initial evidence suggests no direct common-source transmission amongst the studied strains. Our study establishes a comprehensive genome-based phylogenetic methodology for population structure of Nocardia species.

Phylogenomic and synteny analysis of BAHD and SCP/SCPL gene families reveal their evolutionary histories in plant specialized metabolism.

Plant chemical diversity is largely owing to a number of enzymes which catalyse reactions involved in the assembly, and in the subsequent chemical modifications, of the core structures of major classes of plant specialized metabolites. One such reaction is acylation. With this in mind, to study the deep evolutionary history of BAHD and the serine-carboxypeptidase-like (SCPL) acyltransferase genes, we assembled phylogenomic synteny networks based on a large-scale inference analysis of orthologues across whole-genome sequences of 126 species spanning Stramenopiles and Archaeplastida, including Arabidopsis thaliana, tomato (Solanum lycopersicum) and maize (Zea mays). As such, this study combined the study of genomic location with changes in gene sequences. Our analyses revealed that serine-carboxypeptidase (SCP)/serine-carboxypeptidase-like (SCPL) genes had a deeper evolutionary origin than BAHD genes, which expanded massively on the transition to land and with the development of the vascular system. The two gene families additionally display quite distinct patterns of copy number variation across phylogenies as well as differences in cross-phylogenetic syntenic network components. In unlocking the above observations, our analyses demonstrate the possibilities afforded by modern phylogenomic (syntenic) networks, but also highlight their current limitations, as demonstrated by the inability of phylogenetic methods to separate authentic SCPL acyltransferases from standard SCP peptide hydrolases.This article is part of the theme issue 'The evolution of plant metabolism'.

Evolutionary Insights from Association Rule Mining of Co-Occurring Mutations in Influenza Hemagglutinin and Neuraminidase.

Seasonal influenza viruses continuously evolve via antigenic drift. This leads to recurring epidemics, globally significant mortality rates, and the need for annually updated vaccines. Co-occurring mutations in hemagglutinin (HA) and neuraminidase (NA) are suggested to have synergistic interactions where mutations can increase the chances of immune escape and viral fitness. Association rule mining was used to identify temporal relationships of co-occurring HA-NA mutations of influenza virus A/H3N2 and its role in antigenic evolution. A total of 64 clusters were found. These included well-known mutations responsible for antigenic drift, as well as previously undiscovered groups. A majority (41/64) were associated with known antigenic sites, and 38/64 involved mutations across both HA and NA. The emergence and disappearance of N-glycosylation sites in the pattern of N-X-[S/T] were also identified, which are crucial post-translational processes to maintain protein stability and functional balance (e.g., emergence of NA:339ASP and disappearance of HA:187ASP). Our study offers an alternative approach to the existing mutual-information and phylogenetic methods used to identify co-occurring mutations, enabling faster processing of large amounts of data. Our approach can facilitate the prediction of critical mutations given their occurrence in a previous season, facilitating vaccine development for the next flu season and leading to better preparation for future pandemics.

The evolutionary dynamics of grammatical gender in Torricelli languages

Abstract Grammatical gender in New Guinea is an often neglected area in typological research, even though it is extremely diverse. For example, in New Guinea, some languages have grammatical gender systems with two sex-based categories, more than four gender-indexing targets, and no gender marking on nouns, while some languages have grammatical gender systems with much more categories, which are only marginally sex-based. This paper infers the processes of development and change of grammatical gender in Torricelli languages from two perspectives. First, it synthesizes the available data in the existing literature and hypothesizes the evolutionary pathway of gender systems in Torricelli languages. Nineteen Torricelli languages are selected as a representative coverage of the 55 Torricelli languages listed in Glottolog within the limits of the available documentation. These languages are then coded based on 6 presence-absence features relating to gender marking on verbs, adjectives, nouns, numerals, pronouns, and demonstratives. Second, it conducts an analysis with phylogenetic comparative methods to provide a quantitative assessment of the evolutionary possibilities for gender systems in Torricelli languages. The preliminary results show that gender is likely marked at the root of Torricelli languages, with pronouns and verbs being at the core of the system. This is in agreement with trends reflecting the evolution of gender systems in languages across the world.

Multidrug-resistant tuberculosis clusters and transmission in Taiwan: a population-based cohort study.

Multidrug-resistant tuberculosis (MDR-TB) remains a challenge in the TB program of Taiwan, where 0.5% of new cases and 2.1% of previously treated cases were resistant to at least rifampin (RIF) and isoniazid (INH). Since >80% of our MDR-TB are new cases, genotyping of MDR Mycobacterium tuberculosis is implemented to facilitate contact investigation, cluster identification, and outbreak delineation. This is a population-based retrospective cohort study analyzing MDR-TB cases from 2019 to 2022. Whole genome sequencing (WGS) was performed using the Illumina MiSeq and analyzed using the TB Profiler. A single nucleotide polymorphism (SNP) threshold of ≤ 12 and phylogenetic methods were used to identify putative transmission clusters. An outbreak was confirmed using genomic data and epidemiologic links. Of the 297 MDR-TB cases, 246 (82.8%), 45 (15.2%), and 6 (2.0%) were simple MDR, extensively drug-resistant tuberculosis (pre-XDR-TB) and extensively drug-resistant tuberculosis (XDR-TB), respectively. The sublineage 2.2 modern Beijing was the predominant (48.8%) MDR-TB strain in Taiwan. Phylogenetic analysis identified 25.3% isolates in 20 clusters, with cluster sizes ranging from 2 to 13 isolates. Nevertheless, only 2 clusters, one household and one community, were confirmed as outbreaks. In this study, we found that males had a higher risk of MDR-TB transmission compared to females, and those infected with the sublineage 2.1-proto-Beijing genotype isolates were at a higher risk of transmission. Furthermore, 161 (54.2%) isolates harbored compensatory mutations in the rpoC and non-rifampicin resistant determinant region (non-RRDR) of the rpoB gene. MDR-TB strains containing rpoB S450L and other compensatory mutations concurrently were significantly associated with clusters, especially the proto-Beijing genotype strains with the compensatory mutation rpoC E750D or the modern Beijing genotype strains with rpoC D485Y/rpoC E1140D. Routine and continuous surveillance using WGS-based analysis is recommended to warn of risks and delineate transmission clusters of MDR-TB. We proposed the use of compensatory mutations as epidemiological markers of M. tuberculosis to interrupt putative MDR-TB transmission.

Molecular cluster, transmission characteristics, origin and dynamics analysis of HIV-1 CRF59_01B in China: A molecular epidemiology study

PurposeThis study investigated for the HIV-1 CRF59_01B epidemic's spatiotemporal dynamics and its transmission networks in China. MethodsBetween 2007 and 2020, a total of 250 partial pol gene sequences of HIV-1 CRF59_01B were collected from four regions (10 Chinese provinces). Phylogenetic tree construction and cluster identification were then performed. The Bayesian skyline and birth–death susceptible–infected–removed models were employed for the phylodynamic analyses of subtypes and large clusters, respectively. Phylogenetic analyses and trait diffusion of these sequences were performed using Bayesian phylogenetic methods (beast-classic package). Distance-based molecular network analyses were performed to identify putative relationships. ResultsUsing a genetic distance threshold of 1.3 %, We identified 45 clusters that included 62.40 % (156/250) of the sequences. Three clusters (6.67 %, 3/45) had 10 or more sequences, and were considered "large clusters". Six clusters (13.33 %) included sequences from different regions (Southeast, Northeast, Southeast, and Central China). Thirteen clusters (28.89 %) included sequences of men who had sex with men only, three clusters (6.67 %) included sequences of heterosexuals only, and 12 clusters (26.67 %) included sequences of both groups. The substitution rate of CRF59_01B was 1.91 × 10–3 substitutions per site per year [95 % highest posterior density (HPD) interval: 1.39 × 10−3–2.49 × 10−3)], the time to the most recent common ancestor of CRF59_01B was to be 1992.83 (95 % HPD: 1977.97–2002.81). A Bayesian skyline plot revealed that the effective population size of CRF59_01B increased from 2000 to 2015 and remained stable after 2015. The large clusters showed continuous growth from 2013 to 2020. Phylogeographic analysis showed that CRF59_01B B most likely originated in Southeast China, with a posterior probability of 97.44 %, and then spread to other regions. ConclusionsOur study revealed the temporal and geographical origins of HIV-1 CRF59_01B as well as the process of transmission among various regions and risk groups in China, which can help develop targeted HIV prevention strategies.

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. The open source package for the phylogenetic inference framework BEAST2, called ASM, that implements these methods, making them accessible through a user-friendly GUI, is available from https://github.com/rbouckaert/asm/. The open source Python package, called tetres, that provides an interface for these methods enabling their applications beyond BEAST2 can be accessed at https://github.com/bioDS/tetres/.

Recurrent evolution of small body size and loss of the sword ornament in Northern swordtail fish.

Across the tree of life, species have repeatedly evolved similar phenotypes. While well-studied for ecological traits, there is also evidence for recurrent evolution of sexually selected traits. Swordtail fish (Xiphophorus) are a classic model system for studying sexual selection, and female Xiphophorus exhibit strong mate preferences for large male body size and a range of sexually dimorphic ornaments. Interestingly, sexually selected traits have also been lost multiple times in the genus. However, there has been uncertainty over the number of losses of ornamentation and large body size because phylogenetic relationships between species in this group have historically been controversial, partially due to prevalent gene flow. Here, we use whole-genome sequencing approaches to re-examine phylogenetic relationships within a Xiphophorus clade that varies in the presence and absence of sexually selected traits. Using wild-caught individuals, we determine the phylogenetic placement of a small, unornamented species, X. continens, confirming an additional loss of ornamentation and large body size in the clade. With these revised phylogenetic relationships, we analyze evidence for coevolution between body size and other sexually selected traits using phylogenetic comparative methods. These results provide insights into the evolutionary pressures driving the recurrent loss of suites of sexually selected traits.

Compensating for the corolla? Pollen exposure is not associated with pollen collecting hair length.

Secondary pollen-presentation, the relocation of pollen from the anthers to elsewhere on the flower, has evolved multiple times across many plant families. While hypotheses suggest it evolved to promote outcrossing, a byproduct of relocation may be protection of pollen from loss due to abiotic factors. In Campanulaceae pollen is presented on pollen-collecting hairs along the style and the hairs retract over time and release pollen for transfer. Campanulaceae taxa vary in the degree to which pollen is exposed to environmental factors due to variation in the corolla shape and size. We tested the protective function of pollen-collecting hairs by assessing whether there was a tradeoff between the protection provided by the corolla and the pollen-collecting hairs. We used phylogenetic comparative methods to test for associations between pollen-collecting hair length, floral shape and size metrics, and pollen exposure traits across 39 species. We anticipated longer pollen-collecting hairs in taxa with more exposed pollen presentation but found there was no relationship between estimates of pollen exposure and pollen-collecting hair length. However, pollen-collecting hair length scaled allometrically with floral size, and variation in pollen-collecting hairs, as well as most floral traits, was phylogenetically structured. These results indicate that variation in pollen exposure across species does not structure variation in the pollen-collecting hairs, rather hair length scales allometrically and is phylogenetically constrained, therefore pollen-collecting hairs are unlikely to facilitate protection from environmental pressures.

Genomic transfers help to decipher the ancient evolution of filoviruses and interactions with vertebrate hosts.

Although several filoviruses are dangerous human pathogens, there is conflicting evidence regarding their origins and interactions with animal hosts. Here we attempt to improve this understanding using the paleoviral record over a geological time scale, protein structure predictions, tests for evolutionary maintenance, and phylogenetic methods that alleviate sources of bias and error. We found evidence for long branch attraction bias in the L gene tree for filoviruses, and that using codon-specific models and protein structural comparisons of paleoviruses ameliorated conflict and bias. We found evidence for four ancient filoviral groups, each with extant viruses and paleoviruses with open reading frames. Furthermore, we found evidence of repeated transfers of filovirus-like elements to mouse-like rodents. A filovirus-like nucleoprotein ortholog with an open reading frame was detected in three subfamilies of spalacid rodents (present since the Miocene). We provide evidence that purifying selection is acting to maintain amino acids, protein structure and open reading frames in these elements. Our finding of extant viruses nested within phylogenetic clades of paleoviruses informs virus discovery methods and reveals the existence of Lazarus taxa among RNA viruses. Our results resolve a deep conflict in the evolutionary framework for filoviruses and reveal that genomic transfers to vertebrate hosts with potentially functional co-options have been more widespread than previously appreciated.

A Phylogenetic Model of Established and Enabled Biome Shifts.

Where each species actually lives is distinct from where it could potentially survive and persist. This suggests that it may be important to distinguish established from enabled biome affinities when considering how ancestral species moved and evolved among major habitat types. We introduce a new phylogenetic method, called RFBS, to model how anagenetic and cladogenetic events cause established and enabled biome affinities (or, more generally, other discrete realized versus fundamental niche states) to shift over evolutionary timescale. We provide practical guidelines for how to assign established and enabled biome affinity states to extant taxa, using the flowering plant clade Viburnum as a case study. Through a battery of simulation experiments, we show that RFBS performs well, even when we have realistically imperfect knowledge of enabled biome affinities for most analyzed species. We also show that RFBS reliably discerns established from enabled affinities, with similar accuracy to standard competing models that ignore the existence of enabled biome affinities. Lastly, we apply RFBS to Viburnum to infer ancestral biomes throughout the tree and to highlight instances where repeated shifts between established affinities for warm and cold temperate forest biomes were enabled by a stable and slowly-evolving enabled affinity for both temperate biomes.