Genetic consequences of effective and suboptimal dosing with mutagenic drugs in a hamster model of SARS-CoV-2 infection.

Mutagenic antiviral drugs have shown promise against multiple viruses, but concerns have been raised about whether their use might promote the emergence of new and harmful viral variants. Recently, genetic signatures associated with molnupiravir use have been identified in the global SARS-COV-2 population. Here, we examine the consequences of using favipiravir and molnupiravir to treat SARS-CoV-2 infection in a hamster model, comparing viral genome sequence data collected from (1) untreated hamsters, and (2) from hamsters receiving effective and suboptimal doses of treatment. We identify a broadly linear relationship between drug dose and the extent of variation in treated viral populations, with a high proportion of this variation being composed of variants at frequencies of less than 1 per cent, below typical thresholds for variant calling. Treatment with an effective dose of antiviral drug was associated with a gain of between 7 and 10 variants per viral genome relative to drug-free controls: even after a short period of treatment a population founded by a transmitted virus could contain multiple sequence differences to that of the original host. Treatment with a suboptimal dose of drug showed intermediate gains of variants. No dose-dependent signal was identified in the numbers of single-nucleotide variants reaching frequencies in excess of 5 per cent. We did not find evidence to support the emergence of drug resistance or of novel immune phenotypes. Our study suggests that where onward transmission occurs, a short period of treatment with mutagenic drugs may be sufficient to generate a significant increase in the number of viral variants transmitted.

Favipiravir induces HuNoV viral mutagenesis and infectivity loss with clinical improvement in immunocompromised patients.

Chronic human norovirus (HuNoV) infections in immunocompromised patients result in severe disease, yet approved antivirals are lacking. RNA-dependent RNA polymerase (RdRp) inhibitors inducing viral mutagenesis display broad-spectrum in vitro antiviral activity, but clinical efficacy in HuNoV infections is anecdotal and the potential emergence of drug-resistant variants is concerning. Upon favipiravir (and nitazoxanide) treatment of four immunocompromised patients with life-threatening HuNoV infections, viral whole-genome sequencing showed accumulation of favipiravir-induced mutations which coincided with clinical improvement although treatment failed to clear HuNoV. Infection of zebrafish larvae demonstrated drug-associated loss of viral infectivity and favipiravir treatment showed efficacy despite occurrence of RdRp variants potentially causing favipiravir resistance. This indicates that within-host resistance evolution did not reverse loss of viral fitness caused by genome-wide accumulation of sequence changes. This off-label approach supports the use of mutagenic antivirals for treating prolonged RNA viral infections and further informs the debate surrounding their impact on virus evolution.

Persistent low-level variants in a subset of viral genes are highly predictive of poor outcome in immunocompromised patients with cytomegalovirus infection.

BACKGROUND: Human cytomegalovirus is the most common and serious opportunistic infection after solid organ and haematopoietic stem cell transplantation. In this study, we used whole-genome cytomegalovirus data to investigate viral factors associated with the clinical outcome. METHODS: We sequenced cytomegalovirus samples from 16 immunocompromised paediatric patients with persistent viraemia. 8/16 patients died of complications due to cytomegalovirus infection. We also sequenced samples from 35 infected solid organ adult recipients of whom one died with cytomegalovirus infection. RESULTS: We showed that samples from both groups have fixed variants at resistance sites and mixed infections. NGS sequencing also revealed non-fixed variants at resistance sites in most of the patients who died (6/9). A machine learning approach identified 10 genes with non-fixed variants in these patients. These genes formed a viral signature which discriminated patients with cytomegalovirus infection who died from those that survived with high accuracy (AUC=0.96). Lymphocyte numbers for a subset of patients showed no recovery post-transplant in the patients who died. CONCLUSIONS: We hypothesise that the viral signature identified in this study may be a useful biomarker for poor response to antiviral drug treatment and indirectly for poor T cell function, potentially identifying early, those patients requiring non-pharmacological interventions.

Genomic and geographical structure of human cytomegalovirus.

Human cytomegalovirus (CMV) has infected humans since the origin of our species and currently infects most of the world's population. Variability between CMV genomes is the highest of any human herpesvirus, yet large portions of the genome are conserved. Here, we show that the genome encodes 74 regions of relatively high variability each with 2 to 8 alleles. We then identified two patterns in the CMV genome. Conserved parts of the genome and a minority (32) of variable regions show geographic population structure with evidence for African or European clustering, although hybrid strains are present. We find no evidence that geographic segregation has been driven by host immune pressure affecting known antigenic sites. Forty-two variable regions show no geographical structure, with similar allele distributions across different continental populations. These "nongeographical" regions are significantly enriched for genes encoding immunomodulatory functions suggesting a core functional importance. We hypothesize that at least two CMV founder populations account for the geographical differences that are largely seen in the conserved portions of the genome, although the timing of separation and direction of spread between the two are not clear. In contrast, the similar allele frequencies among 42 variable regions of the genome, irrespective of geographical origin, are indicative of a second evolutionary process, namely balancing selection that may preserve properties critical to CMV biological function. Given that genetic differences between CMVs are postulated to alter immunogenicity and potentially function, understanding these two evolutionary processes could contribute important information for the development of globally effective vaccines and the identification of novel drug targets.

Epstein-Barr Virus Reactivation After Paediatric Haematopoietic Stem Cell Transplantation: Risk Factors and Sensitivity Analysis of Mathematical Model.

Epstein-Barr virus (EBV) establishes a lifelong latent infection in healthy humans, kept under immune control by cytotoxic T cells (CTLs). Following paediatric haematopoetic stem cell transplantation (HSCT), a loss of immune surveillance leads to opportunistic outgrowth of EBV-infected cells, resulting in EBV reactivation, which can ultimately progress to post-transplant lymphoproliferative disorder (PTLD). The aims of this study were to identify risk factors for EBV reactivation in children in the first 100 days post-HSCT and to assess the suitability of a previously reported mathematical model to mechanistically model EBV reactivation kinetics in this cohort. Retrospective electronic data were collected from 56 children who underwent HSCT at Great Ormond Street Hospital (GOSH) between 2005 and 2016. Using EBV viral load (VL) measurements from weekly quantitative PCR (qPCR) monitoring post-HSCT, a multivariable Cox proportional hazards (Cox-PH) model was developed to assess time to first EBV reactivation event in the first 100 days post-HSCT. Sensitivity analysis of a previously reported mathematical model was performed to identify key parameters affecting EBV VL. Cox-PH modelling revealed EBV seropositivity of the HSCT recipient and administration of anti-thymocyte globulin (ATG) pre-HSCT to be significantly associated with an increased risk of EBV reactivation in the first 100 days post-HSCT (adjusted hazard ratio (AHR) = 2.32, P = 0.02; AHR = 2.55, P = 0.04). Five parameters were found to affect EBV VL in sensitivity analysis of the previously reported mathematical model. In conclusion, we have assessed the effect of multiple covariates on EBV reactivation in the first 100 days post-HSCT in children and have identified key parameters in a previously reported mechanistic mathematical model that affect EBV VL. Future work will aim to fit this model to patient EBV VLs, develop the model to account for interindividual variability and model the effect of clinically relevant covariates such as rituximab therapy and ATG on EBV VL.

HIV-1 Evolutionary Dynamics under Nonsuppressive Antiretroviral Therapy.

Prolonged virologic failure on 2nd-line protease inhibitor (PI)-based antiretroviral therapy (ART) without emergence of major protease mutations is well recognized and provides an opportunity to study within-host evolution in long-term viremic individuals. Using next-generation sequencing and in silico haplotype reconstruction, we analyzed whole-genome sequences from longitudinal plasma samples of eight chronically infected HIV-1-positive individuals failing 2nd-line regimens from the French National Agency for AIDS and Viral Hepatitis Research (ANRS) 12249 Treatment as Prevention (TasP) trial. On nonsuppressive ART, there were large fluctuations in synonymous and nonsynonymous variant frequencies despite stable viremia. Reconstructed haplotypes provided evidence for selective sweeps during periods of partial adherence, and viral haplotype competition, during periods of low drug exposure. Drug resistance mutations in reverse transcriptase (RT) were used as markers of viral haplotypes in the reservoir, and their distribution over time indicated recombination. We independently observed linkage disequilibrium decay, indicative of recombination. These data highlight dramatic changes in virus population structure that occur during stable viremia under nonsuppressive ART. IMPORTANCE HIV-1 infections are most commonly initiated with a single founder virus and are characterized by extensive inter- and intraparticipant genetic diversity. However, existing literature on HIV-1 intrahost population dynamics is largely limited to untreated infections, predominantly in subtype B-infected individuals. The manuscript characterizes viral population dynamics in long-term viremic treatment-experienced individuals, which has not been previously characterized. These data are particularly relevant for understanding HIV dynamics but can also be applied to other RNA viruses. With this unique data set we propose that the virus is highly unstable, and we have found compelling evidence of HIV-1 within-host viral diversification, recombination, and haplotype competition during nonsuppressive ART.

Altered TMPRSS2 usage by SARS-CoV-2 Omicron impacts infectivity and fusogenicity.

The SARS-CoV-2 Omicron BA.1 variant emerged in 20211 and has multiple mutations in its spike protein2. Here we show that the spike protein of Omicron has a higher affinity for ACE2 compared with Delta, and a marked change in its antigenicity increases Omicron's evasion of therapeutic monoclonal and vaccine-elicited polyclonal neutralizing antibodies after two doses. mRNA vaccination as a third vaccine dose rescues and broadens neutralization. Importantly, the antiviral drugs remdesivir and molnupiravir retain efficacy against Omicron BA.1. Replication was similar for Omicron and Delta virus isolates in human nasal epithelial cultures. However, in lung cells and gut cells, Omicron demonstrated lower replication. Omicron spike protein was less efficiently cleaved compared with Delta. The differences in replication were mapped to the entry efficiency of the virus on the basis of spike-pseudotyped virus assays. The defect in entry of Omicron pseudotyped virus to specific cell types effectively correlated with higher cellular RNA expression of TMPRSS2, and deletion of TMPRSS2 affected Delta entry to a greater extent than Omicron. Furthermore, drug inhibitors targeting specific entry pathways3 demonstrated that the Omicron spike inefficiently uses the cellular protease TMPRSS2, which promotes cell entry through plasma membrane fusion, with greater dependency on cell entry through the endocytic pathway. Consistent with suboptimal S1/S2 cleavage and inability to use TMPRSS2, syncytium formation by the Omicron spike was substantially impaired compared with the Delta spike. The less efficient spike cleavage of Omicron at S1/S2 is associated with a shift in cellular tropism away from TMPRSS2-expressing cells, with implications for altered pathogenesis.

Shrinking Weibel-Palade bodies prevents high platelet recruitment in assays using thrombotic thrombocytopenic purpura plasma.

BACKGROUND: Thrombotic thrombocytopenic purpura (TTP), caused by a genetic or autoimmune-driven lack of ADAMTS-13 activity, leads to high levels of the ultra-large von Willebrand factor (VWF) multimers produced by endothelial cells, causing excess platelet recruitment into forming thrombi, often with mortal consequences. Treatments include plasma infusion or replacement to restore ADAMTS-13 activity, or prevention of platelet recruitment to VWF. OBJECTIVES: We tested a different approach, exploiting the unique cell biology of the endothelium. Upon activation, the VWF released by exocytosis of Weibel-Palade bodies (WPBs), transiently anchored to the cell surface, unfurls as strings into flowing plasma, recruiting platelets. Using plasma from patients with TTP increases platelet recruitment to the surface of cultured endothelial cells under flow. WPBs are uniquely plastic, and shortening WPBs dramatically reduces VWF string lengths and the recruitment of platelets. We wished to test whether the TTP plasma-driven increase in platelet recruitment would be countered by reducing formation of the longest WPBs that release longer strings. METHODS: Endothelial cells grown in flow chambers were treated with fluvastatin, one of 37 drugs shown to shorten WPBs, then activated under flow in the presence of platelets and plasma of either controls or patients with TTP. RESULT: We found that the dramatic increase in platelet recruitment caused by TTP plasma is entirely countered by treatment with fluvastatin, shortening the WPBs. CONCLUSIONS: This potential approach of ameliorating the endothelial contribution to thrombotic risk by intervening far upstream of hemostasis might prove a useful adjunct to more conventional and direct therapies.

SARS-CoV-2 B.1.617.2 Delta variant replication and immune evasion.

The B.1.617.2 (Delta) variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first identified in the state of Maharashtra in late 2020 and spread throughout India, outcompeting pre-existing lineages including B.1.617.1 (Kappa) and B.1.1.7 (Alpha)1. In vitro, B.1.617.2 is sixfold less sensitive to serum neutralizing antibodies from recovered individuals, and eightfold less sensitive to vaccine-elicited antibodies, compared with wild-type Wuhan-1 bearing D614G. Serum neutralizing titres against B.1.617.2 were lower in ChAdOx1 vaccinees than in BNT162b2 vaccinees. B.1.617.2 spike pseudotyped viruses exhibited compromised sensitivity to monoclonal antibodies to the receptor-binding domain and the amino-terminal domain. B.1.617.2 demonstrated higher replication efficiency than B.1.1.7 in both airway organoid and human airway epithelial systems, associated with B.1.617.2 spike being in a predominantly cleaved state compared with B.1.1.7 spike. The B.1.617.2 spike protein was able to mediate highly efficient syncytium formation that was less sensitive to inhibition by neutralizing antibody, compared with that of wild-type spike. We also observed that B.1.617.2 had higher replication and spike-mediated entry than B.1.617.1, potentially explaining the B.1.617.2 dominance. In an analysis of more than 130 SARS-CoV-2-infected health care workers across three centres in India during a period of mixed lineage circulation, we observed reduced ChAdOx1 vaccine effectiveness against B.1.617.2 relative to non-B.1.617.2, with the caveat of possible residual confounding. Compromised vaccine efficacy against the highly fit and immune-evasive B.1.617.2 Delta variant warrants continued infection control measures in the post-vaccination era.