Generation and characterization of interferon-lambda 1-resistant H1N1 influenza A viruses

NS7a of SADS-CoV promotes viral infection via inducing apoptosis to suppress type III interferon production

Swine acute diarrhea syndrome coronavirus (SADS-CoV) is a newly discovered swine coronavirus with potential cross-species transmission risk. Although SADS-CoV-induced host cell apoptosis and innate immunity antagonization has been revealed, underlying signaling pathways remain obscure. Here, we demonstrated that infection of SADS-CoV induced apoptosis in vivo and in vitro, and that viral protein NS7a is mainly responsible for SADS-CoV-induced apoptosis in host cells. Furthermore, we found that NS7a interacted with apoptosis-inducing factor mitochondria associated 1 (AIFM1) to activate caspase-3 via caspase-6 in SADS-CoV-infected cells, and enhanced SADS-CoV replication. Importantly, NS7a suppressed poly(I:C)-induced expression of type III interferon (IFN-λ) via activating caspase-3 to cleave interferon regulatory factor 3 (IRF3), and caspase-3 inhibitor protects piglets against SADS-CoV infection in vivo. These findings reveal how SADS-CoV induced apoptosis to inhibit innate immunity and provide a valuable clue to the development of effective drugs for the clinical control of SADS-CoV infection.IMPORTANCEOver the last 20 years, multiple animal-originated coronaviruses, including severe acute respiratory syndrome coronavirus (SARS-CoV), middle east respiratory syndrome coronavirus (MERS-CoV), and SARS-CoV-2, have caused millions of deaths, seriously jeopardized human health, and hindered social development, indicating that the study of animal-originated coronaviruses with potential for cross-species transmission is particularly important. Bat-originated swine acute diarrhea syndrome coronavirus (SADS-CoV), discovered in 2017, can not only cause fatal diarrhea in piglets, but also infect multiple human cells, with a potential risk of cross-species transmission, but its pathogenesis is unclear. In this study, we demonstrated that NS7a of SADS-CoV suppresses IFN-λ production via apoptosis-inducing factor mitochondria associated 1 (AIFM1)-caspase-6-caspase-3-interferon regulatory factor 3 (IRF3) pathway, and caspase-3 inhibitor (Z-DEVD-FMK) can effectively inhibit SADS-CoV replication and protect infected piglets. Our findings in this study contribute to a better understanding of SADS-CoV-host interactions as a part of the coronaviruses pathogenesis and using apoptosis-inhibitor as a drug as potential therapeutic approaches for prevention and control of SADS-CoV infection.

Genomic Analysis of Influenza A and B Viruses Carrying Baloxavir Resistance-Associated Substitutions Serially Passaged in Human Epithelial Cells

Baloxavir marboxil (baloxavir) is an FDA-approved inhibitor of the influenza virus polymerase acidic (PA) protein. Here, we used next-generation sequencing to compare the genomic mutational profiles of IAV H1N1 and H3N2, and IBV wild type (WT) and mutants (MUT) viruses carrying baloxavir resistance-associated substitutions (H1N1-PA I38L, I38T, and E199D; H3N2-PA I38T; and IBV-PA I38T) during passaging in normal human bronchial epithelial (NHBE) cells. We determined the ratio of nonsynonymous to synonymous nucleotide mutations (dN/dS) and identified the location and type of amino acid (AA) substitutions that occurred at a frequency of ≥30%. We observed that IAV H1N1 WT and MUT viruses remained relatively stable during passaging. While the mutational profiles for IAV H1N1 I38L, I38T, and E199D, and IBV I38T MUTs were relatively similar after each passage compared to the respective WTs, the mutational profile of the IAV H3N2 I38T MUT was significantly different for most genes compared to H3N2 WT. Our work provides insight into how baloxavir resistance-associated substitutions may impact influenza virus evolution in natural settings. Further characterization of the potentially adaptive mutations identified in this study is needed.

Interchangeability of the Assays Used to Assess the Activity of Anti-SARS-CoV-2 Monoclonal Antibodies

The recent global COVID-19 pandemic caused by SARS-CoV-2 lasted for over three years. A key measure in combatting this pandemic involved the measurement of the monoclonal antibody (mAb)-mediated inhibition of binding between the spike receptor-binding domain (RBD) and hACE2 receptor. Potency assessments of therapeutic anti-SARS-CoV-2 mAbs typically include binding or cell-based neutralization assays. We assessed the inhibitory activity of five anti-SARS-CoV-2 mAbs using ELISA, surface plasmon resonance (SPR), and four cell-based neutralization assays using different pseudovirus particles and 293T or A549 cells expressing hACE2 with or without TMPRSS2. We assessed the interchangeability between cell-based and binding assays by applying the Bland-Altman method under certain assumptions. Our data demonstrated that the IC50 [nM] values determined by eight neutralization assays are independent of the cell line, presence of TMPRSS2 enzyme on the cell surface, and pseudovirus backbone used. Moreover, the Bland-Altman analysis showed that the IC50 [nM] and KD [nM] values determined by neutralization/ELISA or by SPR are equivalent and that the anti-spike mAb activity can be attributed to one variable directly related to its tertiary conformational structure conformation, rate dissociation constant Koff. This parameter is independent from the concentrations of the components of the mAb:RBD:hACE2 complexes and can be used for a comparison between the activities of the different mAbs.

Impact of Baloxavir Resistance-Associated Substitutions on Influenza Virus Growth and Drug Susceptibility

Baloxavir marboxil (baloxavir) is a recently FDA-approved influenza virus polymerase acidic (PA) endonuclease inhibitor. Several PA substitutions have been demonstrated to confer reduced susceptibility to baloxavir; however, their impacts on measurements of antiviral drug susceptibility and replication capacity when present as a fraction of the viral population have not been established. We generated recombinant A/California/04/09 (H1N1)-like viruses (IAV) with PA I38L, I38T, or E199D substitutions and B/Victoria/504/2000-like virus (IBV) with PA I38T. These substitutions reduced baloxavir susceptibility by 15.3-, 72.3-, 5.4-, and 54.5-fold, respectively, when tested in normal human bronchial epithelial (NHBE) cells. We then assessed the replication kinetics, polymerase activity, and baloxavir susceptibility of the wild-type:mutant (WT:MUT) virus mixtures in NHBE cells. The percentage of MUT relative to WT virus necessary to detect reduced baloxavir susceptibility in phenotypic assays ranged from 10% (IBV I38T) to 92% (IAV E199D). While I38T did not alter IAV replication kinetics or polymerase activity, IAV PA I38L and E199D MUTs and the IBV PA I38T MUT exhibited reduced replication levels and significantly altered polymerase activity. Differences in replication were detectable when the MUTs comprised ≥90%, ≥90%, or ≥75% of the population, respectively. Droplet digital PCR (ddPCR) and next-generation sequencing (NGS) analyses showed that WT viruses generally outcompeted the respective MUTs after multiple replication cycles and serial passaging in NHBE cells when initial mixtures contained ≥50% of the WT viruses; however, we also identified potential compensatory substitutions (IAV PA D394N and IBV PA E329G) that emerged and appeared to improve the replication capacity of baloxavir-resistant virus in cell culture. IMPORTANCE Baloxavir marboxil, an influenza virus polymerase acidic endonuclease inhibitor, represents a recently approved new class of influenza antivirals. Treatment-emergent resistance to baloxavir has been observed in clinical trials, and the potential spread of resistant variants could diminish baloxavir effectiveness. Here, we report the impact of the proportion of drug-resistant subpopulations on the ability to detect resistance in clinical isolates and the impact of substitutions on viral replication of mixtures containing both drug-sensitive and drug-resistant variants. We also show that ddPCR and NGS methods can be successfully used for detection of resistant subpopulations in clinical isolates and to quantify their relative abundance. Taken together, our data shed light on the potential impact of baloxavir-resistant I38T/L and E199D substitutions on baloxavir susceptibility and other biological properties of influenza virus and the ability to detect resistance in phenotypic and genotypic assays.

Method for screening influenza neutralizing antibodies in crude human plasma and its derivatives using SPR

We applied Surface Plasmon Resonance (SPR) technology to develop a method for potency screening and quantification of anti-influenza antibodies in minimally processed human plasma samples and intravenous immunoglobulin (IGIV) products. We found that specific antibodies in human plasma or IGIV capable of inhibiting binding of influenza hemagglutinin to receptor-analogous glycans do so in concentration-dependent manner. We ranked the inhibitory activity of plasma samples from multiple donors and found a good correlation (r = 0.87) of SPR assay measurements and conventional hemagglutination inhibition (HAI) assay results. This method was also applied to screen for specific anti-influenza antibodies in IGIV lots manufactured pre- and post-2009 H1N1 pandemic. The SPR method was also applied to study binding inhibition of the intact A/California/04/2009 H1N1 and B/Victoria/504/2000 influenza viruses to α2,6 or α2,3-linked synthetic glycans. In contrast to recombinant H1 hemagglutinin, which was found to interact primarily with α2,6-linked terminal sialic acids, intact H1N1 or influenza B virus recognized both types of receptor analogs with different observed dissociation rates and the inhibitory activity of plasma antibodies was dependent on the type of sialic acid link. The SPR method can provide a high-throughput, time-saving and semi-automated alternative to conventional assays such as HAI or microneutralization in situations where screening of large numbers of plasma donations to identify high titer units is needed to product highly potent immunoglobulins.

Modeling an immune response to influenza A virus infection in alveolar epithelial cells

Influenza A viruses (IAV) have been the cause of several influenza pandemics in history and are a significant threat for the next global pandemic. Hospitalized influenza patients often have excess interferon production and a dysregulated immune response to the IAV infection. Obtaining a better understanding of the mechanisms of IAV infection that induce these harmful effects would help drug developers and health professionals create more effective treatments for IAV infection and improve patient outcomes. IAV stimulates viral sensors and receptors expressed by alveolar epithelial cells, like RIG-I and toll-like receptor 3 (TLR3). These two pathways coordinate with one another to induce expression of type III interferons to combat the infection. Presented here is a queuing theory-based model of these pathways that was designed to analyze the timing and amount of interferons produced in response to IAV single stranded RNA and double-stranded RNA detection. The model accurately represents biological data showing the necessary coordination of the RIG-I and TLR3 pathways for effective interferon production. This model can serve as the framework for future studies of IAV infection and identify new targets for potential treatments.

An anti-influenza combined therapy assessed by single cell RNA-sequencing

Influenza makes millions of people ill every year, placing a large burden on the healthcare system and the economy. To develop a treatment against influenza, we combined virucidal sialylated cyclodextrins with interferon lambda and demonstrated, in human airway epithelia, that the two compounds inhibit the replication of a clinical H1N1 strain more efficiently when administered together rather than alone. We investigated the mechanism of action of the combined treatment by single cell RNA-sequencing analysis and found that both the single and combined treatments impair viral replication to different extents across distinct epithelial cell types. We showed that each cell type comprises multiple sub-types, whose proportions are altered by H1N1 infection, and assessed the ability of the treatments to restore them. To the best of our knowledge this is the first study investigating the effectiveness of an antiviral therapy against influenza virus by single cell transcriptomic studies.

When combined with interferon lambda, virucidal sialylated cyclodextrins inhibit the replication of a clinical H1N1 influenza strain in ex vivo human airway epithelia more efficiently than when delivered alone.

Pleiotropic Effects of Influenza H1, H3, and B Baloxavir-Resistant Substitutions on Replication, Sensitivity to Baloxavir, and Interferon Expression

Baloxavir is an anti-influenza endonuclease inhibitor that targets the polymerase acidic (PA) protein of influenza A and B viruses. Our knowledge regarding the pleiotropic effects of baloxavir resistance-associated substitutions is limited. We generated recombinant A/California/04/09 (H1N1)-, A/Hong Kong/218849/2006 (H3N2)-, and B/Victoria/504/2000-like viruses that contained PA substitutions identified in baloxavir clinical trials and surveillance that could potentially be associated with baloxavir resistance. We characterized their susceptibility to baloxavir, impact on polymerase activity, viral growth, and ability to induce interferon (IFN) and IFN-stimulated genes expression in vitro. Four PA substitutions, H1N1 I38L/T, E199D, and B G199R, significantly reduced the sensitivity of the recombinant viruses to baloxavir (14.1-fold). We confirmed our findings by using the luciferase-based ribonucleoprotein minigenome assay and by using virus yield reduction assay in Calu-3 and normal human bronchial epithelial (NHBE) cells. We observed that I38L and E199D resulted in decreased viral replication of the H1N1 wild-type virus (1.4-fold) but the H1N1 I38T and B G199R substitutions did not significantly alter replication capacity in Calu-3 cells. In addition, H1N1 variants with PA I38L/T and E199D induced significantly higher levels of IFNB1 gene expression compared to the wild-type virus (4.2-fold). In contrast, the B variant, G199R, triggered the lowest levels of IFN genes in Calu-3 cells (1.6-fold). Because baloxavir is a novel anti-influenza therapeutic agent, identifying and characterizing substitutions associated with reduced sensitivity to baloxavir, as well as the impact of these substitutions on viral fitness, is paramount to the strategic implementation of this novel countermeasure.

Interferon Lambda Delays the Emergence of Influenza Virus Resistance to Oseltamivir

Influenza viruses are a leading cause of morbidity and mortality worldwide. These air-borne pathogens are able to cross the species barrier, leading to regular seasonal epidemics and sporadic pandemics. Influenza viruses also possess a high genetic variability, which allows for the acquisition of resistance mutations to antivirals. Combination therapies with two or more drugs targeting different mechanisms of viral replication have been considered an advantageous option to not only enhance the effectiveness of the individual treatments, but also reduce the likelihood of resistance emergence. Using an in vitro infection model, we assessed the barrier to viral resistance of a combination therapy with the neuraminidase inhibitor oseltamivir and human interferon lambda against the pandemic H1N1 A/Netherlands/602/2009 (H1N1pdm09) virus. We serially passaged the virus in a cell line derived from human bronchial epithelial cells in the presence or absence of increasing concentrations of oseltamivir alone or oseltamivir plus interferon lambda. While the treatment with oseltamivir alone quickly induced the emergence of antiviral resistance through a single mutation in the neuraminidase gene, the co-administration of interferon lambda delayed the emergence of drug-resistant influenza virus variants. Our results suggest a possible clinical application of interferon lambda in combination with oseltamivir to treat influenza.

Adaptation of influenza B virus by serial passage in human airway epithelial cells

Influenza B viruses cause seasonal epidemics and are a considerable burden to public health. To understand their adaptation capability, we examined the genetic changes that occurred following 15 serial passages of two influenza B viruses, B/Brisbane/60/2008 and B/Victoria/504/2000, in human epithelial cells. Thirteen distinct amino acid mutations were found in the PB1, PA, hemagglutinin (HA), neuraminidase (NA), and M proteins after serial passage in the human lung epithelial cell line, Calu-3, and normal human bronchial epithelial (NHBE) cells. These changes were associated with significantly decreased viral replication levels. Our results demonstrate that adaptation of influenza B viruses for growth in human airway epithelial cells is partially conferred by selection of HA1, NA, and polymerase mutations that regulate receptor specificity, functional compatibility with the HA protein, and polymerase activity, respectively.

Laninamivir-Interferon Lambda 1 Combination Treatment Promotes Resistance by Influenza A Virus More Rapidly than Laninamivir Alone

Each year, 5% to 20% of the population of the United States becomes infected with influenza A virus. Combination therapy with two or more antiviral agents has been considered a potential treatment option for influenza virus infection. However, the clinical results derived from combination treatment with two or more antiviral drugs have been variable. We examined the effectiveness of cotreatment with two distinct classes of anti-influenza drugs, i.e., neuraminidase (NA) inhibitor, laninamivir, and interferon lambda 1 (IFN-λ1), against the emergence of drug-resistant virus variants in vitro We serially passaged pandemic A/California/04/09 [A(H1N1)pdm09] influenza virus in a human lung epithelial cell line (Calu-3) in the presence or absence of increasing concentrations of laninamivir or laninamivir plus IFN-λ1. Surprisingly, laninamivir used in combination with IFN-λ1 promoted the emergence of the E119G NA mutation five passages earlier than laninamivir alone (passage 2 versus passage 7, respectively). Acquisition of this mutation resulted in significantly reduced sensitivity to the NA inhibitors laninamivir (∼284-fold) and zanamivir (∼1,024-fold) and decreased NA enzyme catalytic activity (∼5-fold) compared to the parental virus. Moreover, the E119G NA mutation emerged together with concomitant hemagglutinin (HA) mutations (T197A and D222G), which were selected more rapidly by combination treatment with laninamivir plus IFN-λ1 (passages 2 and 3, respectively) than by laninamivir alone (passage 10). Our results show that treatment with laninamivir alone or in combination with IFN-λ1 can lead to the emergence of drug-resistant influenza virus variants. The addition of IFN-λ1 in combination with laninamivir may promote acquisition of drug resistance more rapidly than treatment with laninamivir alone.

Macrophage Coordination of the Interferon Lambda Immune Response

Lambda interferons (IFN-λs) are a major component of the innate immune defense to viruses, bacteria, and fungi. In human liver, IFN-λ not only drives antiviral responses, but also promotes inflammation and fibrosis in viral and non-viral diseases. Here we demonstrate that macrophages are primary responders to IFN-λ, uniquely positioned to bridge the gap between IFN-λ producing cells and lymphocyte populations that are not intrinsically responsive to IFN-λ. While CD14⁺ monocytes do not express the IFN-λ receptor, IFNLR1, sensitivity is quickly gained upon differentiation to macrophages in vitro. IFN-λ stimulates macrophage cytotoxicity and phagocytosis as well as the secretion of pro-inflammatory cytokines and interferon stimulated genes that mediate immune cell chemotaxis and effector functions. In particular, IFN-λ induced CCR5 and CXCR3 chemokines, stimulating T and NK cell migration, as well as subsequent NK cell cytotoxicity. Using immunofluorescence and cell sorting techniques, we confirmed that human liver macrophages expressing CD14 and CD68 are highly responsive to IFN-λ ex vivo. Together, these data highlight a novel role for macrophages in shaping IFN-λ dependent immune responses both directly through pro-inflammatory activity and indirectly by recruiting and activating IFN-λ unresponsive lymphocytes.

Spatio-temporal dynamics of Host-Virus competition: A model study of influenza A

We present results of a study of the early-time response of the innate immune system to influenza virus infection in an agent-based model (ABM) of epithelial cell layers. We find that the competition between the anti-viral immune response and viral antagonism can lead to viral titers non-monotonic in the initial infection fraction as found in experiments. Our model includes a coarse-grained version of intra-cellular processes and inter-cellular communication via cytokine and virion diffusion. We use ABM to follow the propagation of viral infection in the layer and the increase of the viral load as a function of time for different values of the multiplicity of infection (MOI), the initial number of viruses added per cell. We find that for moderately strong host immune response, the number of infected cells and viral load for a smaller MOI exceeds that for larger MOI, as seen in experiments. We elucidate the mechanism underlying this result as the synergistic action of cytokines secreted by infected cells in controlling viral amplification for larger MOI. We investigate the length and time scales that determine this non-monotonic behavior within the ABM. We study the diffusive spread of virions and cytokines from a single infected cell in an absorbing medium analytically and numerically and deduce the length scale that yields a good estimate of the MOI at which we find non-monotonicity. Detailed computations of the temporal behavior of averaged quantities and spatial measures provide further insights into host-viral interactions and connections to experimental observations.

Identification and quantification of defective virus genomes in high throughput sequencing data using DVG-profiler, a novel post-sequence alignment processing algorithm

Most viruses are known to spontaneously generate defective viral genomes (DVG) due to errors during replication. These DVGs are subgenomic and contain deletions that render them unable to complete a full replication cycle in the absence of a co-infecting, non-defective helper virus. DVGs, especially of the copyback type, frequently observed with paramyxoviruses, have been recognized to be important triggers of the antiviral innate immune response. DVGs have therefore gained interest for their potential to alter the attenuation and immunogenicity of vaccines. To investigate this potential, accurate identification and quantification of DVGs is essential. Conventional methods, such as RT-PCR, are labor intensive and will only detect primer sequence-specific species. High throughput sequencing (HTS) is much better suited for this undertaking. Here, we present an HTS-based algorithm called DVG-profiler to identify and quantify all DVG sequences in an HTS data set generated from a virus preparation. DVG-profiler identifies DVG breakpoints relative to a reference genome and reports the directionality of each segment from within the same read. The specificity and sensitivity of the algorithm was assessed using both in silico data sets as well as HTS data obtained from parainfluenza virus 5, Sendai virus and mumps virus preparations. HTS data from the latter were also compared with conventional RT-PCR data and with data obtained using an alternative algorithm. The data presented here demonstrate the high specificity, sensitivity, and robustness of DVG-profiler. This algorithm was implemented within an open source cloud-based computing environment for analyzing HTS data. DVG-profiler might prove valuable not only in basic virus research but also in monitoring live attenuated vaccines for DVG content and to assure vaccine lot to lot consistency.

Molecular epidemiology survey and characterization of human influenza A viruses circulating among Palestinians in East Jerusalem and the West Bank in 2015

Frequent typing and molecular characterization of influenza A (IAV) strains are crucial for the identification of circulating subtypes and for the selection of the subtypes’ lineages to be included in the annually prepared vaccine cocktail. We investigated IAV sampled from an underrepresented population from Palestine. 200 nasopharyngeal aspirates (NPA) were collected between February and May of 2015 from Palestinians in East Jerusalem and the West Bank suffering from mild to severe symptoms of upper respiratory infections. NPA were screened for the presence of IAV using RT-PCR. Epidemiological data, hemagglutinin (HA) and neuraminidase (NA) gene sequences were analyzed in IAV positive samples. 50 samples tested positive for IAV; 48% of which were identified as A(H1N1)pdm09 and 52% as A(H3N2), respectively. Infection with A(H1N1)pdm09 occurred mainly in April, while A(H3N2) infections were mainly detected in March. Most IAV infections in 6-year-olds and below were attributed to subtype A(H3N2), while A(H1N1)pdm09 was responsible for most infections in adults above 18-year-olds. Analyses of HA and NA amino acid sequences revealed numerous substitutions. Thereafter, and based on the HA analysis, the Palestinian A(H1N1)pdm09 isolates fell into clade 6B, while the A(H3N2) isolates fell into clades 3C.2 and 3C.3, respectively. This study is significant in providing the first insight into the epidemiology and genetic properties of IAV circulating in Palestine. In contrast to international reports for the same season, A(H3N2) was not the dominant subtype as in northern hemisphere, nor was A(H1N1)pdm09 as in WHO reports for the Middle East, however genetic properties of Palestinian A(H3N2) and A(H1N1)pdm09 were in line with global isolates.

The use of plant lectins to regulate H1N1 influenza A virus receptor binding activity

We applied an in vitro selection approach using two different plant lectins that bind to α2,3- or α2,6-linked sialic acids to determine which genetic changes of the A/California/04/09 (H1N1) virus alter hemagglutinin (HA) receptor binding toward α2,3- or α2,6-linked glycans. Consecutive passages of the A/California/04/09 virus with or without lectins in human lung epithelial Calu-3 cells led to development of three HA1 amino acid substitutions, N129D, G155E, and S183P, and one mutation in the neuraminidase (NA), G201E. The S183P mutation significantly increased binding to several α2,6 SA-linked glycans, including YDS, 6'SL(N), and 6-Su-6'SLN, compared to the wild-type virus (↑3.6-fold, P < 0.05). Two other HA1 mutations, N129D and G155E, were sufficient to significantly increase binding to α2,6-linked glycans, 6'SLN and 6-Su-6'SLN, compared to S183P (↑4.1-fold, P < 0.05). These HA1 mutations also increased binding affinity for 3'SLN glycan compared to the wild-type virus as measured by Biacore surface plasmon resonance method. In addition, the HA1 N129D and HA1 G155E substitutions were identified as antigenic mutations. Furthermore, the G201E mutation in NA reduced the NA enzyme activity (↓2.3-fold). These findings demonstrate that the A/California/04/09 (H1N1) virus can acquire enhanced receptor affinity for both α2,3- and α2,6-linked sialic receptors under lectin-induced selective pressure. Such changes in binding affinity are conferred by selection of beneficial HA1 mutations that affect receptor specificity, antigenicity, and/or functional compatibility with the NA protein.