Rapid isolation and resolution immune evasion and viral fitness across contemporary SARS-CoV-2 variants

From late 2020 the world observed the rapid emergence of many distinct SARS-CoV-2 variants. At the same time, pandemic responses coalesced into significant global vaccine roll-out that have now significantly lowered Covid-19 hospital and mortality rates in the developed world. Over this period, we developed a rapid platform (R-20) for viral isolation and characterisation using primary remnant diagnostic swabs. This combined with quarantine testing and genomics surveillance, enabled the rapid isolation and characterisation of all major SARS-CoV-2 variants (all variants of concern and 6 variants of interest) globally with a 4-month period. This platform facilitated viral variant isolation and enabled rapid resolution of variant phenotype by allowing determining end point viral titers from primary nasopharyngeal swabs and through ranking of evasion of neutralising antibodies. In late 2021, when the Delta variant was dominating, Omicron rapidly emerged. Using this platform, we isolated and tested the first cases of this variant within Australia. In this setting we observed Omicron to diverge from other variants at two levels: Firstly, it ranks at the mots evasive to neutralisation antibodies compared to all VOCs and major VUIs. Secondly, it no longer engages TMPRSS2 during the late stages of fusion.

Immunogenicity of BNT162b2 vaccine against the Alpha and Delta variants in immunocompromised patients with systemic inflammatory diseases

ObjectivesThe emergence of strains of SARS-CoV-2 exhibiting increase viral fitness and immune escape potential, such as the Delta variant (B.1.617.2), raises concerns in immunocompromised patients. We aimed to evaluate seroconversion, cross-neutralisation and T-cell responses induced by BNT162b2 in immunocompromised patients with systemic inflammatory diseases.MethodsProspective monocentric study including patients with systemic inflammatory diseases and healthcare immunocompetent workers as controls. Primary endpoints were anti-spike antibodies levels and cross-neutralisation of Alpha and Delta variants after BNT162b2 vaccine. Secondary endpoints were T-cell responses, breakthrough infections and safety.ResultsSixty-four cases and 21 controls not previously infected with SARS-CoV-2 were analysed. Kinetics of anti-spike IgG after BNT162b2 vaccine showed lower and delayed induction in cases, more pronounced with rituximab. Administration of two doses of BNT162b2 generated a neutralising response against Alpha and Delta in 100% of controls, while sera from only one of rituximab-treated patients neutralised Alpha (5%) and none Delta. Other therapeutic regimens induced a partial neutralising activity against Alpha, even lower against Delta. All controls and cases except those treated with methotrexate mounted a SARS-CoV-2 specific T-cell response. Methotrexate abrogated T-cell responses after one dose and dramatically impaired T-cell responses after two doses of BNT162b2. Third dose of vaccine improved immunogenicity in patients with low responses.ConclusionRituximab and methotrexate differentially impact the immunogenicity of BNT162b2, by impairing B-cell and T-cell responses, respectively. Delta fully escapes the humoral response of individuals treated with rituximab. These findings support efforts to improve BNT162b2 immunogenicity in immunocompromised individuals (ClinicalTrials.gov number, NCT04870411).

Structural and computational insights into the SARS-CoV-2 Omicron RBD-ACE2 interaction

Since SARS-CoV-2 Omicron variant (B.1.1.529) was reported in November 2021, it has quickly spread to many countries and outcompeted the globally dominant Delta variant in several countries. The Omicron variant contains the largest number of mutations to date, with 32 mutations located at spike (S) glycoprotein, which raised great concern for its enhanced viral fitness and immune escape[1-4]. In this study, we reported the crystal structure of the receptor binding domain (RBD) of Omicron variant S glycoprotein bound to human ACE2 at a resolution of 2.6 Å. Structural comparison, molecular dynamics simulation and binding free energy calculation collectively identified four key mutations (S477N, G496S, Q498R and N501Y) for the enhanced binding of ACE2 by the Omicron RBD compared to the WT RBD. Representative states of the WT and Omicron RBD-ACE2 systems were identified by Markov State Model, which provides a dynamic explanation for the enhanced binding of Omicron RBD. The effects of the mutations in the RBD for antibody recognition were analyzed, especially for the S371L/S373P/S375F substitutions significantly changing the local conformation of the residing loop to deactivate several class IV neutralizing antibodies.

Full Genome Analysis of Circulating DENV-2 in Senegal Reveals a Regional Diversification Into Separate Clades

To assess the genetic diversity of circulating dengue virus 2 in Senegal in 2018 we performed molecular characterization by complete genome sequencing and performing phylogenetic analysis. Sequenced strains belong to Cosmopolitan genotype of DENV-2 we observed intra-genotype variability leading to a divergence in two clades with differential geographic distribution. We report two variants namely; the “Northern variant” harbouring three nonsynonymous mutations (V1183M, R1405K, P2266T) located respectively on NS2A, NS2B and NS4A and the “Western variant” with two nonsynonymous mutations (V1185E, V3214E) located respectively in the NS2A gene and the NS5 gene. Findings calls for in depth in vitro and functional study to elucidate the impact of observed mutations on viral fitness, spread, epidemiology and disease outcome.

Structural and computational insights into the SARS-CoV-2 Omicron RBD-ACE2 interaction

Since SARS-CoV-2 Omicron variant (B.1.1.529) was reported in November 2021, it has quickly spread to many countries and outcompeted the globally dominant Delta variant in several countries. The Omicron variant contains the largest number of mutations to date, with 32 mutations located at spike (S) glycoprotein, which raised great concern for its enhanced viral fitness and immune escape[1-4]. In this study, we reported the crystal structure of the receptor binding domain (RBD) of Omicron variant S glycoprotein bound to human ACE2 at a resolution of 2.6 angstrom. Structural comparison, molecular dynamics simulation and binding free energy calculation collectively identified four key mutations (S477N, G496S, Q498R and N501Y) for the enhanced binding of ACE2 by the Omicron RBD compared to the WT RBD. Representative states of the WT and Omicron RBD-ACE2 systems were identified by Markov State Model, which provides a dynamic explanation for the enhanced binding of Omicron RBD. The effects of the mutations in the RBD for antibody recognition were analyzed, especially for the S371L/S373P/S375F substitutions significantly changing the local conformation of the residing loop to deactivate several class IV neutralizing antibodies.

SARS-CoV-2 comes up with a perfect recipe for disaster: Each mutation plays a specific role in the evolution of variants of concern

COVID-19 pandemic has extended for close to two years with the continuous emergence of new variants. Mutations in the receptor binding domain (RBD) are of prime importance in dictating the SARS-CoV-2 spike protein function. By studying a series of single, double and triple RBD mutants, we have delineated the individual and collective effects of RBD mutations in a variant of concern (VOC) containing multiple mutations (Gamma variant; K417T/E484K/N501Y) on binding to angiotensin converting enzyme 2 (ACE2) receptor, antibody escape and protein stability. Our results show that each mutation in the VOC serves a distinct function that improves virus fitness landscape supporting its positive selection, even though individual mutations have deleterious effects that make them prone to negative selection. K417T contributes to increased expression, increased stability and escape from class 1 antibodies; however, it has decreased ACE2 binding. E484K contributes to escape from class 2 antibodies; however, it has decreased expression, decreased stability, and decreased ACE2 binding affinity. N501Y increases receptor binding affinity; however, it has decreased stability and decreased expression. But when these mutations come together, the deleterious effects are mitigated in the triple mutant due to the presence of compensatory effects, which improves the chances of selection of mutations together. These results show the implications of presence of multiple mutations on virus evolution and indicate the emergence of future SARS-CoV-2 variants with multiple mutations that enhance viral fitness on different fronts by balancing both positive and negative selection.

A Conserved uORF Regulates APOBEC3G Translation and Is Targeted by HIV-1 Vif Protein to Repress the Antiviral Factor

The HIV-1 Vif protein is essential for viral fitness and pathogenicity. Vif decreases expression of cellular restriction factors APOBEC3G (A3G), A3F, A3D and A3H, which inhibit HIV-1 replication by inducing hypermutation during reverse transcription. Vif counteracts A3G at several levels (transcription, translation, and protein degradation) that altogether reduce the levels of A3G in cells and prevent its incorporation into viral particles. How Vif affects A3G translation remains unclear. Here, we uncovered the importance of a short conserved uORF (upstream ORF) located within two critical stem-loop structures of the 5′ untranslated region (5′-UTR) of A3G mRNA for this process. A3G translation occurs through a combination of leaky scanning and translation re-initiation and the presence of an intact uORF decreases the extent of global A3G translation under normal conditions. Interestingly, the uORF is also absolutely required for Vif-mediated translation inhibition and redirection of A3G mRNA into stress granules. Overall, we discovered that A3G translation is regulated by a small uORF conserved in the human population and that Vif uses this specific feature to repress its translation.

Long-term adaptation following influenza A virus host shifts results in increased within-host viral fitness due to higher replication rates, broader dissemination within the respiratory epithelium and reduced tissue damage

The mechanisms and consequences of genome evolution on viral fitness following host shifts are poorly understood. In addition, viral fitness -the ability of an organism to reproduce and survive- is multifactorial and thus difficult to quantify. Influenza A viruses (IAVs) circulate broadly among wild birds and have jumped into and become endemic in multiple mammalian hosts, including humans, pigs, dogs, seals, and horses. H3N8 equine influenza virus (EIV) is an endemic virus of horses that originated in birds and has been circulating uninterruptedly in equine populations since the early 1960s. Here, we used EIV to quantify changes in infection phenotype associated to viral fitness due to genome-wide changes acquired during long-term adaptation. We performed experimental infections of two mammalian cell lines and equine tracheal explants using the earliest H3N8 EIV isolated (A/equine/Uruguay/63 [EIV/63]), and A/equine/Ohio/2003 (EIV/2003), a monophyletic descendant of EIV/63 isolated 40 years after the emergence of H3N8 EIV. We show that EIV/2003 exhibits increased resistance to interferon, enhanced viral replication, and a more efficient cell-to-cell spread in cells and tissues. Transcriptomics analyses revealed virus-specific responses to each virus, mainly affecting host immunity and inflammation. Image analyses of infected equine respiratory explants showed that despite replicating at higher levels and spreading over larger areas of the respiratory epithelium, EIV/2003 induced milder lesions compared to EIV/63, suggesting that adaptation led to reduced tissue pathogenicity. Our results reveal previously unknown links between virus genotype and the host response to infection, providing new insights on the relationship between virus evolution and fitness.

Compared with SARS-CoV2 wild type's spike protein, the SARS-CoV2 omicron's receptor binding motif has adopted a more SARS-CoV1 and/or bat/civet-like structure.

Our study focuses on free energy calculations of SARS-Cov2 spike protein receptor binding motives (RBMs) from wild type and variants-of-concern with particular emphasis on currently emerging SARS- CoV2 omicron variants of concern (VOC). Our computational free energy analysis underlines the occurrence of positive selection processes that specify omicron host adaption and bring changes on the molecular level into context with clinically relevant observations. Our free energy calculations studies regarding the interaction of omicron's RBM with human ACE2 shows weaker binding to ACE2 than alpha's, delta's, or wild type's RBM. Thus, less virus is predicted to be generated in time per infected cell. Our mutant analyses predict with focus on omicron variants a reduced spike-protein binding to ACE2--receptor protein possibly enhancing viral fitness / transmissibility and resulting in a delayed induction of danger signals as trade-off. Finally, more virus is produced but less per cell accompanied with delayed Covid-19 immunogenicity and pathogenicity. Regarding the latter, more virus is assumed to be required to initiate inflammatory immune responses.

Potential Associations of Mutations within the HIV-1 Env and Gag Genes Conferring Protease Inhibitor (PI) Drug Resistance

An increasing number of patients in Africa are experiencing virological failure on a second-line antiretroviral protease inhibitor (PI)-containing regimen, even without resistance-associated mutations in the protease region, suggesting a potential role of other genes in PI resistance. Here, we investigated the prevalence of mutations associated with Lopinavir/Ritonavir (LPV/r) failure in the Envelope gene and the possible coevolution with mutations within the Gag-protease (gag-PR) region. Env and Gag-PR sequences generated from 24 HIV-1 subtype C infected patients failing an LPV/r inclusive treatment regimen and 344 subtype C drug-naïve isolates downloaded from the Los Alamos Database were analyzed. Fisher’s exact test was used to determine the differences in mutation frequency. Bayesian network probability was applied to determine the relationship between mutations occurring within the env and gag-PR regions and LPV/r treatment. Thirty-five mutations in the env region had significantly higher frequencies in LPV/r-treated patients. A combination of Env and Gag-PR mutations was associated with a potential pathway to LPV/r resistance. While Env mutations were not directly associated with LPV/r resistance, they may exert pressure through the Gag and minor PR mutation pathways. Further investigations using site-directed mutagenesis are needed to determine the impact of Env mutations alone and in combination with Gag-PR mutations on viral fitness and LPV/r efficacy.