SNP rs72552343 in IFNAR1 reduces risk of developing active TB

To understand the function of type I IFN signaling in TB susceptibility in humans, we determined the association of IFNAR1 genetic polymorphisms and the risk of developing active TB disease in Chinese populations. Human IFNAR1 is located on chromosome 21 adjacent to IL10RB and IFNGR2 (Supplementary Fig. ^1a). In this study, six single-nucleotide polymorphisms (SNP) (Supplementary Table ¹) in the gene encoding IFNAR1 were examined for their association with clinical TB. The donor selection criteria and SNP genotyping using the Sequenom MassARRAY platform have been described in our previous study²². Active TB and HC patients were recruited at Shenzhen, a south-eastern city in Guangdong Province. The cohort includes 1533 active TB patients and 1445 healthy controls. Among the 1533 active TB patients, 1432 were diagnosed with pulmonary TB (PTB), 101 with extra-pulmonary TB (ETB). We observed that individuals carrying the SNP (rs72552343) with the deletion of nucleotides TCC in IFNAR1 (TCC/Del), showed significantly decreased risk of developing active TB (multiplicative model: p=0.0002; dominant model: p=0.0003) (Table 1). Further sequencing analysis confirmed that SNP rs72552343 is associated with the in-frame deletion of nucleotides TCC (Supplementary Fig. ^1b).

To validate this genetic association, the distribution of the SNP rs72552343 genotype was further determined in a second TB cohort, which was established independently in Chengdu city in Sichuan province located in south-western China²³. This validation cohort consists of 832 active TB cases (367 pulmonary TB and 465 extra-pulmonary TB) and 1084 controls. The significant association of IFNAR1 SNP rs72552343 (TCC/Del) and reduced risk of TB was replicated in both this validation cohort (p=0.004, multiplicative and dominant models) and the combined discovery and validation cohorts (2365 cases and 2529 controls, p<0.0001, multiplicative and dominant models, Table 1).

SNP rs72552343 in IFNAR1 increases hepatitis B susceptibility

Since type I IFNs are known to mediate anti-viral immunity, we next examined whether SNP rs72552343 influenced disease pathogenesis in viral infection in humans. To this end, the rs72552343 genotype distribution was determined in hepatitis B virus surface antigen (HBsAg)-positive hepatitis B (HepB) patients and HC recruited for a previously published study²⁴. In contrast to mycobacterial infection, IFNAR1 rs72552343 TCC/Del genotype was observed to be associated with significantly increased susceptibility to viral hepatitis (multiplicative model: p=0.01; dominant model: p=0.02, Table 2). Therefore, IFNAR1 signaling can differentially regulate resistance to intracellular bacterial and viral infection in humans.

TCC deletion in IFNAR1 is associated with reduced TB pathology

We next analyzed the association of rs72552343 with the clinical manifestations of TB diseases, sputum smear and mycobacterial culture positivity and lung cavity formation in the combined TB cohort. We found that while the rs72552343 TCC deletion was associated significantly with decreased risk of both pulmonary (PTB vs. HC, OR: 0.47; p<0.0001; dominant model) and extra-pulmonary TB disease (ETB vs. HC, OR: 0.40; p=0.005; dominant model), it did not preferentially influence the development of the either type of the disease (Supplementary Table ³). Similarly, the host protective effect of TCC deletion on pulmonary TB was present in TB patients who were mycobacterial culture and sputum smear positive or negative (Supplementary Tables ⁴ and ⁵). In addition, analysis of high-resolution computed tomography (HRCT) findings revealed that pulmonary TB patients carrying rs72552343 TCC/Del exhibited reduced risk of developing pulmonary cavities compared to their TCC/TCC counterparts (cavity− vs. cavity+, p=0.02; multiplicative and dominant models, Supplementary Table ⁶).

To investigate further the association between SNP rs72552343 and pulmonary pathology in TB patients, HRCT images of pulmonary TB patients enrolled in Shenzhen study were further assessed and scored based on the patterns, profusion and distribution of pulmonary abnormalities including cavity formation²⁵. This analysis revealed that pulmonary TB patients carrying TCC/Del genotype showed significantly lower HRCT scores (less severe TB disease) than those with the common allele (Fig. 2a). When radiographic data collected before and 2 years after treatment were compared in a subset of the TB patients, we found that IFNAR1 TCC deletion was also associated with a more favorable disease outcome following a standard 2HRZE/4HR treatment regimen (Fig. 2b). Altogether, these data establish an association between type I IFN receptor signaling and increased tissue pathology in active pulmonary TB.

Open in a separate window

Fig. 2

TCC deletion in IFNAR1 is associated with decreased risk of developing severe pulmonary pathology. a The HRCT scores of pulmonary TB patients carrying rs72552343 TCC/TCC (n=421) or TCC/Del (n=32) genotypes before the initiation of anti-TB chemotherapy. Circles denote individual patients and bars represent group means. b Paired analysis of HRCT scores before and 2 years after completion of anti-TB chemotherapy in pulmonary TB patients carrying rs72552343 TCC/TCC (n=43) or TCC/Del (n=10). Each paired data set represents an individual patient. Statistical analyses were performed using Student’s t-test. *p<0.05 and **p<0.01

TCC deletion in IFNAR1 impairs type I IFN signaling in PBMCs

To determine whether IFNAR1 SNP rs72552343 plays a role in the cellular response to IFN stimulation, PBMCs from individuals carrying WT (TCC/TCC) or mutant (TCC/Del) IFNAR1 genotypes were stimulated with recombinant human IFN-β, IFN-γ, or both cytokines, and ISG expression analyzed using qRT-PCR. We observed that MX1, IFIT1, OAS1, STAT2, and STAT1, which were highly sensitive to IFN-β induction, were up-regulated to a lesser extent in cells carrying the TCC deletion than in their WT counterparts following IFN-β stimulation (Fig. 3). In contrast, the upregulation of CXCL10 and IRF1 driven predominately by IFN-γ at this time-point was not affected by the IFNAR1 mutation. In addition, there was no significant difference in the level of IFNAR1 or IFNAR2 between genotypes. These findings collectively suggested TCC deletion is associated with defective type I, but not type II, IFN signaling.

Open in a separate window

Fig. 3

TCC/Del genotype PBMCs have reduced response to type I but not type II IFN. PBMC with TCC/TCC (gray) or TCC/Del (red) genotype were stimulated with IFN-β or IFN-γ or both (at 100 U/ml) for 6h and gene expression analyzed using qRT-PCR. The data shown are the median and range (n=6/genotype) of mRNA expression relative to GAPDH. Box plots represent the median (solid line within the box) and 25th and 75th percentile. The whiskers indicate minimal and maximal values

As IFNAR1 signaling can regulate type I IFN production^26,27, we examined whether the IFNAR1 mutation has an impact on the production of type I IFNs. We first compared the levels of IFN-α and IFN-β in the plasma of HC and TB patients carrying TCC/TCC or TCC/Del genotypes. Although IFN-α was not measureable in plasma, a low level of circulating IFN-β was detected in all TB cases irrespective of their genotype (Supplementary Fig. ^2a). Similarly, a comparable induction of IFN-α and IFN-β was observed in the culture supernatants of M. tuberculosis-infected macrophages carrying TCC/TCC or TCC/Del IFNAR1 (Supplementary Fig. ^2b), suggesting that the TCC/Del in IFNAR1 does not appear to impact type I IFN production in TB patients and M. tuberculosis-infected macrophages.

TCC deletion in IFNAR1 decreases type I IFN signal transduction

To establish definitively that the identified IFNAR1 deletion affects type I IFN signaling, mammalian expression vectors encoding human WT IFNAR1 and IFNAR1Del were generated for transfection studies. Together with an IFNAR2 encoding vector and an ISRE-driven luciferase reporter plasmid, the WT or mutant IFNAR1 vectors were transfected into Ifnar1^−/− mouse embryonic fibroblasts (MEFs). Ifnar1^−/− MEFs were used to ensure that the activation of the IFN signaling pathway by human IFNs depends exclusively on the transfected IFNAR1. Because both IFN-β and IFN-α can signal via IFNAR1¹, we also included IFN-α in this series of experiments to determine whether the signaling defect resulting from the TCC deletion is restricted to IFN-β stimulation. Upon stimulation with human IFN-α or IFN-β, WT IFNAR1-transfected cells exhibited significantly stronger luciferase activity compared to their IFNAR1Del-transfected counterparts, whereas IFN-γ induced minimal luciferase activity, similar to levels seen in the untransfected cells (Fig. 4a). Similarly, WT IFNAR1-transfected MEFs expressed significantly higher levels of IFN-inducible Isg15, Ifit1, and Oas1 than the IFNAR1Del-transfected cells when measured at 3 and 6 hr after IFN stimulation (Fig. 4b). The reduced IFN-inducible response in IFNAR1Del-transfected MEFs correlated with similarly impaired STAT1 phosphorylation following IFN-β stimulation (Fig. 4c and d).

Open in a separate window

Fig. 4

TCC deletion in IFNAR1 impairs type I IFN signal transduction. Ifnar1^−/− MEFs were transfected with vectors expressing WT (black) or mutant (red) IFNAR1 together with IFNAR2 and ISRE-driven luciferase reporter plasmids. a Luciferase activity in WT or TCC-deleted IFNAR1-transfected MEFs determined 18h after stimulation with IFN-α, IFN-β, or IFN-α. The data shown are mean±S.D. (triplicate cultures) relative intensity units. White bars: untransfected Ifnar1^−/− MEFs. b mRNA expression of ISGs in WT (black) or TCC-deleted IFNAR1 (red) transfected MEFs at 3 and 6h after stimulation with IFN-α or IFN-β. The data shown are mean±S.D. (triplicate cultures) fold increase over unstimulated controls. c Representative Western blot and d associated densitometric analysis of the levels of phosphorylated and total STAT1 proteins in WT or TCC-deleted cells at the indicated time points (minutes) following IFN-β stimulation. The protein levels were assessed by densitometry and the relative levels of phosphorylated STAT1 are expressed as fold induction with respect to total STAT1 levels. Data shown are representative of 2–3 independent experiments

To confirm the above findings in a human system, we knocked down endogenous IFNAR1 in HEK293 cells using CRISPR/Cas9 gene-editing technology and then transfected the cells with human WT IFNAR1 or IFNAR1Del plasmids together with an ISRE-driven luciferase reporter plasmid. Similar to the studies performed using MEFs, following type I IFN stimulation, HEK293 cells transfected with WT IFNAR1 exhibited significantly stronger luciferase activity compared to IFNAR1Del-transfected counterparts (Fig. 5a). WT IFNAR1-transfected cells also expressed significantly higher levels of IFN-inducible ISG15, IFIT1, and OAS1 compared to the IFNAR1Del-transfected counterparts when measured at 6h after IFN stimulation (Fig. 5b). Altogether, these biochemical studies performed using both murine and human cells provide direct evidence demonstrating that the TCC deletion and the resulting deletion of proline 335 (Pro335) impairs IFNAR1-dependent type I IFN signaling.

Open in a separate window

Fig. 5

TCC deletion in IFNAR1 decreases type I IFN signaling in human cells. HEK 293T cells were transfected with Cas9-expressing vector with or without endogenous IFNAR1 targeting gRNA. The endogenous IFNAR1-silenced (gRNA+) or control (gRNA−) cells were then transfected with WT or mutant IFNAR1 expression vectors together with an ISRE-driven luciferase reporter plasmid. a The luciferase activity was analyzed 18h after IFN-α or IFN-β stimulation. The data shown are mean±S.D. relative luciferase intensity units of triplicate cultures. b The expression of ISGs in WT (black) or TCC-deleted IFNAR1 (red) transfected IFNAR1 knockout HEK293 cells at 6h after stimulation with IFN-β determined using qRT-PCR. The data shown are mean fold increase±S.D. (triplicate cultures) over unstimulated controls. Data shown are representative of 2–3 independent experiments

Ethics Statement

Blood samples were collected and analyzed after written informed consent was obtained from participants and with approval of Shenzhen Third People’s Hospital Ethics Committee (Reference No. 2012–003) and Clinical Trials and Biomedical Ethics Committee of West China Hospital, Sichuan University (Reference No. 198 (2014)).

SNP selection and genotyping

Genomic DNA was prepared from PBMC using QIAamp DNA Blood Minikit (Qiagen, Hilden, Germany). To select SNP candidates in IFNAR1 gene, we used a Position Weight Matrices (PWM)_SCAN algorithm to scan IFNAR1 sequence in Jaspar, UniPROBE, TRANSFAC and PITA databases. For the SNP within the putative binding sites, Sr value was calculated to reflect the change of binding scores between the putative binding site and putative binding protein, and then converted into a p-value based on the FastPval program, the SNP with Sr p-value<0.01 was selected for genotyping⁴⁷, in addition, non-synonymous substitutions in the exon region were also selected. In TB discovery study, three SNPs in the promoter (rs2834191 T>G, rs2843710 C>G, rs17875752 G>T), 1 SNP in exon (rs72552343 TCC>DEL) and 2 SNPs in the intron (rs1012334 A>T, rs1041868 G>A) of IFNAR1 (Supplementary Table ¹) were genotyped in health controls and TB patients using the MassARRAY system (Sequenom, San Diego, CA). The relative height (intensity) of the peaks and the signal-to-noise (SNR) ratio were analyzed using Caller software to call genotypes in real-time. After cluster analysis using Typer software, manual curation of spectra was performed to further validate the outcome. In TB validation and HBV cohort studies, the SNP rs72552343 was genotyped in control and case samples using TaqMan assays (ABI, Carlsbad, CA) by CT Bioscience (Suzhou, China).

HRCT and radiological scoring

HRCT were performed at 10mm section interval (120kV, 50–450 mAs, 1mm slice thickness, 1.5s scanning time) with a window level between 2550 and 40 Hounsfield Units (HU) and window width between 300 and 1600 HU using the Toshiba Aquilion 64 CT Scanner (Toshiba, Tokyo, Japan). HRCT scans were examined by two independent chest radiologists and final conclusions on the findings were reached by consensus. Radio-pathological changes were quantified using a scoring system developed by Ors et al. with the following parameters: (1) micronodule; (2) nodule; (3) consolidation; (4) ground glass opacity; (5) cavity; (6) bronchial lesion; (7) parenchymal bands²⁵. The total score for two lungs is 168.

Antibodies, plasmids, cell lines, and cytokines

The following rabbit primary antibodies (Cell Signaling Technology), anti-STAT1 (clone D1K9Y, 1:1000), anti-phospho-STAT1 (Tyr701, clone D4A7, 1:1000), and β-tubulin (clone 9F3, 1:1000), together with appropriate horseradish peroxidase (HRP)-conjugated secondary antibodies (Promega, 1:5000) were used for Western blotting.

Human IFNAR1 (pEF-BOS-hIFNAR1) and IFNAR2 (pEF-BOS-hIFNAR2) expressing vectors were constructed by cloning hIFNAR1 and hIFNAR2 open reading frame (ORF) into pEF-BOS expression vector respectively. To generate the TCC deletion vector (pEF-BOS-hIFNAR1-DEL), the pEF-BOS-hIFNAR1 vector was subjected to site-directed mutagenesis using Q5 Site-Directed Mutagenesis Kit (New England Biolabs). Interferon Stimulated Response Element (ISRE) firefly luciferase reporter and thymidine kinase-driven renilla luciferase control (pRL-TK) vectors have been described previously⁴⁸. Plasmids pRRL-Cas9-HA-NLS (nuclear localization signal) and pRRL-U6-sgRNA⁴⁹ were kindly provided by Dr Paul Zhou (Shanghai Institute of Pasteur, China). Vector pRRL-Cas9-HA-NLS was derived from the third generation of lentiviral transfer vector and contained human codon optimized fusion gene encoding S. pyogenes Cas9, HA tag, and SV40 NLS.

Ifnar1^−/− MEFs⁵⁰ and HEK 293T cells (ATCC CRL-3216™) were maintained in complete RPMI 1640 medium containing 10% FCS, 100 U/ml penicillin, 100μg/ml streptomycin, 2μM glutamate, and 10mM HEPES. The recombinant human IFN-α-2a, IFN-β and IFN-γ were purchased from PBL Assay Science.

Cell isolation and in vitro cytokine stimulation of PBMCs

PBMCs were isolated by density gradient centrifugation over Ficoll-Hypaque. Pleural fluid mononuclear cells (PFMCs) were obtained by centrifuging 50–200ml of pleural fluid from TB patients with pleural effusions at 300×g for 5min. For ISG expression analysis, PBMCs isolated from the individuals carrying rs72552343 TCC/TCC or TCC/DEL genotypes were cultured in a 96-well plate (2×10⁵ cells/well) with medium only, human IFN-β (100 IU/ml), human IFN-γ (100 IU/ml) or both cytokines for 6h and total RNA isolated for gene expression analysis using qRT-PCR.

Human macrophage culture and M. tuberculosis infection

PBMCs were obtained from healthy donors with either TCC or TCC/Del genotype. Monocytes positively selected from PBMCs using anti-CD14 magnetic beads (Miltenyi Biotec) were cultured in complete RPMI 1640 medium for 7 days. Recombinant human M-CSF (10 ng/ml, PeproTech) was added on days 0, 2, and 4. For infection with mycobacteria, differentiated macrophages seeded in 24-well plates (2×10⁵ cells per well in 0.5ml antibiotic-free medium) were exposed to M. tuberculosis H37Rv at a multiplicity of infection of 10. After 6h of incubation at 37°C, the cells were washed three times with warm antibiotic-free medium to remove extracellular bacteria. The supernatants of macrophage cultures were collected 24h post-infection and the levels of type I IFNs determined using high sensitivity human IFN-β and IFN-α all subtype ELISA kits (PBL Assay Science).

Transfection in MEFs

Ifnar1^−/− MEFs, seeded in 24-well plates (4×10⁴/ml/well) 12h earlier, were co-transfected with plasmid DNA pEF-BOS-hIFNAR1-WT (0.1μg) or pEF-BOS-hIFNAR1-DEL (0.1μg) vector, together with pEF-BOS-hIFNAR2 (0.1μg), ISRE reporter (0.2μg), and pRL-TK control vector (0.1μg) using Lipofectamine 2000 reagent (Life Technology). Twenty-four hours after transfection, the cells were stimulated with 100 IU/ml human IFN-α, IFN-β, or IFN-γ. In some experiments, the transfected cells were lysed in 100μl passive lysis buffer (Promega) after 18h IFN stimulation and the lysates assayed for both the firefly and Renilla luciferase activities using the dual-luciferase reporter assay system (Promega). Firefly luciferase was normalized against Renilla luciferase activity (F/R ratio) to account for variation in cell numbers and viability. Alternatively, transfected and IFN stimulated Ifnar1^−/− MEFs were lysed at indicated time points to extract total RNA and protein for analyzing ISG expression by qRT-PCR and STAT1 phosphorylation by Western blot, respectively.

CRISPR/Cas9 gene editing and transfection

To determine the role of IFNAR1 TCC/Del in IFN signaling in human cells, we developed a transfection system where the constitutively expressed endogenous IFNAR1 in the HEK 293T cells was knocked down using CRISPR/Cas9 gene-editing technology. To ensure CRISPR/Cas9 specifically edited endogenous IFNAR1 gene, but not open reading frame-containing, transfected exogenous IFNAR1 constructs, we designed a gRNA (5′-GTGGGAGGGTCACTTGAACC-3′) using CRISPR Design Tool developed by the laboratory of Feng Zhang, MIT (http://crispr.mit.edu/) to target intron 8 of hINFAR1 and inserted the sequence into pRRL-U6-sgRNA vector (pRRL-U6-sgRNA-intron 8). HEK 293T cells were transfected with Cas9-expressing vector pRRL-Cas9-HA-NLS and pRRL-U6-sgRNA-intron 8 using Lipofectamine 2000 reagent. Some of the endogenous IFNAR1-silenced cells were co-transfected with pEF-BOS-hIFNAR1-WT or pEF-BOS-hIFNAR1-DEL vector, with ISRE reporter and pRL-TK control vectors. After 24h, the cells were stimulated with 100 IU/ml human IFN-α, IFN-β, or IFN-γ and luciferase activities determined at 18h post-stimulation as described above.

Western blotting

Ifnar1^−/− MEFs co-transfected with pEF-BOS-hIFNAR2 and pEF-BOS-hIFNAR1-WT or pEF-BOS-hIFNAR1-DEL vector were lysed at indicated time points after type I IFN stimulation (100 IU/ml) using lysis buffer (10mM Tris, pH 7.5, 150mM NaCl, 1% Triton X-100, 1μM phenylmethylsulfonyl fluoride, 0.2μM sodium orthovanadate, 0.5% Nonidet P-40) supplemented with a cocktail of protease inhibitors (Sigma-Aldrich). The lysates were separated by SDS-PAGE on a 4–15% gradient polyacrylamide gel (Bio-Rad) and then transferred onto a polyvinylidene difluoride membrane. The membranes were sequentially probed with the respective primary antibodies followed by appropriate HRP-conjugated secondary antibodies and then visualized using a CCD camera system (Bio-Rad). The densities of individual protein bands were quantified using the ImageJ software (http://rsb.info.nih.gov/ij/index.html) and the ratio of p-STAT1 to total STAT1 determined. Uncropped Western blotting images are provided in Supplementary Fig. ³.

qRT-PCR

RNA was extracted with Trisure according to the manufacturer’s instructions (Bioline, Australia). Total RNA (2μg) was reverse transcribed using the Tetro cDNA synthesis Kit with random primers according to the manufacturer’s instructions (Bioline). All qRT-PCR were performed using SYBR NoROX master mix (Bioline) on a Roche LightCycler480. Forward and reverse qRT-PCR primers are listed in Supplementary Table ⁷. Data were calculated by the 2^−ΔΔCT method using 18S or GAPDH as the house-keeping genes.

Construction of hIFNAR1-ECD clone and mutagenesis

The coding sequence for the extracellular domain (ECD) of hIFNAR1 was amplified from a clone of full-length hIFNAR1⁵¹ and represents amino acids 28–424 of hIFNAR1 (GenBank accession number CAA42992.1). Amplification of the hIFNAR1-ECD was directed by specific forward (5′-GCCGGATCCAAAAATCTAAAATCTCCTCAA-3′) and reverse (5′-GGCGAATTCTTAT-TTAGAGGTATTTCCTGG-3′) primers; after amplification and digestion of the PCR product with BamHI (5′-end) and EcoRI (3′-end), the fragment was cloned into a modified pFastBac-1 vector permitting expression of recombinant proteins with an N-terminal 6xHis tag as described previously⁵². Upon sequence verification of hIFNAR1-ECD, mutagenesis primers (fwd 5′- GCTTTCCTACTTCCTGTCTTTAACATTAGATC-3′; rev 5′-GATCTAATGTTAAAGA-CAGGAAGTAGGAAAGC-3′) were used to delete P335 from the construct to generate hIFNAR1-ECD-Del) using Pfu DNA polymerase following the manufacturer’s specifications (Promega). This clone was also sequence verified before use.

Recombinant protein expression and purification

For expression in insect cells, the clones of hIFNAR1-ECD and hIFNAR1-ECD-Del in modified pFastBac-1 (above) were used to generate recombinant bacmids and baculoviral stocks using transfection of Spodoptera frugiperda (Sf)-9 incest cells following the protocols in the Bac-to-Bac Expression System manual (LifeTechnologies). Recombinant forms of hIFNAR1-ECD and hIFNAR1-ECD-Del were expressed as soluble proteins and secreted into the culture supernatant of High Five™ (BTI-TN-5B1–4 from Trichoplusia ni purchased from ThermoFisher Scientific, Cat. No. B85502). Recombinant proteins were purified by affinity and size exclusion chromatography following protocols established in our laboratory for protein purification from insect cell culture⁵². All purified proteins were assessed for purity by reducing SDS-PAGE and protein concentration (mg/ml) measured by absorbance at 280nM/extinction co-efficient (per M/cm) (using a molecular weight of 49,000 for hIFNAR1-ECD).

Surface plasmon resonance

SPR was performed to measure and compare the affinity of ReBif (hIFN-β) to recombinant forms of hIFNAR1-ECD and hIFNAR1-ECD-Del. All experiments were carried out in the absence of IFNAR2. All SPR experiments were performed on a ProteOn XPR36 (Bio-rad Labs) using a HTG chip for His-tagged proteins and TBS as the running buffer. Both ligands (hIFNAR1-ECD and hIFNAR1-ECD-Del) were immobilized to the nickel activated chip via their C-terminal His-tags after dilution to 25µg/ml in TBS. ReBif, (hIFN-β) previously dialyzed into TBS, was also diluted in TBS to various concentrations. All data were referenced according to the manufacturer’s instructions (Bio-rad) and analyzed using the Langmuir binding model. Data were considered for inclusion in the analysis only if the Chi₂ value (the measure of error between measured and fitted values) was less than 10% of the R_max as per the manufacturer’s instructions (Bio-rad). K_a (1/Ms), K_d (1/s), and K_D (nM) were calculated by the ProteOn Manager software (Bio-rad) and are represented as mean±S.D. from at least three independent experiments.

Statistical analysis

For genetic association studies, the Hardy–Weinberg Equilibrium for IFNAR1 SNP distribution was analyzed in patients and healthy controls. The Pearson χ² test was used to compare the genotypic and allelic frequencies of SNPs between cases and controls. Unconditional logistic regression adjusted by gender and age was performed to calculate the odd ratios (ORs), 95% confidence intervals (CIs), and corresponding p values using four alternative models (multiplicative, additive, dominant and recessive). For statistical analysis for gene expression and biochemical analysis, one-way analysis of variance (ANOVA) with Newman–Keuls multiple comparison test was used to compare the differences among multiple groups. Student’s t-test was used to compare the difference between two experimental groups. All statistical tests were performed using Prism 6.0 (GraphPad). p<0.05 is considered statistically significant.

This work was supported by a National Health and Medical Research Council (NHMRC) of Australia Grants (APP1051742 to C.G.F., APP1070782 to N.A.d.W., APP1049686 and APP1027020 to P.J.H., and APP1043225 to W.J.B.). Research at the Hudson Institute is supported by the Victorian Government’s Operational Infrastructure Support Program. We acknowledge the support from Natural Science Foundation of China (No. 81471913, 81501714, 81525016); Thirteenth-Fifth Mega-Scientific Projects on “prevention and treatment of AIDS, viral hepatitis and other infectious diseases” (No. 2017ZX10103004); Shenzhen Scientific and Technological Foundation (No. JCYJ20170412151620658, JCYJ201703070950030518, JSGG20160427104724699).