A natural variant in ANP32B impairs influenza virus replication in human cells

Viruses require host factors to support their replication, and genetic variation in such factors can affect susceptibility to infectious disease. Influenza virus replication in human cells relies on ANP32 proteins, which are involved in assembly of replication-competent dimeric influenza virus polymerase (FluPol) complexes. Here, we investigate naturally occurring single nucleotide variants (SNV) in the human Anp32A and Anp32B genes. We note that variant rs182096718 in Anp32B is found at a higher frequency than other variants in either gene. This SNV results in a D130A substitution in ANP32B, which is less able to support FluPol activity than wild-type ANP32B and binds FluPol with lower affinity. Interestingly, ANP32B-D130A exerts a dominant negative effect over wild-type ANP32B and interferes with the functionally redundant paralogue ANP32A. FluPol activity and virus replication are attenuated in CRISPR-edited cells expressing wild-type ANP32A and mutant ANP32B-D130A. We propose a model in which the D130A mutation impairs FluPol dimer formation, thus resulting in compromised replication. We suggest that both homozygous and heterozygous carriers of rs182096718 may have some genetic protection against influenza viruses.


INTRODUCTION
All viruses rely on host factors to support their replication. Genetic variation in such proteins can affect susceptibility to infectious disease [1][2][3][4]. Perhaps the best-known mutation conferring genetic resistance to viral infection is the Δ32 variant of the HIV-1 co-receptor C-C chemokine receptor type 5 (CCR5) [5]. Single nucleotide variants (SNV) that redutce susceptibility to virus infection include a P424A substitution in the filovirus endosomal fusion receptor Niemann-PickC1 (NPC1) [6,7], and a G428A mutation in fucosyltransferase 2 (FUT2) that renders homozygous carriers resistant to norovirus [8,9]. These variants affect viral entry -no SNVs resulting in defects in replication have thus far been described.
Genetic variation in the endosomal transmembrane protein IFITM3 has been associated with increased severity of infection with influenza virus. The single nucleotide polymorphism (SNP) rs12252 leads to a 21-amino acid N-terminal truncation through an alteration of the first splice acceptor site, and is perhaps the best known genetic variant affecting influenza disease severity . Although the association has not always been clear [13,14], a recent meta-analysis including 12 studies and >16 000 subjects confirmed the link between the IFITM3 minority allele and influenza severity [15]. Another important polymorphism in the IFITM3 gene associated with severe influenza is located in the 5′ UTR and leads to transcriptional repression through increased CTCF binding to the promoter [16].

ACCESS
Almost all the host genetic variation relevant to influenza that has been described is associated with exacerbating rather than alleviating influenza virus infection, a bias explained by the higher visibility of severe cases, and the opportunity to perform whole genome sequencing on patients with severe symptoms. Hence, besides a pair of polymorphisms in the non-coding region (NCR) of the Galectin-1 encoding gene LGALS1, which may confer some protection against H7N9 influenza A virus (IAV) in Chinese poultry workers [28,29], no single nucleotide variants protective against influenza virus have thus far been described.
Influenza replication and transcription are carried out in the host cell nucleus by the virus-encoded heterotrimeric RNAdependent RNA polymerase (FluPol) (reviewed in [30,31]), which relies on host factors to support RNA synthesis [32]. Oligomerisation of trimeric FluPol molecules into dimers, trimers or even higher order complexes is now believed to be required for efficient replication [33][34][35][36][37].
Influenza viruses co-opt the host acidic nuclear phosphoproteins of 32 kilodaltons A (ANP32A) and ANP32B to support their replication [38][39][40]. ANP32A and ANP32B are functionally redundant in their support for influenza virus: in human cells lacking either but not both ANP32A (AKO) or ANP32B (BKO), virus proliferation is not impaired. In cells lacking both proteins (dKO), however, FluPol activity and virus replication are completely abrogated [39,40]. ANP32 proteins are expressed in many different human tissues, including fibroblasts and lung [41], and involved in a wide range of cellular processes, including apoptosis and transcriptional regulation (reviewed in [42]). All of the human paralogues -ANP32A, ANP32B and ANP32E -associate with histones [43][44][45][46], placing them at the chromatin within the nucleus.
ANP32A and ANP32B are 249 and 251 amino acids long, respectively, with two major domains: a structured N-terminal leucine-rich repeat region (LRR) and an intrinsically disordered low-complexity acidic region (LCAR) [46][47][48]. Mutational analysis of ANP32 proteins of several species has shown that the identity of amino acid residues 129 and 130 dictates pro-viral function [39,40,49,50]. ANP32 proteins that cannot be co-opted by influenza virus include human ANP32E, chicken ANP32B, and mouse ANP32A. What these orthologues have in common is divergence from the pro-viral dyad 129N-130D: human ANP32E has glutamate at position 129 (129E) rather than asparagine, chicken ANP32B has isoleucine at position 129 (129I) and asparagine at position 130 (130N), and mouse ANP32A has alanine at position 130 (130A), suggesting that these amino acids are important for pro-viral activity of ANP32 proteins.
Recent structural work depicts the ANP32A LRR domain as stabilising a replication-competent dimer of a viral promoter-bound 'replicating' FluPol heterotrimer (FluPol R ) and an 'encapsidating' FluPol E , which is proposed to receive the nascent viral RNA product synthesised by FluPol R [33]. More specifically, the N-terminal half of the LRR interacts with FluPol R , while the C-terminal portion interacts with FluPol E , in effect bridging the dimer. The importance of residues 129 and 130 is explained by direct interaction with basic residues on the PB1 subunit of FluPol E . Interaction of the LCAR domain of ANP32A with FluPol has also been postulated but remains unresolved [33,47,51].
Here we describe natural variation in the human Anp32A and Anp32B genes, demonstrating that the naturally occurring SNV rs182096718 in Anp32B can confer protection from influenza virus replication by compromising the interaction between FluPol and the host factor. We demonstrate that rs182096718 exerts a dominant-negative effect over wildtype ANP32B and moreover interferes with the functionally redundant paralogue ANP32A through paralogue interference [52]. We propose a model explaining why influenza replication is impaired in cells expressing the mutant ANP32B protein.

CRISPR/Cas9 genome editing
Guide RNA GCAC TCTC TCGG TAGT CATTC was designed manually against the protospacer sequence in exon 4 of Anp32B to target DNA endonuclease SpCas9, expressed from Addgene plasmid # 62 988 (PX459), to the target nucleotide. The guide RNA itself was cloned into Addgene plasmid # 80 457 (pmCherry_gRNA). A custom-designed 88 base ssODN (single strand DNA) homology-directed repair template, harbouring the point mutation (cytosine in bold) and a silent PAM mutation (thymidine in italic script) -GAAA AGCC TGGA CCTC TTTA ACTG TGAG GTTACCAA T CTGAATG C CTAC CGAG AGAG TGTC TTCA AGCT CCTG CCCC AGCT TACCTAC -was obtained from Integrated DNA Technologies (IDT). Equal amounts of PX459 and pmCherry_gRNA (total 1 µg), and 1.0 µl 10 µM ssODN template, were transfected by electroporation into approximately 400 000 eHAP cells using the Neon transfection system (Invitrogen). Cells and DNA were mixed into a 10 µl volume of suspension buffer, which was subjected to a single 1200 Volt pulse for a duration of 40 ms. Cells were incubated at 37 °C for 24 h in IMDM growth medium without antibiotics. A fluorescenceactivated cell sorter (FACS) Aria IIIU (BD Biosciences) with an 85 µm nozzle was used to sort cells expressing mCherry (550-650 nm emission) into 96-well plates containing growth medium. Single cells were grown out into monoclonal populations over a period of 10 to 14 days. Total genomic DNA was extracted using the Purelink Genomic DNA Mini Kit (Invitrogen) and amplified by touchdown PCR to generate a 1561-base pair fragment of the edited locus (primers TACC TCTG CCCT CTCA ATCTCT and ACGC ACAC AAAC ACAC ACTATT). PCR products were then incubated at 65 °C for 30 min in the presence of BsmI restriction enzyme (NEB). The resulting DNA fragments were separated by 1.5 % agarose gel electrophoresis. Potentially successfully edited clones were verified by Sanger sequencing (primers TAAA GACC GCTT GATA CCCAGG and TGAGGCTGAGTGGGTAGTGG) and analysed in Geneious prime 2019.

Split luciferase complementation assay
pCAGGS expression plasmids encoding H3N2 Vic/75 PB1-luc1, PB2, PA, and the indicated ANP32-luc2 construct were transfected into ~100 000 293 T cells at a ratio of 1 : 1 : 1 : 1 (15 ng per well). Control conditions contained pCAGGS-luc1 and untagged PB1, or pCAGGS-luc2 and untagged ANP32A, respectively, with all other components remaining constant. Empty pCAGGS plasmid was used to ensure total transfected DNA was equal across conditions. Twenty-four hours after transfection, cells were lysed in 50 µl Renilla lysis buffer (Promega) for 1 h at room temperature with gentle shaking (Gaussia and Renilla luciferase share the same substrate). Bioluminescence generated by Gaussia luciferase was measured using the Renilla luciferase kit (Promega) on a FLUOstar Omega plate reader (BMG Labtech). Normalized luminescence ratios (NLR) were calculated by dividing the signal from the potential interacting partners by the sum of the two controls, as described [51].

Influenza virus infection
A total of 500 000 CRISPR/Cas9-modified monoclonal eHAP cells were infected with ~2500 plaque-forming units (PFU) H3N2 Vic/75 6 : 2 or pH1N1 Eng/195 6 : 2 virus diluted in 200 µl serum-free IMDM for 1 h at 37 °C (MOI 0.005) to allow virus to adsorb and enter the cells. The inoculum was removed and cells were incubated in room-temperature phosphate-buffered saline / HCl at pH 3.0 for 3 min to inactivate residual virus. Cells were incubated at 37 °C in serum-free cell culture medium (IMDM) supplemented with 1 µg ml −1 L-1-tosylamide-2-phenylethyl chloromethyl ketone (TPCK) trypsin (Worthington-Biochemical). Cell supernatants were harvested at indicated time points post-infection. Infectious titres were determined by plaque assay on MDCK cells. Virus infection assays were performed in triplicate on two separate occasions.

Safety/biosecurity
All work with infectious agents was conducted in biosafety level two facilities, approved by the Health and Safety Executive of the United Kingdom and in accordance with local rules, at Imperial College London, United Kingdom.

A relatively common single nucleotide variant encodes ANP32B-D130A
We searched public databases for missense SNVs in the human Anp32A and Anp32B genes that were predicted, on the basis of our previous work, to have functional impact on the pro-influenza virus activity of the proteins. We predominantly used the genome aggregation database gnomAD v2, which holds 125 748 whole exomes and 15 708 whole genomes from unrelated individuals sequenced in disease-specific or population genetic studies [53]. The gnomAD database identifies 54 missense SNVs in Anp32A and 82 in Anp32B (observed/expected score of 0.39 vs 0.68), suggesting that Anp32A may be slightly more conserved. It is not currently possible to obtain specific phenotypic information about individual carriers of a particular SNV from gnomAD.
Four SNVs in Anp32A and five in Anp32B were selected for investigation, taking into account changes in charge, bulk, polarity and putative solvent exposure of the encoded amino acid substitutions [46,48], as well as the previously established importance of positions 129 and 130 for proinfluenza viral activity (Fig. 1a). All the variants in Anp32A are exceedingly rare (MAF <0.001), and no homozygous carriers have been identified. The serine mutants at position 158, S158A and S158T were of interest as this is a putative phosphorylation site [54]. The SNV responsible for a mutant ANP32A-D130N protein (rs771143708) is found in a single heterozygous individual out of 3854 in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. This SNV is not described in the much larger gnomAD database, although it is cited in the NCBI database dbSNP. rs751239680 (ANP32A-R132Q) is found in a single heterozygous individual out of 10 052 in the African cohort of the gnomAD database, rs772530468 (ANP32A-S158T) is present in a single South Asian male in 23 070 in the gnomAD database, as well as two separate homozygous individuals in the Trans-omics for precision medicine (TOPMed) database [55]. Finally, a single African male in 5032 in the gnomAD database is a heterozygous carrier of rs772530468 (ANP32A-S158A) (Fig. 1b).
Most variants in Anp32B were also low-frequency (MAF <0.0001): rs377406514 (ANP32B-L128V) is found in a single African female out of 10 066 in the gnomAD database; rs771977254 (ANP32B-E133Q) occurs in two out of 4900 Ashkenazi Jewish females in the gnomAD database; rs761932651 (ANP32B-L138H) is found in one South Asian male in 23 070 in the gnomAD database; rs770020996 (ANP32B-L142F) occurs in one female in 7544 and one male in 23 072 from the South Asian cohort of the gnomAD database (Fig. 1c).
There was, however, an interesting exception. SNV rs182096718, encoding ANP32B-D130A, was relatively common in the Hispanic / Latino cohort of the gnomAD database, where in a total pool of 35 420 individuals, 1209 heterozygous carriers were identified, as well as 25 homozygous carriers. The D130A substitution in ANP32B is thus present in 3.41 % of the Latino cohort. The global MAF of this variant is 0.0044, much higher than any other SNV in either Anp32A or Anp32B.

Natural variants of ANP32 proteins affect support of FluPol activity
We have previously described a method to measure the capacity of mutated ANP32 proteins to act as pro-viral factors for influenza polymerase activity. This is achieved by exogenous expression of the cloned mutants in human eHAP cells lacking ANP32A and ANP32B (dKO), in which influenza polymerase is unable to function in absence of complementing ANP32 [39,56]. Here we expressed each Compared with wild-type ANP32A (WTA), ANP32A-D130N did not support Eng/195 polymerase activity at all, while R132Q and S158A substitutions had significantly reduced capacity to support FluPol activity (Fig. 2a). An artificial mutant with the phosphomimic S158E also had reduced capacity to support Eng/195 polymerase. In contrast, ANP32A-S158T supported FluPol to an extent similar to wild-type ANP32A. Compared with wild-type ANP32B (WTB), ANP32B-D130A had a deleterious effect on the support for Eng/195 polymerase activity, as did the leucine to valine substitution at position 128 (ANP32B-L128V) (Fig. 2b). ANP32B-E133Q was also significantly less able to support Eng/195 FluPol activity. Substitutions of the leucines at positions 138 and 142 to histidine (L138H) and phenylalanine (L142F), respectively, did not compromise the ability of the mutant ANP32B proteins to support Eng/195 polymerase activity.
In general, the effects of natural variation in ANP32 proteins was similar for Vic/75 polymerase. ANP32A-D130N did not support Vic/75 polymerase activity, compared with WTA (Fig. 2c), but R132Q and S158A substitutions had a smaller effect on Vic/75 FluPol activity than on Eng/195 FluPol activity. The S158T substitution and the phosphomimic S158E were as capable of supporting Vic/75 FluPol activity as WTA. As seen with Eng/195 polymerase, the ANP32B-L128V substitution was unable to rescue Vic/75 polymerase activity, but the D130A mutation had a detrimental effect and resulted in a reduction in Vic/75 polymerase activity (Fig. 2d). The observed differences in FluPol activity were not explained by differences in expression of the FLAG-tagged ANP32 constructs (Western blots accompanying Fig. 2a-d), nor by impaired nuclear localisation of the mutant proteins (Fig. 2e, f).

ANP32 position 130 mutants show impaired binding to FluPol
Using a split luciferase complementation assay previously developed in our laboratory [51] (Fig. 3a) we next assessed whether mutations at amino acid 130 affected interaction between ANP32 and FluPol. We and others have previously shown that murine ANP32A, which naturally harbours 130A, does not support influenza virus polymerase in minireplicon assays. Substituting 130A in human ANP32A greatly reduced FluPol activity, while introducing 130D in mouse ANP32A rescued its capacity to support FluPol activity [39,40]. Indeed, ANP32B KO mice, but not ANP32A KO mice, showed reduced viral loads and mortality when infected with H3N2 or H5N1 influenza A virus [57]. Poor replication of H7N9 viruses in ANP32A knockout mice has also been described however [58]. Interestingly, interaction of murine ANP32A with FluPol was reduced >2-fold compared with mouse ANP32B, which does have pro-influenza virus activity, and with wild-type human ANP32A (Fig. 3b).
The naturally occurring human ANP32 variants ANP32A-D130N and ANP32B-D130A also showed significantly reduced binding to FluPol (Fig. 3C). The luciferase signal indicating binding of ANP32A-D130N to FluPol was reduced >2-fold, relative to wild-type ANP32A, while the signal indicating FluPol interaction of ANP32B-D130A was greatly reduced compared with wild-type ANP32B. These differences in binding affinity were not due to differential expression of the ANP32-Gluc2 constructs (Fig. 3d).

ANP32 position 130 mutants exert dominantnegative effects
To investigate the functional consequence of the ANP32B-D130A variant, we artificially recapitulated heterozygous and homozygous ANP32B-D130A variant genotypes by exogenous expression of the mutated proteins in CRISPRedited eHAP cells that lack ANP32B expression (BKO) [39] (Fig. 4a). Minigenome assays were performed with reconstituted pH1N1 Eng/195 polymerase complemented by transient transfection with increasing amounts of either wild-type ANP32B (homozygous wild-type; blue bars), ANP32B-D130A (homozygous mutant; red bars) or a 1 : 1 ratio of both (heterozygous mutant; purple bars). We found that adding increasing amounts of wild-type ANP32B, essentially mimicking the homozygous wildtype genotype (blue bars), had no effect on FluPol activity compared with FluPol supported by ANP32A alone (grey bar). This is presumably due to the presence of wild-type ANP32A in the BKO cells, since ANP32A and ANP32B serve redundant roles in supporting FluPol [39]. In contrast, polymerase activity decreased significantly when ANP32B-D130A rather than wild-type ANP32B was co-transfected, and this effect was more dramatic the more variant was expressed (red bars). This suggests paralogue interference of ANP32B-D130A over wild-type ANP32A, which would otherwise support polymerase activity effectively on its own (grey bar). Co-transfecting increasing amounts of wildtype and mutant ANP32B in a 1 : 1 ratio, recapitulating the heterozygous mutant genotype (purple bars), also resulted in a significant drop in polymerase activity. This suggests that ANP32B-D130A also exerts a dominant-negative effect over wild-type ANP32B. These observations were unrelated to expression levels of the FLAG-tagged ANP32 constructs (accompanying Western blots). Statistical significance of the differences in Eng/195 polymerase activity between conditions is shown in the accompanying Table. We next investigated whether similar dominant-negative effects were exerted by the rare ANP32A-D130N variant (Fig. 4b). In eHAP cells lacking ANP32A expression (AKO), pH1N1 Eng/195 polymerase activity increased significantly when wild-type ANP32A was provided by transient co-transfection (blue bars), compared with FluPol activity supported by ANP32B alone (grey bar). In contrast, decreased polymerase activity was evident when either ANP32A-D130N alone (red bars) or wild-type and mutant ANP32A in a 1 : 1 ratio (purple bars) were co-transfected, albeit not in comparison with the ANP32B-only condition (grey bar). This suggests that, like ANP32B-D130A, ANP32A-D130N is dominant-negative over wild-type ANP32A, although there is no evidence that it also exerts paralogue interference over wild-type ANP32B. Together these data show that different substitutions of the canonical aspartate at position 130 of both pro-influenza viral human ANP32 proteins can have comparable phenotypic effects and may impact influenza proviral activity even in the heterozygous condition.

IAV polymerase activity and replication are attenuated in ANP32B-D130A mutant cells
To further probe the potential significance of the relatively common ANP32B-D130A variant, we generated a cell line with the wild-type Anp32A / mutant Anp32B (D130A) genotype by CRISPR/Cas9 genome editing, using human codon-optimised SpCas9 and a single-stranded DNA (ssODN) homology-directed repair template. We also generated a cell line lacking the entire pro-viral 128-130 loop (Δ128-130) (Fig. 5a). We selected an unsuccessfully edited clone to serve as a negative control, and carried out Western blotting analysis to ensure endogenous wild-type ANP32A and wild-type or mutant ANP32B were expressed in control and mutant cell lines (Fig. 5b). Using minigenome reporter assays, we found that reconstituted IAV polymerases from Eng/195 (Fig. 5c) or Vic/75 (Fig. 5d) were significantly less active in edited ANP32B-D130A and Δ128-130 cells, compared with control cells. As previously described [39], polymerase activity was not attenuated in BKO cells, again suggesting paralogue interference of mutant ANP32B over wild-type ANP32A. In other words, IAV FluPol activity was higher in cells completely lacking ANP32B expression than in cells expressing the D130A mutant. As described previously [39], no polymerase activity was seen in cells lacking both ANP32A and ANP32B (dKO), nor in the negative control conditions where the polymerase subunit PB2 was lacking (2P).
Next we infected mutant (ANP32B-D130A), ANP32B knockout, or control cells at low MOI (0.005) with 6 : 2 recombinant influenza A viruses containing the internal genes from pH1N1 Eng/195 and H3N2 Vic/75, with the neuraminidase (NA) and haemagglutinin (HA) external genes of the laboratory-adapted H1N1 PR/8 virus. Recombinant rather than wild-type viruses were used to enable us to carry out the work at biosafety level two. Eng/195 6 : 2 virus replication was severely attenuated in mutant D130A cells, in comparison with either control or BKO cells (Fig. 5e). Plaques were below the level of detection (LOD) 16 and 24 h post-infection and still three logs lower than control cells 48 h after infection (P=0.0012). Vic/75 6 : 2 infectious virus production in mutant cells was also lower than in control cells (Fig. 5f) but less affected than Eng/195. Both viruses replicated to higher titres in cells that completely lacked ANP32B (BKO) than in cells expressing the ANP32B-D130A mutant, suggesting paralogue interference of ANP32B-D130A over wild-type ANP32A also in an infectious virus context.

DISCUSSION
Here, we used publicly available databases to perform a biassed screen for naturally occurring single nucleotide variants in human Anp32A and Anp32B, under the hypothesis that some mutations may affect susceptibility to influenza virus infection. One of the mutations, which translates to an aberrant ANP32B-D130A protein, was highly enriched in carriers of Latino descent. Although most carriers of this SNV (>1200) are heterozygous for the mutation, 25 homozygotes have been reported so far in the gnomAD database. Almost all carriers occur in the Hispanic/Latino subpopulation. We generated this naturally occurring homozygous genotype by CRISPR/Cas9 genome editing in low-ploidy human eHAP cells, and found that IAV FluPol activity and virus replication were significantly attenuated in the mutant cell line, compared to wild-type control cells or cells lacking ANP32B. This effect was most pronounced with the 2009 pandemic H1N1 strain Eng/195, which may be due to subtle differences in the interaction surfaces of the FluPol trimer. Importantly, we demonstrate that the mutant ANP32B-D130A protein exerts a dominant-negative effect over wild-type ANP32B, as well as interference over the functionally redundant paralogue ANP32A. This suggests that despite the redundancy in proinfluenza viral function of ANP32A and B, even in heterozygotes there may be some consequence of this genetic variation to outcome of virus infection. We provide the first example of a single nucleotide variant in the coding region of a human gene that might confer some protection against influenza virus infection.
Although we and others have previously reported redundancy in reconstituted systems for ANP32A and B in their support of FluPol activity [39,40], the data presented here imply that deleterious mutations in one ANP32 protein might exert dominant-negative effects on the pro-viral activity of the other, a phenomenon known as paralogue interference. Poor interaction with FluPol of one paralogue might be expected to make no difference when both paralogues are present in the cell, but clearly this does seem to be the case.
Recent structural data suggest that ANP32 proteins assist influenza polymerase by stabilizing a dimer between the replicating FluPol trimer, FluPol R , and a second trimer that is recruited as an encapsidating complex, FluPol E [33]. To explain the dominant-negative effects we observe here, as well as the compromised replicative activity, we propose a model in which the mutant ANP32B-D130A binds competitively to FluPol R , via the unaltered N-terminal portion of the LRR domain, but is unable to recruit FluPol E due to the lack of acidity at position 130 (Fig. 6). Consequently, nascent viral RNA generated by FluPol R is not transferred effectively to FluPol E , and replication is compromised. Furthermore, recruitment of ANP32B-D130A to FluPol R occupies the binding site for proviral wild-type ANP32A or ANP32B proteins, explaining the paralogue interference and dominant-negative effects we observe.
It is important to bear in mind that the outcome of genetic variation will depend on gene dosing in relevant tissues. The genotype-tissue expression (GTEx) portal uses whole genome sequencing and RNA-Seq data to estimate gene expression in 54 non-diseased human tissue types from nearly 1000 individuals [41]. ANP32A and ANP32B are both clearly detectable in healthy lung tissue, at a median transcripts per kilobase million (TPM) of 42 and 201, respectively. There are thus 4.8-fold more ANP32B transcripts than ANP32A transcripts in the lung, suggesting that mutations that impair pro-viral function of ANP32B, such as ANP32B-D130A, might have significant effects on virus replication in the target tissue. Moreover, ANP32B is more stable than ANP32A, judging from half-life measurements of both paralogues [59].
A caveat of the work presented here is the use of human eHAP cells where (primary) respiratory epithelium would have been preferable. We selected eHAP cells for their low ploidy, which renders CRISPR/Cas9 genome editing more straightforward, in the knowledge that more commonly used laboratory cell lines like A549 or HEK-293 derivatives are polyploid [60][61][62]. We and others have, however, shown previously that eHAP cells are a suitable model for influenza virus infection [33,39,[49][50][51]. A drawback of using a low-ploidy cell line is that heterozygous genotypes cannot easily be generated. Thus far we have been unable to trace carriers of the variants but future work using diploid respiratory epithelial cells that reflect the heterozygous genotype could indicate whether true heterozygotes, like the reconstituted cells we used here, also display reduced influenza susceptibility.
In conclusion, we have provided the first example of a single nucleotide variant in the coding region of a human gene that may offer carriers some protection against influenza virus. This work has the potential to inform future intervention.  replicating FluPol (FluPol R , in blue) via the N-terminal region of the LRR domain, while the C-terminal LRR, harbouring the surface-exposed acidic residue 130D, is used to recruit the encapsidating FluPol (FluPol E , in green), via direct interaction with basic residues on the PB1 subunit of FluPol E . In the context of a properly formed FluPol R -ANP32-FluPol E complex the nascent RNA product replicated by FluPol R (in pink) transfers to FluPol E and replication proceeds efficiently. (b) The mutant ANP32B-D130A can still bind efficiently -and therefore competitively -to FluPol R , as the N-terminal portion of the LRR is unaffected, but because there is now an alanine in place of the acidic aspartate at position 130, recruitment of FluPol E is impaired. Furthermore, binding of ANP32B-D130A to FluPol R prevents recruitment of wild-type ANP32A or ANP32B, explaining the paralogue interference and dominantnegative effects, respectively. Thus, in the presence of ANP32B-D130A, the replication-competent FluPol R -ANP32-FluPol E complex is not formed as efficiently, and therefore replication is compromised.