# Type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNA polymerase II CTD

Tim Krischuns, Catherine Isel, Petra Drncova, Maria Lukarska, Alexander Pflug, Sylvain Paisant, Vincent Navratil, Stephen Cusack, Nadia Naffakh

► To cite this version:

Tim Krischuns, Catherine Isel, Petra Drncova, Maria Lukarska, Alexander Pflug, et al., Type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNA polymerase II CTD. PLoS Pathogens, 2022, 18 (5), pp.e1010328. 10.1371/journal.ppat.1010328. pasteur-03840864

HAL Id: pasteur-03840864

https://pasteur.hal.science/pasteur-03840864

Submitted on 6 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# Type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNA polymerase II CTD

Tim Krischuns¹, Catherine Isel¹, Petra Drncova², Maria Lukarska²*, Alexander Pflug², Sylvain Paisant¹, Vincent Navratil³, Stephen Cusack¹*, Nadia Naffakh¹*

1 Institut Pasteur, Université Paris Cité, CNRS UMR3569, Unité Biologie des ARN et Virus Influenza, Paris, France, 2 European Molecular Biology Laboratory, 71 Avenue des Martyrs, Grenoble, France, 3 PRABI, Rhône Alpes Bioinformatics Center, UCBL, Lyon1, Université de Lyon, Lyon, France

<figure><img src="image_2.png" /></figure>

## Abstract

During annual influenza epidemics, influenza B viruses (IBVs) co-circulate with influenza A viruses (IAVs), can become predominant and cause severe morbidity and mortality. Phylogenetic analyses suggest that IAVs (primarily avian viruses) and IBVs (primarily human viruses) have diverged over long time scales. Identifying their common and distinctive features is an effective approach to increase knowledge about the molecular details of influenza infection. The virus-encoded RNA-dependent RNA polymerases (FluPolB and FluPolA) are PB1-PB2-PA heterotrimers that perform transcription and replication of the viral genome in the nucleus of infected cells. Initiation of viral mRNA synthesis requires a direct association of FluPol with the host RNA polymerase II (RNAP II), in particular the repetitive C-terminal domain (CTD) of the major RNAP II subunit, to enable “cap-snatching” whereby 5’-capped oligomers derived from nascent RNAP II transcripts are pirated to prime viral transcription. Here, we present the first high-resolution co-crystal structure of FluPolB bound to a CTD mimicking peptide at a binding site crossing from PA to PB2. By performing structure-based mutagenesis of FluPolB and FluPolA followed by a systematic investigation of FluPol-CTD binding, FluPol activity and viral phenotype, we demonstrate that IBVs and IAVs have evolved distinct binding interfaces to recruit the RNAP II CTD, despite the CTD sequence being highly conserved across host species. We find that the PB2 627 subdomain, a major determinant of FluPol-host cell interactions and IAV host-range, is involved in CTD-binding for IBVs but not for IAVs, and we show that FluPolB and FluPolA bind to the host RNAP II independently of the CTD. Altogether, our results suggest that the CTD-binding modes of IAV and IBV may represent avian- and human-optimized binding modes, respectively, and that their divergent evolution was shaped by the broader interaction network between the FluPol and the host transcriptional machinery.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Krischuns T, Isel C, Drncova P, Lukarska M, Pflug A, Paisant S, et al. (2022) Type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNA polymerase II CTD. PLoS Pathog 18(5): e1010328. https://doi.org/10.1371/journal.ppat.1010328

**Editor:** Anice C. Lowen, Emory University School of Medicine, UNITED STATES

**Received:** February 3, 2022

**Accepted:** April 9, 2022

**Published:** May 23, 2022

**Copyright:** © 2022 Krischuns et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The data is available through the PDB under accessions 7Z42, 7Z43 and 7Z40. The RPB1 sequences from Gallus gallus and Anas platyrhyncos are held on the Zenodo repository: https://doi.org/10.5281/zenodo.6467097 All other relevant data are within the manuscript and its Supporting information files.

**Funding:** This work was funded by the ANR (Agence Nationale de la Recherche, https://anr.fr) grant FluTranscript ANR-18-CE18-0028, held jointly by SC and NN, and grant ANR-10-LABX-62-

<footnote>¤ Current address: Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, California, United States of America</footnote>

<footnote>* cusack@embl.fr (SC); nadia.naffakh@pasteur.fr (NN)</footnote>
IBEID (NN). The funders had no role in study
design, data collection and analysis, decision to
publish, or preparation of the manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

Author summary

During seasonal influenza epidemics, influenza B viruses (IBVs) co-circulate with influ-
enza A viruses (IAVs) and can cause severe outcomes. The influenza polymerase is a key
drug target and it is therefore important to understand the common and distinctive
molecular features of IBV and IAV polymerases. To achieve efficient transcription and
replication in the nucleus of infected cells, influenza polymerases closely cooperate with
the cellular RNA polymerase II (RNAP II) and interact with the repetitive C-terminal
domain (CTD) of its major subunit. Here we gained new insights into the way IBV and
IAV polymerases interact with the CTD of RNAP II. High-resolution structural data was
used to perform structure-based mutagenesis of IBV and IAV polymerases followed by a
systematic investigation of their interaction with RNAP II, transcription/replication activ-
ity and viral phenotype. Strikingly, we found that IBVs and IAVs have evolved distinct
interfaces to interact with the host transcriptional machinery, in particular with the CTD
of RNAP II. We provide evidence that these differences may have evolved as a conse-
quence of the differences in IBV and IAV host range. Our findings are of significant
importance with regard to the development of broad-spectrum antivirals that target the
virus-host interface.

**Introduction**

Influenza viruses are members of the *Orthomyxoviridae* family and are classified into four gen-
era: influenza A, B, C and D viruses. Influenza A viruses (IAVs) and influenza B viruses (IBVs)
are of public health importance, as they co-circulate in humans with a seasonal epidemic pat-
tern and cause a significant morbidity and mortality, especially in the aged or immunocom-
promised population [1]. IBV infections account for an estimated 23% of all influenza cases
[2], can become predominant during annual influenza epidemics, and can cause severe disease
in children [3]. IBVs have received less attention because, unlike IAVs which continuously cir-
culate in a wide range of birds and mammalian species [4], they have no known potential to
cause pandemics. Based on sequence analysis of the viral hemagglutinin, the evolutionary
divergence between IBVs and IAVs was estimated to have occurred about 4000 years ago [5].
The recent identification of IBV-like viruses in non-mammalian vertebrate species suggest
that IBVs and IAVs have actually diverged over much longer time scales [6].

IBVs and IAVs share the same genome organization of eight single-stranded negative RNA segments, and major features of the viral replication cycle such as transcription and replication of the viral genome in the nucleus of infected cells. However their genes have undergone functional divergence, as reflected notably by the lack of intertypic genetic reassortment [7]. To identify common and distinctive features of IBVs and IAVs is an effective approach to improve our understanding of the molecular mechanisms of influenza infection and our ability to fight influenza disease.

The genomic RNA segments of IAVs and IBVs are organized into viral ribonucleoprotein
complexes (vRNPs) [8]. In the vRNP, the 5' and 3' terminal viral RNA sequences are associated
with one copy of the RNA-dependent RNA polymerase complex (FluPol) while the RNA is
covered by multiple copies of the viral nucleoprotein (NP) [9–11]. FluPol is a heterotrimer
composed of PB1 (polymerase basic protein 1), PB2 (polymerase basic protein 2), and PA
(polymerase acidic protein) [12], which replicates and transcribes the viral RNA in the nucleus
of infected host cells. Replication is a primer-independent two-step process, which relies on *de
novo* initiation by FluPol [13,14]. In contrast, viral transcription is primer-dependent and
results in the synthesis of 5' capped and 3' polyadenylated mRNAs, which are translated by the
host translation machinery [15,16]. Polyadenylation is achieved by stuttering of FluPol at a 5'
proximal oligo(U) stretch present on the genomic RNA [17,18]. In contrast to other RNA
virus polymerases, FluPol cannot synthesize 5' cap structures [19]. In a process referred to as
cap-snatching [20], FluPol binds the 5' cap of nascent host RNA polymerase II (RNAP II) tran-
scripts by the PB2 cap-binding domain. Then, the PA endonuclease domain [21] cleaves 10-
15 nts downstream of the 5' cap thereby generating primers that are used by FluPol to initiate
transcription [18,19,22].

To perform cap-snatching, FluPol needs access to nascent capped RNAP II-derived RNAs,
which represents a challenge as host cap structures are rapidly sequestered co-transcriptionally
by the cap-binding complex [23]. The cellular RNAP II consists of 12 subunits [24], and the
largest subunit (RPB1) is characterised by a unique long unstructured C-terminal domain
(CTD) which in mammals consists of 52 repeats of the consensus sequence Tyr-Ser-Pro-Thr-
Ser-Pro-Ser (Y₁S₂P₃T₄S₅P₆S₇). Post-translational modifications of the CTD during the tran-
scription process are controlling the spatiotemporal regulation of RNAP II transcription
[25,26]. FluPol binds specifically to S5 phosphorylated CTD (CTD pS5) [27,28] and it was pro-
posed that it targets RNAP II for cap-snatching in the paused elongation state, of which CTD
pS5 is the hallmark modification [29–31].

Structural studies revealed bipartite CTD binding sites on the FluPol of influenza A, B and C viruses (FluPolA, FluPolB and FluPolc) with notable differences from one type to another [32,33]. However, the original crystal structure data for FluPolB were of insufficient resolution and only one of the CTD binding sites could be modelled, therefore preventing functional studies. In this study, we report the first high-resolution co-crystal structure of FluPolB bound to a CTD pS5 mimicking peptide that allows the modelling of both CTD-binding sites, one exclusively on PA also observed on FluPolA, and another, crossing from PA to PB2, specific for FluPolB. We used these novel data to perform structure-guided mutagenesis of FluPolB and FluPolA, followed by a systematic investigation of cell-based CTD-binding, cell-based polymerase activity and plaque phenotype of recombinant viruses. Our findings demonstrate that type B and type A influenza polymerases have evolved distinct binding interfaces to recruit the RNAP II CTD, which is intriguing as the RNAPII CTD is highly conserved across influenza host species. We find that the PB2 627 subdomain, a major determinant of FluPol-host cell interactions and IAV host-range, is involved in CTD-binding for IBVs but not for IAVs. Finally, we provide evidence for additional FluPol-RNAP II interactions that do not involve the CTD.

**Materials and methods**

Purification, crystallisation, data collection and structure determination of
FluPolB with bound CTD peptide

Influenza B/Memphis/13/2003 polymerase, wild type or with the PA K135A mutation to elimi-
nate endonuclease activity, was expressed and purified as described previously [22].

For crystals enabling high resolution visualisation of CTD binding in site 2B, FluPolB PA mutant K135A at 9 mg ml⁻¹ (35 µM) was mixed with 40 µM of nucleotides 1–13 vRNA 5' end (5'-pAGUAGUAACAAGA-3') and 1.8 mM 28-mer CTD peptide (YSPTpSPS)₄ in a buffer containing 50 mM HEPES pH 7.5, 500 mM NaCl, 5% glycerol, 2 mM TCEP. Hanging drops for crystallisation were set up at 20 °C. Rod-shaped crystals growing up to 700 µm in length appeared one week after set-up in mother liquor containing 100 mM tri-sodium citrate and 13% PEG 3350 with a drop ratio of 0.5 µl + 2 µl protein to well solution. Crystals were cryo-protected with additional 20% glycerol and 1.8 mM CTD peptide in mother liquor and
<table><caption>Table 1. Crystallographic data collection and refinement statistics.</caption><thead><tr><th>Crystal</th><th>B/Memphis polymerase (WT) +3' vRNA 1-18 +5' vRNA 1-13 +13-mer capped RNA 28-mer CTD</th><th>B/Memphis polymerase (PA K135A) +5' vRNA 1-14 28-mer CTD</th><th>A/H7N9 polymerase core 201-PA_PB1_PB2-127 12-mer 5' vRNA pSer5 28-mer CTD peptide</th></tr></thead><tbody><tr><td>wwPDB ID</td><td>7Z42</td><td>7Z43</td><td>7Z4O</td></tr><tr><td colspan="4"><strong>Diffraction Data</strong></td></tr><tr><td>Space group</td><td>P3<sub>2</sub>2<sub>1</sub></td><td>P2<sub>1</sub></td><td>P2<sub>1</sub> 2<sub>1</sub> 2<sub>1</sub></td></tr><tr><td>Cell dimensions (Å)</td><td>a = b = 200.82, c = 257.02</td><td>a = 130.93, b = 202.27, c = 135.73, β = 110.54°</td><td>a = 76.47, b = 144.13, c = 336.20</td></tr><tr><td>Wavelength (Å)</td><td>0.966</td><td>1.254</td><td>0.968</td></tr><tr><td>Beamline (ESRF)</td><td>ID30A1</td><td>ID29</td><td>ID29</td></tr><tr><td>No. Crystals</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Resolution range (last shell) (Å)</td><td>Isotropic 3.56 Å Ellipsoidal 3.12 Å 76.85-3.12 (3.32-3.12)</td><td>Isotropic 2.95 Å Ellipsoidal 2.42 Å 127.1-2.42 (2.69-2.42)</td><td>Ellipsoidal 3.41 Å 49.21-3.41 (3.57-3.41)</td></tr><tr><td>Completeness (last shell) (%)</td><td>Overall 73.1 Ellipsoidal: 91.6 (67.9)</td><td>Overall 61.4 Ellipsoidal: 89.1 (64.0)</td><td>Overall 73.3 Ellipsoidal 94.8 (88.2)</td></tr><tr><td>R-sym (last shell)</td><td>0.128 (1.30)</td><td>0.112 (0.618)</td><td>0.283 (1.873)</td></tr><tr><td>R-meas (last shell)</td><td>0.149 (1.49)</td><td>0.136 (0.769)</td><td>0.306 (1.946)</td></tr><tr><td>I/σI (last shell)</td><td>8.0 (1.7)</td><td>7.7 (1.9)</td><td>7.7 (1.5)</td></tr><tr><td>CC(1/2) (last shell)</td><td>0.994 (0.536)</td><td>0.992 (0.630)</td><td>0.995 (0.696)</td></tr><tr><td>Redundancy (last shell)</td><td>3.9 (3.9)</td><td>3.0 (2.5)</td><td>14.1 (14.6)</td></tr><tr><td colspan="4"><strong>Refinement</strong></td></tr><tr><td>Reflections (total) For refinement: work (free)</td><td>77848 73897 (3951)</td><td>154519 146946 (7573)</td><td>37875 36023 (1852)</td></tr><tr><td>R-work (last shell)</td><td>0.240 (0.369)</td><td>0.223 (0.328)</td><td>0.220 (0.313)</td></tr><tr><td>R-free (last shell)</td><td>0.273 (0.397)</td><td>0.257 (0.257)</td><td>0.273 (0.360)</td></tr><tr><td>No of non-hydrogen atoms</td><td>18464</td><td>35955</td><td>19969</td></tr><tr><td colspan="4"><strong>Validation</strong></td></tr><tr><td>RMS(bonds)</td><td>0.002</td><td>0.003</td><td>0.003</td></tr><tr><td>RMS(angles)</td><td>1.170</td><td>1.198</td><td>1.207</td></tr><tr><td>Average B-factor (Å<sup>2</sup>)</td><td>109.0</td><td>57.89</td><td>86.14</td></tr><tr><td>Ramachandran favored (%)</td><td>92.59</td><td>96.18</td><td>92.65</td></tr><tr><td>Ramachandran outliers (%)</td><td>0.32</td><td>0.09</td><td>0.79</td></tr><tr><td>Molprobity score</td><td>1.50</td><td>1.76</td><td>1.96</td></tr><tr><td>Clash score</td><td>1.20</td><td>2.74</td><td>2.83</td></tr></tbody></table>

https://doi.org/10.1371/journal.ppat.1010328.t001

flash-frozen in liquid nitrogen. Data were collected on ESRF beamline ID29 and integrated with an ellipsoidal mask using AUTOPROC/STARANISO to an anisotropic resolution of 2.42-2.95 Å. The structure was solved using molecular replacement with PHASER [34] using PDB:5FMZ as model [35]. The model was iteratively corrected and refined using COOT [36] and REFMAC5 [37] and quality-controlled using MOLPROBITY [38]. See Table 1 for data collection and refinement statistics.

For crystals enabling simultaneous visualisation of CTD binding in sites 2A and 2B, wild type FluPolB at 11.7 mg ml⁻¹ (45 µM) was mixed with 52 µM of nucleotides 1–14 (5'-pAGUA GUAACAAGAG-3') and 1–18 (5'-UAUACCUCUGCUUCUGCU-3') of respectively the vRNA 5' and 3' ends, 104 µM of capped 13-mer (5'-m⁷GpppAAUCUAUAAUAGC-3'), and 380 µM CTD 28-mer Ser5P peptide, in buffer containing 50 mM HEPES pH 7.5, 500 mM NaCl, 5% glycerol, 2 mM TCEP. Hanging drops for crystallisation were set up at 4° C. Diamond-shaped crystals growing up to 200 µm in size appeared in two to three weeks in mother liquor containing 200 mM di-ammonium phosphate, and 100 mM sodium acetate at pH 4.4,
with a drop ratio of 1 µl + 2 µl protein to well solution. The drops were soaked with 840 µM CTD peptide for 17 days. Crystals were cryo-protected with an additional 30% glycerol and 885 µM peptide in mother liquor and flash-frozen in liquid nitrogen. Data were collected on ESRF beamline ID30A1 (MASSIF) and processed and refined as described above, using PDB:5MSG as model for molecular replacement. See Table 1.

## Structure determination of FluPol<sub>A</sub>(H7N9) core with bound CTD peptide

The core of influenza A/Zhejiang/DTID-ZJU01/2013(H7N9) polymerase comprising PA 201-716, PB1 full-length, PB2 1-127 was expressed and purified from insect cells as described previously [18]. A/H7N9 polymerase core at a concentration of 9 mg/ml was co-crystallised with 60 µm of a 12-mer of the vRNA 5' end (5'-pAGUAGUAACAAG) in sitting drops at 4° C in conditions of 0.1 M Tris pH 7.0, 13% PEG 8K, 0.2 M MgCl<sub>2</sub>, 0.1 M guanidine hydrochloride with drop mixing ratios of 1:2 (protein:well). Crystals grew typically within 4-5 days and diffracted to around 3.5 Å resolution. A four-repeat pS5 CTD mimicking peptide (Tyr-Ser-Pro-Thr-pSer-Pro-Ser)<sub>4</sub> was soaked into existing crystals at a concentration of ~2 mM over a period of 24 h. Data were collected on ESRF beamline ID29 and processed and refined as described above, using previously described apo-H7N9 core structure ([18, 39], PDB:6TU5) as model for molecular replacement. See Table 1 for data collection and refinement statistics. The PDB numbers for the new protein structures are provided in Table 1: 7Z42, 7Z43 and 7Z4O.

## Cells

HEK-293T cells (purchased at ATCC, CRL-3216) were grown in complete Dulbecco's modified Eagle's medium (DMEM, Gibco) supplemented with 10% fetal bovine serum (FBS) and 100 U/mL penicillin and 100 µg/mL streptomycin. MDCK cells (provided by the National Influenza Center, Paris, France) were grown in Modified Eagle's medium (MEM, Gibco) supplemented with 5% FBS, 100 U/mL penicillin and 100 µg/mL streptomycin.

## Plasmids

The reverse genetics plasmids derived from the IAV A/WSN/1933 (WSN) [40] and the IBV B/Brisbane/60/2008 [41] were kindly provided by Pr. G. Brownlee (Sir William Dunn School of Pathology, Oxford, UK) and Pr. D. Perez (College of Veterinary Medicine, University of Georgia), respectively. For polymerase activity assays, a pPolII-Firefly plasmid encoding the Firefly luciferase sequence in negative polarity flanked by the 5' and 3' non-coding regions of either the IAV or IBV NS segment was used. The pRenilla-TK plasmid (Promega) was used as an internal control. The WSN pcDNA3.1-PB2, -PB1, -PA plasmids [32] and B/Memphis/13/2003 (Memphis) pcDNA3.1-PB2, -PB1, -PA and -NP plasmids [42] were described previously. The WSN-NP Open Reading Frame (ORF) was subcloned into the pCI-plasmid. The WSN and Memphis pCI-PB2-G1 and pCI-PA-G1 plasmids used for *Gaussia* luciferase complementation assays were constructed as described previously [43]. The RPB1 ORF was obtained from the Human ORFeome resource (hORFeome v3.1), fused to G2 at the N-terminus by PCR and subcloned into pcDNA3.1 (G2-RPB1). The CTD repeats 4 to 51 were deleted from the G2-RPB1 construct by PCR (G2-RPB1ΔCTD). The pCI-RPB2-G2 construct was kindly provided by Dr. B. Delmas (INRAE, Jouy-en-Josas). The wild-type full-length CTD sequences was fused at the C-terminus to an SV40 NLS by PCR using the G2-RPB1 plasmid as a template. The resulting amplicon was cloned in frame downstream the G2 sequence into the pCI vector (G2-CTD). A sequence in which each CTD serine 5 residue was replaced by an alanine was ordered as synthetic gene (GenScript) and subcloned in place of the wild-type CTD sequence into the G2-CTD construct (G2-CTD-S5A). The pCI-G2-NUP62 plasmid was described
previously [44]. Mutations were introduced by an adapted QuickChange site-directed muta-
genesis (Agilent Technologies) protocol [45]. Primers and plasmid sequences are available
upon request.

Protein complementation and minigenome assays

HEK-293T cells were seeded in 96-well white opaque plates (Greiner Bio-One) the day before transfection. For the split-luciferase complementation assays, cells were co-transfected in technical triplicates with 25 ng plasmid encoding the polymerase subunits PB2, PB1 and PA (either PB2-G1 or PA-G1, respectively) and 100 ng of the G2-tagged targets (CTD, RPB1 or RPB2, respectively) using polyethyleneimine (PEI-max, #24765-1 Polysciences Inc). When indicated, the CDK7 inhibitor BS-181-HC (Tocris Bioscience) was added 24 hours post-transfection (hpt) at a final concentration of 20 µM for 1 h. DMSO 0.2% was used as a control. Cells were lysed 20-24 hpt in *Renilla* lysis buffer (Promega) for 45 min at room temperature under steady shaking (650 rpm) and the *Gaussia princeps* luciferase enzymatic activity was measured on a Centro XS LB960 microplate luminometer (Berthold Technologies, reading time 10 s after injection of 50 µl Renilla luciferase reagent (Promega)). The Normalized Luminescence Ratios (NLRs) were calculated as follows: the luminescence activity (Relative Light Units or RLU) measured in cells co-transfected with the plasmids encoding the PB2/PA-G1 and G2-CTD/RPB1/RPB2 fusion proteins was divided by the sum of the luminescence activities (RLU) measured in control samples co-transfected with either the G2 and PB2/PA-G1 plasmids or the G1 and G2-CTD/RPB1/RPB2 plasmids. For the minigenome assays, cells were co-transfected in technical triplicates with 25 ng of each pcDNA3.1 PB2, PB1, PA, in conjunction with 50, 10 and 5 ng of the pCI-NP, pPolI-Firefly and pTK-Renilla plasmids, respectively. Luciferase activities were measured 20-24 hpt using the Dual-Glo Luciferase Assay system (Promega) according to the manufacturer's instructions.

Antibodies and immunoblots

Total cell lysates were prepared in RIPA cell lysis buffer as described before [46]. Immunoblot membranes were incubated with primary antibodies directed against CTD-pS5 (Active Motif, 3EB), CTD-pS2 (Active Motif, 3E10), *Gaussia princeps* luciferase (New England Biolabs, #E8023) or Tubulin (Sigma-Aldrich, B-5-1-2), and subsequently with the according HRP-tagged secondary antibodies (Jackson Immunoresearch). Membranes were revealed with the ECL2 substrate according to the manufacturer's instructions (Pierce) and chemiluminescence signals were acquired using the ChemiDoc imaging system (Bio-Rad) and analysed with Ima-geLab (Bio-Rad).

Production and characterisation of recombinant viruses

The recombinant viruses were produced by transfection of a co-culture of HEK-293T and MDCK cells as described previously [40,41]. The reverse genetics supernatants were titrated on MDCK cells in a standard plaque assay as described before [47]. Plaque diameters were measured upon staining with crystal violet using Fiji [48].

*In vitro* endonuclease and transcription activity assays

RNA for the activity assays was produced *in vitro* with T7 polymerase. Recombinant polymerases used corresponding to A/little yellow-shouldered bat/Guatemala/060/2010 and B/Memphis/13/2003 were purified as previously described [42]. 23 nt RNA (5'-GAAUCUAUACAUAAAGACCAGGC-3') was capped with vaccinia capping enzyme and 2'-O-methyltransferase
(NEB) and radiolabelled with [α³²P]-GTP. For the endonuclease assay 25 nM FluPolA and 50 nM FluPolB were incubated with 1.2-fold molar excess of v5' RNA (FluPolA: 5'-pAGUAGAAA CAAGGC-3', FluPolB: 5'-pAGUAGUAACAAGAG-3') and the capped RNA in reaction buffer containing 50 mM HEPES pH 7.5, 150 mM NaCl, 5 mM MgCl2, 2 mM tris(2-carboxyethyl) phosphine (TCEP). Transcription reactions were performed with 50 nM FluPolA or FluPolB in the reaction buffer, supplemented with 150 nM v3' template RNA (FluPolA: 5'-AGUUUGC-CUGCUUCUGCU-3', FluPolB: 5'-UAUACCUGCUUCUGCU-3') and 250 µM NTP mix (ThermoFisher). 50 µM CTD peptides were added at concentrations corresponding to at least a 10-fold excess over the K_D of the lowest measured affinity for a two-repeat peptide. Two- and four-repeat phosphoserine 5 (pS5) CTD peptides were purchased from Covalab and six-repeat pS5 CTD peptide was synthesised at the Chemical Biology Core Facility at EMBL Heidelberg. Reactions were incubated at 30 °C for 30 min and quenched with RNA loading dye (formamide, 8 M urea, 0.1% SDS, 0.01% bromophenol blue, 0.01% xylene cyanol), supplemented with 50 mM EDTA and boiled at 95 °C. The reaction products were separated on 20% denaturing acrylamide gel (containing 8 M urea) in Tris-Borate-EDTA (TBE) buffer, exposed on a Storage Phosphor screen and recorded with a Typhoon reader. DECADE marker was used as ladder.

## CTD sequences and alignment

The RPB1 CTD domain amino acid sequences of *Homo sapiens* (NP_000928.1), *Sus scrofa* (XP_020923484.1), *Equus caballus* (XP_014584045.2), *Canis lupus* (XP_038521325.1), and *Mus musculus* (NP_001277997.1) host species were obtained from the RefSeq database [49]. As the predicted RefSeq sequences available for the *Gallus gallus* (XP_040551262) and *Anas platyrhyncos* (XM_038172734) RPB1 subunits were only partial, we designed a targeted protein sequence assembly strategy data based on RNA-seq and/or WGS SRA public data available for these two species. To obtain the *Gallus gallus* RPB1 complete sequence (1969 aa), we first aligned Illumina RNA-seq short reads (ERR2664216) on the human RefSeq curated protein sequence (NP_000928) using DIAMOND algorithm [50], and then used the aligned reads for subsequent Trinity transcript assembly (“-longreads XP_040551262” option to use the partial sequence as a guide) followed by Transdecoder for the ORF prediction [51]. The *Anas platyrhyncos* RPB1 complete sequence (1970 aa) was obtained by aligning Illumina RNA-seq short reads (SRR10176883) and PACBIO long reads (SRR8718129, SRR8718130) on the JACEUL010000271.1 genomic scaffold by using respectively HISAT2 [52] and minimap2 [53] followed by Stringtie2 [54] with the “-mix” option to allow hybrid de novo gene assembly. The RPB1 coding sequences from *Gallus gallus* and *Anas platyrhyncos* are held on the Zenodo repository: https://doi.org/10.5281/zenodo.6467097. The CTD sequences were aligned with SnapGene 6.0 and visualised by Esprit 3.0 [55].

## Results

### Cocrystal structures reveal distinct CTD binding sites in FluPolB and FluPolA

Previous structural studies using a four repeat CTD pS5 peptide mimic (YSPTpSPS)₄ [32] revealed two distinct CTD binding sites on FluPolB, denoted site 1B and site 2B. Site 1B, exclusively on the PA subunit and in which the pS5 phosphate is bound by PA basic residues K631 and R634, is essentially the same as site 1A for bat influenza A polymerase [32] and is thereafter named site 1AB. Site 2B, which extends across the PA-PB2 interface, is unique to FluPolB and distinct from site 2A for FluPolA, which is again exclusive to the PA subunit [32]. However, the original crystal structure data for FluPolB were of insufficient resolution to be able to
construct a model for the CTD peptide in site 2B, nor even to define its directionality. To overcome this limitation, we co-crystallised the four repeat pS5 peptide with influenza B/Memphis/13/2003 polymerase in a different P2₁ crystal form, previously used to obtain a structure with the 5' end of the vRNA [35], and measured anisotropic diffraction data to a resolution of 2.42–2.95 Å (Table 1). The resultant map, which contains two heterotrimers in the asymmetric unit, showed clear electron density in site 2B for both trimers (S1A Fig), into which an unambiguous model for the CTD peptide could be built (Figs 1A and S2A). No significant differences were observed between the two heterotrimers except for small displacements in the poorly defined endonuclease and cap-binding domains, nor in the mode of CTD binding to site 2B, although the density for one of the binding-sites was slightly better due to B-factor differences. Only very weak density for the CTD peptide is observed in site 1AB, perhaps because of competition with a phosphate bound at the position of the phosphoserine. To reconfirm that sites 1B and 2B could be occupied simultaneously, we re-crystallised full promoter-bound FluPolB with the CTD peptide in the original P321 crystal form, but this time with a capped primer and at lower pH. Under these conditions, the extremity of the vRNA 3' end is in the RNA synthesis active site [56]. Anisotropic diffraction data to a resolution of 3.12–3.56 Å was measured and the resultant map showed clear electron density for the CTD peptide bound in both sites 1AB and 2B (Fig 1B and Table 1), as reported previously for this crystal form [32] but with slightly improved resolution. Unexpectedly, the CTD peptides bound in site 1AB and site 2B are orientated such that they cannot be linked by the shortest path, as this would be between both N-termini, which are ~ 17 Å apart, whereas the straight-line distance between the C-ter of site 1AB and N-ter of site 2B is ~ 36 (44) Å. These distances suggest that a minimum of 6, probably 7, heptad repeats would be required to occupy both sites contiguously (Fig 1B, dotted red line). This contrasts with the situation in FluPolA, where the peptide directionality in sites 1AB and 2A allow them to be linked by the shortest path, implying that four heptad repeats is sufficient to occupy both sites (Fig 1C) [32].

Three repeats (designated repeats a, b and c) of the CTD peptide (i.e. Y1aS2aP3aT4apS5a-P6aS7a-Y1bS2bP3bT4bpS5bP6bS7b-Y1cS2cP3cT4cpS5cP6c) are visible in site 2B in both structures, including two well-defined phosphoserines (in bold). The N-terminal part of the CTD peptide (Y1a-S2b) forms a compact structure comprising two successive proline turns stabilised by four intra-peptide hydrogen bonds, with P3a stacked on Y1b and P6a stacked on PA/Y597 (Fig 1D). PB2 R134 partially stacks against the other side of the Y1b sidechain, whose hydroxyl group hydrogen bonds to the main-chain of PB2 I135. The phosphate of pS5a forms a strong salt-bridge with PA R608 as well as hydrogen bonding with S7a. FluPolB-specific PA R608 is in a four-residue insertion (606-GDRV-609) compared to FluPolA, with hydrogen bond interactions from PA D607 and N611 positioning the side-chain of PA Y597 under the CTD peptide. This configuration of residues seems specifically designed to accommodate the compactly folded CTD peptide. Interestingly, recently identified FluPolB-like polymerases from fish and amphibians [57] also possess the four-residue insertion in PA. However, only in the Wuhan spiny-eel influenza virus polymerase, which is remarkably similar to human FluPolB, are all the functional residues Y597, D607, R608, N611 conserved [6,57]. The rest of the CTD peptide (P3b-T4c) has an extended conformation and lies across the PB2 627-domain (Fig 1D). To create the CTD binding surface requires concerted side chain reorientations of PB2 W553, M572 and W575 (Figs 1D and S2A), allowing P6b to pack on W553 and Y1c on M572 and L561, with its hydroxyl group hydrogen bonding to D571. PB2 K556 forms a salt bridge with pS5b.

Most functional studies on CTD are performed with human or avian influenza A polymerase, whereas CTD binding has only been structurally characterised for bat A/little yellow-shouldered bat/Guatemala/060/2010(H17N10) polymerase [32] and C/Johannesburg/1/1966
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Structural analysis of CTD binding to influenza B polymerase. **A.** Overall view of the crystal structure of influenza B polymerase, with bound vRNA 5' hook (violet) and CTD peptide mimic (slate blue spheres) in site 2B. **B.** Overall view of the crystal structure of influenza B polymerase with bound promoter (violet and yellow), capped primer (blue) and CTD peptide mimic (slate blue spheres) bound in sites 1AB and 2B. The polymerase is coloured as in (A). The N and C-termini of the two CTD fragments are marked and the red dotted line shows the shortest connection between them with directionality indicated by the arrow. **C.** Overall view of the crystal structure of bat influenza A polymerase with bound promoter and CTD peptide mimic (slate blue spheres) bound in sites 1AB and 2A ([32], PDB: 5M3H). The colour code is as in (A). The N and C-termini of the two CTD fragments are marked and the red dotted line shows the shortest connection between them with directionality indicated by the arrow. **D.** Details of the interaction between key residues of the influenza B polymerase PA subunit (green), PB2-N (red) and PB2-627 (deep-salmon) with the CTD peptide (slate blue sticks) in site 2B. Three CTD repeats denoted a (orange), b (cyan) and c (dark green) are involved in this interaction. Hydrogen bonds are indicated as dotted green lines. **E.** Sequence alignment of the CTD binding sites in the A/WSN/33 (A0A2Z5U3X0) and B/Memphis/13/2003 (Q5V8X3) polymerase subunits PA and PB2. Protein sequences were obtained from UniProt (https://www.uniprot.org/) and aligned with SnapGene</figcaption></figure>
6.0. Key residues for CTD binding are indicated in bold. Identical, similar and non-similar residues are highlighted in purple, light blue and orange, respectively. Grey boxes indicate residues that form a loop. Residues submitted to mutagenesis in this study are indicated with their numbers above (FluPol<sub>A</sub>) and below (FluPol<sub>B</sub>) the alignment, respectively.

https://doi.org/10.1371/journal.ppat.1010328.g001

[33]. Although sequence alignments and mutational studies strongly suggest that the mode of CTD binding is conserved for all IAV-like polymerase [32], we attempted to confirm this by determining the structure of a CTD mimicking peptide bound to influenza A/Zhejiang/DTID-ZJU01/2013(H7N9) polymerase. Previously, we have reported crystals of the A/H7N9 core (PA 201–716, PB1 full-length, PB2 1–127) in the apo-state, which forms symmetrical dimers as described elsewhere [18,58,59]. Here, we soaked the four-repeat pS5 CTD peptide mimic into co-crystals of H7N9 core with the vRNA 5' hook. The crystals diffracted to a maximum resolution of 3.41 Å (Table 1) and again contain symmetrical dimers of the polymerase core (S2B Fig). We observed clear electron density, not only for the 5' hook, but also for the CTD peptide bound in site 2A (S1B and S2C Figs), essentially identically bound as previously seen for bat influenza A/H17N10 polymerase (S2D Fig). However, there was no CTD peptide bound in site 1AB. The most likely explanation for this is that in the symmetrical dimeric form of influenza A (core only, or full trimer), both polymerases are in the so-called ‘dislocated’ conformation [18] with an open active site. In particular, PA regions 425–452 and 586–618 are rotated by ~20°, compared to the active, monomeric promoter bound state (e.g. A/H3N2 polymerase structure, [58], PDB:6RR7). This particularly affects the position of key site 1AB binding site residues Y445, E449 and F612 (S2D Fig), thus preventing CTD at this site, while not affecting binding to site 2A.

Alignments of PB2 and PA amino acid sequences from the representative influenza A and B viruses used in our study are shown in S3 and S4 Figs, **respectively**. The PB2 and PA sequence are very conserved, with > 94% identity among the A/WSN/1933(H1N1), A/NT/60/1968(H3N2) and A/Zhejiang/DTID-ZJU01/2013(H7N9) viruses (details provided in S3 Fig), and > 99% identity between the B/Memphis/13/2003 and B/Brisbane/60/2008 viruses. Key residues for CTD binding (highlighted in S3 and S4 Figs) are strictly conserved between the two representative type B viruses on one hand, and among the H1N1, H3N2 and H7N9 type A viruses, on the other hand. Most of them are identical or similar in the A/little yellow shoul-dered bat/Guatemala/060/2010 (H17N10) virus compared to the other type A viruses. When compared across type B and type A virus sequences, key residues for CTD binding appear mostly conserved in site 1AB, non-conserved in site 2A and partially conserved in site 2B (Figs 1E and S5).

**The FluPol-CTD interaction can be monitored using a cell-based luciferase complementation assay**

To confirm the structural findings of this study and investigate the distinctive features of CTD binding sites in FluPol<sub>B</sub> and FluPol<sub>A</sub> in the cellular context, we set up a CTD-binding assay using the *Gaussia princeps* luciferase trans-complementing fragments (G1 and G2) [43]. The full-length CTD was fused to G2 and the SV40 nuclear localization signal (G2-CTD). PB2 or PA were fused to G1 at their C-terminus (PB2-G1 and PA-G1, respectively, schematically rep-resented in red in Fig 2A and 2B) as FluPol was shown to retain activity when tagged at these sites [44]. Upon co-expression of G2-CTD and the three polymerase subunits (including PB2-G1 or PA-G1), a luminescence signal resulting from the FluPol-CTD interaction was measured, which was generally higher for FluPol<sub>B</sub> compared to FluPol<sub>A</sub> (Fig 2A and 2B). The interaction signal decreased when PB1 was omitted, in agreement with previous reports that
<figure><img src="image_5.png" /><figcaption>**Fig 2.** *Gaussia princeps* luciferase-based FluPol—CTD binding assay. A-B. G2-CTD was expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA of FluPol~B~ (B/Memphis/13/2003, blue bars) or FluPol~A~ (A/WSN/1933, grey bars). Either PB2 (A, hatched bars) or PA (B, filled bars) were C-terminally tagged in frame with G1. As controls, the wild-type (wt) PB1 was replaced by the catalytic inactive PB1 D444A D445A mutant (i) or was omitted (-). Luciferase activities were measured in cell lysates at 24 hpt. Normalised luciferase ratios (NLRs) were calculated as described in the Materials and Methods section. The data shown are the mean ± SD of at least three independent experiments performed in technical triplicates. **p** ≤ 0.002; ***p** ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test). C-D. The CTD binding of FluPol~A~ mutants PA K289A and R638A was investigated. HEK-293T cells were transfected as described in (A) and (B), respectively. Relative Light Units (RLUs) are expressed as percentages relative to the FluPol~A~ PA wt. The data shown are the mean ± SD of three independent experiments performed in technical triplicates. **p** ≤ 0.002; ***p** ≤ 0.001 (one-way ANOVA; Dunnett's multiple comparisons test). The dotted line labelled "Ctrl" indicates the background signal i.e. the sum of the luminescence activities measured in control samples co-transfected with either the FluPol-G1 and G2 plasmids, or the G1 and G2-CTD plasmids. E. Schematic representation of the CTD constructs used in (F) and (G): the wild-type G2-CTD (wt, top) and the G2-CTD in which all serine 5 residues were replaced with an alanine (S5A, bottom). F-G. The interaction of the wt or the S5A mutated CTD to FluPol~B~ (F) or FluPol~A~ (G) was investigated by transient transfection in HEK-293T cells as described in (A-B). The data shown are the mean ± SD of four independent experiments performed in technical triplicates. \*p* ≤ 0.033; \*\*p* ≤ 0.001 (two-way ANOVA; Sidak's multiple comparisons test). In parallel, cell lysates were analysed by western blot using antibodies specific for the pS5 or pS2 CTD, *G.princeps* luciferase (Gluc) and tubulin. The slow- and fast-migrating bands detected with the pS2 CTD antibody, which likely correspond to the hyperphosphorylated and hypophosphorylated forms of the CTD, respectively, are indicated by a star and a triangle, respectively. The smeared signal in the G2-S5A-CTD samples likely corresponds to the detection of phosphorylation intermediates.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g002

the FluPol-CTD interaction depends on FluPol assembly [28], and it was independent of Flu-
Pol catalytic activity (Fig 2A and 2B, PB1 D444A-D445A mutant [60]). When key CTD-con-
tacting residues of FluPol~A~ were mutated (PA K289A and R638A [32]), the interaction
signal was significantly decreased compared to PA wt (Fig 2C and 2D). To test whether the
FluPol-CTD binding assay reflects the dependency on the phosphorylation of the CTD S5 moiety, all S5 residues of the CTD were mutated to alanine (schematically represented in Fig 2E), which prevented S5 phosphorylation as documented by western blot (Fig 2F and 2G, bottom). Although the wt and S5A CTD showed similar steady-state levels of expression, the binding of the S5A CTD to FluPolB and FluPolA was significantly decreased compared to WT CTD. Consistently, pharmacological inhibition of CDK7, which represents the major kinase for CTD S5 phosphorylation [61], specifically reduced FluPolA/B binding to the CTD but not to the FluPolA interaction partner NUP62 [44] (S6A and S6B Fig). Overall, these data demonstrate that the FluPol-CTD interaction can be accurately monitored in cells using our split-luciferase complementation assay conditions.

## Structure-driven mutagenesis confirms FluPol<sub>B</sub> and FluPol<sub>A</sub> have distinct CTD binding modes on PA

To systematically assess *in vivo* FluPol-CTD binding, we mutated key residues forming the CTD binding sites in the FluPol<sub>B</sub> and/or FluPol<sub>A</sub> co-crystal structures and measured the impact of these mutations on CTD-binding using the split-luciferase complementation assay described above (Fig 2). In parallel, we investigated polymerase activity in a minireplicon assay using a viral-like Firefly luciferase reporter RNA, and we rescued recombinant mutant IBVs and IAVs and measured plaque diameters on reverse genetic supernatants as a read-out for viral growth capacity. The A/WSN/1933 residue nature and numbering is used in the text and figures, except when indicated.

The structure and key residues of the CTD binding site 1AB are conserved between FluPol<sub>B</sub> and FluPol<sub>A</sub> (Figs 1E and 3A). We mutated the pS5 interacting residues PA K635 and R638 to alanines. The mutations did not affect PA accumulation levels (Fig 3B) but significantly decreased *in vivo* binding to the full-length CTD for both FluPol<sub>B</sub> and FluPol<sub>A</sub> (Figs 3C and S7A), which is in line with biochemical data obtained *in vitro* with CTD-mimicking peptides [32]. Consistently, the corresponding recombinant mutant IBVs and IAVs were attenuated or could not be rescued (Fig 3D), and the mutant FluPol<sub>A</sub> activity was generally reduced (Fig 3E). The CTD-binding and minireplicon data were consistent for the PA-K635A-R638A double mutant. However, the single mutants had a stronger impact on FluPol<sub>B</sub> than FluPol<sub>A</sub> in the CTD-binding assay whereas the opposite trend was observed in the minigenome assay (Fig 3C and 3E). The PA-R638A mutation, which decreased FluPol<sub>B</sub> CTD-binding (Fig 3C) and reduced viral growth (Fig 3D), even increased FluPol<sub>B</sub> activity in the minigenome assay (Fig 3E). We and others have previously documented discrepancies between FluPol activity as measured in a minigenome assay and viral growth capacity (e.g. [62,63]). The most likely explanation for our data is that each of the PA-K635A and PA-R638A mutations on its own has no or a mild effect on CTD-binding when FluPol<sub>B</sub> is incorporated into a transiently expressed vRNP.

CTD binding site 2A differs substantially between FluPol<sub>B</sub> and FluPol<sub>A</sub> (Figs 1E and 4A). We introduced mutations at residues PA K289, R454 and S420, which are critical for CTD binding to FluPol<sub>A</sub> and are not conserved in FluPol<sub>B</sub>, and we deleted the PA 550 loop, which buttresses the CTD in FluPol<sub>A</sub> and is considerably shortened in FluPol<sub>B</sub>. These modifications did not affect PA accumulation levels (Fig 4B) and specifically decreased *in vivo* CTD binding of FluPol<sub>A</sub> but not FluPol<sub>B</sub> (Figs 4C and S7B). Consistently, the PA R454A and S420E mutations in the IBV background (IBV numbering: K450A and K416E) did not impair viral growth (Fig 4D) nor did PA K450A affect FluPol<sub>B</sub> polymerase activity (Fig 4E). The PA S420E and PA 550 loop deletion impaired FluPol<sub>B</sub> activity (Fig 4E), indicating that they hinder a function of the polymerase besides CTD binding.
<figure><img src="image_6.png" /><figcaption>**Fig 3.** FluPol<sub>B</sub> and FluPol<sub>A</sub> CTD-binding mode at site 1AB. **A.** Superposition of the similar CTD binding in sites 1AB on the PA subunit for influenza B (B/Memphis/13/2003, green) and bat influenza A (A/little yellow-shouldered bat/Guatemala/060/2010(H17N10), light grey) polymerases with the CTD peptide as a thin tube (respectively slate blue and light grey). Key conserved residues are indicated in their respective colours, as well as the FluPol<sub>B</sub>-specific insertion (PA 608 loop) that is important for part of site 2B. See sequence alignment in Fig 1E. **B.** HEK-293T cells were transfected with the indicated FluPol<sub>A</sub> (A/WSN/1933) and FluPol<sub>B</sub> (B/Memphis/13/2003) site 1AB mutants, which were C-terminally tagged with the G1 fragment. Cells were lysed at 24 hpt and analysed by western blot using antibodies specific for *G.princeps* luciferase (PA-G1) and tubulin. The residue numbering corresponds to FluPol<sub>A</sub> (A/WSN/1933). **C.** *In vivo* CTD binding of the indicated mutants of FluPol<sub>A</sub> (A/WSN/1933, grey bars) and FluPol<sub>B</sub> (B/Memphis/13/2003, blue bars). The G2-tagged CTD was expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA-G1. RLUs are expressed as percentages relative to wt FluPol<sub>A/B</sub>. The data shown are mean ± SD of four independent experiments performed in technical triplicates. The dotted line labelled "Ctrl" indicates the background signal i.e. the sum of the luminescence activities measured in control samples co-transfected with either the FluPol-G1 and G2 plasmids or the G1 and G2-CTD plasmids. ***p* ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test). **D.** Characterisation of recombinant IAV (A/WSN/1933, grey dots) and IBV (B/Brisbane/60/2008, blue dots) viruses. Recombinant viruses with the indicated mutations were generated by reverse genetics as described in the Materials and Methods section. Reverse genetic supernatants were titrated on MDCK cells, stained at 72 hpi and plaque diameters were determined using the Fiji software. Each dot represents the diameter of a viral plaque relative to the mean plaque size of IAV wt or IBV wt recombinant virus. (#) not measurable pinhead-sized plaque diameter; (♀) no viral rescue. **E.** Polymerase activity of CTD-binding site 1AB mutants. FluPol<sub>A</sub> (A/WSN/1933, grey bars) or FluPol<sub>B</sub> (B/Memphis/13/2003, blue bars) was reconstituted in HEK-293T cells by transient transfection of PB2, PB1, PA, NP and a model RNA encoding the Firefly luciferase flanked by the 5' and 3' non-coding regions of the IAV or IBV NS segments, respectively. As an internal control, a RNA-Polymerase II promotor driven Renilla plasmid was used. Luminescence was measured at 24 hpt as described in the Materials and Methods section. Firefly activity was normalised to Renilla activity and is shown as percentages relative to wt FluPol<sub>A/B</sub>. The data shown are the mean ± SD of three independent experiments performed in technical duplicates. The dotted lines indicate the background signals in cells transfected with FluPol<sub>A</sub> or FluPol<sub>B</sub> plasmids minus the PB2<sub>A</sub> or PB2<sub>B</sub> plasmid, respectively. ***p* ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g003

## The PB2 627 domain is involved in CTD binding for FluPol<sub>B</sub> but not FluPol<sub>A</sub>

The key CTD binding residues and the 3D structure of site 2B are often conserved between FluPol<sub>B</sub> and FluPol<sub>A</sub> (Figs 1E and 5A). However, CTD binding at site 2B has never been observed *in vitro* with FluPol<sub>A</sub>, and the inserted PA 608 loop (IBV numbering) which butresses the CTD at the junction between the PA-Cter and PB2 627 domains in FluPol<sub>B</sub> is absent in FluPol<sub>A</sub>. The PA R608A mutation significantly decreased FluPol<sub>B</sub> CTD binding in our cell-based complementation assay (Fig 5B **middle panel**, no counterpart residue in FluPol<sub>A</sub>). Consistently, the corresponding recombinant mutant IBV could not be rescued upon reverse genetics (S8 Fig), and the PA R608 mutant FluPol<sub>B</sub> showed reduced polymerase activity (Fig 5B, **right panel**). We then mutated to alanines the residues W552 and R555 (W553 and R556 according to IBV numbering), which are located on the PB2 627 domain, make contact with the CTD pS5 in the FluPol<sub>B</sub> co-crystal and are conserved between IBVs and IAVs (Figs 1E and 5A). The mutations did not affect PB2 accumulation levels (Fig 5C) and either decreased
<figure><img src="image_7.png" /><figcaption>**Fig 4.** FluPol<sub>B</sub> and FluPol<sub>A</sub> CTD-binding mode at site 2A. A. Superposition of CTD peptide (slate blue tube) bound at site 2A of the PA subunit of FluPol<sub>A</sub> (A/Zhejiang/DTID-ZJU01/2013(H7N9), green) with the equivalent region of FluPol<sub>B</sub> (B/Memphis/13/2003, wheat), showing similarities and differences in CTD interacting residues. See sequence alignment in Fig 1E. B-E. Protein expression (B), *in vivo* CTD binding (C), characterisation of recombinant IAV and IBV viruses (D) and polymerase activity (E) of CTD-binding site 2A mutants. Experiments were performed as described in Fig 3B-3E for FluPol<sub>B</sub> and FluPol<sub>A</sub> site 1AB mutations. C. **p** ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test). D. (#) not measurable pinhead-sized plaque diameter; (‡) no viral rescue, (n.d.) not determined. E. **p** ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g004

(R555A) or increased (W552A) CTD binding of FluPol<sub>B</sub>, whereas they had no effect on FluPol<sub>A</sub> CTD binding (Figs 5D and S7C). We speculate that the W552A mutation (W553A according to FluPol<sub>B</sub> numbering) could make it easier for the FluPol<sub>B</sub> PB2 627 domain to adopt a CTD-bound conformation. To rule out any CTD binding activity on the FluPol<sub>A</sub> PB2 627 domain, we deleted the whole domain as described before [64]. The deletion strongly and specifically decreased CTD binding to FluPol<sub>B</sub> but not to FluPol<sub>A</sub> (Fig 5E). Nevertheless, single amino acid substitutions at residues PB2 W552 and R555 impaired viral growth and polymerase activity of FluPol<sub>B</sub> as well as FluPol<sub>A</sub>, however with weaker effects on FluPol<sub>A</sub> (Figs 1F and 5G). Given the multiple functions attributed to the PB2 627 domain [64], the most likely interpretation of our data is that residues on the PB2 627 domain contribute to the CTD recruitment exclusively for IBVs while they have overlapping CTD-unrelated functions for IBVs and IAVs.

We asked whether this major difference between FluPol<sub>B</sub> and FluPol<sub>A</sub> CTD binding modes results in different levels of transcriptional activation by CTD mimicking peptides *in vitro*. A model has been proposed for FluPol<sub>C</sub> in which the CTD stabilizes a transcription-competent conformation by binding at the interface of PB1, P3 (PA equivalent), and the flexible PB2 C-terminus [33]. Our observations suggest that the same model could apply to FluPol<sub>B</sub> and not to FluPol<sub>A</sub>. Therefore, we tested the impact of pS5 CTD mimicking peptides of varying lengths (two, four, or six YSPTSPS repeats) on FluPol<sub>B</sub> and FluPol<sub>A</sub> *in vitro* transcriptional activity (Fig 6). The FluPol<sub>B</sub> *in vitro* endonuclease activity (Fig 6A, lane 4) and elongation activity (Fig 6A, lane 8) were increased in the presence of the six-repeat pS5 CTD mimicking peptide compared to the mock control (Fig 6A, lanes 1 and 5), and a similar trend was observed with FluPol<sub>A</sub> (Fig 6B, lanes 4 and 8 compared to lanes 1 and 5, respectively). These data complement previous reports that CTD pS5 binding facilitates FluPol<sub>A</sub> and FluPol<sub>C</sub> transcriptional activity [33] and strengthen the hypothesis that the CTD stabilises FluPol in a transcription-competent conformation [20]. However, our finding that the
<figure><img src="image_8.png" /><figcaption>**Fig 5.** FluPol<sub>A</sub> and FluPol<sub>B</sub> CTD-binding mode at site 2B. **A.** Superposition of CTD peptide (slate blue tube) bound in site 2B of FluPol<sub>B</sub> (B/Memphis/13/2003, PA green, PB2-N red, PB2-627 deep salmon) with the equivalent region of FluPol<sub>A</sub> (A/NT/60/1968 (H3N2), [58], PDB: 6RR7, PA light green, PB2 pink), showing similarities and differences in CTD interacting residues. See sequence alignment in Fig 1E. **B.** Protein expression (left), *in vivo* CTD binding (middle), and polymerase activity (right) of FluPol<sub>B</sub> PA R608A. Experiments were performed as described in Fig 3. The data shown are mean ± SD of three independent experiments performed in technical triplicates. \*\*\*p ≤ 0.001 (unpaired t test). **C-D.** Protein expression (C) and *in vivo* CTD binding (D) of CTD-binding site 2B mutants. Experiments were performed as described in Fig 3B and 3C for FluPol<sub>B</sub> and FluPol<sub>A</sub> site 1AB mutations. D. \*\*p ≤ 0.002, \*\*\*p ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test). E. *In vivo* CTD binding of PB2 Δ627 domain deletion mutants was investigated. G2-CTD was expressed by transient transfection in HEK-293T cells together with PB2-G1, PB1 and PA of FluPol<sub>A</sub> (A/WSN/1933, grey bars) or FluPol<sub>B</sub> (B/Memphis/13/2003, blue bars). Luciferase activities were measured in cell lysates at 24 hpt. Normalised luciferase ratios (NLRs) were calculated as described in the Materials and Methods section. \*\*\*p ≤ 0.001 (two-way ANOVA; Sidak's multiple comparisons test). Cell lysates were analysed in parallel by western blot with antibodies specific for the pS5 CTD, *G. princeps* luciferase (PB2-G1) and tubulin. **F-G.** Characterisation of recombinant IAV and IBV viruses (F) and polymerase activity (G) of CTD-binding site 2B mutants. Experiments were performed as described in Fig 3D-E for FluPol<sub>B</sub> and FluPol<sub>A</sub> site 1AB mutations. F. (Ψ) no viral rescue. G. \*\*\*p ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g005

FluPol<sub>A</sub> PB2 627 domain has no CTD binding activity (Fig 5E) indicates that FluPol<sub>A</sub> has evolved a divergent mechanism by which the CTD stabilizes the FluPol transcriptase. It also questions whether bridging of the PB2 627 and PA-Cter domains *per se* is needed for transcriptional activation of FluPols.

## FluPol<sub>B</sub> and FluPol<sub>A</sub> bind to the host RNAP II independently of the CTD

The RNAP II transcriptional machinery is highly conserved across eukaryotes, and the CTD in particular shows almost no sequence differences among vertebrate (mammalian and avian) host species susceptible to IAV or IBV infection (S9 Fig). It is therefore unlikely that differences in the CTD amino acid sequence drove the evolution of divergent IBV and IAV CTD-binding modes. We investigated whether FluPol<sub>A/B</sub> can interact with the two major RNAP II subunits (RPB1, RPB2) independently of the CTD, using the split-gaussia luciferase complementation assay. The G2-RPB1 and RPB2-G2 fusion proteins were co-expressed with
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig 6.** Effect of CTD pS5 peptides on *in vitro* endonuclease and transcription activity of FluPol<sub>B</sub> and FluPol<sub>A</sub>. A-B. CTD pS5 peptides of different lengths (two, four and six repeats) were added to (A) FluPol<sub>B</sub> (B/Memphis/13/2003) or (B) FluPol<sub>A</sub> (bat influenza A (A/little yellow-shouldered bat/ Guatemala/060/2010(H17N10))) *in vitro* activity reactions as described in Materials and Methods. The left four lanes show endonuclease and right four lanes transcription reactions. Lanes 1 and 5: no CTD pS5 peptides. Lanes 2–6, 3–7 and 4–8: addition of two-, four- and six-repeat CTD pS5 peptides, respectively. The different RNA species are labeled. In the FluPol<sub>B</sub> samples, a secondary transcription product larger in size than those expected from the size of the template is detected, most likely resulting from stable hybridization of template and primary product as previously described [18]. Quantification of the reaction products of four independent experiments is shown below (FluPol<sub>B</sub> in blue and FluPol<sub>A</sub> in grey, respectively). The products of the reactions are normalised to the total RNA amount for each reaction and are presented as fractions of the activity of the reaction without peptide. \*p ≤ 0.033, \*\*p ≤ 0.001 (one-way ANOVA; Dunnett's multiple comparisons test).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g006

G1-tagged FluPol (PA-G1) by transient transfection as described above. Both combinations resulted in robust and comparable interaction signals (Fig 7A). Interestingly, in the presence of a truncated RPB1 deleted from the CTD (RPB1ΔCTD), a stable interaction signal with FluPol<sub>A/B</sub> could still be measured, most likely corresponding to a direct interaction between FluPol<sub>A/B</sub> and the core domain of RPB1. RNPA II can form clusters in the nucleus [65] and bridging of the G2-RPB1ΔCTD to G1-tagged FluPol through a wild-type RNAP II complex cannot be formally excluded. However, it is unlikely that an active *G. princeps* luciferase can be reconstituted in these conditions, given the large size of the complexes involved (RNAPII, > 500 kD; FluPol, > 240 kD). Moreover, as the G2-RPB1ΔCTD protein accumulates at substantially higher levels than the endogenous RPB1 protein (Fig 7A), it is likely to compete with the endogenous RPB1 protein for assembly into RNPAII complexes and their subsequent import into the nucleus, and therefore RNAP II nuclear clusters are likely to contain a majority of G2-tagged RNPA II complexes. Mutations in site 1AB which reduced CTD binding to background levels (Fig 3C) had only weak effects on RPB1 and RPB1ΔCTD binding (PA K631A R634 in Fig 7B, PA K635A R638 in Fig 7C). The same was observed with site 2A mutations (PA K450A in Fig 7B, PA R454A in Fig 7C) and site 2B mutations and (PB2 R555A Fig 7B,
<figure><img src="image_11.png" /><figcaption>**Fig 7. FluPol<sub>B</sub> and FluPol<sub>A</sub> binding to RPB1, RPB1ΔCTD and RPB2. A.** Binding of FluPol<sub>A</sub> (A/WSN/1933, grey bars) and FluPol<sub>B</sub> (B/Memphis/13/2003, blue bars) to RPB1, RPB1ΔCTD and RPB2 was evaluated. RPB1, RPB1ΔCTD and RPB2 were tagged with G2 and expressed by transient transfection in HEK-293T cells together with PB2, PB1 and PA-G1. Normalised luciferase activities (NLRs) were calculated as described in the Materials and Methods section. The data shown are the mean ± SD of five independent experiments performed in technical triplicates. \*p ≤ 0.033, \*\*p ≤ 0.002 (two-way ANOVA; Dunnett's multiple comparisons test). Cell lysates were analysed in parallel by western blot with antibodies specific for the pS5 CTD, *G. princeps* luciferase and tubulin. Lanes 1 to 3 and 4 to 6 correspond to the co-transfection of FluPol<sub>A</sub> and FluPol<sub>B</sub>, respectively, with the G2-RPB1, -RPB1ΔCTD and -RPB2 plasmids as indicated. The band detected in lanes 2–3 and 5–6 with the pS5 CTD antibody correspond to the endogenous RPB1 protein. **B–C.** Binding of (B) FluPol<sub>B</sub> (B/Memphis/13/2003) and (C) FluPol<sub>A</sub> (A/WSN/1933) mutants in site 2B (PB2 K556/R555A), site 1 (PA K631A R634A/ PA K635A R638A) and site 2A (PA K450A / PA R454A) to RPB1, RPB1ΔCTD and RPB2 was evaluated as described in (A).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010328.g007

PB2 K556A in Fig 7C). These findings, taken together with the relatively low affinity of FluPol for pS5 CTD peptides [32], suggest that the CTD is not the only interface between FluPol and the host RNAP II, and it may not be essential to connect FluPol to the RNAP II but rather to coordinate FluPol cap-snatching.

**Discussion**

Here we report co-crystal structures of a human FluPol<sub>B</sub> and an avian (isolated from human) FluPol<sub>A</sub> bound to pS5 CTD mimicking peptides. We uncover the conformation and direction-ality of the CTD peptide bound to FluPol<sub>B</sub> at a site that crosses over from the PA-Cter to the PB2 627 domain (site 2B), and has no counterpart on FluPol<sub>A</sub> or FluPol<sub>C</sub>. Two CTD binding sites have been characterised on FluPol<sub>A</sub> (sites 1A and 2A) ([32] and this study) and on Flu-Pol<sub>C</sub> (sites 1C and 2C, distinct from sites 1A and 2A) [33]. On the FluPol<sub>B</sub> co-crystal structure, site 1B is similar to site 1A, whereas site 2B is distinct from site 2A and 2C.

By performing structure-based mutagenesis of FluPol<sub>B</sub> and FluPol<sub>A</sub> followed by a system-atic investigation of FluPol-CTD binding, FluPol transcription/replication activity and viral phenotype, we confirm that CTD binding involves the same key residues at site 1AB for
FluPol<sub>B</sub> and FluPol<sub>A</sub>, but distinct and specific residues at site 2A for FluPol<sub>A</sub> and site 2B for FluPol<sub>B</sub>, respectively. In particular, we demonstrate that the PA 606–609 loop, which buttresses the CTD at the junction between PA and PB2 in the FluPol<sub>B</sub> co-crystal structure and is not conserved in FluPol<sub>A</sub> or FluPol<sub>C</sub>, is an essential component of site 2B.

Our data and others' [32,33] demonstrate that IAVs, IBVs and ICVs have evolved divergent CTD binding modes, and raise questions about the driving force behind this divergent evolution. Large-scale meta-transcriptomic approaches have identified IBV-like and IDV-like viruses in fish and amphibians, suggesting that the influenza viruses of all four genera might be distributed among a much wider range of vertebrate species than recognised so far [6,57]. Phylogenetic analyses, although limited by strong sampling biases across species, indicate that both virus-host co-divergence over long timescales and cross-species transmissions have shaped the evolution of influenza viruses. With one of the two CTD binding sites being conserved between IAVs and IBVs but absent in ICVs, the divergence of the bipartite CTD binding mode apparently matches the evolutionary distance between the three types of influenza viruses [66]. Interestingly however, we demonstrate that, in contrast to what is observed for IBVs and ICVs ([32,33] and this study), the PB2 627 domain is not involved in CTD binding for IAVs. Therefore, from a mechanistic point of view, the CTD-dependent transcriptional activation of FluPol might be closer between IBVs and ICVs than between IBVs and IAVs as a consequence of a distinctive evolutionary pressure exerted on IAVs. The FluPol<sub>A</sub> CTD binding mode presumably reflects an avian-optimised mode and co-evolved with protein interfaces between avian host factors and the PB2 627 domain, known to restrict avian IAV replication in humans (the principal hosts of IBVs and ICVs).

Another example of a functional interaction with RNAP II being achieved through distinct CTD binding is provided by the cellular mRNA capping enzyme (CE). The CEs from *Schizosaccharomyces pombe*, *Candida albicans* and *Mus musculus* were shown to bind directly S5 CTD repeats with very distinct binding interfaces and distinct conformations of the bound CTD [67]. These distantly related species show major differences in the CTD length and sequence [26] which could at least partially account for the divergence in CE-CTD binding modes. In contrast, the CTD is highly conserved among host species susceptible to IAV, IBV and ICV infections (S9 Fig).

There is considerable evidence, however, that the FluPol-CTD interaction is only part of a more complex interaction pattern between the viral and cellular transcription machineries, raising the possibility that interactions between the FluPol and less conserved components of the cellular transcriptional machinery could have indirectly shaped the evolution of distinct CTD binding modes. We observed that a truncated RPB1 subunit, which lacks the CTD, retains partial binding to FluPol (Fig 7). Mass-spectrometry screenings have identified other RNAP II subunits and multiple transcriptional pausing and elongation factors as potential FluPol interaction partners [31]. Host factors involved in transcription such as DDX17 were found to bind FluPol and to determine IAV host-specificity [68]. By analogy, CEs not only bind to the pS5 CTD but also to the transcription pausing DRB Sensitivity-Inducing Factor (DSIF) [67] and make additional direct interactions with the nascent transcript exit site on the body of RNAP II [69]. Likewise, it was shown recently that the integrator complex binds RNAP II in its promotor-proximal paused state through direct interactions with the CTD of RPB1 but also with RPB2, RPB3, and RPB11, and through indirect interaction with the negative elongation factor NELF and DSIF [70]. Intriguingly, FluPol was also found to interact with the DSIF subunit SPT5 [71]. To what extent host-specific features of SPT5 or other cellular factors may have constrained the evolution of CTD-binding sites on FluPol remains to be explored.
We show that the *in vitro* transcriptional activity of FluPol<sub>B</sub> is facilitated by the addition of CTD pS5 mimicking peptides, as reported previously for FluPol<sub>A</sub> and FluPol<sub>C</sub> [33]. The mechanism previously proposed for FluPol<sub>C</sub> [20,33] in which the CTD stabilises FluPol in a transcription-competent conformation by bridging P3 (the PA equivalent for ICVs) and PB2, could possibly apply to FluPol<sub>B</sub> with PA-PB2 bridging occurring at site 2B. Our data show that it does not apply to FluPol<sub>A</sub>, unless another yet unidentified domain of PB2, distinct from the PB2 627 domain, is involved.

As underlined by the different sensitivity of IAV and IBVs to cap-binding inhibitors related to differences in the cap-binding mode of their PB2 subunits [72], a detailed understanding of structural and functional differences between FluPol<sub>A</sub> and FluPol<sub>B</sub> is of significant importance with regard to the development of broad-spectrum antivirals and need to be taken into account when targeting the FluPol-CTD binding interface for antiviral intervention.

**Supporting information**

**S1 Fig. Omit maps for bound CTD peptide structures.** A. CTD peptide bound in site 2B of FluPol<sub>B</sub> (B/Memphis/13/2003) polymerase. Fo-Fc omit map shown at 2.5 σ with clear density for two phosphoserines (pS5a and pS5b). B. CTD peptide bound in site 2A of FluPol<sub>A</sub> (A/Zhejiang/DTID-ZJU01/2013(H7N9)) polymerase. Fo-Fc omit map shown at 2.7 σ with clear density for one phosphoserine.
(TIF)

**S2 Fig. Structural analysis of CTD binding in FluPol<sub>A</sub> and FluPol<sub>B</sub>.** A. Superposition of the PB2-627 domains of CTD-bound (deep salmon with CTD in slate-blue) and unbound (wheat) FluPol<sub>B</sub> (B/Memphis/13/2003) polymerase, showing induced-fit conformational changes of three key hydrophobic residues (PB2 W553, M572, W575). B. Ribbon diagram looking down the 2-fold axis of the symmetrical homodimeric A/H7N9 polymerase core with CTD peptide bound in site 2A. One monomer is coloured with PA (green), PB1 (cyan), PB2 (red) with CTD peptide (slate spheres) and 5' vRNA hook (violet). The second monomer has paler versions of these colours. C. Details of the binding of the CTD peptide (slate blue) in site 2A of FluPol<sub>A</sub> (A/Zhejiang/DTID-ZJU01/2013(H7N9) core) polymerase. D. Comparison of site 1A configuration for CTD bound form of FluPol<sub>A</sub> (bat influenza A (A/little yellow-shouldered bat/Guatemala/060/2010(H17N10) [32], PDB: 5M3H, PA subunit light grey, CTD peptide slate-blue), CTD free, transcription active form of FluPol<sub>A</sub> (A/NT/60/1968 (H3N2), [58], PDB: 6RR7, wheat) and dimeric FluPol<sub>A</sub> (A/Zhejiang/DTID-ZJU01/2013(H7N9) core, light blue, this work). The FluPol<sub>A</sub> (H7N9) polymerase core is the symmetrical dimer with each polymerase in the open, 'dislocated' state [18]. Due to the dislocation, PA regions 425-452 and 586-618 are rotated by ~20°, which particularly effects the position of site 1A binding site residues Y445, E449 and F612. This likely explains the lack of CTD binding observed in site 1A for the H7N9 core, whereas site 2A is undistorted and occupied by CTD (S2B Fig).
(TIF)

**S3 Fig. Sequence alignment of full-length PB2 and PA sequences from the FluPol<sub>A</sub> complexes included in this study.** Protein sequences from A/WSN1933(H1N1), A/NT/60/1968 (H3N2), A/Zhejiang/DTID-ZJU01/2013(H7N9) and A/little yellow-shouldered bat/Guatemala/060/2010(H17N10), were obtained from UniProt (https://www.uniprot.org/), aligned with SnapGene 6.0 and visualized by Esprit 3.0 [55]. Identical and similar residues are indicated in red or yellow, respectively. Key residues for CTD binding are indicated by an orange (site 1AB), black (site 2A) or blue (site 2B) arrowhead. The sequence stretches shown in Fig 1E are highlighted with colored boxes (orange, grey and blue for site 1AB, site 2A and site 2B,
respectively). Residues submitted to mutagenesis or deletion in this study are indicated with their numbers below the alignment. The percentages of sequence identity are provided in a table.

(TIF)

**S4 Fig. Sequence alignment of full-length PB2, -PB1 and PA sequences from the FluPolB complexes included in this study.** Same as in S3 Fig. Memphis: B/Memphis/13/2003. Brisbane: B/Brisbane/60/2008. The percentages of sequence identity are 99.5% (PB2) and 99.3% (PA).

(TIF)

**S5 Fig. Sequence alignment of full-length PB2, PB1 and PA sequences of A/WSN/1933 (A0A2Z5U3X0) and B/Memphis/13/2003 (Q5V8X3).** Protein sequences were obtained from UniProt (https://www.uniprot.org/), aligned with SnapGene 6.0 and visualized by Escript 3.0 [55]. Identical and similar residues are indicated in red or yellow, respectively.

(TIF)

**S6 Fig. FluPol-CTD interaction in the presence of a CDK7 inhibitor.** A. HEK-293T cells were transfected with G2-CTD. At 24 hpt cells were treated for 1 h with 20 µM BS-181-HCl (DMSO final concentration 0.2%). Cell lysates were analysed by western blot with antibodies specific for pS5 or pS2 CTD and tubulin. B. *In vivo* CTD binding of FluPolA (A/WSN/1933) and FluPolB (B/Memphis/13/2003). G2-tagged CTD was expressed by transient transfection in HEK-293T cells together with the viral polymerase subunit PB2, PB1 and PA-G1. At 24 hpt cells were treated for 1 h with 20 µM BS-181-HCl or 0.2% DMSO before cell lysis and measurement of *G. princeps* luciferase activity as described in the Materials and Methods section. As a control, the previously described FluPolA-NUP62 interaction was investigated by co-transfection of G2-NUP62, PB2, PB1 and PA-G1. RLUs are expressed as percentages relative to DMSO treated cells. The data shown are mean ± SD of five independent experiments performed in technical triplicates. ***p ≤ 0.002, **p ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test).

(TIF)

**S7 Fig. In vivo CTD binding of FluPolA and FluPolB mutants.** A-C. *In vivo* CTD binding of the indicated site 1AB (A), site 2A (B) and site 2B (C) mutants of FluPolA (A/WSN/1933, grey hatched bars) and FluPolB (B/Memphis/13/2003, blue hatched bars). The G2-tagged CTD was expressed by transient transfection in HEK-293T cells together with PB2-G1, PB1 and PA. RLUs are expressed as percentages relative to wt FluPolA/B. The data shown are the mean ± SD of at least three independent experiments performed in technical triplicates. **p ≤ 0.002, ***p ≤ 0.001 (two-way ANOVA; Dunnett's multiple comparisons test). (n.d.) not determined.

(TIF)

**S8 Fig. Plaque phenotype of FluPolB PA R608A.** Characterisation of recombinant IBV (B/Brisbane/60/2008) PA R608A mutant virus. Recombinant viruses with the indicated mutations were generated by reverse genetics as described in the Materials and Methods section. Reverse genetic supernatants were titrated on MDCK cells and stained at 72 hpi by crystal violet. The pictures show one representative plaque assay with the indicated ten-fold dilution.

(TIF)

**S9 Fig. Sequence alignment of the RPB1 C-terminal domain (CTD) across species.** The RPB1 CTD sequences of *Sus scrofa* (wild boar), *Equus caballus* (horse), *Homo sapiens* (human), *Canis lupus* (wolf), *Mus musculus* (house mouse), *Gallus gallus* (chicken) and *Anas platyrhynchos* (wild duck) were obtained as described in the Materials and Methods section, aligned with SnapGene
6.0 and visualised by Esprit 3.0 [55]. The CTD repeat numbers are indicated below the sequence alignment. Identical and similar residues are indicated in red or yellow, respectively. (TIF)

**Acknowledgments**

We thank Pr. Daniel Perez (University of Georgia), Dr. Bernard Delmas (INRAE) and Dr. Yves Janin (Institut Pasteur) for sharing material and reagents, and Dr. Yves Jacob for helpful discussions.

**Author Contributions**

**Conceptualization:** Tim Krischuns, Catherine Isel, Stephen Cusack, Nadia Naffakh.

**Formal analysis:** Tim Krischuns, Catherine Isel, Maria Lukarska, Alexander Pflug, Stephen Cusack, Nadia Naffakh.

**Funding acquisition:** Stephen Cusack, Nadia Naffakh.

**Investigation:** Tim Krischuns, Catherine Isel, Petra Drncova, Maria Lukarska, Alexander Pflug, Sylvain Paisant, Vincent Navratil.

**Methodology:** Tim Krischuns, Vincent Navratil, Stephen Cusack, Nadia Naffakh.

**Supervision:** Stephen Cusack, Nadia Naffakh.

**Writing – original draft:** Tim Krischuns, Stephen Cusack, Nadia Naffakh.

**Writing – review & editing:** Tim Krischuns, Catherine Isel, Petra Drncova, Maria Lukarska, Vincent Navratil, Stephen Cusack, Nadia Naffakh.

**References**

1. Krammer F, Smith GJD, Fouchier RAM, Peiris M, Kedzierska K, Doherty PC, et al. Influenza. Nat Rev Dis Primers. 2018; 4:3. https://doi.org/10.1038/s41572-018-0002-y PMID: 29955068

2. Caini S, Kroneman M, Wiegers T, El Guerche-Séblain C, Paget J. Clinical characteristics and severity of influenza infections by virus type, subtype, and lineage: A systematic literature review. Influenza Other Respi Viruses. 2018; 12:780–92. https://doi.org/10.1111/IRV.12575 PMID: 29858537

3. Zaraket H, Hurt AC, Clinch B, Barr I, Lee N. Burden of influenza B virus infection and considerations for clinical management. Antiviral Res. 2021;185. https://doi.org/10.1016/j.antiviral.2020.104970 PMID: 33159999

4. Yoon SW, Webby RJ, Webster RG. Evolution and ecology of influenza A viruses. Curr Top Microbiol Immunol. 2014; 385:359–75. https://doi.org/10.1007/82_2014_396 PMID: 24990620

5. Suzuki Y, Nei M. Origin and evolution of influenza virus hemagglutinin genes. Mol Biol Evol. 2002; 19:501–9. https://doi.org/10.1093/oxfordjournals.molbev.a004105 PMID: 11919291

6. Parry R, Wille M, Turnbull OMH, Geoghegan JL, Holmes EC. Divergent Influenza-Like Viruses of Amphibians and Fish Support an Ancient Evolutionary Association. Viruses. 2020; 12. https://doi.org/10.3390/v12091042 PMID: 32962015

7. Baker SF, Nogales A, Finch C, Tuffy KM, Domm W, Perez DR, et al. Influenza A and B virus intertypic reassortment through compatible viral packaging signals. J Virol. 2014; 88:10778–91. https://doi.org/10.1128/JVI.01440-14 PMID: 25008914

8. Eisfeld AJ, Neumann G, Kawaoka Y. At the centre: influenza A virus ribonucleoproteins. Nat Rev Microbiol. 2015; 13:28–41. https://doi.org/10.1038/nrmicro3367 PMID: 25417656

9. Coloma R, Arranz R, de la Rosa-Trevín JM, Sorzano COS, Munier S, Carlero D, et al. Structural insights into influenza A virus ribonucleoproteins reveal a processive helical track as transcription mechanism. Nat Microbiol. 2020; 5:727–34. https://doi.org/10.1038/s41564-020-0675-3 PMID: 32152587

10. Moeller A, Kirchdoerfer RN, Potter CS, Carragher B, Wilson IA. Organization of the influenza virus replication machinery. Science. 2012; 338:1631–4. https://doi.org/10.1126/science.1227270 PMID: 23180774
11. Arranz R, Coloma R, Chichón FJ, Conesa JJ, Carrascosa JL, Valpuesta JM, et al. The structure of native influenza virion ribonucleoproteins. Science. 2012; 338:1634-7. https://doi.org/10.1126/science.1228172 PMID: 23180776

12. Wandzik JM, Kouba T, Cusack S. Structure and Function of Influenza Polymerase. Cold Spring Harb Perspect Med. 2021; 11. https://doi.org/10.1101/cshperspect.a038372 PMID: 32341065

13. Deng T, Vreede FT, Brownlee GG. Different de novo initiation strategies are used by influenza virus RNA polymerase on its cRNA and viral RNA promoters during viral RNA replication. J Virol. 2006; 80:2337-48. https://doi.org/10.1128/JVI.80.5.2337-2348.2006 PMID: 16474140

14. Oymans J, Te Velthuis AJW. A Mechanism for Priming and Realignment during Influenza A Virus Replication. J Virol. 2018;92. https://doi.org/10.1128/JVI.01773-17 PMID: 29118119

15. Plotch SJ, Bouloy M, Krug RM. Transfer of 5'-terminal cap of globin mRNA to influenza viral complementary RNA during transcription in vitro. Proc Natl Acad Sci U S A. 1979; 76:1618-22. https://doi.org/10.1073/pnas.76.4.1618 PMID: 287003

16. Bouloy M, Plotch SJ, Krug RM. Both the 7-methyl and the 2'-O-methyl groups in the cap of mRNA strongly influence its ability to act as primer for influenza virus RNA transcription. Proc Natl Acad Sci U S A. 1980; 77:3952-6. https://doi.org/10.1073/pnas.77.7.3952 PMID: 6933444

17. Poon LL, Pritlove DC, Fodor E, Brownlee GG. Direct evidence that the poly(A) tail of influenza A virus mRNA is synthesized by reiterative copying of a U track in the virion RNA template. J Virol. 1999; 73:3473-6. https://doi.org/10.1128/JVI.73.4.3473-3476.1999 PMID: 10074205

18. Wandzik JM, Kouba T, Karuppasamy M, Pflug A, Drncova P, Provaznik J, et al. A Structure-Based Model for the Complete Transcription Cycle of Influenza Polymerase. Cell. 2020; 181:877-893.e21. https://doi.org/10.1016/j.cell.2020.03.061 PMID: 32304664

19. Bouloy M, Plotch SJ, Krug RM. Globin mRNAs are primers for the transcription of influenza viral RNA in vitro. Proc Natl Acad Sci U S A. 1978; 75:4886-90. https://doi.org/10.1073/pnas.75.10.4886 PMID: 283399

20. Walker AP, Fodor E. Interplay between Influenza Virus and the Host RNA Polymerase II Transcriptional Machinery. Trends Microbiol. 2019; 27:398-407. https://doi.org/10.1016/j.tim.2018.12.013 PMID: 30642766

21. Dias A, Bouvier D, Crépin T, McCarthy AA, Hart DJ, Baudin F, et al. The cap-snatching endonuclease of influenza virus polymerase resides in the PA subunit. Nature. 2009; 458:914-8. https://doi.org/10.1038/nature07745 PMID: 19194459

22. Kouba T, Drncová P, Cusack S. Structural snapshots of actively transcribing influenza polymerase. Nat Struct Mol Biol. 2019; 26:460-70. https://doi.org/10.1038/s41594-019-0232-z PMID: 31160782

23. Worch R, Niedzwiecka A, Stepinski J, Mazza C, Jankowska-Anyzka M, Darzynkiewicz E, et al. Specificity of recognition of mRNA 5' cap by human nuclear cap-binding complex. RNA. 2005; 11:1355-63. https://doi.org/10.1261/rna.2850705 PMID: 16043498

24. Cramer P, Bushnell DA, Kornberg RD. Structural basis of transcription: RNA polymerase II at 2.8 angstrom resolution. Science. 2001; 292:1863-76. https://doi.org/10.1126/science.1059493 PMID: 11313498

25. Schüller R, Forné I, Straub T, Schreieck A, Texier Y, Shah N, et al. Heptad-Specific Phosphorylation of RNA Polymerase II CTD. Mol Cell. 2016; 61:305-14. https://doi.org/10.1016/j.molcel.2015.12.003 PMID: 26799765

26. Eick D, Geyer M. The RNA polymerase II carboxy-terminal domain (CTD) code. Chem Rev. 2013; 113:8456-90. https://doi.org/10.1021/cr400071f PMID: 23952966

27. Martínez-Alonso M, Hengrung N, Fodor E. RNA-Free and Ribonucleoprotein-Associated Influenza Virus Polymerases Directly Bind the Serine-5-Phosphorylated Carboxyl-Terminal Domain of Host RNA Polymerase II. J Virol. 2016; 90:6014-21. https://doi.org/10.1128/JVI.00494-16 PMID: 27099314

28. Engelhardt OG, Smith M, Fodor E. Association of the influenza A virus RNA-dependent RNA polymerase with cellular RNA polymerase II. J Virol. 2005; 79:5812-8. https://doi.org/10.1128/JVI.79.9.5812-5818.2005 PMID: 15827195

29. Vos SM, Farnung L, Urlaub H, Cramer P. Structure of paused transcription complex Pol II-DSIF-NELF. Nature. 2018; 560:601-6. https://doi.org/10.1038/s41586-018-0442-2 PMID: 30135580

30. Chan AY, Vreede FT, Smith M, Engelhardt OG, Fodor E. Influenza virus inhibits RNA polymerase II elongation. Virology. 2006; 351:210-7. https://doi.org/10.1016/j.virol.2006.03.005 PMID: 16624367

31. Krischuns T, Lukarska M, Naffakh N, Cusack S. Influenza Virus RNA-Dependent RNA Polymerase and the Host Transcriptional Apparatus. Annu Rev Biochem. 2021; 90:321-48. https://doi.org/10.1146/annurev-biochem-072820-100645 PMID: 33770447
32. Lukarska M, Fournier G, Pflug A, Resa-Infante P, Reich S, Naffakh N, et al. Structural basis of an essential interaction between influenza polymerase and Pol II CTD. Nature. 2017; 541:117–21. https://doi.org/10.1038/nature20594 PMID: 28002402

33. Serna Martin I, Hengrung N, Renner M, Sharps J, Martínez-Alonso M, Masiulis S, et al. A Mechanism for the Activation of the Influenza Virus Transcriptase. Mol Cell. 2018; 70:1101–1110.e4. https://doi.org/10.1016/j.molcel.2018.05.011 PMID: 29910112

34. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. Phaser crystallographic software. J Appl Cryst. 2007; 40:658–74. https://doi.org/10.1107/S0021889807021206 PMID: 19461840

35. Thierry E, Guilligay D, Kosinski J, Bock T, Gaudon S, Round A, et al. Influenza Polymerase Can Adopt an Alternative Configuration Involving a Radical Repacking of PB2 Domains. Mol Cell. 2016; 61:125–37. https://doi.org/10.1016/j.molcel.2015.11.016 PMID: 26711008

36. Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr. 2004; 60:2126–32. https://doi.org/10.1107/S0907444904019158 PMID: 15572765

37. Murshudov GN, Vagin AA, Dodson EJ. Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr. 1997; 53:240–55. https://doi.org/10.1107/S0907444996012255 PMID: 15299926

38. Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, et al. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr D Biol Crystallogr. 2010; 66:12–21. https://doi.org/10.1107/S0907444909042073 PMID: 20057044

39. Wandzik JM, Kouba T, Cusack S. Structure and Function of Influenza Polymerase. Cold Spring Harb Perspect Med. 2020; a038372. https://doi.org/10.1016/j.cell.2020.03.061 PMID: 32304664

40. Fodor E, Devenish L, Engelhardt OG, Palese P, Brownlee GG, García-Sastre A. Rescue of influenza A virus from recombinant DNA. J Virol. 1999; 73:9679–82. https://doi.org/10.1128/JVI.73.11.9679-9682. 1999 PMID: 10516084

41. Nogales A, Perez DR, Santos J, Finch C, Martínez-Sobrido L. Reverse Genetics of Influenza B Viruses. Methods Mol Biol. 2017; 1602:205–38. https://doi.org/10.1007/978-1-4939-6964-7_14 PMID: 28508223

42. Reich S, Guilligay D, Pflug A, Malet H, Berger I, Crépin T, et al. Structural insight into cap-snatching and RNA synthesis by influenza polymerase. Nature. 2014; 516:361–6. https://doi.org/10.1038/nature14009 PMID: 25409151

43. Cassonnet P, Rolloy C, Neveu G, Vidalain PO, Chantier T, Pellet J, et al. Benchmarking a luciferase complementation assay for detecting protein complexes. Nat Methods. 2011; 8:990–2. https://doi.org/10.1038/nmeth.1773 PMID: 22127214

44. Munier S, Rolland T, Diot C, Jacob Y, Naffakh N. Exploration of Binary Virus–Host Interactions Using an Infectious Protein Complementation Assay. Mol Cell Proteomics. 2013; 12:2845. https://doi.org/10.1074/mcp.M113.028688 PMID: 23816991

45. Zheng L, Baumann U, Reymond JL. An efficient one-step site-directed and site-saturation mutagenesis protocol. Nucleic Acids Res. 2004; 32. https://doi.org/10.1093/nar/gnh110 PMID: 15304544

46. Krischuns T, Günl F, Henschel L, Binder M, Willemse J, Schloer S, et al. Phosphorylation of TRIM28 Enhances the Expression of IFN-β and Proinflammatory Cytokines During HPAIV Infection of Human Lung Epithelial Cells. Front Immunol. 2018; 9:2229. https://doi.org/10.3389/fimmu.2018.02229 PMID: 30323812

47. Matrosovich M, Matrosovich T, Garten W, Klenk HD. New low-viscosity overlay medium for viral plaque assays. Virol J. 2006; 3. https://doi.org/10.1186/1743-422X-3-63 PMID: 16945126

48. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012; 9:676–82. https://doi.org/10.1038/nmeth.2019 PMID: 22743772

49. O'Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016; 44:D733–45. https://doi.org/10.1093/nar/gkv1189 PMID: 26553804

50. Buchfink B, Reuter K, Drost HG. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat Methods. 2021; 18:366–8. https://doi.org/10.1038/s41592-021-01101-x PMID: 33828273

51. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29:644–52. https://doi.org/10.1038/nbt.1883 PMID: 21572440

52. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37:907–15. https://doi.org/10.1038/s41587-019-0201-4 PMID: 31375807
53. Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34:3094-100.
https://doi.org/10.1093/bioinformatics/bty191 PMID: 29750242

54. Kovaka S, Zimin AV, Pertea GM, Razaghi R, Salzberg SL, Pertea M. Transcriptome assembly from long-read RNA-seq alignments with StringTie2. Genome Biol. 2019; 20. https://doi.org/10.1186/s13059-019-1910-1 PMID: 31842956

55. Robert X, Gouet P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014;42. https://doi.org/10.1093/nar/gku316 PMID: 24753421

56. Reich S, Guilligay D, Cusack S. An in vitro fluorescence based study of initiation of RNA synthesis by influenza B polymerase. Nucleic Acids Res. 2017; 45:3353-68. https://doi.org/10.1093/nar/gkx043
PMID: 28126917

57. Shi M, Lin XD, Chen X, Tian JH, Chen LJ, Li K, et al. The evolutionary history of vertebrate RNA viruses. Nature. 2018; 556:197-202. https://doi.org/10.1038/s41586-018-0012-7 PMID: 29618816

58. Fan H, Walker AP, Carrique L, Keown JR, Serna Martin I, Karia D, et al. Structures of influenza A virus RNA polymerase offer insight into viral genome replication. Nature. 2019; 573:287-90. https://doi.org/10.1038/s41586-019-1530-7 PMID: 31485076

59. Chang S, Sun D, Liang H, Wang J, Li J, Guo L, et al. Cryo-EM structure of influenza virus RNA polymerase complex at 4.3 Å resolution. Mol Cell. 2015; 57:925-35. https://doi.org/10.1016/j.molcel.2014.12.031 PMID: 25620561

60. Biswas SK, Nayak DP. Mutational analysis of the conserved motifs of influenza A virus polymerase basic protein 1. J Virol. 1994; 68:1819-26. https://doi.org/10.1128/JVI.68.3.1819-1826.1994 PMID: 8107244

61. Ali S, Heathcote DA, Kroll SHB, Jogalekar AS, Scheiper B, Patel H, et al. The development of a selective cyclin-dependent kinase inhibitor that shows antitumor activity. Cancer Res. 2009; 69:6208-15.
https://doi.org/10.1158/0008-5472.CAN-09-0301 PMID: 19638587

62. Chen KY, Afonso EDS, Enouf V, Isel C, Naffakh N. Influenza virus polymerase subunits co-evolve to ensure proper levels of dimerization of the heterotrimer. PLoS Pathog. 2019; 15. https://doi.org/10.1371/journal.ppat.1008034 PMID: 31581279

63. Li C, Hatta M, Watanabe S, Neumann G, Kawaoka Y. Compatibility among polymerase subunit proteins is a restricting factor in reassortment between equine H7N7 and human H3N2 influenza viruses. J Virol. 2008; 82:11880-8. https://doi.org/10.1128/JVI.01445-08 PMID: 18815312

64. Nilsson BE, Velthuis AJWt, Fodor E. Role of the PB2 627 Domain in Influenza A Virus Polymerase Function. J Virol. 2017; 91:2467-83. https://doi.org/10.1128/JVI.02467-16 PMID: 28122973

65. Rippe K, Papantonis A. Functional organization of RNA polymerase II in nuclear subcompartments. Curr Opin Cell Biol. 2022; 74:88-96. https://doi.org/10.1016/j.ceb.2022.01.005 PMID: 35217398

66. Wille M, Holmes EC. The Ecology and Evolution of Influenza Viruses. Cold Spring Harb Perspect Med. 2020; 10:1-19. https://doi.org/10.1101/CSHPERSPECT.A038489 PMID: 31871237

67. Doamekpor SK, Sanchez AM, Schwer B, Shuman S, Lima CD. How an mRNA capping enzyme reads distinct RNA polymerase II and Spt5 CTD phosphorylation codes. Genes Dev. 2014; 28:1323-36.
https://doi.org/10.1101/gad.242768.114 PMID: 24939935

68. Bortz E, Westera L, Maamary J, Steel J, Albrecht RA, Manicassamy B, et al. Host- and strain-specific regulation of influenza virus polymerase activity by interacting cellular proteins. MBio. 2011; 2:1-10.
https://doi.org/10.1128/mBio.00151-11 PMID: 21846828

69. Martinez-Rucobo FW, Kohler R, van de Waterbeemd M, Heck AJR, Hemann M, Herzog F, et al. Molec- ular Basis of Transcription-Coupled Pre-mRNA Capping. Mol Cell. 2015; 58:1079-89. https://doi.org/ 10.1016/j.molcel.2015.04.004 PMID: 25959396

70. Fianu I, Chen Y, Dienemann C, Dybkov O, Linden A, Urlaub H, et al. Structural basis of Integrator-medi- ated transcription regulation. Science. 2021; 374:883-7. https://doi.org/10.1126/science.abk0154
PMID: 34762484

71. Bradel-Tretheway BG, Mattiacio JL, Krasnoselsky A, Stevenson C, Purdy D, Dewhurst S, et al. Com- prehensive proteomic analysis of influenza virus polymerase complex reveals a novel association with mitochondrial proteins and RNA polymerase accessory factors. J Virol. 2011; 85:8569-81. https://doi. org/10.1128/JVI.00496-11 PMID: 21715506

72. Pflug A, Gaudon S, Resa-Infante P, Lethier M, Reich S, Schulze WM, et al. Capped RNA primer binding to influenza polymerase and implications for the mechanism of cap-binding inhibitors. Nucleic Acids Res. 2018; 46:956-71. https://doi.org/10.1093/nar/gkx1210 PMID: 29202182