Phylogeny.fr, created in 2008, has been designed to facilitate the execution of phylogenetic workflows, and is nowadays widely used. However, since its development, user needs have evolved, new tools and workflows have been published, and the number of jobs has increased dramatically, thus promoting new practices, which motivated its refactoring. We developed NGPhylogeny.fr to be more flexible in terms of tools and workflows, easily installable, and more scalable. It integrates numerous tools in their latest version (e.g. TNT, FastME, MrBayes, etc.) as well as new ones designed in the last ten years (e.g. PhyML, SMS, FastTree, trimAl, BOOSTER, etc.). These tools cover a large range of usage (sequence searching, multiple sequence alignment, model selection, tree inference and tree drawing) and a large panel of standard methods (distance, parsimony, maximum likelihood and Bayesian). They are integrated in workflows, which have been already configured (‘One click’), can be customized (‘Advanced’), or are built from scratch (‘A la carte’). Workflows are managed and run by an underlying Galaxy workflow system, which makes workflows more scalable in terms of number of jobs and size of data. NGPhylogeny.fr is deployable on any server or personal computer, and is freely accessible at https://ngphylogeny.fr.

The unceasing circulation of SARS-CoV-2 leads to the continuous emergence of novel viral sublineages. Here, we isolate and characterize XBB.1, XBB.1.5, XBB.1.9.1, XBB.1.16.1, EG.5.1.1, EG.5.1.3, XBF, BA.2.86.1 and JN.1 variants, representing >80% of circulating variants in January 2024. The XBB subvariants carry few but recurrent mutations in the spike, whereas BA.2.86.1 and JN.1 harbor >30 additional changes. These variants replicate in IGROV-1 but no longer in Vero E6 and are not markedly fusogenic. They potently infect nasal epithelial cells, with EG.5.1.3 exhibiting the highest fitness. Antivirals remain active. Neutralizing antibody (NAb) responses from vaccinees and BA.1/BA.2-infected individuals are markedly lower compared to BA.1, without major differences between variants. An XBB breakthrough infection enhances NAb responses against both XBB and BA.2.86 variants. JN.1 displays lower affinity to ACE2 and higher immune evasion properties compared to BA.2.86.1. Thus, while distinct, the evolutionary trajectory of these variants combines increased fitness and antibody evasion.

Phylogenetics is nowadays at the center of numerous studies in many fields, ranging from comparative genomics to molecular epidemiology. However, phylogenetic analysis workflows are usually complex and difficult to implement, as they are often composed of many small, reccuring, but important data manipulations steps. Among these, we can find file reformatting, sequence renaming, tree re-rooting, tree comparison, bootstrap support computation, etc. These are often performed by custom scripts or by several heterogeneous tools, which may be error prone, uneasy to maintain and produce results that are challenging to reproduce. For all these reasons, the development and reuse of phylogenetic workflows is often a complex task. We identified many operations that are part of most phylogenetic analyses, and implemented them in a toolkit called Gotree/Goalign. The Gotree/Goalign toolkit implements more than 120 user-friendly commands and an API dedicated to multiple sequence alignment and phylogenetic tree manipulations. It is developed in Go, which makes executables easily installable, integrable in workflow environments, and parallelizable when possible. Moreover, Go is a compiled language, which accelerates computations compared to interpreted languages. This toolkit is freely available on most platforms (Linux, MacOS and Windows) and most architectures (amd64, i386) on GitHub at https://github.com/evolbioinfo/gotree, Bioconda and DockerHub.

Widely applicable, accurate and fast inference methods in phylodynamics are needed to fully profit from the richness of genetic data in uncovering the dynamics of epidemics. Standard methods, including maximum-likelihood and Bayesian approaches, generally rely on complex mathematical formulae and approximations, and do not scale with dataset size. We develop a likelihood-free, simulation-based approach, which combines deep learning with (1) a large set of summary statistics measured on phylogenies or (2) a complete and compact representation of trees, which avoids potential limitations of summary statistics and applies to any phylodynamics model. Our method enables both model selection and estimation of epidemiological parameters from very large phylogenies. We demonstrate its speed and accuracy on simulated data, where it performs better than the state-of-the-art methods. To illustrate its applicability, we assess the dynamics induced by superspreading individuals in an HIV dataset of men-having-sex-with-men in Zurich. Our tool PhyloDeep is available on github.com/evolbioinfo/phylodeep.

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa ^1,2 . It has in the meantime spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of about 32 mutations in the Spike, located mostly in the N-terminal domain (NTD) and the receptor binding domain (RBD), which may enhance viral fitness and allow antibody evasion. Here, we isolated an infectious Omicron virus in Belgium, from a traveller returning from Egypt. We examined its sensitivity to 9 monoclonal antibodies (mAbs) clinically approved or in development ³ , and to antibodies present in 90 sera from COVID-19 vaccine recipients or convalescent individuals. Omicron was totally or partially resistant to neutralization by all mAbs tested. Sera from Pfizer or AstraZeneca vaccine recipients, sampled 5 months after complete vaccination, barely inhibited Omicron. Sera from COVID-19 convalescent patients collected 6 or 12 months post symptoms displayed low or no neutralizing activity against Omicron. Administration of a booster Pfizer dose as well as vaccination of previously infected individuals generated an anti-Omicron neutralizing response, with titers 5 to 31 fold lower against Omicron than against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and to a large extent vaccine-elicited antibodies.

On May 6, 2022, a powerful outbreak of monkeypox virus (MPXV) had been reported outside of Africa, with many continuing new cases being reported around the world. Analysis of mutations among the 2 different lineages present in the 2021 and 2022 outbreaks revealed the presence of G->A mutations occurring in the 5′GpA context, indicative of APOBEC3 cytidine deaminase activity. By using a sensitive polymerase chain reaction (differential DNA denaturation PCR) method allowing differential amplification of AT-rich DNA, we analyzed the level of APOBEC3-induced MPXV editing in infected cells and in patients. We demonstrate that G->A hypermutated MPXV genomes can be recovered experimentally from APOBEC3 transfection followed by MPXV infection. Here, among the 7 human APOBEC3 cytidine deaminases (A3A-A3C, A3DE, A3F–A3H), only APOBEC3F was capable of extensively deaminating cytidine residues in MPXV genomes. Hyperedited genomes were also recovered in ∼42% of analyzed patients. Moreover, we demonstrate that substantial repair of these mutations occurs. Upon selection, corrected G->A mutations escaping drift loss contribute to the MPXV evolution observed in the current epidemic. Stochastic or transient overexpression of the APOBEC3F gene exposes the MPXV genome to a broad spectrum of mutations that may be modeling the mutational landscape after multiple cycles of viral replication.