NF-κB1 Haploinsufficiency Causing Immunodeficiency and EBV-Driven Lymphoproliferation

Purpose NF-κB signaling is critically important for regulation of both innate and adaptive immune responses. While activation of NF-κB has been implicated in malignancies such as leukemia and lymphoma, loss-of-function mutations affecting different NF-κB pathway components have been shown to cause primary immunodeficiency disorders. Recently, haploinsufficiency of NF-κB1 has been described in three families with common variable immunodeficiency (CVID). Methods and Results We studied a patient with recurrent respiratory infections and bacterial parapharyngeal abscess. Immunological investigations revealed normal total B- cell numbers, but hypogammaglobulinemia, decreased frequencies of class-switched B cells and impaired T-cell proliferation. Targeted next-generation sequencing using a custom-designed panel comprising all known PID genes (IUIS 2014 classification) and novel candidate genes identified a novel heterozygous frameshift mutation in the NFKB1 gene leading to a premature stop codon (c.491delG; p.G165A*31). We could show that the mutation leads to reduced phosphorylation of p105 upon stimulation, resulting in decreased protein levels of p50. The further disease course was mainly characterized by two episodes of severe EBV-associated lymphoproliferative disease responsive to rituximab treatment. Due to disease severity, the patient is considered for allogeneic hematopoietic stem cell transplantation. Interestingly, the father carries the same heterozygous NFKB1 mutation and also shows decreased frequencies of memory B cells but has a much milder clinical phenotype, in line with a considerable phenotypic disease heterogeneity. Conclusions Deficiency of NF-κB1 leads to immunodeficiency with a wider phenotypic spectrum of disease manifestation than previously appreciated, including EBV lymphoproliferative diseases as a hitherto unrecognized feature of the disease. Electronic supplementary material The online version of this article (doi:10.1007/s10875-016-0306-1) contains supplementary material, which is available to authorized users.


Introduction
NF-κB signaling is a complex, conserved pathway critical for regulation of both innate and adaptive immune reactions in response to various stimuli. Central players in this pathway in mammals comprise five transcription factors, including REL-A/p65, REL-B, c-REL, NF-κB1/p50, and NF-κB2/p52 [1]. These transcription factors share a Rel homology domain (RHD) in their N-terminal DNA-binding domain responsible for formation of homo-and heterodimers which can bind to DNA to regulate target gene expression [1].
NF-κB signaling cascades are differentiated into two major subsets termed Bcanonical^and Bnon-canonical^pathways. Canonical NF-κB signaling may be activated upon contact with numerous inflammatory stimuli including pattern recognition receptors and cytokine receptors and signals through a complex consisting of IKKα/IKKβ/IKKγ and the formation of p50/p65 heterodimers (reviewed in [2]). By contrast, noncanonical NF-κB signaling is activated through tumor necrosis family receptors upon binding of BAFF, lymphotoxin β, or CD40 ligand, and is predominantly known to regulate immune cell differentiation and also formation of secondary lymphoid organs (reviewed in [3]). Notably, recent studies have demonstrated considerable cross talk between both arms of NF-κB signaling; thus, these signaling cascades are highly interconnected [4].
In view of the central role of NF-κB signaling for multiple processes in both innate and adaptive immunity, it is not surprising that aberrant NF-κB signaling may be implicated in human disease pathogenesis. Thus, while excessive activation of NF-κB has been shown to play a role in hematological malignancies [5], loss-of-function mutations in genes of the canonical and non-canonical NF-κB pathway have been associated with distinct types of primary immunodeficiency disorders including deficiencies of NEMO [6] or IKKβ [7] and BAFF receptor [8], CD40/CD40 ligand [9,10], NIK [11], AID [12], or autosomal dominant mutations affecting NFκB2 [13,14], respectively.
Very recently, haploinsufficiency of NF-κB1 has been identified as a novel genetic etiology of a subtype of common variable immunodeficiency (OMIM616576) [15]. The gene encodes for the p105 protein which is processed to the Bactive^p50 subunit implicated in the canonical NF-κB pathway. As p50 lacks a C-terminal transcription activation domain, p50 can only activate transcription by formation of heterodimers, for example with REL-A, while p50 homodimers act as transcriptional repressors [16].
Here, we identify a novel mutation of NFKB1 in a patient with combined immunodeficiency with impaired B and T cell functions and presentation with severe Epstein-Barr virus (EBV)-associated lymphoproliferation as a hitherto unrecognized clinical disease manifestation.

Patients
All patient material was obtained in accordance with the Declaration of Helsinki. The study was approved by the ethics committee of the Medical University of Vienna.

DNA Isolation and Preparation
Genomic DNA (gDNA) was isolated from EDTA blood using an adapted protocol of the Wizard ® Genomic DNA Purification Kit (Promega). gDNA isolation from buccal swabs was performed using the QIAamp ® DNA Mini Kit (Qiagen), following the spin protocol of the QIAamp ® DNA Mini and Blood Mini Handbook.
For library preparation, gDNA was diluted and then measured on a Qubit 2.0 Fluorometer (Invitrogen/Life Technologies) for a total concentration of 200 ng.

Targeted Exome Sequencing
The patient sample was screened for disease-causing variants by a custom-designed targeted enrichment approach (HaloPlex ™ /Agilent Technologies) followed by nextgeneration sequencing on a HiSeq3000 (Illumina) platform as described previously [17]. In brief, enrichment of the targeted plus 25-bp flanking region was accomplished using the HaloPlex Target Enrichment System (Agilent Technologies Inc., 2013), based on a molecular inversion probe strategy. Library preparation was performed according to the manufacturer's instruction. In brief, 200 ng of gDNA was digested by eight pairs of restriction enzymes, followed by bar code indexing and hybridization to custom-designed capture probes for 16 h at 54°C. Thereafter, the circularized biotinylated target-probe complexes were extracted using magnetic streptavidin beads. The final steps included nick ligation, PCR library amplification, and AmPure XP bead (Beckman Coulter, Inc.) purification prior to qualitative and quantitative assessment of the DNA library using a 2100 Bioanalyzer instrument (Agilent). Next-generation sequencing was performed in a 150-bp paired-end mode using a HiSeq3000 (Illumina) platform.

Data Analysis
The gross data analysis pipeline included adapter trimming of Illumina sequences (Trimmomatic), Burrows-Wheeler Aligner (BWA) for sequence alignment to the human genome 19 (hg19), Indel Realignment on both sequence aliquot and sample level via Genome Analysis Toolkit (GATK; Broad Institute), Base Quality Score Recalibration (GATK), Haplotype Calling (GATK), and Variant Annotation (SnpEFF, GATK). Thereafter, variant filtering included the criteria of being rare (MAF ≤ 0.01), non-synonymous, and within the coding region of the targeted genes. In addition to published data, we assessed the potential relevance of variants by recurrence within ExAC browser (Exome Aggregation Consortium Cambridge) and in our internal dataset comprising of more than 300 sequenced individuals to date. Of note, variants with a VQSLOD score (the log odds of being a true variant versus being false) below 99.9 % of the truth set of a trained Gaussian mixture model can be considered as false positives and are thus not shown herein.

Coverage
The GATK CallableLoci tool was executed in order to assess the proportion of callable bases, as determined by sequencing depth and mapping quality per interrogated position. Hence, targeted genomic regions were assigned different quality categories (pass, no coverage, low coverage, excessive coverage, poor mapping quality) and summarized in a BED file. According to this analysis, 99.76 % of enriched exonic bases were considered callable with a minimum read depth of 4 and minimum mapping quality score of 10.

Variant Validations
Sanger validation was performed according to standard capillary sequencing protocol using a 3130xl Genetic Analyzer instrument.

Flow Cytometry-Based Immunophenotyping and T Cell Proliferation Assays
These analyses were carried out as described previously [18].

Western Blot Analysis
Following 3 h of starvation in IMDM without FCS, EBVtransformed B lymphocytes derived from the affected individual and the father, both harboring the heterozygous NFKB1 c.491delG mutation, as well as from a healthy donor individual, were stimulated with 50 ng/ml PMA plus 1 μg/ml ionomycin (both from Sigma Aldrich) for 30 and 60 min, as indicated. Cells were lysed on ice (lysis buffer 20 mM Tris-HCl, pH 7.5, 150 mM NaCl, 2 mM EDTA, 1 μM Na 3 VO 4 , 50 mM NaF, 1 % Triton X-100; Protease Inhibitor Cocktail) and total protein concentrations determined via Bradford (Biorad). Fifty micrograms per sample were loaded and protein size fractionated using a 12 % polyacrylamide gel. Proteins were transferred onto a PVDF (GE Healthcare) membrane (0.45-μm pore size) overnight at 4°C (200 mA) in accordance with standard Western blotting protocols. P105 and p50 proteins were detected using a rabbit antibody raised against the NF-κB1 N-terminus (no. 3035, Cell Signaling, New England Biolabs); phosphorylated p105 (Ser933) was detected using a monoclonal rabbit antibody (no. 4806, Cell Signaling). GAPDH served as loading control and was detected using a mouse monoclonal antibody (sc-365062; Santa Cruz). Anti-rabbit secondary antibody (Biorad) and antimouse secondary antibody (BD Biosciences) coupled to horseradish-peroxidase were used for signal detection. Immunoreactive bands were visualized with Amersham ECL Select Western Blotting Substrate (GE Healthcare).

Results
The index patient was born to non-consanguineous parents of Caucasian origin. Following an unremarkable first year of life, she started to have recurrent mainly upper respiratory tract infections in early childhood. At 5 years of age, tonsillectomy and adenotomy were performed. In the following years, the patient suffered from recurrent upper airway infections which were mostly self-limiting.
At the age of 15 years, she was referred to our hospital for the first time with fever and problems with swallowing and was found to have parapharyngeal abscess (Fig. 1a) and cervical lymphadenopathy. Blood results showed leukopenia (1.89 G/l), neutropenia (0.06 G/l), thrombocytopenia (minimum 119 G/l), and elevated C-reactive protein levels (CRP; 174 mg/l). No hepatosplenomegaly was apparent. The parapharyngeal abscess was surgically drained; histological evaluation revealed chronic abscess-forming inflammation, microbiologic culture infection with viridans group streptococci and actinomyces. Decreased leukocyte and thrombocyte counts prompted bone marrow aspiration, which revealed predominantly immature granulopoesis, compatible with bacterial infection, but no evidence of myelodysplasia or malignant disease. The patient improved with antibiotic therapy and absolute neutrophil counts normalized concomitantly. Immunological investigations following recovery revealed low absolute numbers of CD19 + B cells with reduced nonswitched and switched memory cells (IgD + /CD27 + ; IgD − / CD27 + ) and low immunoglobulin levels (Fig. 1b, c and Table 1). Absolute numbers of T cells were normal; however, proliferation assays revealed impaired T-cell function in response to tetanus toxoid, CD3, and PMA with normal response to PHA ( Table 1). The patient was started on intravenous and subsequently subcutaneous immunoglobulin substitution. Absolute neutrophil counts and thrombocytes were repeatedly decreased (thrombocytes with a minimum of 32 G/l) with spontaneous normalization; however, anti-granulocyte and anti-thrombocyte antibodies were negative. Few months later, she developed fever, neutropenia (0.1 G/l), elevated CRP (80 mg/l), esophagitis, cervical, axillary and supraclavicular lymphadenopathy, and mild hepatosplenomegaly. Axillary lymph node biopsy showed reactive changes and the patient improved with antibiotic therapy.
In the following years, the patient had three to four uncomplicated infections per year (sinusitis, pharyngitis, urinary tract infection). She complained of intermittent arthralgia; autoantibodies were repeatedly negative (rheumatoid factor (RF), anti-neutrophil cytoplasmatic antibody (ANCA), anti-cardiolipin antibodies, anti-nuclear antibody (ANA)). At the age of 18 years, splenomegaly (21 cm) with multiple hypodense changes, hepatomegaly (20 cm), and generalized lymphadenopathy occurred (not shown). Blood results showed reactivation of Epstein-Barr-virus (EBV) with 1.2 × 10E3 copies/ml. Within a short time, the patient deteriorated with high fever, elevated CRP (50 mg/l), and increasing hepatosplenomegaly and lymphadenopathy (liver 22 cm, spleen 26 cm; Fig. 1d). Further i nve stigations were not suspicious for macrophage-activating syndrome. With rising EBV load (maximum 3.1 × 10E3 copies/ml), EBV-associated lymphoproliferative disease (EBV-LPD) was suspected. The patient was started on corticosteroid treatment (1 mg/kg) and weekly rituximab (375 mg/m 2 ) was administered four times. Prompt clinical improvement was noted, EBV load turned negative, and CD19 + B cells were not detectable for 8 months (Fig. 2a). However with reemergence of peripheral blood B cells, a similar clinical picture with fever, lymphadenopathy, and hepatoslenomegaly reoccured and reactivation of EBV (maximal load 1.2 × 10E2 copies/ml) was detected. The patient was started on cortisone and four courses rituximab, which has led to remission of the disease again. In view of the relatively severe disease course, HLA typing of the patient and family has been initiated to prepare for allogeneic hematopoietic stem cell transplantation. On inquiry of the patient's family, the patient's father declared he had been suffering from frequent non-severe respiratory tract infections but not requiring regular medical attention.

Identification of Underlying Genetic Defect in NFKB1
Given the unclear origin of immunodeficiency, we assessed the patient using our custom-designed next-generation sequencing (NGS)-based panel which covers 419 genes, including all known primary immunodeficiency disease (PID)-causing genes according to the IUIS classification from 2014 [19] and several PID candidate genes in which defects are likely to be causative for monogenic primary immunodeficiencies. Our sequencing approach yielded very high on-target coverage (Fig. 2b). The analysis of this patient did not reveal any homozygous variant. Among few other heterozygous variants   Table 1), we could, however, identify a novel frameshift mutation in NFKB1 (ENST00000226574 c.491delG; p.G165A*31), which resides within the Nterminal RHD and results in a premature stop codon. While the mother and sister were wild type at this position, genetic analysis revealed the same NFKB1 mutation in the patient's father as in our index patient (Fig. 2c), however, associated with a mild clinical phenotype. While his absolute numbers of CD19 + B cells were mildly elevated (0.89 × 10E9/l), switched memory B cells (IgD + /CD27 + 49.9 %; IgD − /CD27 + 1.05 %) as well as IgG (including subclasses IgG1 and IgG3) and IgA were decreased. T cell proliferation assays revealed decreased response to T cell recall antigens such as tetanus toxoid, however not as pronounced as in the index patient, while T cell responses to stimulation with PHA, PMA, or anti-CD3 were normal. Despite a previous immunization, vaccination titers to diphtheria toxin were low, comparable with the index patient's results (Table 1). Western blot analysis revealed detectable but severely decreased levels of p50 in both the index patient (II-1) and her Values in brackets show reference ranges. Abnormal values are printed in italic WBC white blood cells, ALC absolute lymphocyte count, ANC absolute neutrophil count, PMA phorbol-12-myristate-13-acetate, PHA phytohemagglutinin, Ab antibody a Post four cycles of rituximab father (I-2). Upon stimulation with PMA and ionomycin, we observed that phosphorylation of p105 was markedly lower than in a healthy control in line with a signaling defected conferred by the NFKB1 mutation detected in both the index patient and her father (Fig. 2d). No truncated p50 protein was detected (Supplementary Figure 1).

Discussion
Identification of the genetic etiologies underlying primary immunodeficiencies has contributed considerably to our understanding of core signaling processes in the human immune system. In particular, the advent of NGS technologies has ΔG c.491delG  Table 1). The mutation leads to reduced phosphorylation of p105 upon stimulation in both index patient (II-1) and father (I-2), resulting in decreased protein levels of p50 (d) accelerated disease gene discovery for rare monogenetic disorders altogether [20]. While historically large pedigrees were necessary to identify novel causative gene defects and the majority of identified defects were inherited in an autosomal recessive manner (in part because studies in consanguineous pedigrees were used to identify underlying gene defects), more recently, an increasing number of autosomal dominant (AD) PIDs has been discovered, including such that work through AD gain-of-function mutations (e.g., activated PIK3δ syndrome [21][22][23][24]) or haploinsufficiency (e.g., CTLA4 haploinsufficiency syndrome [25,26]).
Recently, haploinsufficiency of the NF-κB subunit p50 has been described in three families [15]. Mutations were located within the RHD domain and lead to altered or absent p50 protein. Clinically, the patients mainly presented with recurrent respiratory tract infections, lymphadenopathy, and chronic-obstructive pulmonary disease, with marked disease severity. Laboratory findings included hypogammaglobulinemia and abnormal response to immunizations in some patients but normal numbers and function of T cells at least in some patients, in line with a diagnosis of common variable immunodeficiency [15].
Here, we describe a young adult patient with a novel frameshift mutation in NFKB1, also residing within the RHD. Clinically, our patient initially had a fairly mild disease course during childhood but then increasingly suffered from upper respiratory tract infections. The immunological investigations, prompted by diagnosis of a parapharyngeal abscess with uncommon microbes, demonstrated hypogammaglobulinemia with paucity of class-switched B cells and impaired generation of antibodies upon vaccination, but also impaired T-cell proliferative responses, suggestive of a picture of combined immunodeficiency. The functional impairment of T cells likely underlies the susceptibility to EBV-induced lymphoproliferation, which is a hitherto unrecognized feature of NF-κB1 haploinsufficiency. Several immunodeficiency disorders associated with EBV lymphoproliferation have been described, including deficiency of ITK, XIAP, STK4, SH2D1A, and CD27 [27][28][29][30]. The coding sequences of these genes were fully covered, and no mutations were detected in our targeted NGS panel (Supplementary Table 2). A defective function of T and NK cells is thought to lead to immune dysregulation and lymphoproliferation [30]. While absolute numbers of T cells were normal in our patient, T-cell proliferation assays showed markedly decreased results possibly contributing to EBV lymphoproliferation. Haploinsufficiency of NF-κB1 has been described so far as a predominant defect of B cells and has been assigned to the group of CVID. In view of the diminished function of T cells in our patient and the severe EBV lymphoproliferative disease, we suggest that haploinsufficiency of NF-κB1 may cause a phenotype of combined immunodeficiency at least in some patients.
Autosomal dominant inherited immunodeficiencies are often associated with incomplete disease penetrance and heterogeneity of clinical and laboratory manifestations. It seems that this also applies to haploinsufficiency of NF-κB1 as our patient's father carries the same NFKB1 mutation as the index patient, however associated with milder clinical and laboratory findings. In view of the disease heterogeneity, therapeutic implications cannot be uniform and need to be adapted to the presentation of the individual patient. In our index patient, disease severity and complications increased over time with two severe episodes of EBV-associated lymphoproliferation within one year. Hence, we are currently considering the patient for allogeneic hematopoietic stem cell transplantation. Comprehensive studies with larger numbers of NFKB1 mutant patients are urgently needed to exhibit the full spectrum of clinical disease manifestations of NF-κB1 haploinsufficiency, identify potential genotype-phenotype correlations, and rationally define criteria for treatment in specific patients and circumstances.

Conclusions
NF-κB1 haploinsufficiency is a primary immunodeficiency resulting in impaired function of B cells and can, at least in some patients, lead to decreased function of T cells as well.
The clinical phenotype of NF-κB1 haploinsufficiency may be more complex with a more heterogeneous clinical disease pattern than previously described, including EBV-associated lymphoproliferation as a previously unrecognized clinical feature of the disease. Future studies will need to delineate the full clinical and immunological disease spectrum and prospectively assess different therapeutic approaches.
Acknowledgments Open access funding provided by Medical University of Vienna. We thank the index patient and her family for their participation in this study. We thank Nina Serwas and Elisabeth Salzer for their help with generating parts of Fig. 1. Next generation sequencing was performed at the Biomedical Sequencing Facility at the CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences and the Medical University of Vienna. The work was supported by the Austrian Science Fund (FWF): project I2250-B28 to K.B.

Compliance with Ethical Standards
Conflict of Interest The authors declare that they have no conflict of interest.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.