Adaptive dating and fast proposals: Revisiting the phylogenetic relaxed clock model

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

Phylogenomics and biogeography of the small carpenter bees (Apidae: Xylocopinae: Ceratina)

Small carpenter bees in the genus Ceratina are behaviourally diverse, species-rich, and cosmopolitan, with over 370 species and a range including all continents except Antarctica. Here, we present the first comprehensive phylogeny of the genus based on ultraconserved element (UCE) phylogenomic data, covering a total of 185 ingroup specimens representing 22 of the 25 current subgenera. Our results support most recognized subgenera as natural groups, but we also highlight several groups in need of taxonomic revision - particularly the nominate subgenus Ceratina sensu stricto - and several clades that likely need to be described as new subgenera. In addition to phylogeny, we explore the evolutionary history of Ceratina through divergence time estimation and biogeographic reconstruction. Our findings suggest that Ceratinini split from its sister tribe Allodapini about 72 million years ago. The common ancestor of Ceratina emerged in the Afrotropical realm approximately 42 million years ago, near the Middle Eocene Climatic Optimum. Multiple subsequent dispersal events led to the present cosmopolitan distribution of Ceratina, with the majority of transitions occurring between the Afrotropics, Indomalaya, and the Palearctic. Additional movements also led to the arrival of Ceratina in Madagascar, Australasia, and a single colonization of the Americas. Dispersal events were asymmetrical overall, with temperate regions primarily acting as destinations for migrations from tropical source regions.

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.

Transition of survival strategies under global climate shifts in the grape family

Faced with environmental changes, plants may either move to track their ancestral niches or evolve to adapt to new niches. Vitaceae, the grape family, has evolved diverse adaptive traits facilitating a global expansion in wide-ranging habitats, making it ideal for investigating transition between move and evolve strategies and exploring the underlying mechanisms. Here we inferred the patterns of biogeographic diversification and trait evolution in Vitaceae based on a robust phylogeny with dense sampling including 495 species (~52% of Vitaceae species). Vitaceae probably originated from Asia-the diversity centre of extant genera and the major source of dispersals. Boundaries of the Eocene, Oligocene and Miocene were identified as turning points in shifting strategies. A significant decrease in move strategy was identified during the Oligocene, followed by increases in move and evolve. After the Miocene, evolve began to dominate, during which increased niche opportunities and key trait innovations played important roles.

Plastid phylogenomics reveals evolutionary relationships in the mycoheterotrophic orchid genus Dipodium and provides insights into plastid gene degeneration

The orchid genus Dipodium R.Br. (Epidendroideae) comprises leafy autotrophic and leafless mycoheterotrophic species, with the latter confined to sect. Dipodium. This study examined plastome degeneration in Dipodium in a phylogenomic and temporal context. Whole plastomes were reconstructed and annotated for 24 Dipodium samples representing 14 species and two putatively new species, encompassing over 80% of species diversity in sect. Dipodium. Phylogenomic analysis based on 68 plastid loci including a broad outgroup sampling across Orchidaceae found that sect. Leopardanthus is the sister lineage to sect. Dipodium. Dipodium ensifolium, the only leafy autotrophic species in sect. Dipodium, was found to be a sister to all leafless, mycoheterotrophic species, supporting a single evolutionary origin of mycoheterotrophy in the genus. Divergence-time estimations found that Dipodium arose ca. 33.3 Ma near the lower boundary of the Oligocene and that crown diversification commenced in the late Miocene, ca. 11.3 Ma. Mycoheterotrophy in the genus was estimated to have evolved in the late Miocene, ca. 7.3 Ma, in sect. Dipodium. The comparative assessment of plastome structure and gene degradation in Dipodium revealed that plastid ndh genes were pseudogenised or physically lost in all Dipodium species, including in leafy autotrophic species of both Dipodium sections. Levels of plastid ndh gene degradation were found to vary among species as well as within species, providing evidence of relaxed selection for retention of the NADH dehydrogenase complex within the genus. Dipodium exhibits an early stage of plastid genome degradation, as all species were found to have retained a full set of functional photosynthesis-related genes and housekeeping genes. This study provides important insights into plastid genome degradation along the transition from autotrophy to mycoheterotrophy in a phylogenomic and temporal context.

Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

Central African dwarf crocodiles found in syntopy are comparably divergent to South American dwarf caimans

Recent molecular taxonomic advancements have expanded our understanding of crocodylian diversity, revealing the existence of previously overlooked species, including the Congo dwarf crocodile (Osteolaemus osborni) in the central Congo Basin rainforests. This study explores the genomic divergence between O. osborni and its better-known relative, the true dwarf crocodile (Osteolaemus tetraspis), shedding light on their evolutionary history. Field research conducted in the northwestern Republic of the Congo uncovered a locality where both species coexist in sympatry/syntopy. Genomic analysis of sympatric individuals reveals a level of divergence comparable to that between ecologically similar South American dwarf caimans (Paleosuchus palpebrosus and Paleosuchus trigonatus), suggesting parallel speciation in the Afrotropics and Neotropics during the Middle to Late Miocene, 10-12 Ma. Comparison of the sympatric and allopatric dwarf crocodiles indicates no gene flow between the analysed sympatric individuals of O. osborni and O. tetraspis. However, a larger sample will be required to answer the question of whether or to what extent these species hybridize. This study emphasizes the need for further research on the biology and conservation status of the Congo dwarf crocodile, highlighting its significance in the unique biodiversity of the Congolian rainforests and thus its potential as a flagship species.

Variational Bayesian phylogenies through matrix representation of tree space

In this article, we study the distance matrix as a representation of a phylogeny by way of hierarchical clustering. By defining a multivariate normal distribution on (a subset of) the entries in a matrix, this allows us to represent a distribution over rooted time trees. Here, we demonstrate tree distributions can be represented accurately this way for a number of published tree distributions. Though such a representation does not map to unique trees, restriction to a subspace, in particular one we call a "cube", makes the representation bijective at the cost of not being able to represent all possible trees. We introduce an algorithm "cubeVB" specifically for cubes and show through well calibrated simulation study that it is possible to recover parameters of interest like tree height and length. Although a cube cannot represent all of tree space, it is a great improvement over a single summary tree, and it opens up exciting new opportunities for scaling up Bayesian phylogenetic inference. We also demonstrate how to use a matrix representation of a tree distribution to get better summary trees than commonly used maximum clade credibility trees. An open source implementation of the cubeVB algorithm is available from https://github.com/rbouckaert/cubevb as the cubevb package for BEAST 2.

Multiple introductions of monkeypox virus to Ireland during the international mpox outbreak, May 2022 to October 2023

BackgroundMpox, caused by monkeypox virus (MPXV), was considered a rare zoonotic disease before May 2022, when a global epidemic of cases in non-endemic countries led to the declaration of a Public Health Emergency of International Concern. Cases of mpox in Ireland, a country without previous mpox reports, could reflect extended local transmission or multiple epidemiological introductions.AimTo elucidate the origins and molecular characteristics of MPXV circulating in Ireland between May 2022 and October 2023.MethodsWhole genome sequencing of MPXV from 75% of all Irish mpox cases (182/242) was performed and compared to sequences retrieved from public databases (n = 3,362). Bayesian approaches were used to infer divergence time between sequences from different subclades and evaluate putative importation events from other countries.ResultsOf 242 detected mpox cases, 99% were males (median age: 35 years; range: 15-60). All 182 analysed genomes were assigned to Clade IIb and, presence of 12 distinguishable subclades suggests multiple introductions into Ireland. Estimation of time to divergence of subclades further supports the hypothesis for multiple importation events from numerous countries, indicative of extended and sustained international spread of mpox. Further analysis of sequences revealed that 92% of nucleotide mutations were from cytosine to thymine (or from guanine to adenine), leading to a high number of non-synonymous mutations across subclades; mutations associated with tecovirimat resistance were not observed.ConclusionWe provide insights into the international transmission dynamics supporting multiple introductions of MPXV into Ireland. Such information supported the implementation of evidence-informed public health control measures.

Worldwide forest surveys reveal forty-three new species in Phytophthora major Clade 2 with fundamental implications for the evolution and biogeography of the genus and global plant biosecurity

During 25 surveys of global Phytophthora diversity, conducted between 1998 and 2020, 43 new species were detected in natural ecosystems and, occasionally, in nurseries and outplantings in Europe, Southeast and East Asia and the Americas. Based on a multigene phylogeny of nine nuclear and four mitochondrial gene regions they were assigned to five of the six known subclades, 2a-c, e and f, of Phytophthora major Clade 2 and the new subclade 2g. The evolutionary history of the Clade appears to have involved the pre-Gondwanan divergence of three extant subclades, 2c, 2e and 2f, all having disjunct natural distributions on separate continents and comprising species with a soilborne and aquatic lifestyle and, in addition, a few partially aerial species in Clade 2c; and the post-Gondwanan evolution of subclades 2a and 2g in Southeast/East Asia and 2b in South America, respectively, from their common ancestor. Species in Clade 2g are soilborne whereas Clade 2b comprises both soil-inhabiting and aerial species. Clade 2a has evolved further towards an aerial lifestyle comprising only species which are predominantly or partially airborne. Based on high nuclear heterozygosity levels ca. 38 % of the taxa in Clades 2a and 2b could be some form of hybrid, and the hybridity may be favoured by an A1/A2 breeding system and an aerial life style. Circumstantial evidence suggests the now 93 described species and informally designated taxa in Clade 2 result from both allopatric non-adaptive and sympatric adaptive radiations. They represent most morphological and physiological characters, breeding systems, lifestyles and forms of host specialism found across the Phytophthora clades as a whole, demonstrating the strong biological cohesiveness of the genus. The finding of 43 previously unknown species from a single Phytophthora clade highlight a critical lack of information on the scale of the unknown pathogen threats to forests and natural ecosystems, underlining the risk of basing plant biosecurity protocols mainly on lists of named organisms. More surveys in natural ecosystems of yet unsurveyed regions in Africa, Asia, Central and South America are needed to unveil the full diversity of the clade and the factors driving diversity, speciation and adaptation in Phytophthora. Taxonomic novelties: New species: Phytophthora amamensis T. Jung, K. Kageyama, H. Masuya & S. Uematsu, Phytophthora angustata T. Jung, L. Garcia, B. Mendieta-Araica, & Y. Balci, Phytophthora balkanensis I. Milenković, Ž. Tomić, T. Jung & M. Horta Jung, Phytophthora borneensis T. Jung, A. Durán, M. Tarigan & M. Horta Jung, Phytophthora calidophila T. Jung, Y. Balci, L. Garcia & B. Mendieta-Araica, Phytophthora catenulata T. Jung, T.-T. Chang, N.M. Chi & M. Horta Jung, Phytophthora celeris T. Jung, L. Oliveira, M. Tarigan & I. Milenković, Phytophthora curvata T. Jung, A. Hieno, H. Masuya & M. Horta Jung, Phytophthora distorta T. Jung, A. Durán, E. Sanfuentes von Stowasser & M. Horta Jung, Phytophthora excentrica T. Jung, S. Uematsu, K. Kageyama & C.M. Brasier, Phytophthora falcata T. Jung, K. Kageyama, S. Uematsu & M. Horta Jung, Phytophthora fansipanensis T. Jung, N.M. Chi, T. Corcobado & C.M. Brasier, Phytophthora frigidophila T. Jung, Y. Balci, K. Broders & I. Milenković, Phytophthora furcata T. Jung, N.M. Chi, I. Milenković & M. Horta Jung, Phytophthora inclinata N.M. Chi, T. Jung, M. Horta Jung & I. Milenković, Phytophthora indonesiensis T. Jung, M. Tarigan, L. Oliveira & I. Milenković, Phytophthora japonensis T. Jung, A. Hieno, H. Masuya & J.F. Webber, Phytophthora limosa T. Corcobado, T. Majek, M. Ferreira & T. Jung, Phytophthora macroglobulosa H.-C. Zeng, H.-H. Ho, F.-C. Zheng & T. Jung, Phytophthora montana T. Jung, Y. Balci, K. Broders & M. Horta Jung, Phytophthora multipapillata T. Jung, M. Tarigan, I. Milenković & M. Horta Jung, Phytophthora multiplex T. Jung, Y. Balci, K. Broders & M. Horta Jung, Phytophthora nimia T. Jung, H. Masuya, A. Hieno & C.M. Brasier, Phytophthora oblonga T. Jung, S. Uematsu, K. Kageyama & C.M. Brasier, Phytophthora obovoidea T. Jung, Y. Balci, L. Garcia & B. Mendieta-Araica, Phytophthora obturata T. Jung, N.M. Chi, I. Milenković & M. Horta Jung, Phytophthora penetrans T. Jung, Y. Balci, K. Broders & I. Milenković, Phytophthora platani T. Jung, A. Pérez-Sierra, S.O. Cacciola & M. Horta Jung, Phytophthora proliferata T. Jung, N.M. Chi, I. Milenković & M. Horta Jung, Phytophthora pseudocapensis T. Jung, T.-T. Chang, I. Milenković & M. Horta Jung, Phytophthora pseudocitrophthora T. Jung, S.O. Cacciola, J. Bakonyi & M. Horta Jung, Phytophthora pseudofrigida T. Jung, A. Durán, M. Tarigan & M. Horta Jung, Phytophthora pseudoccultans T. Jung, T.-T. Chang, I. Milenković & M. Horta Jung, Phytophthora pyriformis T. Jung, Y. Balci, K.D. Boders & M. Horta Jung, Phytophthora sumatera T. Jung, M. Tarigan, M. Junaid & A. Durán, Phytophthora transposita T. Jung, K. Kageyama, C.M. Brasier & H. Masuya, Phytophthora vacuola T. Jung, H. Masuya, K. Kageyama & J.F. Webber, Phytophthora valdiviana T. Jung, E. Sanfuentes von Stowasser, A. Durán & M. Horta Jung, Phytophthora variepedicellata T. Jung, Y. Balci, K. Broders & I. Milenković, Phytophthora vietnamensis T. Jung, N.M. Chi, I. Milenković & M. Horta Jung, Phytophthora ×australasiatica T. Jung, N.M. Chi, M. Tarigan & M. Horta Jung, Phytophthora ×lusitanica T. Jung, M. Horta Jung, C. Maia & I. Milenković, Phytophthora ×taiwanensis T. Jung, T.-T. Chang, H.-S. Fu & M. Horta Jung. Citation: Jung T, Milenković I, Balci Y, Janoušek J, Kudláček T, Nagy ZÁ, Baharuddin B, Bakonyi J, Broders KD, Cacciola SO, Chang T-T, Chi NM, Corcobado T, Cravador A, Đorđević B, Durán A, Ferreira M, Fu C-H, Garcia L, Hieno A, Ho H-H, Hong C, Junaid M, Kageyama K, Kuswinanti T, Maia C, Májek T, Masuya H, Magnano di San Lio G, Mendieta-Araica B, Nasri N, Oliveira LSS, Pane A, Pérez-Sierra A, Rosmana A, Sanfuentes von Stowasser E, Scanu B, Singh R, Stanivuković Z, Tarigan M, Thu PQ, Tomić Z, Tomšovský M, Uematsu S, Webber JF, Zeng H-C, Zheng F-C, Brasier CM, Horta Jung M (2024). Worldwide forest surveys reveal forty-three new species in Phytophthora major Clade 2 with fundamental implications for the evolution and biogeography of the genus and global plant biosecurity. Studies in Mycology 107: 251-388. doi: 10.3114/sim.2024.107.04.

HetMM: A Michaelis-Menten model for non-homogeneous enzyme mixtures

The Michaelis-Menten model requires its reaction velocities to come from a preparation of homogeneous enzymes, with identical or near-identical catalytic activities. However, this condition is not always met. We introduce a kinetic model that relaxes this requirement, by assuming there are an unknown number of enzyme species drawn from a probability distribution whose standard deviation is estimated. Through simulation studies, we demonstrate the method accurately discriminates between homogeneous and heterogeneous data, even with moderate levels of experimental error. We applied this model to three homogeneous and three heterogeneous biological systems, showing that the standard and heterogeneous models outperform respectively. Lastly, we show that heterogeneity is not readily distinguished from negatively cooperative binding under the Hill model. These two distinct attributes-inequality in catalytic ability and interference between binding sites-yield similar Michaelis-Menten curves that are not readily resolved without further experimentation. Our user-friendly software package allows homogeneity testing and parameter estimation.

Young evolutionary origins of dioecy in the genus Asparagus

PREMISE

Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking.

METHODS

We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature.

RESULTS

We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers.

CONCLUSIONS

Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

Enzymic recognition of amino acids drove the evolution of primordial genetic codes

How genetic information gained its exquisite control over chemical processes needed to build living cells remains an enigma. Today, the aminoacyl-tRNA synthetases (AARS) execute the genetic codes in all living systems. But how did the AARS that emerged over three billion years ago as low-specificity, protozymic forms then spawn the full range of highly-specific enzymes that distinguish between 22 diverse amino acids? A phylogenetic reconstruction of extant AARS genes, enhanced by analysing modular acquisitions, reveals six AARS with distinct bacterial, archaeal, eukaryotic, or organellar clades, resulting in a total of 36 families of AARS catalytic domains. Small structural modules that differentiate one AARS family from another played pivotal roles in discriminating between amino acid side chains, thereby expanding the genetic code and refining its precision. The resulting model shows a tendency for less elaborate enzymes, with simpler catalytic domains, to activate amino acids that were not synthesised until later in the evolution of the code. The most probable evolutionary route for an emergent amino acid type to establish a place in the code was by recruiting older, less specific AARS, rather than adapting contemporary lineages. This process, retrofunctionalisation, differs from previously described mechanisms through which amino acids would enter the code.

The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Unraveling phylogenetic relationships and species boundaries in the arid adapted Gerbillus rodents (Muridae: Gerbillinae) by RAD-seq data

Gerbillus is one of the most speciose genera among rodents, with ca. 51 recognized species. Previous attempts to reconstruct the evolutionary history of Gerbillus mainly relied on the mitochondrial cyt-b marker as a source of phylogenetic information. In this study, we utilize RAD-seq genomic data from 37 specimens representing 11 species to reconstruct the phylogenetic tree for Gerbillus, applying concatenation and coalescence methods. We identified four highly supported clades corresponding to the traditionally recognized subgenera: Dipodillus, Gerbillus, Hendecapleura and Monodia. Only two uncertain branches were detected in the resulting trees, with one leading to diversification of the main lineages in the genus, recognized by quartet sampling analysis as uncertain due to possible introgression. We also examined species boundaries for four pairs of sister taxa, including potentially new species from Morocco, using SNAPP. The results strongly supported a speciation model in which all taxa are treated as separate species. The dating analyses confirmed the Plio-Pleistocene diversification of the genus, with the uncertain branch coinciding with the beginning of aridification of the Sahara at the the Plio-Pleistocene boundary. This study aligns well with the earlier analyses based on the cyt-b marker, reaffirming its suitability as an adequate marker for estimating genetic diversity in Gerbillus.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

Shift in symbiotic lifestyle as the major process shaping the evolution of pea crabs (Decapoda: Brachyura: Pinnotheroidea)

The pea crabs, superfamily Pinnotheroidea, are exceptional among brachyuran crabs in their diverse symbiotic associations involving both inquilinism and protective symbiosis. While this group presents a rare opportunity for evolutionary comparative study of host switching and morphological evolution in marine macroinvertebrates, previous phylogenetic studies have been focused on systematics. Here, we reconstructed the most extensive phylogeny of Pinnotheroidea based on two mitochondrial and six nuclear markers, with the aim of elucidating the host switching pathways and the correlation between symbiotic lifestyles and selected morphological adaptations. Ancestral state reconstruction of host association revealed a monophyletic origin of symbiosis in the form of inquilinism. Subsequent shifts in microhabitat preference for burrows or worm tubes, and the move to protective symbiosis, primarily in the switch to mollusc endosymbiosis, contributed to radiation in Pinnotheridae. Further parallel colonisations of echinoderms and tunicates occurred but did not lead to extensive diversification, except in the Clypeasterophilus + Dissodactylus lineage, which experienced a unique switch to echinoderm ectosymbiosis. The evolution of the third maxillipeds, carapace shape and ambulatory pereiopods suggests a rather strong coupling with the symbiotic lifestyle (whether inquilinism or protective symbiosis). Phenotypic diversity of these characters was higher among species engaged in protective symbiosis, with convergence in form (or function) among those sharing the same host affiliation. Species having different host affiliations or symbiotic lifestyles might also exhibit convergence in the form of the three morphological traits, suggesting a common adaptive value of the specialisations. Pinnotherid crabs overall exhibited a lower trait diversity than the also symbiotic palaemonid shrimps with comparable species diversity. This may plausibly be attributed to differences in potential for morphological modification to serve additional functions among the traits analysed in the two groups, the less frequent host switching and the less diverse host affiliations, and thus a less complicated evolutionary history in pinnotherids.

A comprehensive genus-level phylogeny and biogeographical history of the Lythraceae based on whole plastome sequences

BACKGROUND AND AIMS

The Lythraceae are a mainly subtropical to tropical family of the order Myrtales with 28 currently accepted genera and approximately 600 species. There is currently no well-supported phylogenetic and biogeographical hypothesis of the Lythraceae incorporating all currently accepted genera, which we sought to provide.

METHODS

Plastomes of representative species of 18 distinct Lythraceae genera were sequenced and annotated. Together with existing sequences, plastomes of all 28 currently accepted genera in the Lythraceae were brought together for the first time. The plastomes were aligned and a Bayesian phylogenetic hypothesis was produced. We then conducted a time-calibrated Bayesian analysis and a biogeographical analysis.

KEY RESULTS

Plastome-based Bayesian and maximum-likelihood phylogenetic trees are generally congruent with recent nuclear phylogenomic data and resolve two deeply branching major clades in the Lythraceae. One major clade concentrates shrubby and arboreal South American and African genera that inhabit seasonally dry environments, with larger, often winged seeds, adapted to dispersal by the wind. The second major clade concentrates North American, Asian, African and several near-cosmopolitan herbaceous, shrubby and arboreal genera, often inhabiting humid or aquatic environments, with smaller seeds possessing structures that facilitate dispersal by water.

CONCLUSIONS

We hypothesize that the Lythraceae dispersed early in the Late Cretaceous from South American to North American continents, with subsequent expansion in the Late Cretaceous of a North American lineage through Laurasia to Africa via a boreotropical route. Two later expansions of South American clades to Africa in the Palaeocene and Eocene, respectively, are also hypothesized. Transoceanic dispersal in the family is possibly facilitated by adaptations to aquatic environments that are common to many extant genera of the Lythraceae, where long-distance dispersal and vicariance may be invoked to explain several remarkable disjunct distributions in Lythraceae clades.

Detecting Episodic Evolution through Bayesian Inference of Molecular Clock Models

Molecular evolutionary rate variation is a key aspect of the evolution of many organisms that can be modeled using molecular clock models. For example, fixed local clocks revealed the role of episodic evolution in the emergence of SARS-CoV-2 variants of concern. Like all statistical models, however, the reliability of such inferences is contingent on an assessment of statistical evidence. We present a novel Bayesian phylogenetic approach for detecting episodic evolution. It consists of computing Bayes factors, as the ratio of posterior and prior odds of evolutionary rate increases, effectively quantifying support for the effect size. We conducted an extensive simulation study to illustrate the power of this method and benchmarked it to formal model comparison of a range of molecular clock models using (log) marginal likelihood estimation, and to inference under a random local clock model. Quantifying support for the effect size has higher sensitivity than formal model testing and is straight-forward to compute, because it only needs samples from the posterior and prior distribution. However, formal model testing has the advantage of accommodating a wide range molecular clock models. We also assessed the ability of an automated approach, known as the random local clock, where branches under episodic evolution may be detected without their a priori definition. In an empirical analysis of a data set of SARS-CoV-2 genomes, we find "very strong" evidence for episodic evolution. Our results provide guidelines and practical methods for Bayesian detection of episodic evolution, as well as avenues for further research into this phenomenon.

Integrative species delimitation helps to find the hidden diversity of the leaf-litter frog Ischnocnema manezinho (Garcia, 1996) (Anura, Brachycephalidae), endemic to the southern Atlantic Forest

Background

The delimitation of cryptic species is a challenge for biodiversity conservation. Anurans show high cryptic diversity levels, and molecular species delimitation methods could help identify putative new species. Additionally, species delimitation approaches can provide important results for cryptic species conservation, with integrative methods adding robustness to results. Ischnocnema manezinho was described from Santa Catarina Island (SCI), southern Brazil. More recently, some inventories indicated continental populations supposedly similar in morphology to it. If these records are confirmed as I. manezinho, it would likely change its endangered status on National Red List, removing the species from conservation agendas. We investigated the threatened frog Ischnocnema manezinho, to evaluate if the continental populations belong to this species or if they form an undescribed species complex.

Methods

We used coalescent, distance, and allele-sharing-based species delimitation methods and integrative analyses of morphometric and bioacoustics traits to test evolutionary independence between I. manezinho from SCI, Arvoredo Island, and continental populations.

Results

Ischnocnema manezinho is restricted to Santa Catarina Island, while the five remaining lineages should be further investigated through a taxonomic review. Our results point to a small geographic range of Ischnocnema manezinho. Additionally, the species occurs in isolated fragments of forest in SCI surrounded by expanding urban areas, confirming its status as Endangered. Thus, the protection and monitoring of I. manezinho and the taxonomic description of the continental and Arvoredo Island candidate species should be priorities.

Extensive Cryptic Diversity and Ecological Associations Uncovered among Mexican and Global Collections of Naegleria and Vermamoeba Species by 18S Ribosomal DNA, Internal Transcribed Spacer, and Cytochrome Oxidase Subunit I Sequence Analysis

Free-living amoebae (FLA) are phagocytic protists that play crucial roles in microbial communities as significant microbial grazers. However, our current knowledge of their diversity, ecology, and population genetic structures is marginal due to the shallow and biased sampling of ecosystems and the use of few, poorly resolving molecular markers. Thirty-two FLA were isolated from soil and water samples collected across representative ecosystems of the State of Morelos in Central Mexico, including the drinking water distribution system (DWDS) from the state capital. We classified our isolates as members of Acanthamoeba, Vermamoeba, Naegleria, and Tetramitus by 18S ribosomal DNA (rDNA) sequencing. Vermamoeba isolates were recovered exclusively from the DWDS samples. In contrast, Naegleria strains displayed a broad distribution in soil and water samples across the natural ecosystems. We used a combination of phylogenetic and population genetic analyses of internal transcribed spacer (ITS) and cytochrome oxidase subunit I (COI) sequences from our isolates and a comprehensive set of reference sequences to analyze the currently known diversity of Naegleria spp. Significant associations were uncovered between the most prevalent lineages of Naegleria and Vermamoeba and broad ecological and geographical variables at regional and global levels. The population structure and cryptic diversity within the Naegleria galeacystis-Naegleria americana and Vermamoeba vermiformis species complexes were thoroughly analyzed. Our results prove that the genus Vermamoeba, which was previously thought to consist of only one species, actually encompasses at least seven widely distributed species, as indicated by consistent evidence from Bayesian phylogenetics, two species-delimitation programs, and population genetics analyses. IMPORTANCE Our study sheds new light on the population genetic structure of V. vermiformis and diverse Naegleria species. Using improved molecular markers and advanced analytical approaches, we discovered that N. americana, previously considered a single species, actually contains multiple distinct lineages, as revealed by COI sequencing. These lineages are highly differentiated, with little gene flow between them. Our findings demonstrate that the genus Vermamoeba holds multiple cryptic species, requiring a significant taxonomic revision in light of multilocus sequence analyses. These results advance our understanding of the ecology, molecular systematics, and biogeography of these genera and species complexes at both regional and global scales. This study has significant implications for diagnosing amoebal infections and evaluating health risks associated with FLA in domestic and recreational waters.

Is it inappropriate to ask for your age? Evaluating parameter impact on tree dating in a challenging clade (Macroscelidea)

Sengis (order Macroscelidea) are small mammals endemic to Africa. The taxonomy and phylogeny of sengis has been difficult to resolve due to a lack of clear morphological apomorphies. Molecular phylogenies have already significantly revised sengi systematics, but until now no molecular phylogeny has included all 20 extant species. In addition, the age of origin of the sengi crown clade and the divergence age of its two extant families remain unclear. Two recently published studies based on different datasets and age-calibration parameters (DNA type, outgroup selection, fossil calibration points) proposed highly different divergent age estimates and evolutionary scenarios. We obtained nuclear and mitochondrial DNA from mainly museum specimens using target enrichment of single-stranded DNA libraries to generate the first phylogeny of all extant macroscelidean species. We then explored the effects of different parameters (type of DNA, ratio of ingroup to outgroup sampling, number and type of fossil calibration points) and their resulting impacts on age estimates for the origin and initial diversification of Macroscelidea. We show that, even after correcting for substitution saturation, both using mitochondrial DNA in conjunction with nuclear DNA or alone results in much older ages and different branch lengths than when using nuclear DNA alone. We further show that the former effect can be attributed to insufficient amounts of nuclear data. If multiple calibration points are included, the age of the sengi crown group fossil prior has minimal impact on the estimated time frame of sengi evolution. In contrast, the inclusion or exclusion of outgroup fossil priors has a major effect on the resulting node ages. We also find that a reduced sampling of ingroup species does not significantly affect overall age estimates and that terminal specific substitution rates can serve as a means to evaluate the biological likeliness of the produced temporal estimates. Our study demonstrates how commonly varied parameters in temporal calibration of phylogenies affect age estimates. Dated phylogenies should therefore always be seen in the context of the dataset which was used to produce them.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

Genomic Analysis Revealed the International and Domestic Transmission of Carbapenem-Resistant Klebsiella pneumoniae in Chinese Pediatric Patients

Carbapenem-resistant Klebsiella pneumoniae (CRKP) poses a severe threat to public health worldwide. Based on the genomic analysis of 198 CRKP isolates collected at Shanghai Children's Medical Center over the last 8 years (2013 to 2021), we reported the clinical risk, genetic diversity, and prevalence of antimicrobial resistance (AMR) of CRKP in pediatric patients at the genomic level. We found that the blaNDM genes were the predominant carbapenemase genes, followed by blaKPC-2 and blaIMP. All of the carbapenemases were disseminated mainly by four main types of plasmids, among which one plasmid was associated with a higher risk of bloodstream infections. Notably, we tracked disease outbreaks caused by recent introductions of ST14 CRKP from southeast Asia or western countries, and we reported frequent, repetitive introductions of ST11 from other domestic hospitals that were associated interhospital movement of the patients. The cocirculation of K. pneumoniae and AMR plasmids in hospitals highlights the importance of genome sequencing for monitoring and controlling CRKP infections. IMPORTANCE Carbapenem-resistant Klebsiella pneumoniae (CRKP) infection in pediatric patients differs from that in adults patients in terms of both genetic and phenotypic features, which remain to be elucidated. We present a summary of prevalent CRKP isolates from Chinese pediatric patients over 8 years, demonstrating the prevalence and clinical importance of New Delhi metallo-β-lactamase genes in pediatric patients, mainly describing the genomic features of two predominant CRKP clones (ST11 and ST14) in Chinese children, and identifying four carbapenemase-encoding plasmids that contribute to the transmission of most carbapenemase genes in hospitals. Overall, our research provides valuable information about the international and domestic transmission of CRKP isolates that are prevalent in Chinese children and shows the urgent need for genome sequencing-based surveillance systems for monitoring the transmission of CRKP.

The importance of the Andes in the evolutionary radiation of Sigmodontinae (Rodentia, Cricetidae), the most diverse group of mammals in the Neotropics

The Andean mountains stand out for their striking species richness and endemicity that characterize many emblematic Neotropical clades distributed in or around these mountains. The radiation of the Sigmodontinae subfamily, the most diversified mammalian group in the Neotropics, has been historically related to Andean orogenesis. We aim to evaluate this interplay between geological processes and biological responses through the diversification dynamics, the biogeographical history, and the range evolution of the subfamily. For these, we built the most comprehensive phylogeny and gathered 14,836 occurrences for the subfamily. We identified one shift in the speciation rate in the genus Akodon, which suffered their Andean radiation after the arrival of non-Andean ancestors. Our biogeographic analyses show multiple dispersal paths throughout the evolution that allowed this subfamily to colonize all Neotropics. The Northern Andes and Central-Southern Andes were the most important sources of diversity. In addition, the Central-Southern Andes were the most relevant sink, receiving the highest number of lineages. The Andean region exhibited higher speciation and turnover rates than non-Andean regions. Thus, our results support the crucial role of the Andean Mountains in the Sigmodontinae radiation, acting as a "macroevolutionary cradle" and "species attractor" for several sigmodontine lineages at different times, and as a "species pump" becoming the biogeographic source of multiple widely distributed neotropical lineages. Then, complex macroevolutionary dynamics would explain these rodents' high extant Andean diversity and their wide distribution in the Neotropics.

Dated Phylogeny of Astragalus Section Stereothrix (Fabaceae) and Allied Taxa in the Hypoglottis Clade

The Astragalus subgenus Hypoglottis Bunge, which consists of several sections, is one of the taxonomically most complicated groups in the genus. The Astragalus section Stereothrix Bunge belongs to this subgenus and is a significant element of the Irano-Turanian floristic region. A molecular phylogenetic analysis of this section and its closely related taxa using nuclear ribosomal DNA internal transcribed spacers (ITS) and external transcribed spacer (ETS) regions as well as plastid matK sequences were conducted. Parsimony analyses and Bayesian phylogenetic inference revealed that the section is not monophyletic in its current form, as some taxa belonging to closely related sections such as Hypoglottidei DC. and the Malacothrix Bunge group within the sect. Stereothrix render it paraphyletic. Moreover, species groups belonging to sect. Stereothrix are placed in different clades within the phylogenetic tree of subgenus Hypoglottis, which indicates polyphyly, i.e., multiple independent origins of taxa placed in the sect. Stereothrix. Molecular dating of the group estimated an age of 3.62 (1.73-5.62) My for this assemblage with the major diversification events happening during the last 2 My. Many species groups separated only within the last 0.5 to 1 My. Based on morphological and molecular data, we discuss the phylogenetic relationships of the groups and synonymy of species. In addition, the included taxa of sect. Hypoglottidei are not monophyletic and include species belonging to sects. Hololeuce, Koelziana, Malacothrix, Onobrychoideae, and Ornithodpodium group within the sect. Stereothrix taxa. We conclude that only an analysis including all groups and nearly all species of the sections within the Hypoglottis clade can finally result in an new evolutionary-based system for these taxa.

No phylogenomic support for a Cenozoic origin of the "living fossil" Isoetes

PREMISE

The isoetalean lineage has a rich fossil record that extends to the Devonian, but the age of the living clade is unclear. Recent results indicate that it is young, from the Cenozoic, whereas earlier work based on less data from a denser taxon sampling yielded Mesozoic median ages.

METHODS

We analyzed node ages in Isoetes using two genomic data sets (plastome and nuclear ribosomal cistron), three clock models implemented in MrBayes (ILN, WN, and TK02 models), and a conservative approach to calibration.

RESULTS

While topological results were consistently resolved in Isoetes estimated crown group ages range from the latest Paleozoic (mid-Permian) to the Mesozoic depending on data type and clock model. The oldest estimates were retrieved using the autocorrelated TK02 clock model. An (early) Cenozoic age was only obtained under one specific condition (plastome data analyzed with the uncorrelated ILN clock model). That same plastome data set also yielded the oldest (mid-Permian) age estimate when analyzed with the autocorrelated TK02 clock model. Adding the highly divergent, recently established sister species Isoetes wormaldii to the data set approximately doubled the average median node depth to the Isoetes crown group.

CONCLUSIONS

There is no consistent support for a Cenozoic origin of the living clade Isoetes. We obtained seemingly well-founded, yet strongly deviating results depending on data type and clock model. The single most important future improvement is probably to add calibration points, which requires an improved understanding of the isoetalean fossil record or alternative bases for calibration.

Population structure and evolutionary history of the greater cane rat (Thryonomys swinderianus) from the Guinean Forests of West Africa

Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (F_ST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.

SIEVE: joint inference of single-nucleotide variants and cell phylogeny from single-cell DNA sequencing data

We present SIEVE, a statistical method for the joint inference of somatic variants and cell phylogeny under the finite-sites assumption from single-cell DNA sequencing. SIEVE leverages raw read counts for all nucleotides and corrects the acquisition bias of branch lengths. In our simulations, SIEVE outperforms other methods in phylogenetic reconstruction and variant calling accuracy, especially in the inference of homozygous variants. Applying SIEVE to three datasets, one for triple-negative breast (TNBC), and two for colorectal cancer (CRC), we find that double mutant genotypes are rare in CRC but unexpectedly frequent in the TNBC samples.

Tracing the international arrivals of SARS-CoV-2 Omicron variants after Aotearoa New Zealand reopened its border

In the second quarter of 2022, there was a global surge of emergent SARS-CoV-2 lineages that had a distinct growth advantage over then-dominant Omicron BA.1 and BA.2 lineages. By generating 10,403 Omicron genomes, we show that Aotearoa New Zealand observed an influx of these immune-evasive variants (BA.2.12.1, BA.4, and BA.5) through the border. This is explained by the return to significant levels of international travel following the border's reopening in March 2022. We estimate one Omicron transmission event from the border to the community for every ~5,000 passenger arrivals at the current levels of travel and restriction. Although most of these introductions did not instigate any detected onward transmission, a small minority triggered large outbreaks. Genomic surveillance at the border provides a lens on the rate at which new variants might gain a foothold and trigger new waves of infection.

A Multigene Phylogeny of Native American Hawkweeds (Hieracium Subgen. Chionoracium, Cichorieae, Asteraceae): Origin, Speciation Patterns, and Migration Routes

Native American hawkweeds are mainly mountainous species that are distributed all over the New World. They are severely understudied with respect to their origin, colonization of the vast distribution area, and species relationships. Here, we attempt to reconstruct the evolutionary history of the group by applying seven molecular markers (plastid, nuclear ribosomal and low-copy genes). Phylogenetic analyses revealed that Chionoracium is a subgenus of the mainly Eurasian genus Hieracium, which originated from eastern European hawkweeds about 1.58-2.24 million years ago. Plastid DNA suggested a single origin of all Chionoracium species. They colonized the New World via Beringia and formed several distinct lineages in North America. Via one Central American lineage, the group colonized South America and radiated into more than a hundred species within about 0.8 million years, long after the closure of the Isthmus of Panama and the most recent uplift of the Andes. Despite some incongruences shown by different markers, most of them revealed the same crown groups of closely related taxa, which were, however, largely in conflict with traditional sectional classifications. We provide a basic framework for further elucidation of speciation patterns. A thorough taxonomic revision of Hieracium subgen. Chionoracium is recommended.

Nucleotide Substitutions during Speciation may Explain Substitution Rate Variation

Although molecular mechanisms associated with the generation of mutations are highly conserved across taxa, there is widespread variation in mutation rates between evolutionary lineages. When phylogenies are reconstructed based on nucleotide sequences, such variation is typically accounted for by the assumption of a relaxed molecular clock, which is a statistical distribution of mutation rates without much underlying biological mechanism. Here, we propose that variation in accumulated mutations may be partly explained by an elevated mutation rate during speciation. Using simulations, we show how shifting mutations from branches to speciation events impacts inference of branching times in phylogenetic reconstruction. Furthermore, the resulting nucleotide alignments are better described by a relaxed than by a strict molecular clock. Thus, elevated mutation rates during speciation potentially explain part of the variation in substitution rates that is observed across the tree of life. [Molecular clock; phylogenetic reconstruction; speciation; substitution rate variation.].

Quantitatively defining species boundaries with more efficiency and more biological realism

We introduce a widely applicable species delimitation method based on the multispecies coalescent model that is more efficient and more biologically realistic than existing methods. We extend a threshold-based method to allow the ancestral speciation rate to vary through time as a smooth piecewise function. Furthermore, we introduce the cutting-edge proposal kernels of StarBeast3 to this model, thus enabling rapid species delimitation on large molecular datasets and allowing the use of relaxed molecular clock models. We validate these methods with genomic sequence data and SNP data, and show they are more efficient than existing methods at achieving parameter convergence during Bayesian MCMC. Lastly, we apply these methods to two datasets (Hemidactylus and Galagidae) and find inconsistencies with the published literature. Our methods are powerful for rapid quantitative testing of species boundaries in large multilocus datasets and are implemented as an open source BEAST 2 package called SPEEDEMON.

Introducing SPEEDEMON, a package for BEAST 2 that better defines species boundaries based on molecular data demonstrated on gecko and loris datasets.

Genomic epidemiology of Delta SARS-CoV-2 during transition from elimination to suppression in Aotearoa New Zealand

New Zealand's COVID-19 elimination strategy heavily relied on the use of genomics to inform contact tracing, linking cases to the border and to clusters during community outbreaks. In August 2021, New Zealand entered its second nationwide lockdown after the detection of a single community case with no immediately apparent epidemiological link to the border. This incursion resulted in the largest outbreak seen in New Zealand caused by the Delta Variant of Concern. Here we generated 3806 high quality SARS-CoV-2 genomes from cases reported in New Zealand between 17 August and 1 December 2021, representing 43% of reported cases. We detected wide geographical spread coupled with undetected community transmission, characterised by the apparent extinction and reappearance of genomically linked clusters. We also identified the emergence, and near replacement, of genomes possessing a 10-nucleotide frameshift deletion that caused the likely truncation of accessory protein ORF7a. By early October, New Zealand moved from an elimination strategy to a suppression strategy and the role of genomics changed markedly from being used to track and trace, towards population-level surveillance.

An Efficient Coalescent Epoch Model for Bayesian Phylogenetic Inference

We present a two-headed approach called Bayesian Integrated Coalescent Epoch PlotS (BICEPS) for efficient inference of coalescent epoch models. Firstly, we integrate out population size parameters, and secondly, we introduce a set of more powerful Markov chain Monte Carlo (MCMC) proposals for flexing and stretching trees. Even though population sizes are integrated out and not explicitly sampled through MCMC, we are still able to generate samples from the population size posteriors. This allows demographic reconstruction through time and estimating the timing and magnitude of population bottlenecks and full population histories. Altogether, BICEPS can be considered a more muscular version of the popular Bayesian skyline model. We demonstrate its power and correctness by a well-calibrated simulation study. Furthermore, we demonstrate with an application to SARS-CoV-2 genomic data that some analyses that have trouble converging with the traditional Bayesian skyline prior and standard MCMC proposals can do well with the BICEPS approach. BICEPS is available as open-source package for BEAST 2 under GPL license and has a user-friendly graphical user interface.[Bayesian phylogenetics; BEAST 2; BICEPS; coalescent model.].

Multiple Recent Colonizations of the Australian Region by the Chydorus sphaericus Group (Crustacea: Cladocera)

Biotic introductions are an ongoing disruption for many ecosystems. For passively dispersed freshwater zooplankton, transcontinental introductions have been common but are poorly studied in the southern hemisphere. Here we assess the hypothesis of recent introduction for populations of the Chydorus sphaericus group (Crustacea: Cladocera) in Australia. We analyzed 254 sequences (63 original sequences) from the cytochrome oxidase I region of mitochondrial DNA of Chydorus sp., which included global representation. Three Australian populations were connected with separate clades in the northern hemisphere, suggesting multiple colonization events for Australia. The timescale of the divergences was consistent with recent (Quaternary) dispersal. As Australian populations are exposed to migrating birds from the northern hemisphere, both avian and anthropogenic sources are candidates for dispersal vectors. We concluded that recent cross-hemisphere dispersal in the Chydorus sphaericus group is more common than previously believed.

The Emergence of SARS-CoV-2 Variants of Concern Is Driven by Acceleration of the Substitution Rate

The ongoing SARS-CoV-2 pandemic has seen an unprecedented amount of rapidly generated genome data. These data have revealed the emergence of lineages with mutations associated to transmissibility and antigenicity, known as variants of concern (VOCs). A striking aspect of VOCs is that many of them involve an unusually large number of defining mutations. Current phylogenetic estimates of the substitution rate of SARS-CoV-2 suggest that its genome accrues around two mutations per month. However, VOCs can have 15 or more defining mutations and it is hypothesized that they emerged over the course of a few months, implying that they must have evolved faster for a period of time. We analyzed genome sequence data from the GISAID database to assess whether the emergence of VOCs can be attributed to changes in the substitution rate of the virus and whether this pattern can be detected at a phylogenetic level using genome data. We fit a range of molecular clock models and assessed their statistical performance. Our analyses indicate that the emergence of VOCs is driven by an episodic increase in the substitution rate of around 4-fold the background phylogenetic rate estimate that may have lasted several weeks or months. These results underscore the importance of monitoring the molecular evolution of the virus as a means of understanding the circumstances under which VOCs may emerge.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and across the tree of life.

OUP accepted manuscript

As genomic sequence data become increasingly available, inferring the phylogeny of the species as that of concatenated genomic data can be enticing. However, this approach makes for a biased estimator of branch lengths and substitution rates and an inconsistent estimator of tree topology. Bayesian multispecies coalescent (MSC) methods address these issues. This is achieved by constraining a set of gene trees within a species tree and jointly inferring both under a Bayesian framework. However, this approach comes at the cost of increased computational demand. Here, we introduce StarBeast3—a software package for efficient Bayesian inference under the MSC model via Markov chain Monte Carlo. We gain efficiency by introducing cutting-edge proposal kernels and adaptive operators, and StarBeast3 is particularly efficient when a relaxed clock model is applied. Furthermore, gene-tree inference is parallelized, allowing the software to scale with the size of the problem. We validated our software and benchmarked its performance using three real and two synthetic data sets. Our results indicate that StarBeast3 is up to one-and-a-half orders of magnitude faster than StarBeast2, and therefore more than two orders faster than *BEAST, depending on the data set and on the parameter, and can achieve convergence on large data sets with hundreds of genes. StarBeast3 is open-source and is easy to set up with a friendly graphical user interface. [Adaptive; Bayesian inference; BEAST 2; effective population sizes; high performance; multispecies coalescent; parallelization; phylogenetics.]