A novel E198K substitution in the PA gene of influenza A virus with reduced susceptibility to baloxavir acid

High throughput profiling identified PA-L106R amino acid substitution in A(H1N1)pdm09 influenza virus that confers reduced susceptibility to baloxavir in vitro

Baloxavir acid (BXA) is a pan-influenza antiviral that targets the cap-dependent endonuclease of the polymerase acidic (PA) protein required for viral mRNA synthesis. To gain a comprehensive understanding on the molecular changes associated with reduced susceptibility to BXA and their fitness profile, we performed a deep mutational scanning at the PA endonuclease domain of an A(H1N1)pdm09 virus. The recombinant virus libraries were serially passaged in vitro under increasing concentrations of BXA followed by next-generation sequencing to monitor PA amino acid substitutions with increased detection frequencies. Enriched PA amino acid changes were each introduced into a recombinant A(H1N1)pdm09 virus to validate their effect on BXA susceptibility and viral replication fitness in vitro. The I38T/M substitutions known to confer reduced susceptibility to BXA were invariably detected from recombinant virus libraries within 5 serial passages. In addition, we identified a novel L106R substitution that emerged in the third passage and conferred greater than 10-fold reduced susceptibility to BXA. PA-L106 is highly conserved among seasonal influenza A and B viruses. Compared to the wild-type virus, the L106R substitution resulted in reduced polymerase activity and a minor reduction of the peak viral load, suggesting the amino acid change may result in moderate fitness loss. Our results support the use of deep mutational scanning as a practical tool to elucidate genotype-phenotype relationships, including mapping amino acid substitutions with reduced susceptibility to antivirals.

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.

Development and validation of an ultra-performance liquid chromatography-tandem mass spectrometry method for the determination of baloxavir in rat plasma and its application to pharmacokinetic studies

In this study, an ultra-performance liquid chromatography-tandem mass spectrometry method was established for the development and validation of baloxavir acid (BXA) concentrations and the active ingredients of the antiviral drug baloxavir marboxil (BXM). Further, the method was applied to study the pharmacokinetics of BXA. BXA was determined by the electrospray ionization multiple reaction monitoring positive ion mode, and the mass-to-charge ratios (m/z) of BXA and internal standard baloxavir-d4 were 484.2 → 247.2 and 488.1 → 247.2. An Oasis max online column (2.1 × 20 mm, 30 μm) was used with 1% formic acid in water (A) and 2% formic acid in acetonitrile (B) as mobile phases at a flow rate of 0.5 mL·min-1 for chromatographic separation. The linearity was good in the range of 3-200 ng·mL-1 (r = 0.9994), with 3.00 ng·mL-1 lower limit of quantification. The relative standard deviation of the inter-assay precision was less than or equal to 6.51%, and the accuracy was in the range of 91.28%-104.29%. This method is suitable for the determination of BXA and for performing pharmacokinetic studies in clinical research.

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.