CpG-creating mutations are costly in many human viruses

Structural and functional characterization of the extended-diKH domain from the antiviral endoribonuclease KHNYN

Zinc finger antiviral protein (ZAP) binds CpG dinucleotides in viral RNA and targets them for decay. ZAP interacts with several cofactors to form the ZAP antiviral system, including KHNYN, a multidomain endoribonuclease required for ZAP-mediated RNA decay. However, it is unclear how the individual domains in KHNYN contribute to its activity. Here, we demonstrate that the KHNYN amino-terminal extended-diKH (ex-diKH) domain is required for antiviral activity and present its crystal structure. The structure belongs to a rare group of KH-containing domains, characterized by a noncanonical arrangement between two type 1 KH modules, with an additional helical bundle. N4BP1 is a KHNYN paralog with an ex-diKH domain that functionally complements the KHNYN ex-diKH domain. Interestingly, the ex-diKH domain structure is present in N4BP1-like proteins in lancelets, which are basal chordates, indicating that it is evolutionarily ancient. While many KH domains demonstrate RNA binding activity, biolayer interferometry and electrophoretic mobility shift assays indicate that the KHNYN ex-diKH domain does not bind RNA. Furthermore, residues required for canonical KH domains to bind RNA are not required for KHNYN antiviral activity. By contrast, an inter-KH domain cleft in KHNYN is a potential protein-protein interaction site, and mutations that eliminate arginine salt bridges at the edge of this cleft decrease KHNYN antiviral activity. This suggests that this domain could be a binding site for an unknown KHNYN cofactor.

Genomic Epidemiology of SARS-CoV-2 in Urban Settings in Senegal

We used whole genome sequencing to identify and analyze mutations in SARS-CoV-2 in urban settings during the deadliest wave of the COVID-19 epidemic-from March to April 2021-in Senegal. Nasopharyngeal samples testing positive for SARS-CoV-2 were sequenced on the Illumina NovaSeq 6000 sequencing system using the COVIDSeq protocol. A total of 291 genotypable consensus genome sequences were obtained. Phylogenetic analyses grouped the genomes into 16 distinct PANGOLIN lineages. The major lineage was B.1.1.420, despite circulation of the Alpha variant of concern (VOC). A total of 1125 different SNPs, identified relative to the Wuhan reference genome, were detected. These included 13 SNPs in non-coding regions. An average density of 37.2 SNPs per 1000 nucleotides was found, with the highest density occurring in ORF10. This analysis allowed, for the first time, the detection of a Senegalese SARS-CoV-2 strain belonging to the P.1.14 (GR/20J, Gamma V3) sublineage of the Brazilian P.1 lineage (or Gamma VOC). Overall, our results highlight substantial SARS-CoV-2 diversification in Senegal during the study period.

Assessing in vivo mutation frequencies and creating a high-resolution genome-wide map of fitness costs of Hepatitis C virus

Like many viruses, Hepatitis C Virus (HCV) has a high mutation rate, which helps the virus adapt quickly, but mutations come with fitness costs. Fitness costs can be studied by different approaches, such as experimental or frequency-based approaches. The frequency-based approach is particularly useful to estimate in vivo fitness costs, but this approach works best with deep sequencing data from many hosts are. In this study, we applied the frequency-based approach to a large dataset of 195 patients and estimated the fitness costs of mutations at 7957 sites along the HCV genome. We used beta regression and random forest models to better understand how different factors influenced fitness costs. Our results revealed that costs of nonsynonymous mutations were three times higher than those of synonymous mutations, and mutations at nucleotides A or T had higher costs than those at C or G. Genome location had a modest effect, with lower costs for mutations in HVR1 and higher costs for mutations in Core and NS5B. Resistance mutations were, on average, costlier than other mutations. Our results show that in vivo fitness costs of mutations can be site and virus specific, reinforcing the utility of constructing in vivo fitness cost maps of viral genomes.

S-farnesylation is essential for antiviral activity of the long ZAP isoform against RNA viruses with diverse replication strategies

The zinc finger antiviral protein (ZAP) is a broad inhibitor of virus replication. Its best-characterized function is to bind CpG dinucleotides present in viral RNAs and, through the recruitment of TRIM25, KHNYN and other cofactors, target them for degradation or prevent their translation. The long and short isoforms of ZAP (ZAP-L and ZAP-S) have different intracellular localization and it is unclear how this regulates their antiviral activity against viruses with different sites of replication. Using ZAP-sensitive and ZAP-insensitive human immunodeficiency virus type I (HIV-1), which transcribe the viral RNA in the nucleus and assemble virions at the plasma membrane, we show that the catalytically inactive poly-ADP-ribose polymerase (PARP) domain in ZAP-L is essential for CpG-specific viral restriction. Mutation of a crucial cysteine in the C-terminal CaaX box that mediates S-farnesylation and, to a lesser extent, the residues in place of the catalytic site triad within the PARP domain, disrupted the activity of ZAP-L. Addition of the CaaX box to ZAP-S partly restored antiviral activity, explaining why ZAP-S lacks antiviral activity for CpG-enriched HIV-1 despite conservation of the RNA-binding domain. Confocal microscopy confirmed the CaaX motif mediated localization of ZAP-L to vesicular structures and enhanced physical association with intracellular membranes. Importantly, the PARP domain and CaaX box together jointly modulate the interaction between ZAP-L and its cofactors TRIM25 and KHNYN, implying that its proper subcellular localisation is required to establish an antiviral complex. The essential contribution of the PARP domain and CaaX box to ZAP-L antiviral activity was further confirmed by inhibition of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication, which replicates in double-membrane vesicles derived from the endoplasmic reticulum. Thus, compartmentalization of ZAP-L on intracellular membranes provides an essential effector function in ZAP-L-mediated antiviral activity against divergent viruses with different subcellular replication sites.

Targeted Restriction of Viral Gene Expression and Replication by the ZAP Antiviral System

The zinc finger antiviral protein (ZAP) restricts the replication of a broad range of RNA and DNA viruses. ZAP directly binds viral RNA, targeting it for degradation and inhibiting its translation. While the full scope of RNA determinants involved in mediating selective ZAP activity are unclear, ZAP binds CpG dinucleotides, dictating at least part of its target specificity. ZAP interacts with many cellular proteins, although only a few have been demonstrated to be essential for its antiviral activity, including the 3'-5' exoribonuclease exosome complex, TRIM25, and KHNYN. In addition to inhibiting viral gene expression, ZAP also directly and indirectly targets a subset of cellular messenger RNAs to regulate the innate immune response. Overall, ZAP protects a cell from viral infection by restricting viral replication and regulating cellular gene expression. Further understanding of the ZAP antiviral system may allow for novel viral vaccine and anticancer therapy development. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Comparative Analysis of Within-Host Mutation Patterns and Diversity of Hepatitis C Virus Subtypes 1a, 1b, and 3a

Understanding within-host evolution is critical for predicting viral evolutionary outcomes, yet such studies are currently lacking due to difficulty involving human subjects. Hepatitis C virus (HCV) is an RNA virus with high mutation rates. Its complex evolutionary dynamics and extensive genetic diversity are demonstrated in over 67 known subtypes. In this study, we analyzed within-host mutation frequency patterns of three HCV subtypes, using a large number of samples obtained from treatment-naïve participants by next-generation sequencing. We report that overall mutation frequency patterns are similar among subtypes, yet subtype 3a consistently had lower mutation frequencies and nucleotide diversity, while subtype 1a had the highest. We found that about 50% of genomic sites are highly conserved across subtypes, which are likely under strong purifying selection. We also compared within-host and between-host selective pressures, which revealed that Hyper Variable Region 1 within hosts was under positive selection, but was under slightly negative selection between hosts, which indicates that many mutations created within hosts are removed during the transmission bottleneck. Examining the natural prevalence of known resistance-associated variants showed their consistent existence in the treatment-naïve participants. These results provide insights into the differences and similarities among HCV subtypes that may be used to develop and improve HCV therapies.

The Mutation Profile of SARS-CoV-2 Is Primarily Shaped by the Host Antiviral Defense

Understanding SARS-CoV-2 evolution is a fundamental effort in coping with the COVID-19 pandemic. The virus genomes have been broadly evolving due to the high number of infected hosts world-wide. Mutagenesis and selection are two inter-dependent mechanisms of virus diversification. However, which mechanisms contribute to the mutation profiles of SARS-CoV-2 remain under-explored. Here, we delineate the contribution of mutagenesis and selection to the genome diversity of SARS-CoV-2 isolates. We generated a comprehensive phylogenetic tree with representative genomes. Instead of counting mutations relative to the reference genome, we identified each mutation event at the nodes of the phylogenetic tree. With this approach, we obtained the mutation events that are independent of each other and generated the mutation profile of SARS-CoV-2 genomes. The results suggest that the heterogeneous mutation patterns are mainly reflections of host (i) antiviral mechanisms that are achieved through APOBEC, ADAR, and ZAP proteins, and (ii) probable adaptation against reactive oxygen species.

Biased Mutation and Selection in RNA Viruses

RNA viruses are responsible for some of the worst pandemics known to mankind, including outbreaks of Influenza, Ebola, and COVID-19. One major challenge in tackling RNA viruses is the fact they are extremely genetically diverse. Nevertheless, they share common features that include their dependence on host cells for replication, and high mutation rates. We set out to search for shared evolutionary characteristics that may aid in gaining a broader understanding of RNA virus evolution, and constructed a phylogeny-based data set spanning thousands of sequences from diverse single-stranded RNA viruses of animals. Strikingly, we found that the vast majority of these viruses have a skewed nucleotide composition, manifested as adenine rich (A-rich) coding sequences. In order to test whether A-richness is driven by selection or by biased mutation processes, we harnessed the effects of incomplete purifying selection at the tips of virus phylogenies. Our results revealed consistent mutational biases toward U rather than A in genomes of all viruses. In +ssRNA viruses, we found that this bias is compensated by selection against U and selection for A, which leads to A-rich genomes. In -ssRNA viruses, the genomic mutational bias toward U on the negative strand manifests as A-rich coding sequences, on the positive strand. We investigated possible reasons for the advantage of A-rich sequences including weakened RNA secondary structures, codon usage bias, and selection for a particular amino acid composition, and conclude that host immune pressures may have led to similar biases in coding sequence composition across very divergent RNA viruses.

Host-directed editing of the SARS-CoV-2 genome

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.