Identifying along which lineages shifts in diversification rates occur is a central goal of comparative phylogenetics; these shifts may coincide with key evolutionary events such as the development of novel morphological characters, the acquisition of adaptive traits, polyploidization or other structural genomic changes, or dispersal to a new habitat and subsequent increase in environmental niche space. However, while multiple methods now exist to estimate diversification rates and identify shifts using phylogenetic topologies, the appropriate use and accuracy of these methods are hotly debated. Here we test whether five Bayesian methods-Bayesian Analysis of Macroevolutionary Mixtures (BAMM), two implementations of the Lineage-Specific Birth-Death-Shift model (LSBDS and PESTO), the approximate Multi-Type Birth-Death model (MTBD; implemented in BEAST2), and the Cladogenetic Diversification Rate Shift model (ClaDS2)-produce comparable results. We apply each of these methods to a set of 65 empirical time-calibrated phylogenies and compare inferences of speciation rate, extinction rate, and net diversification rate. We find that the five methods often infer different speciation, extinction, and net-diversification rates. Consequently, these different estimates may lead to different interpretations of the macroevolutionary dynamics. The different estimates can be attributed to fundamental differences among the compared models. Therefore, the inference of shifts in diversification rates is strongly method dependent. We advise biologists to apply multiple methods to test the robustness of the conclusions or to carefully select the method based on the validity of the underlying model assumptions to their particular empirical system.

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species\' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny [30]. Here, we introduce a new family of models - the so-called \"Phylogenetic Integrated Velocity\" (PIV) models - that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the U.S.A. indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Infectious bronchitis virus (IBV) is prevalent in Pakistan causing enormous economic losses. To date no clear data are available on circulating genotypes and phylogeographic spread of the virus. Hence current study assessed these parameters for all available IBV Pakistani isolates, based on the 9 new sequences, with respect to other Asian and non-Asian countries. Results indicated that all Pakistani isolates belonged to genotype I (GI), with more than half of them (16/27) belonging to the GI-24 lineage, against which no vaccine is available. Three possible introduction events of the GI-13 IBV lineage into Pakistan, based on the estimated IBV population using isolates from this study, were observed possibly from Afghanistan, China, and/or Egypt. These events were further analyzed on the S1 amino acid level which showed unique alterations (S250H, T270K, and Q298S) in 1 isolate (IBV4, GI-13) when compared to GI-1 lineage. Both GI-1 and GI-13 Pakistani strains showed close homology with homologous vaccine strains that are used in Pakistan. For GI-24 strains, none of the used vaccines showed substantial homology, necessitating the need for further exploration of this lineage and vaccine design. In addition, our findings highlight the importance of genomic surveillance to support phylogeographical studies on IBV in genotyping and molecular epidemiology.

Jingmen virus (JMV) associated with ticks and vertebrates have been found to be related to human disease. We obtained the genome of a Jingmen tick virus (JMTV) strain from Rhipicephalus microplus in Guizhou province and compared the genomes of seven JMV species associated with ticks and vertebrates to understand the evolutionary relationships. The topology of the phylogenetic tree of segment 1 and segment 3 is similar, and segment 2 and segment 4 formed two different topologies, with the main differences being between Alongshan virus (ALSV), Takachi virus, Yanggou tick virus and Pteropus lylei jingmen virus (PLJV), and the possibility of genetic reassortment among these viruses. Moreover, we detected recombination within JMTV and between PLJV and ALSV. The genetic reassortment and recombination that occurs during cross-species transmission of these JMV associated with ticks and vertebrates not only complicates their evolutionary relationships, but also raises the risk of these viruses to humans.

Several hantaviruses result in zoonotic infections of significant public health concern, causing hemorrhagic fever with renal syndrome (HFRS) or hantavirus cardiopulmonary syndrome (HCPS) in the Old and New World, respectively. Given a 35% case fatality rate, disease-causing New World hantaviruses require a greater understanding of their biology, genetic diversity, and geographical distribution. Juquitiba hantaviruses have been identified in Oligoryzomys nigripes in Brazil, Paraguay, and Uruguay. Brazil has reported the most HCPS cases associated with this virus. We used a multiplexed, amplicon-based PCR strategy to screen and deep-sequence the virus harbored within lung tissues collected from Oligoryzomys species during rodent field collections in southern (Itapúa) and western (Boquerón) Paraguay. No Juquitiba-like hantaviruses were identified in Boquerón. Herein, we report the full-length S and M segments of the Juquitiba hantaviruses identified in Paraguay from O. nigripes. We also report the phylogenetic relationships of the Juquitiba hantaviruses in rodents collected from Itapúa with those previously collected in Canindeyú. We showed, using the TN93 nucleotide substitution model, the coalescent (constant-size) population tree model, and Bayesian inference implemented in the Bayesian evolutionary analysis by sampling trees (BEAST) framework, that the Juquitiba virus lineage in Itapúa is distinct from that in Canindeyú. Our spatiotemporal analysis showed significantly different time to the most recent ancestor (TMRA) estimates between the M and S segments, but a common geographic origin. Our estimates suggest the additional geographic diversity of the Juquitiba virus within the Interior Atlantic Forest and highlight the need for more extensive sampling across this biome.

H5Nx highly pathogenic avian influenza (HPAI) viruses of clade 2.3.4.4 have caused outbreaks in Europe among wild and domestic birds since 2016 and were introduced to North America via wild migratory birds in December 2021. We examined the spatiotemporal extent of HPAI viruses across continents and characterized ecological and environmental predictors of virus spread between geographic regions by constructing a Bayesian phylodynamic generalized linear model (phylodynamic-GLM). The findings demonstrate localized epidemics of H5Nx throughout Europe in the first several years of the epizootic, followed by a singular branching point where H5N1 viruses were introduced to North America, likely via stopover locations throughout the North Atlantic. Once in the United States (US), H5Nx viruses spread at a greater rate between US-based regions as compared to prior spread in Europe. We established that geographic proximity is a predictor of virus spread between regions, implying that intercontinental transport across the Atlantic Ocean is relatively rare. An increase in mean ambient temperature over time was predictive of reduced H5Nx virus spread, which may reflect the effect of climate change on declines in host species abundance, decreased persistence of the virus in the environment, or changes in migratory patterns due to ecological alterations. Our data provide new knowledge about the spread and directionality of H5Nx virus dispersal in Europe and the US during an actively evolving intercontinental outbreak, including predictors of virus movement between regions, which will contribute to surveillance and mitigation strategies as the outbreak unfolds, and in future instances of uncontained avian spread of HPAI viruses.

Biological phenotypes are products of complex evolutionary processes in which selective forces influence multiple biological trait measurements in unknown ways. Phylogenetic comparative methods seek to disentangle these relationships across the evolutionary history of a group of organisms. Unfortunately, most existing methods fail to accommodate high-dimensional data with dozens or even thousands of observations per taxon. Phylogenetic factor analysis offers a solution to the challenge of dimensionality. However, scientists seeking to employ this modeling framework confront numerous modeling and implementation decisions, the details of which pose computational and replicability challenges.We develop new inference techniques that increase both the computational efficiency and modeling flexibility of phylogenetic factor analysis. To facilitate adoption of these new methods, we present a practical analysis plan that guides researchers through the web of complex modeling decisions. We codify this analysis plan in an automated pipeline that distills the potentially overwhelming array of decisions into a small handful of (typically binary) choices.We demonstrate the utility of these methods and analysis plan in four real-world problems of varying scales. Specifically, we study floral phenotype and pollination in columbines, domestication in industrial yeast, life history in mammals, and brain morphology in New World monkeys.General and impactful community employment of these methods requires a data scientific analysis plan that balances flexibility, speed and ease of use, while minimizing model and algorithm tuning. Even in the presence of non-trivial phylogenetic model constraints, we show that one may analytically address latent factor uncertainty in a way that (a) aids model flexibility, (b) accelerates computation (by as much as 500-fold) and (c) decreases required tuning. These efforts coalesce to create an accessible Bayesian approach to high-dimensional phylogenetic comparative methods on large trees.

Botryosphaeriales (Dothideomycetes, Ascomycota) occur in a wide range of habitats as endophytes, saprobes, and pathogens. The order Botryosphaeriales has not been subjected to evaluation since 2019 by Phillips and co-authors using phylogenetic and evolutionary analyses. Subsequently, many studies introduced novel taxa into the order and revised several families separately. In addition, no ancestral character studies have been conducted for this order. Therefore, in this study, we re-evaluated the character evolution and taxonomic placements of Botryosphaeriales species based on ancestral character evolution, divergence time estimation, and phylogenetic relationships, including all the novel taxa that have been introduced so far. Maximum likelihood, maximum parsimony, and Bayesian inference analyses were conducted on a combined LSU and ITS sequence alignment. Ancestral state reconstruction was carried out for conidial colour, septation, and nutritional mode. Divergence times estimates revealed that Botryosphaeriales originated around 109 Mya in the early epoch of the Cretaceous period. All six families in Botryosphaeriales evolved in the late epoch of the Cretaceous period (66-100 Mya), during which Angiosperms also appeared, rapidly diversified and became dominant on land. Families of Botryosphaeriales diversified during the Paleogene and Neogene periods in the Cenozoic era. The order comprises the families Aplosporellaceae, Botryosphaeriaceae, Melanopsaceae, Phyllostictaceae, Planistromellaceae and Saccharataceae. Furthermore, current study assessed two hypotheses; the first one being \"All Botryosphaeriales species originated as endophytes and then switched into saprobes when their hosts died or into pathogens when their hosts were under stress\"; the second hypothesis states that \"There is a link between the conidial colour and nutritional mode in botryosphaerialean taxa\". Ancestral state reconstruction and nutritional mode analyses revealed a pathogenic/saprobic nutritional mode as the ancestral character. However, we could not provide strong evidence for the first hypothesis mainly due to the significantly low number of studies reporting the endophytic botryosphaerialean taxa. Results also showed that hyaline and aseptate conidia were ancestral characters in Botryosphaeriales and supported the relationship between conidial pigmentation and the pathogenicity of Botryosphaeriales species.

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.