Distinct evolution of SARS-CoV-2 Omicron XBB and BA.2.86/JN.1 lineages combining increased fitness and antibody evasion.

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.

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor.

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Integron cassettes integrate into bacterial genomes via widespread non-classical attG sites.

Integrons are genetic elements involved in bacterial adaptation which capture, shuffle and express genes encoding adaptive functions embedded in cassettes. These events are governed by the integron integrase through site-specific recombination between attC and attI integron sites. Using computational and molecular genetic approaches, here we demonstrate that the integrase also catalyses cassette integration into bacterial genomes outside of its known att sites. Once integrated, these cassettes can be expressed if located near bacterial promoters and can be excised at the integration point or outside, inducing chromosomal modifications in the latter case. Analysis of more than 5 × 105 independent integration events revealed a very large genomic integration landscape. We identified consensus recombination sequences, named attG sites, which differ greatly in sequence and structure from classical att sites. These results unveil an alternative route for dissemination of adaptive functions in bacteria and expand the role of integrons in bacterial evolution.

Distinct evolution of SARS-CoV-2 Omicron XBB and BA.2.86/JN.1 lineages combining increased fitness and antibody evasion.

The unceasing circulation of SARS-CoV-2 leads to the continuous emergence of novel viral sublineages. Here, we isolated and characterized 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 replicated in IGROV-1 but no longer in Vero E6 and were not markedly fusogenic. They potently infected nasal epithelial cells, with EG.5.1.3 exhibiting the highest fitness. Antivirals remained active. Neutralizing antibody (NAb) responses from vaccinees and BA.1/BA.2-infected individuals were markedly lower compared to BA.1, without major differences between variants. An XBB breakthrough infection enhanced NAb responses against both XBB and BA.2.86 variants. JN.1 displayed 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.

Robustness of Felsenstein's Versus Transfer Bootstrap Supports With Respect to Taxon Sampling.

The bootstrap method is based on resampling sequence alignments and re-estimating trees. Felsenstein's bootstrap proportions (FBP) are the most common approach to assess the reliability and robustness of sequence-based phylogenies. However, when increasing taxon sampling (i.e., the number of sequences) to hundreds or thousands of taxa, FBP tend to return low support for deep branches. The transfer bootstrap expectation (TBE) has been recently suggested as an alternative to FBP. TBE is measured using a continuous transfer index in [0,1] for each bootstrap tree, instead of the binary {0,1} index used in FBP to measure the presence/absence of the branch of interest. TBE has been shown to yield higher and more informative supports while inducing a very low number of falsely supported branches. Nonetheless, it has been argued that TBE must be used with care due to sampling issues, especially in datasets with a high number of closely related taxa. In this study, we conduct multiple experiments by varying taxon sampling and comparing FBP and TBE support values on different phylogenetic depths, using empirical datasets. Our results show that the main critique of TBE stands in extreme cases with shallow branches and highly unbalanced sampling among clades, but that TBE is still robust in most cases, while FBP is inescapably negatively impacted by high taxon sampling. We suggest guidelines and good practices in TBE (and FBP) computing and interpretation.

BioConvert: a comprehensive format converter for life sciences.

Bioinformatics is a field known for the numerous standards and formats that have been developed over the years. This plethora of formats, sometimes complementary, and often redundant, poses many challenges to bioinformatics data analysts. They constantly need to find the best tool to convert their data into the suitable format, which is often a complex, technical and time consuming task. Moreover, these small yet important tasks are often difficult to make reproducible. To overcome these difficulties, we initiated BioConvert, a collaborative project to facilitate the conversion of life science data from one format to another. BioConvert aggregates existing software within a single framework and complemented them with original code when needed. It provides a common interface to make the user experience more streamlined instead of having to learn tens of them. Currently, BioConvert supports about 50 formats and 100 direct conversions in areas such as alignment, sequencing, phylogeny, and variant calling. In addition to being useful for end-users, BioConvert can also be utilized by developers as a universal benchmarking framework for evaluating and comparing numerous conversion tools. Additionally, we provide a web server implementing an online user-friendly interface to BioConvert, hence allowing direct use for the community.

Reprohackathons: promoting reproducibility in bioinformatics through training.

MOTIVATION: The reproducibility crisis has highlighted the importance of improving the way bioinformatics data analyses are implemented, executed, and shared. To address this, various tools such as content versioning systems, workflow management systems, and software environment management systems have been developed. While these tools are becoming more widely used, there is still much work to be done to increase their adoption. The most effective way to ensure reproducibility becomes a standard part of most bioinformatics data analysis projects is to integrate it into the curriculum of bioinformatics Master's programs. RESULTS: In this article, we present the Reprohackathon, a Master's course that we have been running for the last 3 years at Université Paris-Saclay (France), and that has been attended by a total of 123 students. The course is divided into two parts. The first part includes lessons on the challenges related to reproducibility, content versioning systems, container management, and workflow systems. In the second part, students work on a data analysis project for 3-4 months, reanalyzing data from a previously published study. The Reprohackaton has taught us many valuable lessons, such as the fact that implementing reproducible analyses is a complex and challenging task that requires significant effort. However, providing in-depth teaching of the concepts and the tools during a Master's degree program greatly improves students' understanding and abilities in this area.

APOBEC3F Is a Mutational Driver of the Human Monkeypox Virus Identified in the 2022 Outbreak.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Developing and reusing bioinformatics data analysis pipelines using scientific workflow systems.

Data analysis pipelines are now established as an effective means for specifying and executing bioinformatics data analysis and experiments. While scripting languages, particularly Python, R and notebooks, are popular and sufficient for developing small-scale pipelines that are often intended for a single user, it is now widely recognized that they are by no means enough to support the development of large-scale, shareable, maintainable and reusable pipelines capable of handling large volumes of data and running on high performance computing clusters. This review outlines the key requirements for building large-scale data pipelines and provides a mapping of existing solutions that fulfill them. We then highlight the benefits of using scientific workflow systems to get modular, reproducible and reusable bioinformatics data analysis pipelines. We finally discuss current workflow reuse practices based on an empirical study we performed on a large collection of workflows.

Enhanced neutralization escape to therapeutic monoclonal antibodies by SARS-CoV-2 omicron sub-lineages.

The landscape of SARS-CoV-2 variants dramatically diversified with the simultaneous appearance of multiple subvariants originating from BA.2, BA.4, and BA.5 Omicron sub-lineages. They harbor a specific set of mutations in the spike that can make them more evasive to therapeutic monoclonal antibodies. In this study, we compared the neutralizing potential of monoclonal antibodies against the Omicron BA.2.75.2, BQ.1, BQ.1.1, and XBB variants, with a pre-Omicron Delta variant as a reference. Sotrovimab retains some activity against BA.2.75.2, BQ.1, and XBB as it did against BA.2/BA.5, but is less active against BQ.1.1. Within the Evusheld/AZD7442 cocktail, Cilgavimab lost all activity against all subvariants studied, resulting in loss of Evusheld activity. Finally, Bebtelovimab, while still active against BA.2.75, also lost all neutralizing activity against BQ.1, BQ.1.1, and XBB variants.

In-Depth Characterization of Full-Length Archived Viral Genomes after Nine Years of Posttreatment HIV Control.

In the search for control of human immunodeficiency virus type 1 (HIV-1) infection without antiretroviral therapy, posttreatment controllers (PTCs) are models of HIV remission. To better understand their mechanisms of control, we characterized the HIV blood reservoirs of 8 PTCs (median of 9.4 years after treatment interruption) in comparison with those of 13 natural HIV infection controllers (HICs) (median of 18 years of infection) and with those of individuals receiving efficient antiretroviral therapy initiated during either primary HIV infection (PHIs; n = 8) or chronic HIV infection (CHIs; n = 6). This characterization was performed with single-genome amplification and deep sequencing. The proviral diversity, which reflects the history of past viral replication, was lower in the PTCs, PHIs, and aviremic HICs than in the blipper HICs and CHIs. The proportions of intact and defective proviruses among the proviral pool in PTCs were not significantly different from those of other groups. When looking at the quantities of proviruses per million peripheral blood mononuclear cells (PBMCs), they had similar amounts of intact proviruses as other groups but smaller amounts of defective proviruses than CHIs, suggesting a role of these forms in HIV pathogenesis. Two HICs but none of the PTCs harbored only proviruses with deletion in nef; these attenuated strains could contribute to viral control in these participants. We show, for the first time, the presence of intact proviruses and low viral diversity in PTCs long after treatment interruption, as well as the absence of evolution of the proviral quasispecies in subsequent samples. This reflects low residual replication over time. Further data are necessary to confirm these results. IMPORTANCE Most people living with HIV need antiretroviral therapy to control their infection and experience viral relapse in case of treatment interruption, because of viral reservoir (proviruses) persistence. Knowing that proviruses are very diverse and most of them are defective in treated individuals, we aimed to characterize the HIV blood reservoirs of posttreatment controllers (PTCs), rare models of drug-free remission, in comparison with spontaneous controllers and treated individuals. At a median time of 9 years after treatment interruption, which is unprecedented in the literature, we showed that the proportions and quantities of intact proviruses were similar between PTCs and other individuals. Unlike 2/7 spontaneous controllers who harbored only nef-deleted proviruses, which are attenuated strains, which could contribute to their control, no such case was observed in PTCs. Furthermore, PTCs displayed low viral genetic diversity and no evolution of their reservoirs, indicating very low residual replication, despite the presence of intact proviruses.

Deciphering RNA G-quadruplex function during the early steps of HIV-1 infection.

G-quadruplexes (G4s) are four-stranded nucleic acid structures formed by the stacking of G-tetrads. Here we investigated their formation and function during HIV-1 infection. Using bioinformatics and biophysics analyses we first searched for evolutionary conserved G4-forming sequences in HIV-1 genome. We identified 10 G4s with conservation rates higher than those of HIV-1 regulatory sequences such as RRE and TAR. We then used porphyrin-based G4-binders to probe the formation of the G4s during infection of human cells by native HIV-1. The G4-binders efficiently inhibited HIV-1 infectivity, which is attributed to the formation of G4 structures during HIV-1 replication. Using a qRT-PCR approach, we showed that the formation of viral G4s occurs during the first 2 h post-infection and their stabilization by the G4-binders prevents initiation of reverse transcription. We also used a G4-RNA pull-down approach, based on a G4-specific biotinylated probe, to allow the direct detection and identification of viral G4-RNA in infected cells. Most of the detected G4-RNAs contain crucial regulatory elements such as the PPT and cPPT sequences as well as the U3 region. Hence, these G4s would function in the early stages of infection when the viral RNA genome is being processed for the reverse transcription step.

Prehospital Lyophilized Plasma Transfusion for Trauma-Induced Coagulopathy in Patients at Risk for Hemorrhagic Shock: A Randomized Clinical Trial.

Importance: Blood transfusion is a mainstay of therapy for trauma-induced coagulopathy, but the optimal modalities for plasma transfusion in the prehospital setting remain to be defined. Objective: To determine whether lyophilized plasma transfusion can reduce the incidence of trauma-induced coagulopathy compared with standard care consisting of normal saline infusion. Design, Setting, and Participants: This randomized clinical trial was performed at multiple centers in France involving prehospital medical teams. Participants included 150 adults with trauma who were at risk for hemorrhagic shock and associated coagulopathy between April 1, 2016, and September 30, 2019, with a 28-day follow-up. Data were analyzed from November 1, 2019, to July 1, 2020. Intervention: Patients were randomized in a 1:1 ratio to receive either plasma or standard care with normal saline infusion (control). Main Outcomes and Measures: The primary outcome was the international normalized ratio (INR) on arrival at the hospital. Secondary outcomes included the need for massive transfusion and 30-day survival. As a safety outcome, prespecified adverse events included thrombosis, transfusion-related acute lung injury, and transfusion-associated circulatory overload. Results: Among 150 randomized patients, 134 were included in the analysis (median age, 34 [IQR, 26-49] years; 110 men [82.1%]), with 68 in the plasma group and 66 in the control group. Median INR values were 1.21 (IQR, 1.12-1.49) in the plasma group and 1.20 (IQR, 1.10-1.39) in the control group (median difference, -0.01 [IQR, -0.09 to 0.08]; P = .88). The groups did not differ significantly in the need for massive transfusion (7 [10.3%] vs 4 [6.1%]; relative risk, 1.78 [95% CI, 0.42-8.68]; P = .37) or 30-day survival (hazard ratio for death, 1.07 [95% CI, 0.44-2.61]; P = .89). In the full intention-to-treat population (n = 150), the groups did not differ in the rates of any of the prespecified adverse events. Conclusions and Relevance: In this randomized clinical trial including severely injured patients at risk for hemorrhagic shock and associated coagulopathy, prehospital transfusion of lyophilized plasma was not associated with significant differences in INR values vs standard care with normal saline infusion. Nevertheless, these findings show that lyophilized plasma transfusion is a feasible and safe procedure for this patient population. Trial Registration: ClinicalTrials.gov Identifier: NCT02736812.

Considerable escape of SARS-CoV-2 Omicron to antibody neutralization.

The SARS-CoV-2 Omicron variant was first identified in November 2021 in Botswana and South Africa1-3. It has since spread to many countries and is expected to rapidly become dominant worldwide. The lineage is characterized by the presence of around 32 mutations in spike-located mostly in the N-terminal domain and the receptor-binding domain-that may enhance viral fitness and enable antibody evasion. Here we isolated an infectious Omicron virus in Belgium from a traveller returning from Egypt. We examined its sensitivity to nine monoclonal antibodies that have been clinically approved or are in development4, and to antibodies present in 115 serum samples from COVID-19 vaccine recipients or individuals who have recovered from COVID-19. Omicron was completely or partially resistant to neutralization by all monoclonal antibodies tested. Sera from recipients of the Pfizer or AstraZeneca vaccine, sampled five months after complete vaccination, barely inhibited Omicron. Sera from COVID-19-convalescent patients collected 6 or 12 months after 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 titres 6-fold to 23-fold lower against Omicron compared with those against Delta. Thus, Omicron escapes most therapeutic monoclonal antibodies and, to a large extent, vaccine-elicited antibodies. However, Omicron is neutralized by antibodies generated by a booster vaccine dose.

Correction to 'Gotree/Goalign: toolkit and Go API to facilitate the development of phylogenetic workflows'.

[This corrects the article DOI: 10.1093/nargab/lqab075.].

Gotree/Goalign: toolkit and Go API to facilitate the development of phylogenetic workflows.

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.

Maternal and neonatal outcomes related to Zika virus in pregnant women in Southern Vietnam: An epidemiological and virological prospective analysis.

BACKGROUND: In 2016-2017, 68 women in Southern Vietnam had RT-PCR confirmed Zika virus (ZIKV) infection during pregnancy. We report here the outcomes of the pregnancies and the virological analyses related to this outbreak. METHODS: We collected clinical and epidemiological information from the women who were enrolled in the study. Medical records related to the pregnancy in 2016-2017 were retrieved for those who were not able to be enrolled in the study. Children born to women with ZIKV infection during pregnancy were also enrolled. Serum samples were evaluated for presence of ZIKV antibodies. Phylogenetic analyses were performed on Zika virus genomes sequenced from the 2016-2017 serum samples. FINDINGS: Of the 68 pregnancies, 58 were livebirths and 10 were medically terminated. Four of the medical records from cases of fetal demise were able to be retrieved, of which one was consistent with congenital ZIKV infection. Of the 58 women with a livebirth, 21 participated in the follow-up investigation. All but two women had serologic evidence of ZIKV infection. Of the 21 children included in the study (mean age: 30.3 months), 3 had microcephaly at birth. No other clinical abnormalities were reported and no differences in neurodevelopment were observed compared to a control group. Phylogenetic analysis revealed a clade within the ZIKV Asian lineage and branch at the root of samples from the 2013-2014 French Polynesian outbreak. The prM S139N mutation was not observed. INTERPRETATION: We have been able to demonstrate a clade within the ZIKV Asian lineage implicated in adverse pregnancy outcomes in Southern Vietnam. FUNDING: INCEPTION project (PIA/ANR-16-CONV-0005) and a grant received from BNP Paribas Simplidon.

Automatic external defibrillator provided by unmanned aerial vehicle (drone) in Greater Paris: A real world-based simulation.

AIM: To reduce the delay in defibrillation of out-of-hospital cardiac arrest (OHCA) patients, recent publications have shown that drones equipped with an automatic external defibrillator (AED) appear to be effective in sparsely populated areas. To study the effectiveness of AED-drones in high-density urban areas, we developed an algorithm based on emergency dispatch parameters for the rate and detection speed of cardiac arrests and technical and meteorological parameters. METHODS: We ran a numerical simulation to compare the actual time required by the Basic Life Support team (BLSt) for OHCA patients in Greater Paris in 2017 to the time required by an AED-drone. Endpoints were the proportion of patients with "AED-drone first" and the defibrillation time gained. We built an open-source website (https://airborne-aed.org/) to allow modelling by modifying one or more parameters and to help other teams model their own OHCA data. RESULTS: Of 3014 OHCA patients, 72.2 ±â€¯0.7% were in the "no drone flight" group, 25.8 ±â€¯0.2% in the "AED-drone first" group, and 2.1 ±â€¯0.2% in the "BLSt-drone first" group. When a drone flight was authorized, it arrived an average 190 s before BLSt in 93% of cases. The possibility of flying the drone during the aeronautical night improved the results of the "AED-drone first" group the most (+60%). CONCLUSIONS: In our very high-density urban model, at most 26% of OHCA patients received an AED from an AED-drone before BLSt. The flexible parameters of our website model allows evaluation of the impact of each choice and concrete implementation of the AED-drone.

COVID-Align: accurate online alignment of hCoV-19 genomes using a profile HMM.

MOTIVATION: The first cases of the COVID-19 pandemic emerged in December 2019. Until the end of February 2020, the number of available genomes was below 1000 and their multiple alignment was easily achieved using standard approaches. Subsequently, the availability of genomes has grown dramatically. Moreover, some genomes are of low quality with sequencing/assembly errors, making accurate re-alignment of all genomes nearly impossible on a daily basis. A more efficient, yet accurate approach was clearly required to pursue all subsequent bioinformatics analyses of this crucial data. RESULTS: hCoV-19 genomes are highly conserved, with very few indels and no recombination. This makes the profile HMM approach particularly well suited to align new genomes, add them to an existing alignment and filter problematic ones. Using a core of ∼2500 high quality genomes, we estimated a profile using HMMER, and implemented this profile in COVID-Align, a user-friendly interface to be used online or as standalone via Docker. The alignment of 1000 genomes requires ∼50 minutes on our cluster. Moreover, COVID-Align provides summary statistics, which can be used to determine the sequencing quality and evolutionary novelty of input genomes (e.g. number of new mutations and indels). AVAILABILITY AND IMPLEMENTATION: https://covalign.pasteur.cloud, hub.docker.com/r/evolbioinfo/covid-align. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.