Coupling high-throughput genetics with phylogenetic information reveals an epistatic interaction on the influenza A virus M segment

Deep mutational scanning of influenza A virus neuraminidase facilitates the identification of drug resistance mutations in vivo

IMPORTANCE

NA is a crucial surface antigen and drug target of influenza A virus. A comprehensive understanding of NA's mutational effect and drug resistance profiles in vivo is essential for comprehending the evolutionary constraints and making informed choices regarding drug selection to combat resistance in clinical settings. In the current study, we established an efficient deep mutational screening system in mouse lung tissues and systematically evaluated the fitness effect and drug resistance to three neuraminidase inhibitors of NA single-nucleotide mutations. The fitness of NA mutants is generally correlated with a natural mutation in the database. The fitness of NA mutants is influenced by biophysical factors such as protein stability, complex formation, and the immune response triggered by viral infection. In addition to confirming previously reported drug-resistant mutations, novel mutations were identified. Interestingly, we identified an allosteric drug-resistance mutation that is not located within the drug-binding pocket but potentially affects drug binding by interfering with NA tetramerization. The dual assessments performed in this study provide a more accurate assessment of the evolutionary potential of drug-resistant mutations and offer guidance for the rational selection of antiviral drugs.

Antigenic evolution of human influenza H3N2 neuraminidase is constrained by charge balancing

As one of the main influenza antigens, neuraminidase (NA) in H3N2 virus has evolved extensively for more than 50 years due to continuous immune pressure. While NA has recently emerged as an effective vaccine target, biophysical constraints on the antigenic evolution of NA remain largely elusive. Here, we apply combinatorial mutagenesis and next-generation sequencing to characterize the local fitness landscape in an antigenic region of NA in six different human H3N2 strains that were isolated around 10 years apart. The local fitness landscape correlates well among strains and the pairwise epistasis is highly conserved. Our analysis further demonstrates that local net charge governs the pairwise epistasis in this antigenic region. In addition, we show that residue coevolution in this antigenic region is correlated with the pairwise epistasis between charge states. Overall, this study demonstrates the importance of quantifying epistasis and the underlying biophysical constraint for building a model of influenza evolution.

Applications of Deep Mutational Scanning in Virology

Several recently developed high-throughput techniques have changed the field of molecular virology. For example, proteomics studies reveal complete interactomes of a viral protein, genome-wide CRISPR knockout and activation screens probe the importance of every single human gene in aiding or fighting a virus, and ChIP-seq experiments reveal genome-wide epigenetic changes in response to infection. Deep mutational scanning is a relatively novel form of protein science which allows the in-depth functional analysis of every nucleotide within a viral gene or genome, revealing regions of importance, flexibility, and mutational potential. In this review, we discuss the application of this technique to RNA viruses including members of the Flaviviridae family, Influenza A Virus and Severe Acute Respiratory Syndrome Coronavirus 2. We also briefly discuss the reverse genetics systems which allow for analysis of viral replication cycles, next-generation sequencing technologies and the bioinformatics tools that facilitate this research.

Adaptation of Oxford Nanopore technology for hepatitis C whole genome sequencing and identification of within-host viral variants

Background

Hepatitis C (HCV) and many other RNA viruses exist as rapidly mutating quasi-species populations in a single infected host. High throughput characterization of full genome, within-host variants is still not possible despite advances in next generation sequencing. This limitation constrains viral genomic studies that depend on accurate identification of hemi-genome or whole genome, within-host variants, especially those occurring at low frequencies. With the advent of third generation long read sequencing technologies, including Oxford Nanopore Technology (ONT) and PacBio platforms, this problem is potentially surmountable. ONT is particularly attractive in this regard due to the portable nature of the MinION sequencer, which makes real-time sequencing in remote and resource-limited locations possible. However, this technology (termed here ‘nanopore sequencing’) has a comparatively high technical error rate. The present study aimed to assess the utility, accuracy and cost-effectiveness of nanopore sequencing for HCV genomes. We also introduce a new bioinformatics tool (Nano-Q) to differentiate within-host variants from nanopore sequencing.

Results

The Nanopore platform, when the coverage exceeded 300 reads, generated comparable consensus sequences to Illumina sequencing. Using HCV Envelope plasmids (~ 1800 nt) mixed in known proportions, the capacity of nanopore sequencing to reliably identify variants with an abundance as low as 0.1% was demonstrated, provided the autologous reference sequence was available to identify the matching reads. Successful pooling and nanopore sequencing of 52 samples from patients with HCV infection demonstrated its cost effectiveness (AUD$ 43 per sample with nanopore sequencing versus $100 with paired-end short read technology). The Nano-Q tool successfully separated between-host sequences, including those from the same subtype, by bulk sorting and phylogenetic clustering without an autologous reference sequence (using only a subtype-specific generic reference). The pipeline also identified within-host viral variants and their abundance when the parameters were appropriately adjusted.

Conclusion

Cost effective HCV whole genome sequencing and within-host variant identification without haplotype reconstruction are potential advantages of nanopore sequencing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07460-1.

MaveDB: an open-source platform to distribute and interpret data from multiplexed assays of variant effect

Multiplex assays of variant effect (MAVEs), such as deep mutational scans and massively parallel reporter assays, test thousands of sequence variants in a single experiment. Despite the importance of MAVE data for basic and clinical research, there is no standard resource for their discovery and distribution. Here, we present MaveDB ( https://www.mavedb.org ), a public repository for large-scale measurements of sequence variant impact, designed for interoperability with applications to interpret these datasets. We also describe the first such application, MaveVis, which retrieves, visualizes, and contextualizes variant effect maps. Together, the database and applications will empower the community to mine these powerful datasets.

Contrasting selective patterns across the segmented genome of bluetongue virus in a global reassortment hotspot

For segmented viruses, rapid genomic and phenotypic changes can occur through the process of reassortment, whereby co-infecting strains exchange entire segments creating novel progeny virus genotypes. However, for many viruses with segmented genomes, this process and its effect on transmission dynamics remain poorly understood. Here, we assessed the consequences of reassortment for selection on viral diversity through time using bluetongue virus (BTV), a segmented arbovirus that is the causative agent of a major disease of ruminants. We analysed ninety-two BTV genomes isolated across four decades from India, where BTV diversity, and thus opportunities for reassortment, are among the highest in the world. Our results point to frequent reassortment and segment turnover, some of which appear to be driven by selective sweeps and serial hitchhiking. Particularly, we found evidence for a recent selective sweep affecting segment 5 and its encoded NS1 protein that has allowed a single variant to essentially invade the full range of BTV genomic backgrounds and serotypes currently circulating in India. In contrast, diversifying selection was found to play an important role in maintaining genetic diversity in genes encoding outer surface proteins involved in virus interactions (VP2 and VP5, encoded by segments 2 and 6, respectively). Our results support the role of reassortment in driving rapid phenotypic change in segmented viruses and generate testable hypotheses for in vitro experiments aiming at understanding the specific mechanisms underlying differences in fitness and selection across viral genomes.

Deep Mutational Scan of the Highly Conserved Influenza A Virus M1 Matrix Protein Reveals Substantial Intrinsic Mutational Tolerance

Influenza A virus matrix protein M1 is involved in multiple stages of the viral infectious cycle. Despite its functional importance, our present understanding of this essential viral protein is limited. The roles of a small subset of specific amino acids have been reported, but a more comprehensive understanding of the relationship between M1 sequence, structure, and virus fitness remains elusive. In this study, we used deep mutational scanning to measure the effect of every amino acid substitution in M1 on viral replication in cell culture. The map of amino acid mutational tolerance we have generated allows us to identify sites that are functionally constrained in cell culture as well as sites that are less constrained. Several sites that exhibit low tolerance to mutation have been found to be critical for M1 function and production of viable virions. Surprisingly, significant portions of the M1 sequence, especially in the C-terminal domain, whose structure is undetermined, were found to be highly tolerant of amino acid variation, despite having extremely low levels of sequence diversity among natural influenza virus strains. This unexpected discrepancy indicates that not all sites in M1 that exhibit high sequence conservation in nature are under strong constraint during selection for viral replication in cell culture.IMPORTANCE The M1 matrix protein is critical for many stages of the influenza virus infection cycle. Currently, we have an incomplete understanding of this highly conserved protein's function and structure. Key regions of M1, particularly in the C terminus of the protein, remain poorly characterized. In this study, we used deep mutational scanning to determine the extent of M1's tolerance to mutation. Surprisingly, nearly two-thirds of the M1 sequence exhibits a high tolerance for substitutions, contrary to the extremely low sequence diversity observed across naturally occurring M1 isolates. Sites with low mutational tolerance were also identified, suggesting that they likely play critical functional roles and are under selective pressure. These results reveal the intrinsic mutational tolerance throughout M1 and shape future inquiries probing the functions of this essential influenza A virus protein.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

Mutation and Epistasis in Influenza Virus Evolution

Influenza remains a persistent public health challenge, because the rapid evolution of influenza viruses has led to marginal vaccine efficacy, antiviral resistance, and the annual emergence of novel strains. This evolvability is driven, in part, by the virus's capacity to generate diversity through mutation and reassortment. Because many new traits require multiple mutations and mutations are frequently combined by reassortment, epistatic interactions between mutations play an important role in influenza virus evolution. While mutation and epistasis are fundamental to the adaptability of influenza viruses, they also constrain the evolutionary process in important ways. Here, we review recent work on mutational effects and epistasis in influenza viruses.

Genome-wide identification of interferon-sensitive mutations enables influenza vaccine design

In conventional attenuated viral vaccines, immunogenicity is often suboptimal. Here we present a systematic approach for vaccine development that eliminates interferon (IFN)-modulating functions genome-wide while maintaining virus replication fitness. We applied a quantitative high-throughput genomics system to influenza A virus that simultaneously measured the replication fitness and IFN sensitivity of mutations across the entire genome. By incorporating eight IFN-sensitive mutations, we generated a hyper-interferon-sensitive (HIS) virus as a vaccine candidate. HIS virus is highly attenuated in IFN-competent hosts but able to induce transient IFN responses, elicits robust humoral and cellular immune responses, and provides protection against homologous and heterologous viral challenges. Our approach, which attenuates the virus and promotes immune responses concurrently, is broadly applicable for vaccine development against other pathogens.

Mapping the Evolutionary Potential of RNA Viruses

The deterministic force of natural selection and stochastic influence of drift shape RNA virus evolution. New deep-sequencing and microfluidics technologies allow us to quantify the effect of mutations and trace the evolution of viral populations with single-genome and single-nucleotide resolution. Such experiments can reveal the topography of the genotype-fitness landscapes that shape the path of viral evolution. By combining historical analyses, like phylogenetic approaches, with high-throughput and high-resolution evolutionary experiments, we can observe parallel patterns of evolution that drive important phenotypic transitions. These developments provide a framework for quantifying and anticipating potential evolutionary events. Here, we examine emerging technologies that can map the selective landscapes of viruses, focusing on their application to pathogenic viruses. We identify areas where these technologies can bolster our ability to study the evolution of viruses and to anticipate and possibly intervene in evolutionary events and prevent viral disease.

High-Throughput Fitness Profiling of Zika Virus E Protein Reveals Different Roles for Glycosylation during Infection of Mammalian and Mosquito Cells

Zika virus (ZIKV) infection causes Guillain-Barré syndrome and severe birth defects. ZIKV envelope (E) protein is the major viral protein involved in cell receptor binding and entry and is therefore considered one of the major determinants in ZIKV pathogenesis. Here we report a gene-wide mapping of functional residues of ZIKV E protein using a mutant library, with changes covering every nucleotide position. By comparing the replication fitness of every viral mutant between mosquito and human cells, we identified that mutations affecting glycosylation display the most divergence. By characterizing individual mutants, we show that ablation of glycosylation selectively benefits ZIKV infection of mosquito cells by enhancing cell entry, whereas it either has little impact on ZIKV infection on certain human cells or leads to decreased infection through the entry factor DC-SIGN. In conclusion, we define the roles of individual residues of ZIKV envelope protein, which contribute to ZIKV replication fitness in human and mosquito cells.

Comprehensive Mapping of HIV-1 Escape from a Broadly Neutralizing Antibody

Precisely defining how viral mutations affect HIV's sensitivity to antibodies is vital to develop and evaluate vaccines and antibody immunotherapeutics. Despite great effort, a full map of escape mutants has not been delineated for an anti-HIV antibody. We describe a massively parallel experimental approach to quantify how all single amino acid mutations to HIV Envelope (Env) affect neutralizing antibody sensitivity in the context of replication-competent virus. We apply this approach to PGT151, a broadly neutralizing antibody recognizing a combination of Env residues and glycans. We confirm sites previously defined by structural and functional studies and reveal additional sites of escape, such as positively charged mutations in the antibody-Env interface. Evaluating the effect of each amino acid at each site lends insight into biochemical mechanisms of escape throughout the epitope, highlighting roles for charge-charge repulsions. Thus, comprehensively mapping HIV antibody escape gives a quantitative, mutation-level view of Env evasion of neutralization.

Deep mutational scanning identifies sites in influenza nucleoprotein that affect viral inhibition by MxA

The innate-immune restriction factor MxA inhibits influenza replication by targeting the viral nucleoprotein (NP). Human influenza virus is more resistant than avian influenza virus to inhibition by human MxA, and prior work has compared human and avian viral strains to identify amino-acid differences in NP that affect sensitivity to MxA. However, this strategy is limited to identifying sites in NP where mutations that affect MxA sensitivity have fixed during the small number of documented zoonotic transmissions of influenza to humans. Here we use an unbiased deep mutational scanning approach to quantify how all single amino-acid mutations to NP affect MxA sensitivity in the context of replication-competent virus. We both identify new sites in NP where mutations affect MxA resistance and re-identify mutations known to have increased MxA resistance during historical adaptations of influenza to humans. Most of the sites where mutations have the greatest effect are almost completely conserved across all influenza A viruses, and the amino acids at these sites confer relatively high resistance to MxA. These sites cluster in regions of NP that appear to be important for its recognition by MxA. Overall, our work systematically identifies the sites in influenza nucleoprotein where mutations affect sensitivity to MxA. We also demonstrate a powerful new strategy for identifying regions of viral proteins that affect inhibition by host factors.

During viral infection, human cells express proteins that can restrict virus replication. However, in many cases it remains unclear what determines the sensitivity of a given viral strain to a particular restriction factor. Here we use a high-throughput approach to measure how all amino-acid mutations to the nucleoprotein of influenza virus affect restriction by the human protein MxA. We find several dozen sites where mutations substantially affect the sensitivity of influenza virus to MxA. While a few of these sites are known to have fixed mutations during past adaptations of influenza virus to humans, most of the sites are broadly conserved across all influenza strains and have never previously been described as affecting MxA resistance. Our results therefore show that the known historical evolution of influenza has only involved substitutions at a small fraction of the sites where mutations can in principle affect MxA resistance. We suggest that this is because many sites are already broadly fixed at amino acids that confer high resistance.

Complete mapping of viral escape from neutralizing antibodies

Identifying viral mutations that confer escape from antibodies is crucial for understanding the interplay between immunity and viral evolution. We describe a high-throughput approach to quantify the selection that monoclonal antibodies exert on all single amino-acid mutations to a viral protein. This approach, mutational antigenic profiling, involves creating all replication-competent protein variants of a virus, selecting with antibody, and using deep sequencing to identify enriched mutations. We use mutational antigenic profiling to comprehensively identify mutations that enable influenza virus to escape four monoclonal antibodies targeting hemagglutinin, and validate key findings with neutralization assays. We find remarkable mutation-level idiosyncrasy in antibody escape: for instance, at a single residue targeted by two antibodies, some mutations escape both antibodies while other mutations escape only one or the other. Because mutational antigenic profiling rapidly maps all mutations selected by an antibody, it is useful for elucidating immune specificities and interpreting the antigenic consequences of viral genetic variation.

Limited naturally occurring escape in broadly neutralizing antibody epitopes in hepatitis C glycoprotein E2 and constrained sequence usage in acute infection

Broadly neutralizing antibodies have been associated with spontaneous clearance of the hepatitis C infection as well as viral persistence by immune escape. Further study of neutralizing antibody epitopes is needed to unravel pathways of resistance to virus neutralization, and to identify conserved regions for vaccine design. All reported broadly neutralizing antibody (BNAb) epitopes in the HCV Envelope (E2) glycoprotein were identified. The critical contact residues of these epitopes were mapped onto the linear E2 sequence. All publicly available E2 sequences were then downloaded and the contact residues within the BNAb epitopes were assessed for the level of conservation, as well as the frequency of occurrence of experimentally-proven resistance mutations. Epitopes were also compared between two sequence datasets obtained from samples collected at well-defined time points from acute (<180days) and chronic (>180days) infections, to identify any significant differences in residue usage. The contact residues for all BNAbs were contained within 3 linear regions of the E2 protein sequence. An analysis of 1749 full length E2 sequences from public databases showed that only 10 out of 29 experimentally-proven resistance mutations were present at a frequency >5%. Comparison of subtype 1a viral sequences obtained from samples collected during acute or chronic infection revealed significant differences at positions 610 and 655 with changes in residue (p<0.05), and at position 422 (p<0.001) with a significant difference in variability (entropy). The majority of experimentally-described escape variants do not occur frequently in nature. The observed differences between acute and chronically isolated sequences suggest constraints on residue usage early in infection.

Viral factors in influenza pandemic risk assessment

The threat of an influenza A virus pandemic stems from continual virus spillovers from reservoir species, a tiny fraction of which spark sustained transmission in humans. To date, no pandemic emergence of a new influenza strain has been preceded by detection of a closely related precursor in an animal or human. Nonetheless, influenza surveillance efforts are expanding, prompting a need for tools to assess the pandemic risk posed by a detected virus. The goal would be to use genetic sequence and/or biological assays of viral traits to identify those non-human influenza viruses with the greatest risk of evolving into pandemic threats, and/or to understand drivers of such evolution, to prioritize pandemic prevention or response measures. We describe such efforts, identify progress and ongoing challenges, and discuss three specific traits of influenza viruses (hemagglutinin receptor binding specificity, hemagglutinin pH of activation, and polymerase complex efficiency) that contribute to pandemic risk.

Accurate Measurement of the Effects of All Amino-Acid Mutations on Influenza Hemagglutinin

Influenza genes evolve mostly via point mutations, and so knowing the effect of every amino-acid mutation provides information about evolutionary paths available to the virus. We and others have combined high-throughput mutagenesis with deep sequencing to estimate the effects of large numbers of mutations to influenza genes. However, these measurements have suffered from substantial experimental noise due to a variety of technical problems, the most prominent of which is bottlenecking during the generation of mutant viruses from plasmids. Here we describe advances that ameliorate these problems, enabling us to measure with greatly improved accuracy and reproducibility the effects of all amino-acid mutations to an H1 influenza hemagglutinin on viral replication in cell culture. The largest improvements come from using a helper virus to reduce bottlenecks when generating viruses from plasmids. Our measurements confirm at much higher resolution the results of previous studies suggesting that antigenic sites on the globular head of hemagglutinin are highly tolerant of mutations. We also show that other regions of hemagglutinin—including the stalk epitopes targeted by broadly neutralizing antibodies—have a much lower inherent capacity to tolerate point mutations. The ability to accurately measure the effects of all influenza mutations should enhance efforts to understand and predict viral evolution.