BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis

Estimating the mean in the space of ranked phylogenetic trees

MOTIVATION

Reconstructing evolutionary histories of biological entities, such as genes, cells, organisms, populations, and species, from phenotypic and molecular sequencing data is central to many biological, palaeontological, and biomedical disciplines. Typically, due to uncertainties and incompleteness in data, the true evolutionary history (phylogeny) is challenging to estimate. Statistical modelling approaches address this problem by introducing and studying probability distributions over all possible evolutionary histories, but can also introduce uncertainties due to misspecification. In practice, computational methods are deployed to learn those distributions typically by sampling them. This approach, however, is fundamentally challenging as it requires designing and implementing various statistical methods over a space of phylogenetic trees (or treespace). Although the problem of developing statistics over a treespace has received substantial attention in the literature and numerous breakthroughs have been made, it remains largely unsolved. The challenge of solving this problem is 2-fold: a treespace has nontrivial often counter-intuitive geometry implying that much of classical Euclidean statistics does not immediately apply; many parametrizations of treespace with promising statistical properties are computationally hard, so they cannot be used in data analyses. As a result, there is no single conventional method for estimating even the most fundamental statistics over any treespace, such as mean and variance, and various heuristics are used in practice. Despite the existence of numerous tree summary methods to approximate means of probability distributions over a treespace based on its geometry, and the theoretical promise of this idea, none of the attempts resulted in a practical method for summarizing tree samples.

RESULTS

In this paper, we present a tree summary method along with useful properties of our chosen treespace while focusing on its impact on phylogenetic analyses of real datasets. We perform an extensive benchmark study and demonstrate that our method outperforms currently most popular methods with respect to a number of important 'quality' statistics. Further, we apply our method to three empirical datasets ranging from cancer evolution to linguistics and find novel insights into corresponding evolutionary problems in all of them. We hence conclude that this treespace is a promising candidate to serve as a foundation for developing statistics over phylogenetic trees analytically, as well as new computational tools for evolutionary data analyses.

AVAILABILITY AND IMPLEMENTATION

An implementation is available at https://github.com/bioDS/Centroid-Code.

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and ∼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Phylogenetic climatic niche evolution and diversification of the Neurergus species (Salamandridae) in the Irano-Anatolian biodiversity hotspot

This study explores how climate variables influenced the evolution and diversification of Neurergus newts within the Irano-Anatolian biodiversity hotspot. We use a dated phylogenetic tree and climatic niche models to analyze their evolutionary history and ecological preferences. Using genetic data from nuclear (KIAA) and mitochondrial (16s and 12s) genes, we estimate divergence times and identify four major Neurergus clades. The initial speciation event occurred approximately 11.3 million years ago, coinciding with the uplift of the Zagros and Anatolian mountains. This geological transformation isolated newt populations, likely triggering the first speciation event. By integrating potential geographic distribution with climate variables, we reconstruct ancestral niche occupancy profiles. This highlights the critical roles of temperature and precipitation in shaping Neurergus habitat preferences and distribution. We observe both phylogenetic niche conservatism and divergence, with niche divergence playing a dominant role in diversification. This research emphasizes the complex interplay of geography, climate, and ecology in speciation and the vulnerability of isolated mountain newt populations to environmental changes.

Application of molecular epidemiology in revealing HIV-1 transmission network and recombination patterns in Tianjin, China

Using a comprehensive molecular epidemiological approach, we characterized the transmission dynamics of HIV-1 among the MSM population in Tianjin, China. Our findings revealed that 38.56% (386/1001) of individuals clustered across 109 molecular transmission clusters (TCs), with MSM aged 50 and below being the group most commonly transmitting HIV-1. Among the identified TCs, CRF01_AE predominated, followed by CRF07_BC. Notably, CRF07_BC demonstrated a higher propensity for forming large clusters compared to CRF01_AE. Birth-death skyline analyses of the two largest clusters indicated that the HIV/AIDS transmission may be at a critical point, nearly all had Re approximately 1 by now. A retrospective analysis revealed that the rapid expansion of these large clusters was primarily driven by the introduction of viruses in 2021, highlighting the crucial importance of continuous molecular surveillance in identifying newly emerging high-risk transmission chains and adapting measures to address evolving epidemic dynamics. Furthermore, we detected the transmission of drug-resistant mutations (DRMs) within the TCs, particularly in the CRF07_BC clusters (K103N, Y181C, and K101E) and CRF01_AE clusters (P225H and K219R), emphasizing the importance of monitoring to support the continued efficacy of first-line therapies and pre-exposure prophylaxis (PrEP). Recombination analyses indicated that complex recombinant patterns, associated with increased amino acid variability, could confer adaptive traits to the viruses, potentially providing a competitive advantage in certain host populations or regions. Our study highlights the potential of integrating molecular epidemiological and phylodynamic approaches to inform targeted interventions.

Population structure of Dugong dugon across the Indo-Pacific revealed by historical mitogenomes

Sirenia, an iconic marine taxon with a tropical and subtropical worldwide distribution, face an uncertain future. All species are designated 'Vulnerable' to extinction by the IUCN. Nonetheless, a comprehensive understanding of geographic structuring across the global range is lacking, impeding our ability to highlight particularly vulnerable populations for conservation priority. Here, we use ancient DNA to investigate dugong (Dugong dugon) population structure, analysing 56 mitogenomes from specimens comprising the known historical range. Our results reveal geographically structured and distinct monophyletic clades characterized by contrasting evolutionary histories. We observe deep-rooted and divergent lineages in the East (Indo-Pacific) and obtain new evidence for the relatively recent dispersal of dugongs into the western Indian Ocean. All populations are significantly differentiated from each other with western populations having approximately 10-fold lower levels of genetic variation than eastern Indo-Pacific populations. Additionally, we find a significant temporal loss of genetic diversity in western Indian Ocean dugongs since the mid-twentieth century, as well as a decline in population size beginning approximately 1000 years ago. Our results add to the growing body of evidence that dugong populations are becoming ever more susceptible to ongoing human action and global climate change.

Molecular characterisation of Amblyomma integrum circulating in southern India

Amblyomma integrum is a large gooseberry sized longirostrate tick (when fully repleted) found in India and Sri Lanka. In Kerala (India), this tick is commonly found in the forest and its fringe areas frequently infesting deer and hence it is locally known as "maan chellu / maanunny" (deer tick). In the present study, molecular characterisation and phylogenetic analysis of A. integrum collected from the area grazed by the sambar deer (Rusa unicolor) of Kerala, south India was performed using three molecular markers viz., the mitochondrial cytochrome c oxidase subunit 1 (COI), mitochondrial 16S ribosomal RNA, and nuclear 18S ribosomal RNA genes. Cytochrome c oxidase subunit 1 (COI) gene showed better resolving ability for elucidating the evolutionary relationship of A. integrum and identified two distinct clades, viz., A and B. The Tamil Nadu isolates of south India and Marayoor isolate 1 (from Idukki district of Kerala bordering with Tamil Nadu) belonged to clade A. Majority of Wayanad isolates from Kerala, occupied clade B. The intraspecific genetic distance among the A. integrum species ranged from 0.00 to 13.34%. Between clades A and B, the genetic distance observed was 11.49%. The clade B isolates were genetically close to A. geoemydae (GD: 1.22%). Morphological variations between the clades included darker exoskeletal coloration in clade A and distinct differences in the shape of basis capitulum. Further analysis using Assemble Species by Automatic Partitioning (ASAP) and Generalized Mixed Yule Coalescent (GMYC) provided additional insights. Assemble Species by Automatic Partitioning (ASAP) identified 26 Molecular Operational Taxonomic Units (MOTUs) at a threshold distance of 5.38%, supporting the species partition of A. integrum clade B. Generalized Mixed Yule Coalescent (GMYC) analysis retained the same species complex (A. integrum-geoemydae Complex) inferred from the ASAP analyses. It could be inferred from the present study that the A. integrum clades A and B could be two different putative pseudocryptic species.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Evolutionary patterns of variations in chromosome counts and genome sizes show positive correlations with taxonomic diversity in tropical gingers

PREMISE

Cytogenetic traits such as an organism's chromosome number and genome size are taxonomically critical as they are instrumental in defining angiosperm diversity. Variations in these traits can be traced to evolutionary processes such as polyploidization, although geographic variations across cytogenetic traits remain underexplored. In the pantropical monocot family Zingiberaceae (~1500 species), cytogenetic traits have been well documented; however, the role of these traits in shaping taxonomic diversity and biogeographic patterns of gingers is not known.

METHODS

A time-calibrated Bayesian phylogenetic tree was constructed for 290 taxa covering three of the four subfamilies in Zingiberaceae. We tested models of chromosome number and genome size evolution within the family and whether lineage age, taxonomic diversity, and distributional range explain the variations in the cytogenetic traits. Tests were carried out at two taxonomic ranks: within Zingiberaceae and within genus Hedychium using correlations, generalized linear models and phylogenetic least square models.

RESULTS

The most frequent changes in chromosome number within Zingiberaceae were noted to be demi-polyploidization and polyploidization (~57% of the time), followed by ascending dysploidy (~27%). The subfamily Zingiberoideae showed descending dysploidy at its base, while Alpinioideae showed polyploidization at its internal nodes. Although chromosome counts and genome sizes did not corroborate with each other, suggesting that they are not equivalent; higher chromosome number variations and higher genome size variations were associated with higher taxonomic diversity and wider biogeographic distribution.

CONCLUSIONS

Within Zingiberaceae, multiple incidences of polyploidization were discovered, and cytogenetic events appear to have reduced the genome sizes and increased taxonomic diversity, distributional ranges and invasiveness.

Establishing species boundaries in Bornean geckos

Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus)

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

Niche breadth and divergence in sympatric cryptic coral species (Pocillopora spp.) across habitats within reefs and among algal symbionts

While the presence of morphologically cryptic species is increasingly recognized, we still lack a useful understanding of what causes and maintains co-occurring cryptic species and its consequences for the ecology, evolution, and conservation of communities. We sampled 724 Pocillopora corals from five habitat zones (the fringing reef, back reef, and fore reef at 5, 10, and 20 m) at four sites around the island of Moorea, French Polynesia. Using validated genetic markers, we identified six sympatric species of Pocillopora, most of which cannot be reliably identified based on morphology: P. meandrina (42.9%), P. tuahiniensis (25.1%), P. verrucosa (12.2%), P. acuta (10.4%), P. grandis (7.73%), and P. cf. effusa (2.76%). For 423 colonies (58% of the genetically identified hosts), we also used psbA ncr or ITS2 markers to identify symbiont species (Symbiodiniaceae). The relative abundance of Pocillopora species differed across habitats within the reef. Sister taxa P. verrucosa and P. tuahiniensis had similar niche breadths and hosted the same specialist symbiont species (mostly Cladocopium pacificum) but the former was more common in the back reef and the latter more common deeper on the fore reef. In contrast, sister taxa P. meandrina and P. grandis had the highest niche breadths and overlaps and tended to host the same specialist symbiont species (mostly C. latusorum). Pocillopora acuta had the narrowest niche breadth and hosted the generalist, and more thermally tolerant, Durusdinium gynnii. Overall, there was a positive correlation between reef habitat niche breadth and symbiont niche breadth-Pocillopora species with a broader habitat niche also had a broader symbiont niche. Our results show how fine-scale variation within reefs plays an important role in the generation and coexistence of cryptic species. The results also have important implications for how niche differences affect community resilience, and for the success of coral restoration practices, in ways not previously appreciated.

Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae)

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

Phylogeography of Parnassius citrinarius Based on Mitochondrial DNA Reveals Large Differences in Genetic Structure between the Eastern and Western Japan

The Japanese Archipelago hosts a rich butterfly fauna, and elucidating the genetic structures of multiple species is necessary to clarify their formation processes. This study aimed to reveal the genetic structure and distribution formation process of Parnassius citrinarius, which is widely distributed across the Japanese Archipelago from Hokkaido to Shikoku, through phylogeographic analysis based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequence. Thirty haplotypes were revealed from 311 individuals from 47 sites, indicating significant differences in the genetic structures between the eastern and western parts of the Japanese Archipelago. In Eastern Japan, multiple genetic clusters were found, with some sites harboring two clusters. The divergence times among populations in Eastern Japan were relatively recent, and no genetic differentiation was observed between regions, including between Hokkaido and Honshu, which are separated by a narrow strait. In contrast, in Western Japan, including Shikoku, unique genetic clusters were observed in each region. The phylogenetic relationships among populations were regionally clustered, and the divergence times were relatively ancient. The distribution and genetic structure of P. citrinarius in the Japanese Archipelago have been significantly influenced by temperature fluctuations and the presence of geographical barriers during the Pleistocene glacial-interglacial cycles, including the potential formation of refugia in Western Japan.

Selection against domestication alleles in introduced rabbit populations

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

One hundred years of influenza A evolution

Leveraging the simplicity of nucleotide mismatch distributions, we provide an intuitive window into the evolution of the human influenza A 'nonstructural' (NS) gene segment. In an analysis suggested by the eminent Danish biologist Freddy B. Christiansen, we illustrate the existence of a continuous genetic "backbone" of influenza A NS sequences, steadily increasing in nucleotide distance to the 1918 root over more than a century. The 2009 influenza A/H1N1 pandemic represents a clear departure from this enduring genetic backbone. Utilizing nucleotide distance maps and phylogenetic analyses, we illustrate remaining uncertainties regarding the origin of the 2009 pandemic, highlighting the complexity of influenza evolution. The NS segment is interesting precisely because it experiences less pervasive positive selection, and departs less strongly from neutral evolution than e.g. the HA antigen. Consequently, sudden deviations from neutral diversification can indicate changes in other genes via the hitchhiking effect. Our approach employs two measures based on nucleotide mismatch counts to analyze the evolutionary dynamics of the NS gene segment. The rooted Hamming map of distances between a reference sequence and all other sequences over time, and the unrooted temporal Hamming distribution which captures the distribution of genotypic distances between simultaneously circulating viruses, thereby revealing patterns of nucleotide diversity and epi-evolutionary dynamics.

Shared alleles and genetic structures in different Thai domestic cat breeds: the possible influence of common racial origins

Over hundreds of years, cats have been domesticated and selectively bred, resulting in numerous pedigreed breeds expedited by recent cat shows and breeding associations. Concerns have been raised about the limited breeding options and the genetic implications of inbreeding, indicating challenges in maintaining genetic diversity and accurate identification in purebred cats. In this study, genetic variability and structure were examined in 5 Thai domestic cat breeds using 15 microsatellite markers and mitochondrial DNA (mtDNA) D-loop sequencing. In total, 184 samples representing the Wichien Maat (WCM), Suphalak (SL), Khao-Manee (KM), Korat (KR), and Konja (KJ) breeds were analyzed. High genetic diversity (Ho and He > 0.5) was observed in all breeds, and mtDNA analysis revealed two primary haplogroups (A and B) that were shared among all domestic cat breeds in Thailand and globally. However, minor differences were observed between Thai domestic cat breeds based on clustering analyses, in which a distinct genetic structure was observed in the WCM breed. This suggests that allele fixation for distinctive morphological traits has occurred in Thai domestic cat breeds that emerged in isolated regions with shared racial origins. Analysis of relationships among individuals within the breed revealed high identification efficiency in Thai domestic cat breeds (P(ID)sibs < 10-4). Additionally, diverse and effective individual identification can be ensured by optimizing marker efficiency by using only nine loci. This comprehensive genetic characterization provides valuable insights into conservation strategies and breeding practices for Thai domestic cat breeds.

Investigating linguistic and genetic shifts in East Indian tribal groups

South Asia is home to almost a quarter of the world's total population and is home to significant ethnolinguistic diversity. Previous studies of linguistic and genetic affiliations of Indian populations suggest that the formation of these distinct groups was a protracted and complex phenomenon involving multiple waves of migration, cultural assimilation, and genetic admixture. The evolutionary processes of migration, mixing and merging of populations thus impact the culture and linguistic diversity of different groups, some of which may retain their linguistic affinities despite genetic admixture with other groups, or vice versa. Our study examines the relationship of genetic and linguistic affinities between Austroasiatic and Indo-European speakers in adjacent geographical regions of Eastern India. We analyzed 224 mitogenomes and 0.65 million SNP genotypes from 40 unrelated individuals belonging to the Bathudi, Bhumij, Ho, and Mahali ethnic groups from the Eastern Indian state of Odisha. These four groups are speakers of Austroasiatic languages who have adopted elements from Indo-European languages spoken in neighbouring regions. Our results suggest that these groups have the greatest maternal genetic affinity with other Austroasiatic-speaking groups in India. Allele frequency-based analyses, genome-wide SNPs, haplotype-based methods and IBD sharing further support the genetic similarity of these East Indian groups to Austroasiatic speakers of South Asia rather than regional populations speaking Indo-European and Dravidian languages. Our study shows that these populations experienced linguistic mixing, likely due to industrialization and modernization that brought them into close cultural contact with neighbouring Indo-European-speaking groups. However, linguistic change in these groups is not reflected in genetic mixing in these populations, as they appear to maintain strict genetic boundaries while simultaneously experiencing cultural mixing.

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.

Phylogeny and disease associations of a widespread and ancient intestinal bacteriophage lineage

Viruses are core components of the human microbiome, impacting health through interactions with gut bacteria and the immune system. Most human microbiome viruses are bacteriophages, which exclusively infect bacteria. Until recently, most gut virome studies focused on low taxonomic resolution (e.g., viral operational taxonomic units), hampering population-level analyses. We previously identified an expansive and widespread bacteriophage lineage in inhabitants of Amsterdam, the Netherlands. Here, we study their biodiversity and evolution in various human populations. Based on a phylogeny using sequences from six viral genome databases, we propose the Candidatus order Heliusvirales. We identify heliusviruses in 82% of 5441 individuals across 39 studies, and in nine metagenomes from humans that lived in Europe and North America between 1000 and 5000 years ago. We show that a large lineage started to diversify when Homo sapiens first appeared some 300,000 years ago. Ancient peoples and modern hunter-gatherers have distinct Ca. Heliusvirales populations with lower richness than modern urbanized people. Urbanized people suffering from type 1 and type 2 diabetes, as well as inflammatory bowel disease, have higher Ca. Heliusvirales richness than healthy controls. We thus conclude that these ancient core members of the human gut virome have thrived with increasingly westernized lifestyles.

Convergent Adaptation of True Crabs (Decapoda: Brachyura) to a Gradient of Terrestrial Environments

For much of terrestrial biodiversity, the evolutionary pathways of adaptation from marine ancestors are poorly understood and have usually been viewed as a binary trait. True crabs, the decapod crustacean infraorder Brachyura, comprise over 7600 species representing a striking diversity of morphology and ecology, including repeated adaptation to non-marine habitats. Here, we reconstruct the evolutionary history of Brachyura using new and published sequences of 10 genes for 344 tips spanning 88 of 109 brachyuran families. Using 36 newly vetted fossil calibrations, we infer that brachyurans most likely diverged in the Triassic, with family-level splits in the late Cretaceous and early Paleogene. By contrast, the root age is underestimated with automated sampling of 328 fossil occurrences explicitly incorporated into the tree prior, suggesting such models are a poor fit under heterogeneous fossil preservation. We apply recently defined trait-by-environment associations to classify a gradient of transitions from marine to terrestrial lifestyles. We estimate that crabs left the marine environment at least 7 and up to 17 times convergently, and returned to the sea from non-marine environments at least twice. Although the most highly terrestrial- and many freshwater-adapted crabs are concentrated in Thoracotremata, Bayesian threshold models of ancestral state reconstruction fail to identify shifts to higher terrestrial grades due to the degree of underlying change required. Lineages throughout our tree inhabit intertidal and marginal marine environments, corroborating the inference that the early stages of terrestrial adaptation have a lower threshold to evolve. Our framework and extensive new fossil and natural history datasets will enable future comparisons of non-marine adaptation at the morphological and molecular level. Crabs provide an important window into the early processes of adaptation to novel environments, and different degrees of evolutionary constraint that might help predict these pathways. [Brachyura; convergent evolution; crustaceans; divergence times; fossil calibration; molecular phylogeny; terrestrialization; threshold model.].

Insights from the first phylogenomic analysis of flat wasps (Hymenoptera, Bethylidae) reveal two new subfamilies

Despite significant advances in alpha level taxonomy in the past few decades, the higher-level phylogeny of flat wasps (Hymenoptera, Bethylidae) remains poorly explored. Herein we provide the first phylogenomic assessment of the family based on data from ultraconserved elements for 96 species in 61 genera of the family, with material from 29 countries and all biogeographic regions. Dataset cutoffs including ultraconserved element loci recovered for 50 and 70% of terminals resulted in matrices with 1513 and 451 loci, which were analysed in both parsimony and maximum likelihood frameworks. We also provide the first analyses of divergence dating for the family based on the calibration of 12 nodes. All analyses recovered the Bethylidae as a monophyletic group and estimated the origin of the family at around 143 Ma. The results suggest that all extant subfamilies had already diversified by the Late Cretaceous. All topologies suggest that Glenosema and Chilepyris form a clade separate from other Scleroderminae; owing to the morphological distinctiveness of each genus, we propose that they are accommodated in two new subfamilies, Glenoseminae subf. nov. and Chilepyrinae subf. nov. The monotypic sclerodermine genus Galodoxa was consistently recovered within Epyrinae and is transferred to the latter subfamily.

Phylogenetic and Comparative Genomics Study of Papilionidae Based on Mitochondrial Genomes

Most species of Papilionidae are large and beautiful ornamental butterflies. They are recognized as model organisms in ecology, evolutionary biology, genetics, and conservation biology but present numerous unresolved phylogenetic problems. Complete mitochondrial genomes (mitogenomes) have been widely used in phylogenetic studies of butterflies, but mitogenome knowledge within the family Papilionidae is limited, and its phylogeny is far from resolved. In this study, we first report the mitogenome of Byasa confusa from the subfamily Papilioninae of Papilionidae. The mitogenome of B. confusa is 15,135 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and an AT-rich control region (CR), closely mirroring the genomic structure observed in related butterfly species. Comparative analysis of 77 Papilionidae mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, Ka/Ks, and relative synonymous codon usage (RSCU) are all consistent with that of other reported butterfly mitogenomes. We conducted phylogenetic analyses using maximum-likelihood (ML) and Bayesian-inference (BI) methods, with 77 Papilionidae species as ingroups and two species of Nymphalidae and Lycaenidae as outgroups. The phylogenetic analysis indicated that B. confusa were clustered within Byasa. The phylogenetic trees show the monophyly of the subfamily Papilioninae and the tribes Leptocircini, Papilionini, and Troidini. The data supported the following relationships in tribe level on Papilioninae: (((Troidini + Papilionini) + Teinopalpini) + Leptocircini). The divergence time analysis suggests that Papilionidae originated in the late Creataceous. Overall, utilizing the largest number of Papilionidae mitogenomes sequenced to date, with the current first exploration in a phylogenetic analysis on Papilionidae (including four subfamilies), this study comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of the Papilionidae family.

Global Update on Measles Molecular Epidemiology

Molecular surveillance of circulating measles variants serves as a line of evidence for the absence of endemic circulation and provides a means to track chains of transmission. Molecular surveillance for measles (genotyping) is based on the sequence of 450 nucleotides at the end of the nucleoprotein coding region (N450) of the measles genome. Genotyping was established in 1998 and, with over 50,000 sequence submissions to the Measles Nucleotide Surveillance database, has proven to be an effective resource for countries attempting to trace pathways of transmission. This review summarizes the tools used for the molecular surveillance of measles and describes the challenge posed by the decreased number of circulating measles genotypes. The Global Measles and Rubella Laboratory Network addressed this challenge through the development of new tools such as named strains and distinct sequence identifiers that analyze the diversity within the currently circulating genotypes. The advantages and limitations of these approaches are discussed, together with the need to generate additional sequence data including whole genome sequences to ensure the continued utility of strain surveillance for measles.

Waves of Colonization and Gene Flow in a Great Speciator

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

Stomatal dynamics in Alloteropsis semialata arise from the evolving interplay between photosynthetic physiology, stomatal size and biochemistry

C4 plants are expected to have faster stomatal movements than C3 species because they tend to have smaller guard cells. However, little is known about how the evolution of C4 photosynthesis influences stomatal dynamics in relation to guard cell size and environmental factors. We studied photosynthetically diverse populations of the grass Alloteropsis semialata, showing that the origin of C4 photosynthesis in this species was associated with a shortening of stomatal guard and subsidiary cells. However, for a given cell size, C4 and C3-C4 intermediate individuals had similar or slower light-induced stomatal opening speeds than C3 individuals. Conversely, when exposed to decreasing light, stomata in C4 plants closed as fast as those in non-C4 plants. Polyploid formation in some C4 plants led to larger stomatal cells and was associated with slower stomatal opening. Conversely, diversification of C4 diploid plants into wetter environments was associated with an acceleration of stomatal opening. Overall, there was significant relationship between light-saturated photosynthesis and stomatal opening speed in the C4 plants, implying that photosynthetic energy production was limiting for stomatal opening. Stomatal dynamics in this wild grass therefore arise from the evolving interplay between photosynthetic physiology and the size and biochemical function of stomatal complexes.

An 8000 years old genome reveals the Neolithic origin of the zoonosis Brucella melitensis

Brucella melitensis is a major livestock bacterial pathogen and zoonosis, causing disease and infection-related abortions in small ruminants and humans. A considerable burden to animal-based economies today, the presence of Brucella in Neolithic pastoral communities has been hypothesised but we lack direct genomic evidence thus far. We report a 3.45X B. melitensis genome preserved in an ~8000 year old sheep specimen from Menteşe Höyük, Northwest Türkiye, demonstrating that the pathogen had evolved and was circulating in Neolithic livestock. The genome is basal with respect to all known B. melitensis and allows the calibration of the B. melitensis speciation time from the primarily cattle-infecting B. abortus to approximately 9800 years Before Present (BP), coinciding with a period of consolidation and dispersal of livestock economies. We use the basal genome to timestamp evolutionary events in B. melitensis, including pseudogenization events linked to erythritol response, the supposed determinant of the pathogen's placental tropism in goats and sheep. Our data suggest that the development of herd management and multi-species livestock economies in the 11th-9th millennium BP drove speciation and host adaptation of this zoonotic pathogen.

Naming the other cousin: A new goldie barb (Cyprinidae: Smiliogastrininae) from the northeast escarpment in South Africa, with proposed taxonomic rearrangement of the goldie barb group in southern Africa

A growing body of evidence indicates that the global diversity of freshwater fishes has not been fully documented. Studies of freshwater fishes that were previously thought to be morphologically variable have revealed the existence of deeply divergent lineages, with many distinct species. In southern Africa a number of Enteromius species exhibit either exceedingly wide or divided distribution patterns that should be rare for freshwater fishes with limited dispersal opportunities between river systems. One such species is the sidespot barb, Enteromius neefi. As currently defined, E. neefi has a disjunct distribution that is divided between rivers in the northeast escarpment in South Africa and Eswatini, and tributaries of the Upper Zambezi in Zambia and southern Congo in the Democratic Republic of Congo, with a large geographic gap between these two populations. With the use of molecular and morphological methods, the level of divergence between the two populations was examined, and a new species was described from the Steelpoort River in the Limpopo River system of South Africa. Findings from this study provide further evidence for a number of taxonomic problems within the goldie barbs of southern Africa, and some taxonomic rearrangements are proposed for this group.

Spatial and temporal transmission dynamics of respiratory syncytial virus in New Zealand before and after the COVID-19 pandemic

Human respiratory syncytial virus (RSV) is a major cause of acute respiratory infection. In 2020, RSV was effectively eliminated from the community in New Zealand due to non-pharmaceutical interventions (NPI) used to control the spread of COVID-19. However, in April 2021, following a brief quarantine-free travel agreement with Australia, there was a large-scale nationwide outbreak of RSV that led to reported cases more than five times higher, and hospitalisations more than three times higher, than the typical seasonal pattern. In this study, we generated 1,471 viral genomes of both RSV-A and RSV-B sampled between 2015 and 2022 from across New Zealand. Using a phylodynamics approach, we used these data to better understand RSV transmission patterns in New Zealand prior to 2020, and how RSV became re-established in the community following the relaxation of COVID-19 restrictions. We found that in 2021, there was a large epidemic of RSV in New Zealand that affected a broader age group range compared to the usual pattern of RSV infections. This epidemic was due to an increase in RSV importations, leading to several large genomic clusters of both RSV-A ON1 and RSV-B BA9 genotypes in New Zealand. However, while a number of importations were detected, there was also a major reduction in RSV genetic diversity compared to pre-pandemic seasonal outbreaks. These genomic clusters were temporally associated with the increase of migration in 2021 due to quarantine-free travel from Australia at the time. The closest genetic relatives to the New Zealand RSV genomes, when sampled, were viral genomes sampled in Australia during a large, off-season summer outbreak several months prior, rather than cryptic lineages that were sustained but not detected in New Zealand. These data reveal the impact of NPI used during the COVID-19 pandemic on other respiratory infections and highlight the important insights that can be gained from viral genomes.

Early mutational signatures and transmissibility of SARS-CoV-2 Gamma and Lambda variants in Chile

Genomic surveillance (GS) programmes were crucial in identifying and quantifying the mutating patterns of SARS-CoV-2 during the COVID-19 pandemic. In this work, we develop a Bayesian framework to quantify the relative transmissibility of different variants tailored for regions with limited GS. We use it to study the relative transmissibility of SARS-CoV-2 variants in Chile. Among the 3443 SARS-CoV-2 genomes collected between January and June 2021, where sampling was designed to be representative, the Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were predominant. We found that Lambda and Gamma variants' reproduction numbers were 5% (95% CI: [1%, 14%]) and 16% (95% CI: [11%, 21%]) larger than Alpha's, respectively. Besides, we observed a systematic mutation enrichment in the Spike gene for all circulating variants, which strongly correlated with variants' transmissibility during the studied period (r = 0.93, p-value = 0.025). We also characterised the mutational signatures of local samples and their evolution over time and with the progress of vaccination, comparing them with those of samples collected in other regions worldwide. Altogether, our work provides a reliable method for quantifying variant transmissibility under subsampling and emphasises the importance of continuous genomic surveillance.

Characterizing the Phan Rang Sheep: A First Look at the Y Chromosome, Mitochondrial DNA, and Morphometrics

The Phan Rang sheep, considered the sole indigenous breed of Vietnam, are primarily concentrated in the two central provinces of Ninh Thuan and Binh Thuan, with Ninh Thuan accounting for more than 90% of the country's sheep population. These provinces are known for their high temperatures and frequent droughts. The long-standing presence of the Phan Rang sheep in these regions suggests their potential resilience to heat stress-a trait of increasing interest in the face of global climate change. Despite the breed's significance, a critical knowledge gap hinders conservation and breeding programs. To address this, our study employed a two-pronged approach. First, we collected body conformational data to aid in breed identification. Second, we analyzed mitochondrial DNA (D-loop) and Y chromosome markers (SRY and SRYM18) to elucidate the maternal and paternal lineages. Among the 68 Phan Rang sheep analyzed for their D-loop, 19 belonged to mitochondrial haplogroup A, while 49 belonged to haplogroup B. The haplogroups can be subdivided into 16 unique haplotypes. All 19 rams surveyed for their paternal lineages belonged to haplotypes H5 and H6. These findings strongly support the hypothesis of dual origins for the Phan Rang sheep. This study presents the first genetic data for the Phan Rang breed, providing crucial insights for future research and conservation efforts.

Phylogenomic perspectives on speciation and reproductive isolation in a North American biodiversity hotspot: an example using California sages (Salvia subgenus Audibertia: Lamiaceae)

BACKGROUND AND AIMS

The California Floristic Province (CA-FP) is the most species-rich region of North America north of Mexico. One of several proposed hypotheses explaining the exceptional diversity of the region is that the CA-FP harbours myriad recently diverged lineages with nascent reproductive barriers. Salvia subgenus Audibertia is a conspicuous element of the CA-FP, with multiple sympatric and compatible species.

METHODS

Using 305 nuclear loci and both organellar genomes, we reconstruct species trees, examine genomic discordance, conduct divergence-time estimation, and analyse contemporaneous patterns of gene flow and mechanical reproductive isolation.

KEY RESULTS

Despite strong genomic discordance, an underlying bifurcating tree is supported. Organellar genomes capture additional introgression events not detected in the nuclear genome. Most interfertility is found within clades, indicating that reproductive barriers arise with increasing genetic divergence. Species are generally not mechanically isolated, suggesting that it is unlikely to be the primary factor leading to reproductive isolation.

CONCLUSIONS

Rapid, recent speciation with some interspecific gene flow in conjunction with the onset of a Mediterranean-like climate is the underlying cause of extant diversity in Salvia subgenus Audibertia. Speciation has largely not been facilitated by gene flow. Its signal in the nuclear genome seems to mostly be erased by backcrossing, but organellar genomes each capture different instances of historical gene flow, probably characteristic of many CA-FP lineages. Mechanical reproductive isolation appears to be only part of a mosaic of factors limiting gene flow.

A new semi-fossorial thescelosaurine dinosaur from the Cenomanian-age Mussentuchit Member of the Cedar Mountain Formation, Utah

Thescelosaurines are a group of early diverging, ornithischian dinosaurs notable for their conservative bauplans and mosaic of primitive features. Although abundant within the latest Cretaceous ecosystems of North America, their record is poor to absent in earlier assemblages, leaving a large gap in our understanding of their evolution, origins, and ecological roles. Here we report a new small bodied thescelosaurine-Fona herzogae gen. et sp. nov.-from the Mussentuchit Member of the Cedar Mountain Formation, Utah, USA. Fona herzogae is represented by multiple individuals, representing one of the most comprehensive skeletal assemblages of a small bodied, early diverging ornithischian described from North America to date. Phylogenetic analysis recovers Fona as the earliest member of Thescelosaurinae, minimally containing Oryctodromeus, and all three species of Thescelosaurus, revealing the clade was well-established in North America by as early as the Cenomanian, and distinct from, yet continental cohabitants with, their sister clade, Orodrominae. To date, orodromines and thescelosaurines have not been found together within a single North American ecosystem, suggesting different habitat preferences or competitive exclusion. Osteological observations reveal extensive intraspecific variation across cranial and postcranial elements, and a number of anatomical similarities with Oryctodromeus, suggesting a shared semi-fossorial lifestyle.

Compound osteoderms preserved in amber reveal the oldest known skink

Scincidae is one of the most species-rich and cosmopolitan clades of squamate reptiles. Abundant disarticulated fossil material has also been attributed to this group, however, no complete pre-Cenozoic crown-scincid specimens have been found. A specimen in Burmite (99 MYA) is the first fossil that can be unambiguously referred to this clade. Our analyses place it as nested within extant skinks, supported by the presence of compound osteoderms formed by articulated small ostedermites. The specimen has a combination of dorsal and ventral compound osteoderms and overlapping cycloid scales that is limited to skinks. We propose that this type of osteoderm evolved as a response to an increased overlap of scales, and to reduced stiffness of the dermal armour. Compound osteoderms could be a key innovation that facilitated diversification in this megadiverse family.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

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.

Comparing Phylogeographies to Reveal Incompatible Geographical Histories within Genomes

Modern phylogeography aims at reconstructing the geographic movement of organisms based on their genomic sequences and spatial information. Phylogeographic approaches are often applied to pathogen sequences and therefore tend to neglect the possibility of recombination, which decouples the evolutionary and geographic histories of different parts of the genome. Genomic regions of recombining or reassorting pathogens often originate and evolve at different times and locations, which characterize their unique spatial histories. Measuring the extent of these differences requires new methods to compare geographic information on phylogenetic trees reconstructed from different parts of the genome. Here we develop for the first time a set of measures of phylogeographic incompatibility, aimed at detecting differences between geographical histories in terms of distances between phylogeographies. We study the effect of varying demography and recombination on phylogeographic incompatibilities using coalescent simulations. We further apply these measures to the evolutionary history of human and livestock pathogens, either reassorting or recombining, such as the Victoria and Yamagata lineages of influenza B and the O/Ind-2001 foot-and-mouth disease virus strain. Our results reveal diverse geographical paths of migration that characterize the origins and evolutionary histories of different viral genes and genomic segments. These incompatibility measures can be applied to any phylogeography, and more generally to any phylogeny where each tip has been assigned either a continuous or discrete "trait" independent of the sequence. We illustrate this flexibility with an analysis of the interplay between the phylogeography and phylolinguistics of Uralic-speaking human populations, hinting at patrilinear language transmission.

Population structure and history of Mycobacterium bovis European 3 clonal complex reveal transmission across ecological corridors of unrecognized importance in Portugal

Mycobacterium bovis causes animal tuberculosis in livestock and wildlife, with an impact on animal health and production, wildlife management, and public health. In this work, we sampled a multi-host tuberculosis community from the official hotspot risk area of Portugal over 16 years, generating the largest available data set in the country. Using phylogenetic and ecological modeling, we aimed to reconstruct the history of circulating lineages across the livestock-wildlife interface to inform intervention and the implementation of genomic surveillance within the official eradication plan. We find evidence for the co-circulation of M. bovis European 1 (Eu1), Eu2, and Eu3 clonal complexes, with Eu3 providing sufficient temporal signal for further phylogenetic investigation. The Eu3 most recent common ancestor (bovine) was dated in the 1990s, subsequently transitioning to wildlife (red deer and wild boar). Isolate clustering based on sample metadata was used to inform phylogenetic inference, unravelng frequent transmission between two clusters that represent an ecological corridor of previously unrecognized importance in Portugal. The latter was associated with transmission at the livestock-wildlife interface toward locations with higher temperature and precipitation, lower agriculture and road density, and lower host densities. This is the first analysis of M. bovis Eu3 complex in Iberia, shedding light on background ecological factors underlying long-term transmission and informing where efforts could be focused within the larger hotspot risk area of Portugal.

IMPORTANCE

Efforts to strengthen surveillance and control of animal tuberculosis (TB) are ongoing worlwide. Here, we developed an eco-phylodynamic framework based on discrete phylogenetic approaches informed by M. bovis whole-genome sequence data representing a multi-host transmission system at the livestock-wildlife interface, within a rich ecological landscape in Portugal, to understand transmission processes and translate this knowledge into disease management benefits. We find evidence for the co-circulation of several M. bovis clades, with frequent transmission of the Eu3 lineage among cattle and wildlife populations. Most transition events between different ecological settings took place toward host, climate and land use gradients, underscoring animal TB expansion and a potential corridor of unrecognized importance for M. bovis maintenance. Results stress that animal TB is an established wildlife disease without ecological barriers, showing that control measures in place are insufficient to prevent long-distance transmission and spillover across multi-host communities, demanding new interventions targeting livestock-wildlife interactions.

Genetic Insights into the Historical Attribution of Variety Names of Sweet Chestnut (Castanea sativa Mill.) in Northern Italy

The sweet chestnut (Castanea sativa Mill.) is subject to the progressive disappearance of its traditional chestnut groves. In the northern part of Italy, where distribution of the sweet chestnut is fragmented, many local varieties continue to be identified mostly by oral tradition. We characterised by SSRs eleven historically recognised varieties of sweet chestnut in the area surrounding Lake Como, with the goal of giving a genetic basis to the traditional classification. We performed classical analysis about differentiation and used Bayesian approaches to detect population structure and to reconstruct demography. The results revealed that historical and genetic classifications are loosely linked when chestnut fruits are just "castagne", that is, normal fruits, but increasingly overlap where "marroni" (the most prized fruits) are concerned. Bayesian classification allowed us to identify a homogeneous gene cluster not recognised in the traditional assessment of the varieties and to reconstruct possible routes used for the propagation of sweet chestnut. We also reconstructed ancestral relationships between the different gene pools involved and dated ancestral lineages whose results fit with palynological data. We suggest that conservation strategies based on a genetic evaluation of the resource should also rely on traditional cultural heritage, which could reveal new sources of germplasm.

Exploring selection signatures in the divergence and evolution of lipid droplet (LD) associated genes in major oilseed crops

BACKGROUND

Oil bodies or lipid droplets (LDs) in the cytosol are the subcellular storage compartments of seeds and the sites of lipid metabolism providing energy to the germinating seeds. Major LD-associated proteins are lipoxygenases, phospholipaseD, oleosins, TAG-lipases, steroleosins, caleosins and SEIPINs; involved in facilitating germination and enhancing peroxidation resulting in off-flavours. However, how natural selection is balancing contradictory processes in lipid-rich seeds remains evasive. The present study was aimed at the prediction of selection signatures among orthologous clades in major oilseeds and the correlation of selection effect with gene expression.

RESULTS

The LD-associated genes from the major oil-bearing crops were analyzed to predict natural selection signatures in phylogenetically close-knit ortholog clusters to understand adaptive evolution. Positive selection was the major force driving the evolution and diversification of orthologs in a lineage-specific manner. Significant positive selection effects were found in 94 genes particularly in oleosin and TAG-lipases, purifying with excess of non-synonymous substitution in 44 genes while 35 genes were neutral to selection effects. No significant selection impact was noticed in Brassicaceae as against LOX genes of oil palm. A heavy load of deleterious mutations affecting selection signatures was detected in T-lineage oleosins and LOX genes of Arachis hypogaea. The T-lineage oleosin genes were involved in mainly anther, tapetum and anther wall morphogenesis. In Ricinus communis and Sesamum indicum > 85% of PLD genes were under selection whereas selection pressures were low in Brassica juncea and Helianthus annuus. Steroleosin, caleosin and SEIPINs with large roles in lipid droplet organization expressed mostly in seeds and were under considerable positive selection pressures. Expression divergence was evident among paralogs and homeologs with one gene attaining functional superiority compared to the other. The LOX gene Glyma.13g347500 associated with off-flavor was not expressed during germination, rather its paralog Glyma.13g347600 showed expression in Glycine max. PLD-α genes were expressed on all the tissues except the seed,δ genes in seed and meristem while β and γ genes expressed in the leaf.

CONCLUSIONS

The genes involved in seed germination and lipid metabolism were under strong positive selection, although species differences were discernable. The present study identifies suitable candidate genes enhancing seed oil content and germination wherein directional selection can become more fruitful.

Disentangling the drivers of urban bird diversity in the non-breeding season: A general synthesis

Current knowledge about the impacts of urbanisation on bird assemblages is based on evidence from studies partly or wholly undertaken in the breeding season. In comparison, the non-breeding season remains little studied, despite the fact that winter conditions at higher latitudes are changing more rapidly than other seasons. During the non-breeding season, cities may attract or retain bird species because they offer milder conditions or better feeding opportunities than surrounding habitats. However, the range of climatic, ecological and anthropogenic mechanisms shaping different facets of urban bird diversity in the non-breeding season are poorly understood. We explored these mechanisms using structural equation modelling to assess how urbanisation affects the taxonomic, phylogenetic and functional diversity of avian assemblages sampled worldwide in the non-breeding season. We found that minimum temperature, elevation, urban area and city age played a critical role in determining taxonomic diversity while a range of factors-including productivity, precipitation, elevation, distance to coasts and rivers, socio-economic (as a proxy of human facilitation) and road density-each contributed to patterns of phylogenetic and functional diversity. The structure and function of urban bird assemblages appear to be predominantly shaped by temperature, productivity and city age, with effects of these factors differing across seasons. Our results underline the importance of considering multiple hypotheses, including seasonal effects, when evaluating the impacts of urbanisation on biodiversity.

Archival mitogenomes identify invasion by the Batrachochytrium dendrobatidis CAPE lineage caused an African amphibian extinction in the wild

Outbreaks of emerging infectious diseases are influenced by local biotic and abiotic factors, with host declines occurring when conditions favour the pathogen. Deterioration in the population of the micro-endemic Tanzanian Kihansi spray toad (Nectophrynoides asperginis) occurred after the construction of a hydropower dam, implicating habitat modification in this species decline. Population recovery followed habitat augmentation; however, a subsequent outbreak of chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd) led to the spray toad's extinction in the wild. We show using spatiotemporal surveillance and mitogenome assembly of Bd from archived toad mortalities that the outbreak was caused by invasion of the BdCAPE lineage and not the panzootic lineage BdGPL. Molecular dating reveals an emergence of BdCAPE across southern Africa overlapping with the timing of the spray toad's extinction. That our post-outbreak surveillance of co-occurring amphibian species in the Udzungwa Mountains shows widespread infection by BdCAPE yet no signs of ill-health or decline suggests these other species can tolerate Bd when environments are stable. We conclude that, despite transient success in mitigating the impact caused by dams' construction, invasion by BdCAPE caused the ultimate die-off that led to the extinction of the Kihansi spray toad.

Time-calibrated phylogenetic and chromosomal mobilome analyses of Staphylococcus aureus CC398 reveal geographical and host-related evolution

An international collection of Staphylococcus aureus of clonal complex (CC) 398 from diverse hosts spanning all continents and a 30 year-period is studied based on whole-genome sequencing (WGS) data. The collection consists of publicly available genomic data from 2994 strains and 134 recently sequenced Swiss methicillin-resistant S. aureus (MRSA) CC398 strains. A time-calibrated phylogeny reveals the presence of distinct phylogroups present in Asia, North and South America and Europe. European MRSA diverged from methicillin-susceptible S. aureus (MSSA) at the beginning of the 1950s. Two major European phylogroups (EP4 and EP5), which diverged approximately 1974, are the main drivers of MRSA CC398 spread in Europe. Within EP5, an emergent MRSA lineage spreading among the European horse population (EP5-Leq) diverged approximately 1996 from the pig lineage (EP5-Lpg), and also contains human-related strains. EP5-Leq is characterized by staphylococcal cassette chromosome mec (SCCmec) IVa and spa type t011 (CC398-IVa-t011), and EP5-Lpg by CC398-SCCmecVc-t011. The lineage-specific antibiotic resistance and virulence gene patterns are mostly mediated by the acquisition of mobile genetic elements like SCCmec, S. aureus Genomic Islands (SaGIs), prophages and transposons. Different combinations of virulence factors are present on S. aureus pathogenicity islands (SaPIs), and novel antimicrobial resistance gene containing elements are associated with certain lineages expanding in Europe. This WGS-based analysis reveals the actual evolutionary trajectory and epidemiological trend of the international MRSA CC398 population considering host, temporal, geographical and molecular factors. It provides a baseline for global WGS-based One-Health studies of adaptive evolution of MRSA CC398 as well as for local outbreak investigations.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

CDC48 in plants and its emerging function in plant immunity

Protein homeostasis, namely the balance between protein synthesis and degradation, must be finely controlled to ensure cell survival, notably through the ubiquitin-proteasome system (UPS). In all species, including plants, homeostasis is disrupted by biotic and abiotic stresses. A key player in the maintenance of protein balance, the protein CDC48, shows emerging functions in plants, particularly in response to biotic stress. In this review on CDC48 in plants, we detail its highly conserved structure, describe a gene expansion that is only present in Viridiplantae, discuss its various functions and regulations, and finally highlight its recruitment, still not clear, during the plant immune response.

Phylogeography of Morella nana: The Wumeng Mountains as a natural geographical isolation boundary on the Yunnan-Guizhou Plateau

The Yunnan-Guizhou Plateau (YGP) is characterized by the distinctive isolated habitat of the limestone Karst Islands and features the Wumeng Mountains, which divide the YGP into the two Plateaus of Yunnan and Guizhou. This study aimed to assess the effects of geographic isolation and past climate fluctuation on the distribution of flora in the YGP. To achieve this, we carried out the phylogeographical pattern and genetic structure based on chloroplast and nuclear ribosomal DNA sequence in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modeling for Morella nana, an important wild plant resource and endemic to the YGP once considered a vulnerable species. The results suggest that the genetic and chlorotype network structures of M. nana are divided into at least two groups: cpDNA chlorotype H2 (or dominant nrDNA haplotypes h1 and h2), distributed primarily to the east of the Wumeng Mountains, and cpDNA chlorotypes H1 and H3-H10 (or dominant nrDNA haplotype h2 and h3), distributed to the west of the Wumeng Mountains. A deep genetic split was noted within the two groups to reach 25 steps, especially for the cpDNA fragment variation. This east-west divergence reveals the existence of a natural geographical isolation boundary in the form of the Wumeng Mountains, and supports the existence of at least two glacial refuges during the Quaternary glacial period, along with two genetic diversity center, and at least two large geographic protection units for the important species of M. nana. This study indicates that the phylogeographical pattern of M. nana can be attributed to geographic/environmental isolation caused by the Wumeng Mountains and climate fluctuation during the last glacial maximum, and proposes an effective strategy to protecting this important plant resource.

Evolutionary déjà vu? A case of convergent evolution in an ant-plant association

Obligatory ant-plant symbioses often appear to be single evolutionary shifts within particular ant lineages; however, convergence can be revealed once natural history observations are complemented with molecular phylogenetics. Here, we describe a remarkable example of convergent evolution in an ant-plant symbiotic system. Exclusively arboreal, Myrmelachista species can be generalized opportunists nesting in several plant species or obligately symbiotic, live-stem nesters of a narrow set of plant species. Instances of specialization within Myrmelachista are known from northern South America and throughout Middle America. In Middle America, a diverse radiation of specialists occupies understory treelets of lowland rainforests. The morphological and behavioural uniformity of specialists suggests that they form a monophyletic assemblage, diversifying after a single origin of specialization. Using ultraconserved element phylogenomics and ancestral state reconstructions, we show that shifts from opportunistic to obligately symbiotic evolved independently in South and Middle America. Furthermore, our analyses support a remarkable case of convergence within the Middle American radiation, with two independently evolved specialist clades, arising nearly simultaneously from putative opportunistic ancestors during the late Pliocene. This repeated evolution of a complex phenotype suggests similar mechanisms behind trait shifts from opportunists to specialists, generating further questions about the selective forces driving specialization.

Rapid identification and subsequent contextualization of an outbreak of methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit using nanopore sequencing

Outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) are well described in the neonatal intensive care unit (NICU) setting. Genomics has revolutionized the investigation of such outbreaks; however, to date, this has largely been completed retrospectively and has typically relied on short-read platforms. In 2022, our laboratory established a prospective genomic surveillance system using Oxford Nanopore Technologies sequencing for rapid outbreak detection. Herein, using this system, we describe the detection and control of an outbreak of sequence-type (ST)97 MRSA in our NICU. The outbreak was identified 13 days after the first MRSA-positive culture and at a point where there were only two known cases. Ward screening rapidly defined the extent of the outbreak, with six other infants found to be colonized. There was minimal transmission once the outbreak had been detected and appropriate infection control measures had been instituted; only two further ST97 cases were detected, along with three unrelated non-ST97 MRSA cases. To contextualize the outbreak, core-genome single-nucleotide variants were identified for phylogenetic analysis after de novo assembly of nanopore data. Comparisons with global (n=45) and national surveillance (n=35) ST97 genomes revealed the stepwise evolution of methicillin resistance within this ST97 subset. A distinct cluster comprising nine of the ten ST97-IVa genomes from the NICU was identified, with strains from 2020 to 2022 national surveillance serving as outgroups to this cluster. One ST97-IVa genome presumed to be part of the outbreak formed an outgroup and was retrospectively excluded. A second phylogeny was created using Illumina sequencing, which considerably reduced the branch lengths of the NICU isolates on the phylogenetic tree. However, the overall tree topology and conclusions were unchanged, with the exception of the NICU outbreak cluster, where differences in branch lengths were observed. This analysis demonstrated the ability of a nanopore-only prospective genomic surveillance system to rapidly identify and contextualize an outbreak of MRSA in a NICU.