Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods

A step in the deep evolution of Alvinellidae (Annelida: Polychaeta): A phylogenomic comparative approach based on transcriptomes

The Alvinellidae are a family of worms that are endemic to deep-sea hydrothermal vents in the Pacific and Indian Oceans. These annelid worms, a sister group to the Ampharetidae, occupy a wide range of thermal habitats. The family includes the most thermotolerant marine animals described to date such as the Pompeii worm Alvinella pompejana, and other species living at much lower temperatures such as Paralvinella grasslei or Paralvinella pandorae. The phylogeny of this family has not been studied extensively. It is, however, a complex case where molecular phylogenies have given conflicting results, especially concerning the monophyletic or polyphyletic nature of the genus Paralvinella. We carried out a comprehensive study of the phylogeny of this family using the best molecular data currently available from RNAseq datasets, leading to several hundred orthologous transcripts for 11 of the 14 species currently described or in description. The results obtained by the two most popular phylogenetic inference models (using either gene concatenation with maximum likelihood, or a coalescent-based model from gene trees) were compared using a series of ampharetid and terebellid outgroups. Our study shows that the global phylogenetic signal favors the hypothesis of paraphyly for the Paralvinella genus, with P. pandorae being sister species of the other Alvinellidae. However, a high number of gene trees also supports the hypothesis of alternative trees in which the monophyly of the Paralvinella genus, as initially proposed by Desbruyères and Laubier, is valid with the species P. pandorae and Paralvinella unidentata being sister species. According to molecular dating, the radiation of the Alvinellidae was rapid and took place in a short period of time between 80 and 110 million years ago. This is reflected at the genomic scale by gene trees equally separated between different phylogenetic hypothesis, showing high rates of incomplete lineage sorting between the first lineages of the Alvinellidae and probable gene transfers. Although different genomic regions seem to have different phylogenetic stories in the early step of the alvinellid radiation, our study supports the view that the two P. pandorae species can be grouped into a separate genus (possibly Nautalvinella) and that the Miralvinella subgenus, defined by Desbruyères and Laubier, is not valid anymore.

A multivalent mRNA-LNP vaccine protects against Clostridioides difficile infection

Clostridioides difficile infection (CDI) is an urgent public health threat with limited preventative options. In this work, we developed a messenger RNA (mRNA)-lipid nanoparticle (LNP) vaccine targeting C. difficile toxins and virulence factors. This multivalent vaccine elicited robust and long-lived systemic and mucosal antigen-specific humoral and cellular immune responses across animal models, independent of changes to the intestinal microbiota. Vaccination protected mice from lethal CDI in both primary and recurrent infection models, and inclusion of non-toxin cellular and spore antigens improved decolonization of toxigenic C. difficile from the gastrointestinal tract. Our studies demonstrate mRNA-LNP vaccine technology as a promising platform for the development of novel C. difficile therapeutics with potential for limiting acute disease and promoting bacterial decolonization.

Farmed fur animals harbour viruses with zoonotic spillover potential

Animals such as raccoon dogs, mink and muskrats are farmed for fur and are sometimes used as food or medicinal products1,2, yet they are also potential reservoirs of emerging pathogens3. Here we performed single-sample metatranscriptomic sequencing of internal tissues from 461 individual fur animals that were found dead due to disease. We characterized 125 virus species, including 36 that were novel and 39 at potentially high risk of cross-species transmission, including zoonotic spillover. Notably, we identified seven species of coronaviruses, expanding their known host range, and documented the cross-species transmission of a novel canine respiratory coronavirus to raccoon dogs and of bat HKU5-like coronaviruses to mink, present at a high abundance in lung tissues. Three subtypes of influenza A virus-H1N2, H5N6 and H6N2-were detected in the lungs of guinea pig, mink and muskrat, respectively. Multiple known zoonotic viruses, such as Japanese encephalitis virus and mammalian orthoreovirus4,5, were detected in guinea pigs. Raccoon dogs and mink carried the highest number of potentially high-risk viruses, while viruses from the Coronaviridae, Paramyxoviridae and Sedoreoviridae families commonly infected multiple hosts. These data also reveal potential virus transmission between farmed animals and wild animals, and from humans to farmed animals, indicating that fur farming represents an important transmission hub for viral zoonoses.

Evolutionary origin and population diversity of a cryptic hybrid pathogen

Cryptic fungal pathogens pose disease management challenges due to their morphological resemblance to known pathogens. Here, we investigated the genomes and phenotypes of 53 globally distributed isolates of Aspergillus section Nidulantes fungi and found 30 clinical isolates-including four isolated from COVID-19 patients-were A. latus, a cryptic pathogen that originated via allodiploid hybridization. Notably, all A. latus isolates were misidentified. A. latus hybrids likely originated via a single hybridization event during the Miocene and harbor substantial genetic diversity. Transcriptome profiling of a clinical isolate revealed that both parental subgenomes are actively expressed and respond to environmental stimuli. Characterizing infection-relevant traits-such as drug resistance and growth under oxidative stress-revealed distinct phenotypic profiles among A. latus hybrids compared to parental and closely related species. Moreover, we identified four features that could aid A. latus taxonomic identification. Together, these findings deepen our understanding of the origin of cryptic pathogens.

Untangling Zebrafish Genetic Annotation: Addressing Complexities and Nomenclature Issues in Orthologous Evaluation of TCOF1 and NOLC1

Treacher Collins syndrome (TCS) is a genetic disorder affecting facial development, primarily caused by mutations in the TCOF1 gene. TCOF1, along with NOLC1, play important roles in ribosomal RNA transcription and processing. Previously, a zebrafish model of TCS successfully recapitulated the main characteristics of the syndrome by knocking down the expression of a gene on chromosome 13 (coding for Uniprot ID B8JIY2), which was identified as the TCOF1 orthologue. However, database updates renamed this gene as nolc1 and the zebrafish database (ZFIN) identified a different gene on chromosome 14 as the TCOF1 orthologue (coding for Uniprot ID E7F9D9). NOLC1 and TCOF1 are large proteins with unstructured regions and repetitive sequences that complicate alignments and comparisons. Also, the additional whole genome duplication of teleosts sets further difficulty. In this study, we present evidence that endorses that NOLC1 and TCOF1 are paralogs, and that the zebrafish gene on chromosome 14 is a low-complexity LisH domain-containing factor that displays homology to NOLC1 but lacks essential sequence features to accomplish TCOF1 nucleolar functions. Our analysis also supports the idea that zebrafish, as has been suggested for other non-tetrapod vertebrates, lack the TCOF1 gene that is associated with tripartite nucleolus. Using BLAST searches in a group of teleost genomes, we identified fish-specific sequences similar to E7F9D9 zebrafish protein. We propose naming them "LisH-containing Low Complexity Proteins" (LLCP). Interestingly, the gene on chromosome 13 (nolc1) displays the sequence features, developmental expression patterns, and phenotypic impact of depletion that are characteristic of TCOF1 functions. These findings suggest that in teleost fish, the nucleolar functions described for both NOLC1 and TCOF1 mediated by their repeated motifs, are carried out by a single gene, nolc1. Our study, which is mainly based on computational tools available as free web-based algorithms, could help to solve similar conflicts regarding gene orthology in zebrafish.

A Late Devonian coelacanth reconfigures actinistian phylogeny, disparity, and evolutionary dynamics

The living coelacanth Latimeria (Sarcopterygii: Actinistia) is an iconic, so-called 'living fossil' within one of the most apparently morphologically conservative vertebrate groups. We describe a new, 3-D preserved coelacanth from the Late Devonian Gogo Formation in Western Australia. We assemble a comprehensive analysis of the group to assess the phylogeny, evolutionary rates, and morphological disparity of all coelacanths. We reveal a major shift in morphological disparity between Devonian and post-Devonian coelacanths. The newly described fossil fish fills a critical transitional stage in coelacanth disparity and evolution. Since the mid-Cretaceous, discrete character changes (representing major morphological innovations) have essentially ceased, while meristic and continuous characters have continued to evolve within coelacanths. Considering a range of putative environmental drivers, tectonic activity best explains variation in the rates of coelacanth evolution.

In Silico Predicting the Presence of the S100B Motif in Edible Plants and Detecting Its Immunoreactive Materials: Perspectives for Functional Foods, Dietary Supplements and Phytotherapies

The protein S100B is a part of the S100 protein family, which consists of at least 25 calcium-binding proteins. S100B is highly conserved across different species, supporting important biological functions. The protein was shown to play a role in gut microbiota eubiosis and is secreted in human breast milk, suggesting a physiological trophic function in newborn development. This study explores the possible presence of the S100B motif in plant genomes, and of S100B-like immunoreactive material in different plant extracts, opening up potential botanical uses for dietary supplementation. To explore the presence of the S100B motif in plants, a bioinformatic workflow was used. In addition, the immunoreactivity of S100B from vegetable and fruit samples was tested using an ELISA assay. The S100B motif was expected in silico in the genome of different edible plants belonging to the Viridiplantae clade, such as Durio zibethinus or Malus domestica and other medicinal species. S100B-like immunoreactive material was also detected in samples from fruits or leaves. The finding of S100B-like molecules in plants sheds new light on their role in phylogenesis and in the food chain. This study lays the foundation to elucidate the possible beneficial effects of plants or derivatives containing the S100B-like principle and their potential use in nutraceuticals.

Mapping Transmission Dynamics and Drug Resistance Surveillance in the Cyprus HIV-1 Epidemic (2017-2021)

The human immunodeficiency virus type 1 (HIV-1) epidemic has been a major public health threat on a global scale since the early 1980s. Despite the introduction of combination antiretroviral therapy (cART), the incidence of new HIV-1 infections continues to rise in some regions around the world. Thus, with the continuous transmission of HIV-1 and the lack of a cure, it is imperative for molecular epidemiological studies to be performed, to monitor the infection and ultimately be able to control the spread of this virus. This work provides a comprehensive molecular epidemiological analysis of the HIV-1 infection in Cyprus, through examining 305 HIV-1 sequences collected between 9 March 2017 and 14 October 2021. Employing advanced statistical and bioinformatic techniques, the research delved deeply into understanding the transmission dynamics of the HIV-1 epidemic in Cyprus, as well as the monitoring of HIV-1's genetic diversity and the surveillance of transmitted drug resistance. The characterization of Cyprus's HIV-1 epidemic revealed a diverse landscape, comprising 21 HIV-1 group M pure subtypes and circulating recombinant forms (CRFs), alongside numerous uncharacterized recombinant strains. Subtypes A1 and B emerged as the most prevalent strains, followed by CRF02_AG. The findings of this study also revealed high levels of transmitted drug resistance (TDR) patterns, raising concerns for the efficacy of cART. The demographic profiles of individuals involved in HIV-1 transmission underscored the disproportionate burden borne by young to middle-aged Cypriot males, particularly those in the MSM community, who reported contracting the virus in Cyprus. An assessment of the spatiotemporal evolutionary dynamics illustrated the global interconnectedness of HIV-1 transmission networks, implicating five continents in the dissemination of strains within Cyprus: Europe, Africa, Asia, North America, and Oceania. Overall, this study advances the comprehension of the HIV-1 epidemic in Cyprus and highlights the importance of understanding HIV-1's transmission dynamics through continuous surveillance efforts. Furthermore, this work emphasizes the critical role of state-of-the-art bioinformatics analyses in addressing the challenges posed by HIV-1 transmission globally, laying the groundwork for public health interventions aimed at curbing its spread and improving patient outcomes.

The evolutionary and transmission dynamics of HIV-1 CRF08_BC

HIV-1 CRF08_BC is a significant subtype in China, though its origin and spread remain incompletely understood. Previous studies using partial genomic data have provided insights but lack comprehensive analysis. Here, we investigate the early evolutionary and spatiotemporal dynamics of HIV-1 CRF08_BC in China and Myanmar using near-complete genome sequences. We analyzed 28 near-complete HIV-1 CRF08_BC genomes from China and Myanmar (1997-2013). Phylogenetic, molecular clock, and Bayesian discrete trait analyses were performed to infer the virus's origin, spread, and associated risk groups. Based on Bayesian time-scaled inference with the best-fitting combination of models determined by marginal likelihood estimation (MLE), we inferred the time to the most recent common ancestor (TMRCA) and evolutionary rate of HIV-1 CRF08_BC to be at 3 October 1991 (95% HPD: 22 February1989-27 November 1993) and 2.30 × 10-3 substitutions per site per year (95% HPD: 1.96 × 10-3-2.63 × 10-3), respectively. Our analysis suggests that HIV-1 CRF08_BC originated in Yunnan Province, China, among injecting drug users, and subsequently spread to other regions. This study provides valuable insights into the early dynamics of HIV-1 CRF08_BC through combined genomic and epidemiological data, which may inform effective prevention and mitigation efforts. However, the limited genomic data influenced the extent of our findings, and challenges in collecting accurate risk group information during surveillance were evident.

Movement of the A-strain maize streak virus in and out of Madagascar

The maize streak virus belongs in the genus Mastrevirus, in the family Geminiviridae. The A-strain of the virus (MSV-A) is recognised as the principal causative agent of the most severe manifestation of maize streak disease (MSD). This disease continues to be a persistent limitation on maize output across sub-Saharan Africa and the nearby Indian Ocean islands. Irrespective of the causes behind the spread of MSV-A, we can determine the paths and speeds with which MSV-A spreads by analysing MSV genome sequence data along with information on when and where samples were taken. This information is valuable for identifying the geographical origins of viral strains that cause sporadic MSD epidemics in specific places and the geographical regions where viruses remain in reservoirs and contribute to prolonged epidemics during outbreaks. Our aim is to utilise these analyses to estimate the timing and origin of the MSV-A that arrived on the island of Madagascar in the Indian Ocean. Specifically, we employ model-based phylogeographic analyses on 524 complete MSV-A genome sequences, which consist of 56 newly obtained genomes from infected maize plants collected in Madagascar. These studies allow us to reconstruct the most likely paths of MSV-A to Madagascar. We found strong evidence for the existence of at least four separate movements of MSV-A variants from East and southern Africa to Madagascar. These movements took place between roughly 1979 (with a 95% highest probability density interval [HPD] ranging from 1976 to 1982) and 2003 (with a 95% HPD ranging from 2002 to 2003). While we inferred that MSV-A variants are spreading at an average rate of 38.9 km/year (with a 95% highest posterior density interval of 34.0-44.4) across their geographical range. Since their arrival in Madagascar, MSV-A variants have been migrating at an average rate of 47.6 km/year (with a 95% highest posterior density interval of 36.05-61.70). Human influences are likely significant contributors to both sporadic long-range movements of MSV-A between mainland Africa and Madagascar, as well as shorter to medium range movements within the island.

Random-Effects Substitution Models for Phylogenetics via Scalable Gradient Approximations

Phylogenetic and discrete-trait evolutionary inference depend heavily on an appropriate characterization of the underlying character substitution process. In this paper, we present random-effects substitution models that extend common continuous-time Markov chain models into a richer class of processes capable of capturing a wider variety of substitution dynamics. As these random-effects substitution models often require many more parameters than their usual counterparts, inference can be both statistically and computationally challenging. Thus, we also propose an efficient approach to compute an approximation to the gradient of the data likelihood with respect to all unknown substitution model parameters. We demonstrate that this approximate gradient enables scaling of sampling-based inference, namely Bayesian inference via Hamiltonian Monte Carlo, under random-effects substitution models across large trees and state-spaces. Applied to a dataset of 583 SARS-CoV-2 sequences, an HKY model with random-effects shows strong signals of nonreversibility in the substitution process, and posterior predictive model checks clearly show that it is a more adequate model than a reversible model. When analyzing the pattern of phylogeographic spread of 1441 influenza A virus (H3N2) sequences between 14 regions, a random-effects phylogeographic substitution model infers that air travel volume adequately predicts almost all dispersal rates. A random-effects state-dependent substitution model reveals no evidence for an effect of arboreality on the swimming mode in the tree frog subfamily Hylinae. Simulations reveal that random-effects substitution models can accommodate both negligible and radical departures from the underlying base substitution model. We show that our gradient-based inference approach is over an order of magnitude more time efficient than conventional approaches.

GTRpmix: A Linked General Time-Reversible Model for Profile Mixture Models

Profile mixture models capture distinct biochemical constraints on the amino acid substitution process at different sites in proteins. These models feature a mixture of time-reversible models with a common matrix of exchangeabilities and distinct sets of equilibrium amino acid frequencies known as profiles. Combining the exchangeability matrix with each profile generates the matrix of instantaneous rates of amino acid exchange for that profile. Currently, empirically estimated exchangeability matrices (e.g. the LG matrix) are widely used for phylogenetic inference under profile mixture models. However, these were estimated using a single profile and are unlikely optimal for profile mixture models. Here, we describe the GTRpmix model that allows maximum likelihood estimation of a common exchangeability matrix under any profile mixture model. We show that exchangeability matrices estimated under profile mixture models differ from the LG matrix, dramatically improving model fit and topological estimation accuracy for empirical test cases. Because the GTRpmix model is computationally expensive, we provide two exchangeability matrices estimated from large concatenated phylogenomic-supermatrices to be used for phylogenetic analyses. One, called Eukaryotic Linked Mixture (ELM), is designed for phylogenetic analysis of proteins encoded by nuclear genomes of eukaryotes, and the other, Eukaryotic and Archaeal Linked mixture (EAL), for reconstructing relationships between eukaryotes and Archaea. These matrices, combined with profile mixture models, fit data better and have improved topology estimation relative to the LG matrix combined with the same mixture models. Starting with version 2.3.1, IQ-TREE2 allows users to estimate linked exchangeabilities (i.e. amino acid exchange rates) under profile mixture models.

Molecular basis of one-step methyl anthranilate biosynthesis in grapes, sweet orange, and maize

Plants synthesize an array of volatile compounds, many of which serve ecological roles in attracting pollinators, deterring herbivores, and communicating with their surroundings. Methyl anthranilate (MeAA) is an anti-herbivory defensive volatile responsible for grape aroma that is emitted by several agriculturally relevant plants, including citrus, grapes, and maize. Unlike maize, which uses a one-step anthranilate methyltransferase (AAMT), grapes have been thought to use a two-step pathway for MeAA biosynthesis. By mining available transcriptomics data, we identified two AAMTs in Vitis vinifera (wine grape), as well as one ortholog in "Concord" grape. Many angiosperms methylate the plant hormone salicylic acid (SA) to produce methyl salicylate, which acts as a plant-to-plant communication molecule. Because the Citrus sinensis (sweet orange) SA methyltransferase can methylate both anthranilate (AA) and SA, we used this enzyme to examine the molecular basis of AA activity by introducing rational mutations, which identified several active site residues that increase activity with AA. Reversing this approach, we introduced mutations that imparted activity with SA in the maize AAMT, which uncovered different active site residues from those in the citrus enzyme. Sequence and phylogenetic analysis revealed that one of the Vitis AAMTs shares an ancestor with jasmonic acid methyltransferases, similar to the AAMT from strawberry (Frageria sp.). Collectively, these data demonstrate the molecular mechanisms underpinning AA activity across methyltransferases and identify one-step enzymes by which grapes synthesize MeAA.

Intrahost evolution leading to distinct lineages in the upper and lower respiratory tracts during SARS-CoV-2 prolonged infection

Accumulating evidence points to persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in immunocompromised individuals as a source of novel lineages. While intrahost evolution of the virus in chronically infected patients has previously been reported, existing knowledge is primarily based on samples from the nasopharynx. In this study, we investigate the intrahost evolution and genetic diversity that accumulated during a prolonged SARS-CoV-2 infection with the Omicron BF.7 sublineage, which is estimated to have persisted for >1 year in an immunosuppressed patient. Based on the sequencing of eight samples collected at six time points, we identified 87 intrahost single-nucleotide variants, 2 indels, and a 362-bp deletion. Our analysis revealed distinct viral genotypes in the nasopharyngeal (NP), endotracheal aspirate, and bronchoalveolar lavage samples. This suggests that NP samples may not offer a comprehensive representation of the overall intrahost viral diversity. Our findings not only demonstrate that the Omicron BF.7 sublineage can further diverge from its already exceptionally mutated state but also highlight that patients chronically infected with SARS-CoV-2 can develop genetically specific viral populations across distinct anatomic compartments. This provides novel insights into the intricate nature of viral diversity and evolution dynamics in persistent infections.

Analysis of fine-scale phylogeny of Burkholderia pseudomallei in relation to regional geography and drug susceptibility in Thailand

Melioidosis caused by Burkholderia pseudomallei (Bp) is a public health threat. Genomic-epidemiology research on this deadly disease is scarce. We investigated whole-genome sequences of Bp isolates in relation to environmental source and drug susceptibility. In total, 563 Bp isolates were collected from 11 Northeast Thai provinces during the period 2004-2021. Patients (n = 530 isolates), infected animals (n = 8), and environmental sources (n = 25) provided samples. Phylogenetic analysis revealed genetic diversity among the Bp isolates, including numerous well-supported clusters of varying sizes. Through in-depth analysis of 38 monophyletic clades (MCs), we found eleven associated with province of origin (p-value < 0.001). Closely related clusters (CRCs) within MCs resembled MLST-identified "sequence types" (STs). We found 102 known and 52 novel STs. ST-70 was the most prevalent in this area (n = 78; 13.85%). Sample type (human/environmental) and sampling time intervals were not correlated with genetic distance among clonal Bp isolates. Some members of 12 CRCs had acquired resistance to co-trimoxazole and one against amoxicillin-clavulanic acid. Within Northeast Thailand, there is an association between Bp genotype and geographical origin.

phyloBARCODER: A Web Tool for Phylogenetic Classification of Eukaryote Metabarcodes Using Custom Reference Databases

We developed phyloBARCODER (https://github.com/jun-inoue/phyloBARCODER), a new web tool that can identify short DNA sequences to the species level using metabarcoding. phyloBARCODER estimates phylogenetic trees based on the uploaded anonymous DNA sequences and reference sequences from databases. Without such phylogenetic contexts, alternative, similarity-based methods independently identify species names and anonymous sequences of the same group by pairwise comparisons between queries and database sequences, with the caveat that they must match exactly or very closely. By putting metabarcoding sequences into a phylogenetic context, phyloBARCODER accurately identifies (i) species or classification of query sequences and (ii) anonymous sequences associated with the same species or even with populations of query sequences, with clear and accurate explanations. Version 1 of phyloBARCODER stores a database comprising all eukaryotic mitochondrial gene sequences. Moreover, by uploading their own databases, phyloBARCODER users can conduct species identification specialized for sequences obtained from a local geographic region or those of nonmitochondrial genes, e.g. ITS or rbcL.

Functional trait mismatch between native and introduced bee pollinators servicing a global fruit crop

BACKGROUND

Understanding connections between biodiversity and ecosystem services can be enhanced by shifting focus from species richness to functional trait-based approaches, that when paired with comparative phylogenetic methods can provide even deeper insights. We investigated the functional ecology and phylogenetic diversity of pollination services provided by hymenopteran insects visiting apple flowers in orchards surrounded by either 'natural' or 'disturbed' landscapes in New South Wales, Australia. We assessed whether morphological and behavioural traits (hairiness, body size, glossa length, pollen load purity, and probability of loose pollen) exhibited non-random phylogenetic patterns. Then, explored whether bees, the primary pollinators in this system, filled unique or overlapping functional entities (FEs). For each landscape, we calculated phylogenetic diversity and used FEs to assess functional richness, evenness, and diversion.

RESULTS

A phylogenomic matrix based on ultraconserved elements (UCEs; 1,382,620 bp from 1,969 loci) was used to infer a fully-resolved and well-supported maximum likelihood phylogeny for 48 hymenopteran morphospecies. There was no significant difference in species richness between landscape categories. Pollinator communities at natural sites had higher phylogenetic complexity (X = 2.37) and functional divergence (x̄ = 0.74 ± 0.02 s.e.) than disturbed sites (X = 1.65 and x̄ = 0.6 ± 0.01 s.e.). Hairiness showed significant phylogenetic clustering (K = 0.94), whereas body size, glossa length, and loose pollen showed weaker non-random phylogenetic patterns (K between 0.3-0.5). Pollen load purity showed no association with phylogeny. The assemblage of 17 bee morphospecies comprised nine FEs: eight FEs consisted of native bees with three containing 65% of all native bee taxa. The introduced honey bee (Apis mellifera) occupied a unique FE, likely due to its different evolutionary history. Both landscape types supported six FEs each with three overlapping: two native bee FEs and the honey bee FE.

CONCLUSIONS

Bee hairiness was the only functional trait to exhibit demonstrable phylogenetic signal. Despite differences in species richness, and functional and phylogenetic diversity between orchard landscape types, both maintained equal bee FE numbers. While no native bee taxon was analogous to the honey bee FE, four native bee FEs shared the same hairiness level as honey bees. Health threats to honey bee populations in Australia will likely disrupt pollination services to apple, and other pollination-dependent food crops, given the low level of functional redundancy within the investigated pollinator assemblages.

Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

A new Paleogene fossil and a new dataset for waterfowl (Aves: Anseriformes) clarify phylogeny, ecological evolution, and avian evolution at the K-Pg Boundary

Despite making up one of the most ecologically diverse groups of living birds, comprising soaring, diving and giant flightless taxa, the evolutionary relationships and ecological evolution of Anseriformes (waterfowl) remain unresolved. Although Anseriformes have a comparatively rich, global Cretaceous and Paleogene fossil record, morphological datasets for this group that include extinct taxa report conflicting relationships for all known extinct taxa. Correct placement of extinct taxa is necessary to understand whether ancestral anseriform feeding ecology was more terrestrial or one of a set of diverse aquatic ecologies and to better understand avian evolution around the K-T boundary. Here, we present a new morphological dataset for Anseriformes that includes more extant and extinct taxa than any previous anseriform-focused dataset and describe a new anseriform species from the early Eocene Green River Formation of North America. The new taxon has a mediolaterally narrow bill which is rarely found in previously described anseriform fossils. The matrix created to assess the placement of this taxon comprises 41 taxa and 719 discrete morphological characters describing skeletal morphology, musculature, syringeal morphology, ecology, and behavior. We additionally combine the morphological dataset with published sequences using Bayesian methods and perform ancestral state reconstruction for select morphological, ecological and behavioral characters. We recover the new Eocene taxon as the sister taxon to (Anseranatidae+Anatidae) across all analyses, and find that the new taxon represents a novel ecology within known Anseriformes and the Green River taxa. Results provide insight into avian evolution during and following the K-Pg mass extinction and indicate that Anseriformes were likely ancestrally aquatic herbivores with rhamphothecal lamellae..

MAST: Phylogenetic Inference with Mixtures Across Sites and Trees

Hundreds or thousands of loci are now routinely used in modern phylogenomic studies. Concatenation approaches to tree inference assume that there is a single topology for the entire dataset, but different loci may have different evolutionary histories due to incomplete lineage sorting (ILS), introgression, and/or horizontal gene transfer; even single loci may not be treelike due to recombination. To overcome this shortcoming, we introduce an implementation of a multi-tree mixture model that we call mixtures across sites and trees (MAST). This model extends a prior implementation by Boussau et al. (2009) by allowing users to estimate the weight of each of a set of pre-specified bifurcating trees in a single alignment. The MAST model allows each tree to have its own weight, topology, branch lengths, substitution model, nucleotide or amino acid frequencies, and model of rate heterogeneity across sites. We implemented the MAST model in a maximum-likelihood framework in the popular phylogenetic software, IQ-TREE. Simulations show that we can accurately recover the true model parameters, including branch lengths and tree weights for a given set of tree topologies, under a wide range of biologically realistic scenarios. We also show that we can use standard statistical inference approaches to reject a single-tree model when data are simulated under multiple trees (and vice versa). We applied the MAST model to multiple primate datasets and found that it can recover the signal of ILS in the Great Apes, as well as the asymmetry in minor trees caused by introgression among several macaque species. When applied to a dataset of 4 Platyrrhine species for which standard concatenated maximum likelihood (ML) and gene tree approaches disagree, we observe that MAST gives the highest weight (i.e., the largest proportion of sites) to the tree also supported by gene tree approaches. These results suggest that the MAST model is able to analyze a concatenated alignment using ML while avoiding some of the biases that come with assuming there is only a single tree. We discuss how the MAST model can be extended in the future.

Bayesian Phylogenetic Analysis on Multi-Core Compute Architectures: Implementation and Evaluation of BEAGLE in RevBayes With MPI

Phylogenies are central to many research areas in biology and commonly estimated using likelihood-based methods. Unfortunately, any likelihood-based method, including Bayesian inference, can be restrictively slow for large datasets-with many taxa and/or many sites in the sequence alignment-or complex substitutions models. The primary limiting factor when using large datasets and/or complex models in probabilistic phylogenetic analyses is the likelihood calculation, which dominates the total computation time. To address this bottleneck, we incorporated the high-performance phylogenetic library BEAGLE into RevBayes, which enables multi-threading on multi-core CPUs and GPUs, as well as hardware specific vectorized instructions for faster likelihood calculations. Our new implementation of RevBayes+BEAGLE retains the flexibility and dynamic nature that users expect from vanilla RevBayes. In addition, we implemented native parallelization within RevBayes without an external library using the message passing interface (MPI); RevBayes+MPI. We evaluated our new implementation of RevBayes+BEAGLE using multi-threading on CPUs and 2 different powerful GPUs (NVidia Titan V and NVIDIA A100) against our native implementation of RevBayes+MPI. We found good improvements in speedup when multiple cores were used, with up to 20-fold speedup when using multiple CPU cores and over 90-fold speedup when using multiple GPU cores. The improvement depended on the data type used, DNA or amino acids, and the size of the alignment, but less on the size of the tree. We additionally investigated the cost of rescaling partial likelihoods to avoid numerical underflow and showed that unnecessarily frequent and inefficient rescaling can increase runtimes up to 4-fold. Finally, we presented and compared a new approach to store partial likelihoods on branches instead of nodes that can speed up computations up to 1.7 times but comes at twice the memory requirements.

Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

CMAPLE: Efficient Phylogenetic Inference in the Pandemic Era

We have recently introduced MAPLE (MAximum Parsimonious Likelihood Estimation), a new pandemic-scale phylogenetic inference method exclusively designed for genomic epidemiology. In response to the need for enhancing MAPLE's performance and scalability, here we present two key components: (i) CMAPLE software, a highly optimized C++ reimplementation of MAPLE with many new features and advancements, and (ii) CMAPLE library, a suite of application programming interfaces to facilitate the integration of the CMAPLE algorithm into existing phylogenetic inference packages. Notably, we have successfully integrated CMAPLE into the widely used IQ-TREE 2 software, enabling its rapid adoption in the scientific community. These advancements serve as a vital step toward better preparedness for future pandemics, offering researchers powerful tools for large-scale pathogen genomic analysis.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Dispersal history of SARS-CoV-2 in Galicia, Spain

The dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron-BA.1 variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

Human Y chromosome haplogroup L1-M22 traces Neolithic expansion in West Asia and supports the Elamite and Dravidian connection

West and South Asian populations profoundly influenced Eurasian genetic and cultural diversity. We investigate the genetic history of the Y chromosome haplogroup L1-M22, which, while prevalent in these regions, lacks in-depth study. Robust Bayesian analyses of 165 high-coverage Y chromosomes favor a West Asian origin for L1-M22 ∼20.6 thousand years ago (kya). Moreover, this haplogroup parallels the genome-wide genetic ancestry of hunter-gatherers from the Iranian Plateau and the Caucasus. We characterized two L1-M22 harboring population groups during the Early Holocene. One expanded with the West Asian Neolithic transition. The other moved to South Asia ∼8-6 kya but showed no expansion. This group likely participated in the spread of Dravidian languages. These South Asian L1-M22 lineages expanded ∼4-3 kya, coinciding with the Steppe ancestry introduction. Our findings advance the current understanding of Eurasian historical dynamics, emphasizing L1-M22's West Asian origin, associated population movements, and possible linguistic impacts.

The trunk replaces the longer mandible as the main feeding organ in elephant evolution

The long-trunked elephantids underwent a significant evolutionary stage characterized by an exceptionally elongated mandible. The initial elongation and subsequent regression of the long mandible, along with its co-evolution with the trunk, present an intriguing issue that remains incompletely understood. Through comparative functional and eco-morphological investigations, as well as feeding preference analysis, we reconstructed the feeding behavior of major groups of longirostrine elephantiforms. In the Platybelodon clade, the rapid evolutionary changes observed in the narial region, strongly correlated with mandible and tusk characteristics, suggest a crucial evolutionary transition where feeding function shifted from the mandible to the trunk, allowing proboscideans to expand their niches to more open regions. This functional shift further resulted in elephantids relying solely on their trunks for feeding. Our research provides insights into how unique environmental pressures shape the extreme evolution of organs, particularly in large mammals that developed various peculiar adaptations during the late Cenozoic global cooling trends.

STMND1 is a phylogenetically ancient stathmin which localizes to motile cilia and exhibits nuclear translocation that is inhibited when soluble tubulin concentration increases

Stathmins are small, unstructured proteins that bind tubulin dimers and are implicated in several human diseases, but whose function remains unknown. We characterized a new stathmin, STMND1 (Stathmin Domain Containing 1) as the human representative of an ancient subfamily. STMND1 features a N-terminal myristoylated and palmitoylated motif which directs it to membranes and a tubulin-binding stathmin-like domain (SLD) that contains an internal nuclear localization signal. Biochemistry and proximity labeling showed that STMND1 binds tubulin, and live imaging showed that tubulin binding inhibits translocation from cellular membranes to the nucleus. STMND1 is highly expressed in multiciliated epithelial cells, where it localizes to motile cilia. Overexpression in a model system increased the length of primary cilia. Our study suggests that the most ancient stathmins have cilium-related functions that involve sensing soluble tubulin.

Long-term evolution of Streptococcus mitis and Streptococcus pneumoniae leads to higher genetic diversity within rather than between human populations

Evaluation of the apportionment of genetic diversity of human bacterial commensals within and between human populations is an important step in the characterization of their evolutionary potential. Recent studies showed a correlation between the genomic diversity of human commensal strains and that of their host, but the strength of this correlation and of the geographic structure among human populations is a matter of debate. Here, we studied the genomic diversity and evolution of the phylogenetically related oro-nasopharyngeal healthy-carriage Streptococcus mitis and Streptococcus pneumoniae, whose lifestyles range from stricter commensalism to high pathogenic potential. A total of 119 S. mitis genomes showed higher within- and among-host variation than 810 S. pneumoniae genomes in European, East Asian and African populations. Summary statistics of the site-frequency spectrum for synonymous and non-synonymous variation and ABC modelling showed this difference to be due to higher ancestral bacterial population effective size (Ne) in S. mitis, whose genomic variation has been maintained close to mutation-drift equilibrium across (at least many) generations, whereas S. pneumoniae has been expanding from a smaller ancestral bacterial population. Strikingly, both species show limited differentiation among human populations. As genetic differentiation is inversely proportional to the product of effective population size and migration rate (Nem), we argue that large Ne have led to similar differentiation patterns, even if m is very low for S. mitis. We conclude that more diversity within than among human populations and limited population differentiation must be common features of the human microbiome due to large Ne.

Nature of epigenetic aging from a single-cell perspective

Age-related changes in DNA methylation (DNAm) form the basis of the most robust predictors of age-epigenetic clocks-but a clear mechanistic understanding of exactly which aspects of aging are quantified by these clocks is lacking. Here, to clarify the nature of epigenetic aging, we juxtapose the dynamics of tissue and single-cell DNAm in mice. We compare these changes during early development with those observed during adult aging in mice, and corroborate our analyses with a single-cell RNA sequencing analysis within the same multiomics dataset. We show that epigenetic aging involves co-regulated changes as well as a major stochastic component, and this is consistent with transcriptional patterns. We further support the finding of stochastic epigenetic aging by direct tissue and single-cell DNAm analyses and modeling of aging DNAm trajectories with a stochastic process akin to radiocarbon decay. Finally, we describe a single-cell algorithm for the identification of co-regulated and stochastic CpG clusters showing consistent transcriptomic coordination patterns. Together, our analyses increase our understanding of the basis of epigenetic clocks and highlight potential opportunities for targeting aging and evaluating longevity interventions.

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.

Analysis of Genetic Structure and Genetic Diversity in Japanese Grey-Headed Lapwing Population Using mtDNA

The grey-headed lapwing (Vanellus cinereus) is a wading species in East Asia. However, examples of regional population dynamics and genetic research are limited. To reconsider the natural history and current status of the grey-headed lapwing in Japan, we analyzed the genetic diversity of the Japanese grey-headed lapwing population. We collected 77 grey-headed lapwing samples from 12 locations across Japan during the breeding season and three individuals during the wintering season and extracted DNA; 496-bp sequences of the ND2, which form part of the mitochondrial DNA, were determined for genetic analysis of the population. Consequently, 10 haplotypes were detected in 80 individuals, and 67 individuals, 84% of the total, shared two haplotypes, namely Vc1 and Vc2. Furthermore, the results showed that the prevalence of Vc1 was higher mainly in northern Japan, while that of Vc2 was higher mainly in southern Japan. Genetic diversity analysis showed that the overall haplotype diversity in Japan was 0.617, which is not particularly low. The sequence of Vc1 was exactly the same as that of grey-headed lapwing in China. Our study revealed the genetic structure of the grey-headed lapwing, suggesting that as the grey-headed lapwing expanded its distribution area into southern Japan, many Vc2-positive individuals migrated southward, resulting in a higher detection rate of Vc2 in southern Japan.

Complete Mitochondrial Genomes of Nannostomus Pencilfish: Genome Characterization and Phylogenetic Analysis

Although the pencilfish is a globally popular economic fish in the aquarium market, its taxonomic classification could be further refined. In order to understand the taxonomy of species of the genus Nannostomus (Characiformes, Lebiasinidae) and their phylogenetic position within the order Characiformes, in this study, we characterized mitochondrial genomes (mitogenomes) from four Nannostomus species for the first time. The four mitogenomes exhibited the typical circular structure, with overall sizes varying from 16,661 bp to 16,690 bp. They contained 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 1 control region (CR). Nucleotide composition analysis suggested that the mitochondrial sequences were biased toward A and T. Bayesian inference and maximum likelihood analyses based on PCGs support the family Lebiasinidae classification, described using four Nannostomus species, clustering together with Lebiasina multimaculata from the same family. The results of this study support the current taxonomic classification of the family Lebiasinidae. Phylogenetic analysis also suggested that gene rearrangement would not significantly impact the phylogenetic relationships within the order Characiformes. These results might provide new data regarding the phylogeny and classification of the order Characiformes, thus providing a theoretical basis for the economic development of aquarium fish markets.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Is Over-parameterization a Problem for Profile Mixture Models?

Biochemical constraints on the admissible amino acids at specific sites in proteins lead to heterogeneity of the amino acid substitution process over sites in alignments. It is well known that phylogenetic models of protein sequence evolution that do not account for site heterogeneity are prone to long-branch attraction (LBA) artifacts. Profile mixture models were developed to model heterogeneity of preferred amino acids at sites via a finite distribution of site classes each with a distinct set of equilibrium amino acid frequencies. However, it is unknown whether the large number of parameters in such models associated with the many amino acid frequency vectors can adversely affect tree topology estimates because of over-parameterization. Here, we demonstrate theoretically that for long sequences, over-parameterization does not create problems for estimation with profile mixture models. Under mild conditions, tree, amino acid frequencies, and other model parameters converge to true values as sequence length increases, even when there are large numbers of components in the frequency profile distributions. Because large sample theory does not necessarily imply good behavior for shorter alignments we explore the performance of these models with short alignments simulated with tree topologies that are prone to LBA artifacts. We find that over-parameterization is not a problem for complex profile mixture models even when there are many amino acid frequency vectors. In fact, simple models with few site classes behave poorly. Interestingly, we also found that misspecification of the amino acid frequency vectors does not lead to increased LBA artifacts as long as the estimated cumulative distribution function of the amino acid frequencies at sites adequately approximates the true one. In contrast, misspecification of the amino acid exchangeability rates can severely negatively affect parameter estimation. Finally, we explore the effects of including in the profile mixture model an additional "F-class" representing the overall frequencies of amino acids in the data set. Surprisingly, the F-class does not help parameter estimation significantly and can decrease the probability of correct tree estimation, depending on the scenario, even though it tends to improve likelihood scores.

A leading-edge scenario in the phylogeography and evolutionary history of East Asian insular Taxus in Taiwan and the Philippines

The cool temperate origin of gymnosperm Taxus species in East Asia is specifically diverse and widespread. Certain lineages have managed to extend their distribution further south to subtropical and tropical islands such as Taiwan and the Philippines. To address questions including whether these insular lineages, recently identified as T. phytonii, have become genetically distinct from each other and from their continental relatives, and when and how they colonized their residing islands, we sampled over 11 populations, covering 179 Taxus individuals from Taiwan and the Philippines. Using four cpDNA and one nuclear marker, we showed in population genetic and genealogical analyses that the two insular lineages were genetically distinct from each other and also from other continental Taxus and that they represented each other's closest relative. Estimated with the coalescent-based multi-type tree (MTT) analyses, we inferred an origin of Taiwanese T. phytonii more ancient than 2.49 Mya and that of Philippine T. phytonii more ancient than 1.08 Mya. In addition, the divergence demographic history revealed by both MTT and isolation with migration (IM) analyses indicated the presence of recent post-split migrations from a continental taxon, T. mairei, to Taiwanese T. phytonii, as well as from Taiwanese T. phytonii to Philippine T. phytonii. Overall, this study suggests Taiwan as a stepping stone through which the temperate-origin yew trees can extend their distributions to tropical regions such as the Philippines.

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Ancestral sequence reconstruction of Mic60 reveals a residue signature supporting respiration in yeast

In eukaryotes, the essential process of cellular respiration takes place in the cristae of mitochondria. The protein Mic60 is known to stabilize crista junctions; however, how the C-terminal Mitofilin domain of Mic60 mediates cristae-supported respiration remains elusive. Here, we used ancestral sequence reconstruction to generate Mitofilin ancestors up to and including the last opisthokont common ancestor (LOCA). We found that yeast-lineage derived Mitofilin ancestors as far back as the LOCA rescue respiration. By comparing Mitofilin ancestors with different respiratory phenotypes, we identify four residues that explain the difference between respiration functional yeast- and non-functional animal-derived common Mitofilin ancestors. Our results imply that Mitofilin-supported respiration in yeast stems from a conserved mechanism, and provide a foundation for investigating the divergence of candidate crista junction interactions present during the emergence of eukaryotes.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Comparison of the performance of multiple whole-genome sequence-based tools for the identification of Bacillus cereus sensu stricto biovar Thuringiensis

The Bacillus cereus sensu stricto (s.s.) species comprises strains of biovar Thuringiensis (Bt) known for their bioinsecticidal activity, as well as strains with foodborne pathogenic potential. Bt strains are identified (i) based on the production of insecticidal crystal proteins, also known as Bt toxins, or (ii) based on the presence of cry, cyt, and vip genes, which encode Bt toxins. Multiple bioinformatics tools have been developed for the detection of crystal protein-encoding genes based on whole-genome sequencing (WGS) data. However, the performance of these tools is yet to be evaluated using phenotypic data. Thus, the goal of this study was to assess the performance of four bioinformatics tools for the detection of crystal protein-encoding genes. The accuracy of sequence-based identification of Bt was determined in reference to phenotypic microscope-based screening for the production of crystal proteins. A total of 58 diverse B. cereus sensu lato strains isolated from clinical, food, environmental, and commercial biopesticide products underwent WGS. Isolates were examined for crystal protein production using phase contrast microscopy. Crystal protein-encoding genes were detected using BtToxin_Digger, BTyper3, IDOPS (identification of pesticidal sequences), and Cry_processor. Out of 58 isolates, the phenotypic production of crystal proteins was confirmed for 18 isolates. Specificity and sensitivity of Bt identification based on sequences were 0.85 and 0.94 for BtToxin_Digger, 0.97 and 0.89 for BTyper3, 0.95 and 0.94 for IDOPS, and 0.88 and 1.00 for Cry_processor, respectively. Cry_processor predicted crystal protein production with the highest specificity, and BtToxin_Digger and IDOPS predicted crystal protein production with the highest sensitivity. Three out of four tested bioinformatics tools performed well overall, with IDOPS achieving high sensitivity and specificity (>0.90).IMPORTANCEStrains of Bacillus cereus sensu stricto (s.s.) biovar Thuringiensis (Bt) are used as organic biopesticides. Bt is differentiated from the foodborne pathogen Bacillus cereus s.s. by the production of insecticidal crystal proteins. Thus, reliable genomic identification of biovar Thuringiensis is necessary to ensure food safety and facilitate risk assessment. This study assessed the accuracy of whole-genome sequencing (WGS)-based identification of Bt compared to phenotypic microscopy-based screening for crystal protein production. Multiple bioinformatics tools were compared to assess their performance in predicting crystal protein production. Among them, identification of pesticidal sequences performed best overall at WGS-based Bt identification.

Quantitative tests of albendazole resistance in beta-tubulin mutants

Benzimidazole (BZ) anthelmintics are among the most important treatments for parasitic nematode infections in the developing world. Widespread BZ resistance in veterinary parasites and emerging resistance in human parasites raise major concerns for the continued use of BZs. Knowledge of the mechanisms of resistance is necessary to make informed treatment decisions and circumvent resistance. Benzimidazole resistance has traditionally been associated with mutations and natural variants in the C. elegans beta-tubulin gene ben-1 and orthologs in parasitic species. However, variants in ben-1 alone do not explain the differences in BZ responses across parasite populations. Here, we examine the roles of five C. elegans beta-tubulin genes (tbb-1, mec-7, tbb-4, ben-1, and tbb-6) to identify the role each gene plays in BZ response. We generated C. elegans strains with a loss of each beta-tubulin gene, as well as strains with a loss of tbb-1, mec-7, tbb-4, or tbb-6 in a genetic background that also lacks ben-1 to test beta-tubulin redundancy in BZ response. We found that only the individual loss of ben-1 conferred a substantial level of BZ resistance, although the loss of tbb-1 was found to confer a small benefit in the presence of albendazole (ABZ). The loss of ben-1 was found to confer an almost complete rescue of animal development in the presence of 30 μM ABZ, likely explaining why no additive effects caused by the loss of a second beta-tubulin were observed. We demonstrate that ben-1 is the only beta-tubulin gene in C. elegans where loss confers substantial BZ resistance.

Punctuational evolution is pervasive in distal site metastatic colonization

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

The genomic evolutionary dynamics and global circulation patterns of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection in young children and the second leading cause of infant death worldwide. While global circulation has been extensively studied for respiratory viruses such as seasonal influenza, and more recently also in great detail for SARS-CoV-2, a lack of global multi-annual sampling of complete RSV genomes limits our understanding of RSV molecular epidemiology. Here, we capitalise on the genomic surveillance by the INFORM-RSV study and apply phylodynamic approaches to uncover how selection and neutral epidemiological processes shape RSV diversity. Using complete viral genome sequences, we show similar patterns of site-specific diversifying selection among RSVA and RSVB and recover the imprint of non-neutral epidemic processes on their genealogies. Using a phylogeographic approach, we provide evidence for air travel governing the global patterns of RSVA and RSVB spread, which results in a considerable degree of phylogenetic mixing across countries. Our findings highlight the potential of systematic global RSV genomic surveillance for transforming our understanding of global RSV spread.

Evolutionary genetics of pulmonary anatomical adaptations in deep-diving cetaceans

BACKGROUND

Cetaceans, having experienced prolonged adaptation to aquatic environments, have undergone evolutionary changes in their respiratory systems. This process of evolution has resulted in the emergence of distinctive phenotypic traits, notably the abundance of elastic fibers and thickened alveolar walls in their lungs, which may facilitate alveolar collapse during diving. This structure helps selective exchange of oxygen and carbon dioxide, while minimizing nitrogen exchange, thereby reducing the risk of DCS. Nevertheless, the scientific inquiry into the mechanisms through which these unique phenotypic characteristics govern the diving behavior of marine mammals, including cetaceans, remains unresolved.

RESULTS

This study entails an evolutionary analysis of 42 genes associated with pulmonary fibrosis across 45 mammalian species. Twenty-one genes in cetaceans exhibited accelerated evolution, featuring specific amino acid substitutions in 14 of them. Primarily linked to the development of the respiratory system and lung morphological construction, these genes play a crucial role. Moreover, among marine mammals, we identified eight genes undergoing positive selection, and the evolutionary rates of three genes significantly correlated with diving depth. Specifically, the SFTPC gene exhibited convergent amino acid substitutions. Through in vitro cellular experiments, we illustrated that convergent amino acid site mutations in SFTPC contribute positively to pulmonary fibrosis in marine mammals, and the presence of this phenotype can induce deep alveolar collapse during diving, thereby reducing the risk of DCS during diving.

CONCLUSIONS

The study unveils pivotal genetic signals in cetaceans and other marine mammals, arising through evolution. These genetic signals may influence lung characteristics in marine mammals and have been linked to a reduced risk of developing DCS. Moreover, the research serves as a valuable reference for delving deeper into human diving physiology.

Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation

Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

Jurassic shuotheriids show earliest dental diversification of mammaliaforms

Shuotheriids are Jurassic mammaliaforms that possess pseudotribosphenic teeth in which a pseudotalonid is anterior to the trigonid in the lower molar, contrasting with the tribosphenic pattern of therian mammals (placentals, marsupials and kin) in which the talonid is posterior to the trigonid1-4. The origin of the pseudotribosphenic teeth remains unclear, obscuring our perception of shuotheriid affinities and the early evolution of mammaliaforms1,5-9. Here we report a new Jurassic shuotheriid represented by two skeletal specimens. Their complete pseudotribosphenic dentitions allow reidentification of dental structures using serial homology and the tooth occlusal relationship. Contrary to the conventional view1,2,6,10,11, our findings show that dental structures of shuotheriids can be homologized to those of docodontans and partly support homologous statements for some dental structures between docodontans and other mammaliaforms6,12. The phylogenetic analysis based on new evidence removes shuotheriids from the tribosphenic ausktribosphenids (including monotremes) and clusters them with docodontans to form a new clade, Docodontiformes, that is characterized by pseudotribosphenic features. In the phylogeny, docodontiforms and 'holotherians' (Kuehneotherium, monotremes and therians)13 evolve independently from a Morganucodon-like ancestor with triconodont molars by labio-lingual widening their posterior teeth for more efficient food processing. The pseudotribosphenic pattern passed a cusp semitriangulation stage9, whereas the tribosphenic pattern and its precursor went through a stage of cusp triangulation. The two different processes resulted in complex tooth structures and occlusal patterns that elucidate the earliest diversification of mammaliaforms.

Targeting Fks1 proteins for novel antifungal drug discovery

Fungal infections are a major threat to human health. The limited availability of antifungal drugs, the emergence of drug resistance, and a growing susceptible population highlight the critical need for novel antifungal agents. The enzymes involved in fungal cell wall synthesis offer potential targets for antifungal drug development. Recent studies have enhanced our focus on the enzyme Fks1, which synthesizes β-1,3-glucan, a critical component of the cell wall. These studies provide a deeper understanding of Fks1's function in cell wall biosynthesis, pathogenicity, structural biology, evolutionary conservation across fungi, and interaction with current antifungal drugs. Here, we discuss the role of Fks1 in the survival and adaptation of fungi, guided by insights from evolutionary and structural analyses. Furthermore, we delve into the dynamics of Fks1 modulation with novel antifungal strategies and assess its potential as an antifungal drug target.

Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.