Comparison of Th1/Th2 cytokine profiles between primary and secondary haemophagocytic lymphohistiocytosis

Background Haemophagocytic lymphohistiocytosis (HLH) is a life-threatening disorder of immune regulation, and HLH patients with mutations in genes including PRF1, UNC13D, STX11, STXBP2, SH2D1A, XIAP, and ITK were reported to be primary HLH. Due to the different treatment options, the differentiation between primary and secondary HLH is critical. Our previous studies have showed that a Th1/Th2 cytokine profile is diagnostic for HLH, yet the cytokine profiles between primary and secondary HLH have not been compared. The aim of the study was to test whether the Th1/Th2 cytokine profile could be used as a tool to differentiate between primary and secondary HLH. Methods A total of 45 hospitalized Chinese children with HLH during the period of February 2010 through September 2012 were enrolled in the study. Fifty healthy children were enrolled as controls. Primary HLH related genes were sequenced using genomic DNA samples. The Th1/Th2 cytokine levels including interferon-γ (IFN-γ), tumor necrosis factor-alpha (TNF-α), interleukin (IL)-10, IL-6, IL-4 and IL-2 were quantitatively determined by cytometric bead assay techniques. Results Primary HLH group (n = 4) included one patient with biallelic heterozygous mutations in PRF1 gene, and three patients with hemizygous mutation in SH2D1A gene. Based on the available genetic data, the other 41 patients were classified into the secondary HLH group. When compared the cytokine levels between the two groups, IL-4 level in primary-HLH was significantly lower than that in secondary HLH (P = 0.025), while IFN-γ level in primary HLH had a tendency of statistically lower than that in secondary HLH (P = 0.051). Area under receiver operating characteristic (ROC) curves of IL-4 and IFN-γ, IL-10, TNF-α, IL-2, and IL-6 levels were 0.841, 0.799, 0.506, 0.494, 0.457, and 0.250, respectively. ROC curves showed that 1.7 pg/ml of IL-4 had sensitivity and specificity for differentiation between primary and secondary HLH as 70.7 and 100.0 %, while 433.9 pg/ml of IFN-γ had sensitivity and specificity as 51.2 and 100.0 %, respectively. Conclusions HLH patients with lower IL-4 and IFN-γ levels have higher possibility to be primary HLH. The cytokine profile may be used as an additional tool for the quick differential diagnosis between primary and secondary HLH. Electronic supplementary material The online version of this article (doi:10.1186/s13052-016-0262-7) contains supplementary material, which is available to authorized users.


Background
Haemophagocytic lymphohistiocytosis (HLH) is a lifethreatening disorder of immune regulation, characterized by a highly stimulated, but ineffective immune response to antigens, which results in cytokine storm and inflammatory reaction [1]. HLH is not a single entity, but a clinical syndrome that can be encountered in association with various underlying diseases leading to similar characteristic clinical and laboratory presentations. Briefly, the diagnosis of HLH requires either a genetic diagnosis, or fulfillment of 5 out of 8 clinical criteria including fever, splenomegaly, bicytopenia, either hypofibrinoginemia or hypertriglyceridemia, hyperferritinemia, elevated soluble interleukin-2 receptor (sCD25 or sIL2R), impaired natural killer (NK) cell cytolytic function, or the observation of hemophagocytosis in bone marrow, spleen, or lymph nodes [2]. Based on the etiology, HLH can be classified into two types, primary HLH (with genetic mutations), and secondary HLH (no known mutations). Epstein-Barr virus (EBV) infections seem to be very common in HLH patients, especially in Asian countries [3][4][5][6][7]. Beside EBV infection, cytomegalovirus (CMV) infection is another common cause in HLH patients [8,9]. Clinically, prompt and accurate diagnosis is critical to initiate definitive therapy. Unfortunately, the clinical features are identical in primary and secondary HLH, and both forms are often triggered by infections, so it is difficult to distinguish between these two types [8]. Furthermore, the diagnosis of HLH based on the current combination of clinical, laboratory and immunological criteria is challenging because all of the criteria are not specifically diagnostic for any HLH subtypes.
Primary HLH can be further divided into familial hemophagocytic lymphohistiocytosis (FHL) and immune deficiencies associated HLH [10]. Five genetic defects have been identified in FHL. A potential gene locus (FHL1) has been reported to be associated with HLH, but the specific gene involved has not yet been identified [11]. The first described FHL related gene is PRF1 gene (FHL2) [12], a gene encodes perforin protein. The next identified cause of FHL (FHL3) is UNC13D gene [13], which encodes Munc13-4 protein. Later, mutations of STX11 gene and STXBP2 gene are found to be the causes of FHL4 [14] and FHL5 [15,16], respectively. Immune deficiencies associated HLH occur with significant frequency in X-linked lymphoproliferative syndrome type 1 (XLP1) [17], XLP type 2 (XLP2) [18], and IL-2inducible T cell kinase deficiency-associated lymphoproliferation [19], which are characterized by mutations in SH2D1A, XIAP, and ITK genes, respectively.
Hypercytokinemia is a hallmark of HLH. Henter JI et al. suggested that hypercytokinemia may be caused by a genetic defect in cytokine regulation as FHL patients showed high cytokine levels [20]. Using the rapid cytometric bead array (CBA) technique, a specific cytokine profile (significant increase of IFN-γ and IL-10, combined with a slightly increased level of IL-6) for childhood HLH was described by our group in 2008 [21]. In our clinical practice, this cytokine profile was helpful for the early diagnosis of HLH and for the differentiation from other disease entities [22,23]. However, one limitation of our previous cytokine assay was that we did not perform genetic sequencing of HLH involved genes simultaneously. In this study, we tried to sequence PRF1, UNC13D, STX11, STXBP2, SH2D1A, XIAP, and ITK genes, together with cytokine determination for HLH patients to test whether cytokine profile could be used as a tool to distinguish between primary HLH and secondary HLH patients.

Patients
The diagnoses of all our patients were made based on the HLH-2004 criteria [2]. A cohort of 45 consecutively hospitalized HLH patients were enrolled in the study, which included 27 males and 18 females with a male-to-female ratio of 1.5:1, and a mean age of 3.7 years with a range of 8 days through 12.3 years. Patients who had been treated by steroids before referral to our hospital were excluded. A total of 50 unrelated healthy individuals (30 males, and 20 females) matched for age, and race were recruited as healthy controls. All the patients and healthy controls were from our hospital between February 2010 and September 2012. No HLH patient showed signs of mucocutaneous albinism. The study protocol was reviewed and approved by the Ethics Committee at Children's Hospital Zhejiang University School of Medicine and written informed consents were obtained from all participants' parents or guardians before this study.

Genetic analysis
Genomic DNA was isolated from peripheral blood using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany). The coding regions and flanking intronic sequences of PRF1, UNC13D, STX11, STXBP2, SH2D1A, XIAP, and ITK genes, and the deep intronic sequences of intron 1 in UNC13D gene were amplified by a polymerase chain reaction (PCR) machine (BIO-RAD) in our laboratory, and the PCR products were sequenced by a DNA sequencer in Invitrogen Company (Shanghai, China). We analyzed PRF1, UNC13D, STX11, STXBP2, and ITK genes in all 45 patients and 50 controls, and all the male patients (n = 27) and male controls (n = 30) were sequenced in their SH2D1A and XIAP genes. The mutations were validated by re-sequencing an independent PCR-generated amplicon from the subjects. The variants were named according to the Human Genome Variation Society and journal requirements.

Cytokine determination
All 45 HLH patients and 50 controls were tested for the six cytokines including IL-2, IL-4, IL-6, IL-10, TNF-α, and IFN-γ. Peripheral blood samples were collected, transferred to a serum separating tube and centrifuged at 1000 g at 20°C for 20 min after clotting. The serum was carefully harvested, and the determination of the cytokines was performed immediately, or if the situation was not so urgent, the aliquot was temporarily stored at 2°C to 8°C until analysis (usually within 12 h). Concentrations of the six cytokines aforementioned were quantitatively determined using the CBA Human Th1/Th2 Cytokine Kit II (BD Biosciences, San Jose, California) as described previously [21]. The minimum and maximum limits of detection for all six cytokines were 1 and 5000 pg/ml, respectively.

Degranulation assays with flow cytometry
Degranulation assay were tested in 36 HLH patients and 43 controls. The assay was performed as previously described with modifications [24]. Peripheral blood mononuclear cells (PBMC) were isolated by Ficoll gradient centrifugation, rested in Iscove's modified Dulbecco medium (IMDM, Invitrogen, Carlsbad, California, USA) supplemented with 10 % fetal bovine serum (FBS, Sijiqing, Hangzhou, China) at 37°C in a humidified atmosphere of 5 % CO 2 for 2 h, and washed in 1 × PBS (phosphate buffer saline). K562 cells were obtained from the American Type Culture Collection (ATCC, Rockefeller, Maryland, USA) and cultured in RPMI1640 medium (Invitrogen, Carlsbad, California, USA) supplemented with 10 % FBS. All the antibodies including CD3-FITC, CD8-PerCP, CD56-APC, CD107a-PE (H4A3, IgG1), and isotype controls were purchased from Becton Dickinson (San Jose, CA, USA). 5 × 10 5 PBMCs were mixed with 5 × 10 4 sensitive target cells (K562) as stimulants in 12 flat bottom well plates and incubated for 2 h at 37°C. An equal volume of RPMI-1640 culture medium to replace the volume of K562 cells were used as negative control. Then the cells were stained with 10μl of CD3-FITC, 10μl of CD8-PerCP, 2.5μl of CD56-APC, and 10μl of CD107a-PE (isotype IgG1-PE as negative control) for 30 min at 4°C in the dark. After two washes with 1 × PBS of the above samples, flow cytometric analyses were performed by a FACSCalibur cytometer (Becton-Dickinson, San Jose, CA, USA). CD3 -CD56 + and CD3 + CD8 + T cells were gated as NK cells and Cytotoxic T cells (CTLs), respectively. Data were acquired with CellQuest software (BD Bioscience).

Clinical data
Information at diagnosis of age, sex, with fever or not, hemoglobin levels, platelets numbers, white blood cell counts, percentage of neutrophils, total neutrophil counts, percentage of lymphocytes, total lymphocyte counts, triglyceride levels, fibrinogen levels, LDH levels, ferritin levels, sIL2R levels, and sIL2R/ferritin ratio were collected. Quantitative real-time PCR was used to detect EBV-DNA and CMV-DNA copies in sera.

Statistical analysis
Serum concentrations of individual cytokines, clinical data were compared between groups using the Mann-Whitney U test. A chi-square test was used to assess ratio differences between groups. Receiver operating characteristic (ROC) curves were derived from the cytokine levels for all HLH patients. In a ROC curve, the sensitivity and specificity of all six cytokines were calculated for the differentiation between primary and secondary HLH. All statistical analyses were performed using SPSS 17.0 software (SPSS Inc, Chicago, Illinois). A two-sided P-value <0.05 was considered to be statistically significant.

Genetic results
The genetic findings of the 45 HLH patients were shown in Table 1, Additional file 1: Figure S1, Additional file 2: Figure S2, Additional file 3: Figure S3, Additional file 4: Figure S4, Additional file 5: Figure S5, and Additional file 6: Figure S6. All variants classified as pathogenic were not detected in the controls, while those classified as single nucleotide polymorphisms (SNPs) were found in both HLH patients and healthy controls. One patient with biallelic heterozygous mutations in PRF1 gene (P2), and three patients with hemizygous mutation in SH2D1A gene (P16, P17 and P26) were categorized into primary HLH group (n = 4). Based on the available genetic data, the other 41 patients were classified into the secondary HLH group, including nine HLH patients with single heterozygous mutation (mutations in patients P6, P10, and P45 were already reported in literature, while mutations in patients P1, P3, P4, P5, P7, and P11 were not reported before), three HLH patients with only SNPs (P8, P9, and P38), and the remaining 29 HLH patients without any SNPs.
Database investigation and in silico prediction of mutations in HLH patients were shown in Table 2. The in silico prediction results between PolyPhen-2 and Sorting Intolerant From Tolerant (SIFT) were consistent with each other, except p.Arg1087Trp of UNC13D gene (PolyPhen-2 predicted this change was BENIGN with a score of 0.002, while SIFT indicated that this change would AFFECT PROTEIN FUNCTION with a score of 0.00).
When we compared the levels of the cytokines between the primary and secondary HLH groups, the IL-4 level in primary-HLH was significantly lower than that in secondary HLH (P = 0.025), with mean (95 % CI) in primary and secondary HLH groups showed as 1.1 (0.6-1.6) and 2.3 (1.9-2.7), respectively. Additionally, IFN-γ level in primary HLH had a tendency of statistically lower than that in secondary HLH (P = 0.051), with mean (95 % CI) in primary and secondary HLH groups showed as 106.6 (minus-393.9) pg/ml and 905.7 (530.7-1280.6) pg/ml, respectively. The 95 % CI gap of IL-4 between primary and secondary HLH groups was 1.6-1.9 pg/ml, while the gap of IFN-γ between the two groups was 393.9-530.7 pg/ml. The levels of the remaining four cytokines including IL-2, IL-6, IL-10, and TNF-α were not significantly different between primary and secondary HLH groups (P = 0.78, P = 0.102, P = 0.968, and P = 0.968, respectively). The area under the ROC curve, referred to as the AUC, is an appropriate measure for describing the overall accuracy of a diagnostic test, and higher AUC value mean better diagnostic value. AUCs of IL-4, IFN-γ, IL-10, TNF-α, IL-2, and IL-6 levels were calculated to be 0.841, 0.799, 0.506, 0.494, 0.457, and 0.250, respectively, indicating that levels of IL-4 and IFN-γ may be used as additional tools for the quick differential diagnosis between primary and secondary HLH (Fig. 2). ROC curves of IL-4 levels showed that 1.6 pg/ml, 1.7 pg/ml, 1.8 pg/ml, and 1.9 pg/ml had sensitivity and specificity as 73.2 and 75.0 %, 70.7 and 100.0 %, 68.3 and 100.0 %, 63.4 and 100 %, respectively. Furthermore, ROC curves of IFN-γ levels showed that 360.6 pg/ ml, 433.9 pg/ml, and 539.7 pg/ml had sensitivity and specificity as 51.2 and 75.0 %, 51.2 and 100.0 %, 48.8 and 100 %, respectively. Taken the results of IL-4 and IFN-γ together, we proposed that HLH patients with IL-4 below 1.7 pg/ml and IFN-γ below 433.9 pg/ml had a higher chance to be primary HLH.

Degranulation and clinical data
Information of surface CD107a level in resting NK cells, age at diagnosis, sex, with fever or not, hemoglobin levels, platelet counts, white blood cell counts, percentage of neutrophils, absolute neutrophil counts, percentage of lymphocytes, absolute lymphocyte counts, triglyceride levels, fibrinogen levels, LDH levels, ferritin levels, sIL2R levels, sIL2R/ferritin ratio, EBV-DNA, and CMV-DNA copies were shown in Table 3. The results showed that, except the percentage of neutrophils (P = 0.012) and the percentage of lymphocytes (P = 0.012), there was no significant difference between primary HLH and secondary HLH groups in any other factors.

Discussion
Inheritance of HLH is known to be autosomal recessive. Based on the data from published studies, primary HLH usually includes patients with homozygous, hemizygous, or compound heterozygous mutations, as mutations of these types affect protein function in a known way. In our HLH patients, no homozygous mutation was detected, and all were hemizygous or heterozygous variants. Similar to our results, Zhizhuo H et al. investigated 67 Chinese HLH patients and identified eight patients with variants in primary HLH related genes, included one patient with a hemizygous mutation of XIAP gene, three patients with compound heterozygous mutations, and four patients with single heterozygous mutation [27]. How does the single heterozygous mutation affect its coding protein's function is complicated and remains to be elucidated. Based on the recent study of Spessott WA et al. in 2015 [28], some single heterozygous mutation in special site can cause primary HLH, yet the dominant negative effect can not be generalized, and people should do further study to prove the related protein function is damaged if they want to assign some HLH patients with monoallelic mutation as primary HLH. Another possibility is that the risk of having missed pathogenic variants is high, and possible additional variants, such as deep intronic variants similar to mutations in deep intron 1 of UNC13D gene [29][30][31], might have been missed for some cases.
In our previous study, we found that HLH patients presented a specific cytokine profile of highly increased levels of IFN-γ and IL-10, and a moderately increased level of IL-6 [23]. In this study, the levels of IL-6, IL-10, and IFN-γ in both primary and secondary HLH groups were significantly higher (P < 0.05) than those of controls, and no patient was overlapped with our previous study, which further confirmed the usefulness of this cytokine profile for the diagnosis of HLH.
When we compared the levels of the cytokines between the primary and secondary HLH groups, the IL-4 level in primary-HLH was significantly lower than that in secondary HLH (P = 0.025), with the gap of IL-4 between the two groups was 1.6-1.9 pg/ml, and IFN-γ level in primary HLH had a tendency of statistically lower than that in secondary HLH (P = 0.051), with the gap of IFN-γ between the two groups was 393.9-530.7 pg/ml. AUCs of six cytokines indicated that levels of IL-4 and IFN-γ may be used as additional tools for the quick differential diagnosis between primary and secondary HLH. Taken the results of IL-4 and IFN-γ together, we propose that HLH patients with IL-4 below 1.7 pg/ml and IFN-γ below 433.9 pg/ml have a higher possibility to be primary HLH.
There have been studies trying to find other discriminators to differentiate primary from secondary cases of HLH. According to the report by Bryceson et al. [26], degranulation assay has a high sensitivity and specificity rate for discrimination between two types of HLH (type1 includes FHL3, FHL4, and FHL5, while type2 includes FHL2, XLP1, XLP2, and secondary HLH cases). In this cohort, all primary HLH cases belong to type2 HLH including FHL2 and XLP1, and it is not surprising that we can not find difference of degranulation between primary HLH and secondary HLH groups. Yasumi et al. proposed that the percentage of total lymphocytes, serum levels of LDH, and the sIL2R/ferritin ratio could differentiate familial from secondary HLH [32]. In this study, when we compare to the secondary HLH patients, primary HLH patients have higher percentage of lymphocytes (mean level 75.3 % vs 48.0 %, P = 0.012), and lower percentage of neutrophils (12.8 % vs 40.7 %, P = 0.012), which is consistent with Yasumi et al's findings.
Our study has some limitations. First, the risk of having missed pathogenic variants is high, and possible additional variants, such as deep intronic variants, regulatory variants, and complex structural variants, might have been missed for some cases, which would have great impact on the grouping strategy and then influence the interpretation of the results. Second, as this is a single-center study, and the results are based on data from small number of samples and incomplete genetic sequencing results, a multi-center study was required to be performed to validate the results.