• Synergistic effects were observed in the granule mediated lymphocyte cytotoxicity.

  • Digenic pathogenesis contributed to the development of hemophagocytic lymphohistiocytosis.

Several molecules (LYST, AP3, RAB27A, STX11, STXBP2, MUNC13-4, and PRF1) have been associated with the function of cytotoxic lymphocytes. Biallelic defects in all of these molecules have been associated with familial hemophagocytic lymphohistiocytosis (FHL). We retrospectively reviewed the genetic and immunology test results from 2701 patients with a clinically suspected diagnosis of hemophagocytic lymphohistiocytosis and found 28 patients with single heterozygous mutations in 2 FHL-associated genes. Of these patients, 21 had mutations within PRF1 and a degranulation gene, and 7 were found to have mutations within 2 genes involved in the degranulation pathway. In patients with combination defects involving 2 genes in the degranulation pathway, CD107a degranulation was decreased, comparable to patients with biallelic mutations in one of the genes in the degranulation pathway. This suggests a potential digenic mode of inheritance of FHL as a result of a synergistic function effect within genes involved in cytotoxic lymphocyte degranulation.

Familial hemophagocytic lymphohistiocytosis (FHL) is an autosomal recessive immune disorder with defective lymphocyte granule–mediated cytotoxicity.1,2  Several genes have been implicated in FHL. FHL2, -3, -4, and -5 and Griscelli syndrome involve PRF1, UNC13D (MUNC13-4), STX11, STXBP2, and RAB27A, respectively.3-6  All of these genes, except for PRF1, play a role in the degranulation of cytotoxic lymphocytes, which affects the precise steps of docking, priming, and fusion of the cytotoxic lymphocytes to the target cell.7,8  Perforin is then released from granules to assist in the penetration and delivery of granzymes into the cytosol of the target cell, where the granzymes cleave key substrates to initiate apoptotic cell death.3,7  An uncompromised and coordinated function of all of these molecules is essential for normal lymphocyte cytotoxic activity. Although it is well known that FHL is inherited in an autosomal recessive manner, symptomatic heterozygous mutation carriers have been presented in several reports.9,10  In these cases, additional synergistic defects in the molecules of the cytotoxic pathway may contribute to the final development of hemophagocytic lymphohistiocytosis, which suggests a digenic inheritance of FHL.

In digenic inheritance, the combination of concurrent partial defects in 2 genes within the same pathway gives rise to a clinical phenotype, whereas a heterozygous state in either of the genes alone results in a less severe phenotype or none at all. Examples of a digenic mode of inheritance have been described for an array of disorders, including retinitis pigmentosa, holoprosencephaly, deafness, epidermolysis bullosa, Hirschsprung disease, insulin resistance, and polycystic kidney disease.11,12  To investigate potential digenic inheritance for FHL, we reviewed patients with clinical FHL and identified those with heterozygous variants in 2 FHL-associated genes; we then examined the age of onset, defective degranulation, perforin levels, and natural killer (NK) cell function activity compared with patients with one (heterozygous) or 2 sequence variants in a specific gene. We present data suggestive of synergistic heterozygosity in FHL between 2 genes involved in cytotoxic lymphocyte degranulation.

This was a retrospective chart review, approved by the Cincinnati Children’s Hospital Medical Center institutional review board, which included 2701 patients referred by their physicians for genetic testing because of clinically suspected FHL. This cohort included affected individuals with homozygous or compound heterozygous mutations in a single gene, single heterozygous individuals with a single mutation in one gene, and double heterozygous (digenic) individuals with 2 heterozygous mutations in 2 separate genes. We analyzed clinical information, clinical immunology testing results, and/or genetic testing results that were available for these patients.13-15 

Interpretation of sequence variants was based on the algorithm described by the Human Genome Variation Society (http://www.hgvs.org/mutnomen/recs.html) and included the use of predictions with Alamut 2.2.2 (Interactive Biosoftware, Rouen, France) for SIFT (Sorting Intolerant From Tolerant),16  PolyPhen-217  and the Grantham Scale.18  Alamut also integrates databases for variant frequencies (dbSNP, 1000 Genomes Project, and Exome Variant Server) as well as splice site prediction algorithms (Human Splicing Finder) and the Human Gene Mutation Database for reported mutations (data not shown).

In this study, a total of 2701 patients with clinically suspected FHL underwent genetic and immunologic testing. Among these patients, 28 patients (P1 to P28) were found to be heterozygous for either a known mutation or a likely pathogenic variant in 2 distinct cytotoxic pathway genes (PRF1, MUNC13-4, STXBP2, STX11, and RAB27A). Genetic variants and clinical information from these patients were used to examine the possibility of digenic inheritance in FHL (Table 1). Fifteen patients were found to have heterozygous variants in PRF1 and MUNC13-4, 6 in PRF1 and STXBP2, 4 in MUNC13-4 and STXBP2, 1 in MUNC13-4 and STX11, 1 in STXBP2 and STX11, and 1 in STXBP2 and RAB27A.

Table 1

Double heterozygous patients with FHL

Patient IDAge at diagnosis, yPRF1MUNC13-4STXBP2STX11RAB27A
P1 0.25 c.1310 C>T (p.A437V) c.169 G>T (p.E57X) ND ND ND 
P2 0.75 c.272 C>T (p.A91V) c.2709+6 G>T ND ND ND 
P3 0.92 c.992 C>T (p.S331L) c.1232 G>A (p.R411Q) NMI NMI NMI 
P4 2.25 c.272 C>T (p.A91V) c.227 C>T (p.T76M) ND ND ND 
P5 c.272 C>T (p.A91V) c.869 C>T (p.S290L) ND ND ND 
P6 c.272 C>T (p.A91V) c.2243 C>T (p.A748V) ND NMI ND 
P7 c.1229 G>A (p.R410Q) c.1036 G>A (p.D346N) ND NMI NMI 
P8 c.272 C>T (p.A91V) p.3160 A>G (p.I1054V) ND ND ND 
P9 c.10 C>T (p.R4C) c.3232 G>C (p.A1078P) NMI NMI ND 
P10 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI ND ND 
P11 10 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI NMI NMI 
P12 12 c.50 delT c.1579 C>T (p.R527W) NMI NMI NMI 
P13 13 c.445 G>A c.2896 C>T (p.R966W) ND ND ND 
P14 13 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI NMI NMI 
P15 28 c.272 C>T (p.A91V) c.182 A>G (p.Y61C) NMI NMI NMI 
P16 c.272 C>T (p.A91V) NMI c.1034 C>T (p.T345M) NMI NMI 
P17 10 c.272 C>T (p.A91V) ND c.1034 C>T (p.T345M) ND ND 
P18 16 655 T>A (Y219N) NMI c.1034 C>T (p.T345M) NMI ND 
P19 21 c.272 C>T (p.A91V) NMI c.1586 G>C (p.R529P) NMI NMI 
P20 24 c.50 delT NMI c.1459 G>A (p.V487M) NMI NMI 
P21 24 c.272 C>T (p.A91V) NMI c.795−4 C>T NMI NMI 
P22 0.167 NMI c.2896 C>T (p.R966W) c.911 C>T (p.T304M) NMI NMI 
P23 0.417 NMI c.1389+1 G>A c.1782*12 G>A NMI ND 
P24 0.667 NMI c.2828 A>G (p.N943S) c.1782*12 G>A NMI ND 
P25 NMI c.2828 A>G (p.N943S) c.715 C>T (p.P239S) NMI NMI 
P26 0.167 NMI c.2030 T>C (p.I677T) NMI c.221 C>T (p.T74M) ND 
P27 14 NMI NMI c.568 C>T (p.R190C) c.9 C>A (p.D3E) NMI 
P28 NMI NMI c.1034 C>T (p.T345M) NMI c.295 T>G (p.F99V) 
Patient IDAge at diagnosis, yPRF1MUNC13-4STXBP2STX11RAB27A
P1 0.25 c.1310 C>T (p.A437V) c.169 G>T (p.E57X) ND ND ND 
P2 0.75 c.272 C>T (p.A91V) c.2709+6 G>T ND ND ND 
P3 0.92 c.992 C>T (p.S331L) c.1232 G>A (p.R411Q) NMI NMI NMI 
P4 2.25 c.272 C>T (p.A91V) c.227 C>T (p.T76M) ND ND ND 
P5 c.272 C>T (p.A91V) c.869 C>T (p.S290L) ND ND ND 
P6 c.272 C>T (p.A91V) c.2243 C>T (p.A748V) ND NMI ND 
P7 c.1229 G>A (p.R410Q) c.1036 G>A (p.D346N) ND NMI NMI 
P8 c.272 C>T (p.A91V) p.3160 A>G (p.I1054V) ND ND ND 
P9 c.10 C>T (p.R4C) c.3232 G>C (p.A1078P) NMI NMI ND 
P10 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI ND ND 
P11 10 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI NMI NMI 
P12 12 c.50 delT c.1579 C>T (p.R527W) NMI NMI NMI 
P13 13 c.445 G>A c.2896 C>T (p.R966W) ND ND ND 
P14 13 c.272 C>T (p.A91V) c.2896 C>T (p.R966W) NMI NMI NMI 
P15 28 c.272 C>T (p.A91V) c.182 A>G (p.Y61C) NMI NMI NMI 
P16 c.272 C>T (p.A91V) NMI c.1034 C>T (p.T345M) NMI NMI 
P17 10 c.272 C>T (p.A91V) ND c.1034 C>T (p.T345M) ND ND 
P18 16 655 T>A (Y219N) NMI c.1034 C>T (p.T345M) NMI ND 
P19 21 c.272 C>T (p.A91V) NMI c.1586 G>C (p.R529P) NMI NMI 
P20 24 c.50 delT NMI c.1459 G>A (p.V487M) NMI NMI 
P21 24 c.272 C>T (p.A91V) NMI c.795−4 C>T NMI NMI 
P22 0.167 NMI c.2896 C>T (p.R966W) c.911 C>T (p.T304M) NMI NMI 
P23 0.417 NMI c.1389+1 G>A c.1782*12 G>A NMI ND 
P24 0.667 NMI c.2828 A>G (p.N943S) c.1782*12 G>A NMI ND 
P25 NMI c.2828 A>G (p.N943S) c.715 C>T (p.P239S) NMI NMI 
P26 0.167 NMI c.2030 T>C (p.I677T) NMI c.221 C>T (p.T74M) ND 
P27 14 NMI NMI c.568 C>T (p.R190C) c.9 C>A (p.D3E) NMI 
P28 NMI NMI c.1034 C>T (p.T345M) NMI c.295 T>G (p.F99V) 

ND, not done; NMI, no mutation identified.

*

Denotes the variant is located at the 3'UTR of the gene.

Patients were separated into 2 groups: heterozygotes with variants in 2 genes involved in degranulation (MUNC13-4, STXBP2, STX11, and RAB27A [denoted as Deg/Deg]; 7 total) and heterozygotes with variants in PRF1 and one degranulation gene (denoted as PRF1/Deg; 21 total). Controls included a group of affected patients who were homozygous or compound heterozygous with PRF1 (120 patients identified), MUNC13-4 (72), or STXBP2 (33) mutations and a group of patients who were single heterozygous for single mutations within PRF1 (163), MUNC13-4 (119), or STXBP2 (29).

Deg/Deg patients had an earlier age of onset (≤24 months) for 71.4% (5 of 7) of patients (Figure 1A). In contrast, PRF1/Deg patients were observed to have an early age of onset (≤24 months) for only 14.3% (3 of 21) of patients (Figure 1A). The age of onset for Deg/Deg patients (≤24 months) was similar to that of patients with biallelic mutations in any one of these homozygous or compound heterozygous genes: 83.3% (100 of 120) of patients with PRF1, 83.3% (60 of 72) patients with MUNC13-4, and 69.7% (23 of 33) patients with STXBP2 (Figure 1A). In contrast, the single heterozygous patients were observed to have a later age of onset (Figure 1A). Thus, early age of onset in patients with 2 genetic defects within the degranulation pathway supports a potential digenic mode of inheritance for FHL.

Figure 1

Comparison of age of onset, CD107a degranulation, and perforin expression. Percentage of patients with (A) early age of onset (≤24 months), (B) degranulation based on CD107a expression, and (C) perforin expression levels were compared among the groups. Patients were grouped as Affected (homozygous or compound heterozygous), Single heterozygotes, or Double heterozygotes and then further separated on the basis of the genes affected, as shown by the color bar and the dotted lines between groups.

Figure 1

Comparison of age of onset, CD107a degranulation, and perforin expression. Percentage of patients with (A) early age of onset (≤24 months), (B) degranulation based on CD107a expression, and (C) perforin expression levels were compared among the groups. Patients were grouped as Affected (homozygous or compound heterozygous), Single heterozygotes, or Double heterozygotes and then further separated on the basis of the genes affected, as shown by the color bar and the dotted lines between groups.

Close modal

Upon degranulation, NK cells express CD107a on their surface. Thus, decreased measure of CD107a expression is indicative of defective degranulation.8,19  To determine which individuals displayed defective degranulation, CD107a expression data from immunologic testing was analyzed. CD107a expression, perforin expression, and NK function results were available for only a small portion of these patients. From the limited number of patients, we found 75% (3 of 4) of patients with Deg/Deg double heterozygotes had defective CD107a degranulation. This finding is similar to the defective degranulation seen in individuals with biallelic MUNC13-4 (91.7%; 11 of 12) and STXBP2 (92.3%; 12 of 13) mutations. In contrast, PRF1/Deg double heterozygotes had normal degranulation (0%; 0 of 8) (Figure 1C). As expected, defective NK function was observed in all patient groups (data not shown).

Patients with defects in PRF1 had decreased perforin levels. Numerically, 100% (41 of 41) of PRF1-affected, 63.6% (42 of 66) of PRF1 heterozygotes, and 72.7% (8 of 11) of PRF1/Deg patients had decreased perforin expression (Figure 1C). Interestingly, only 20% (1 of 5) of Deg/Deg double heterozygotes, 16% (4 of 25) of MUNC13-4–affected, 30.8% (4 of 13) of STXBP2-affected, 33.3% (14 of 42) of MUNC13-4 heterozygotes, and 40% (6 of 15) of STXBP2 heterozygotes had low perforin in any of these categories (Figure 1C). Patients without PRF1 mutations had increased or normal levels of perforin (data not shown).

Synergistic function was hypothesized for the genes in the degranulation pathway because of the known highly precise and coordinated functions that lead to the secretion of granzyme B to target cells. Findings from this retrospective study suggest that heterozygous defects in 2 degranulation pathway genes have a synergistic deleterious effect and thus imply that a digenic mode of inheritance can lead to FHL. Patients with mutations in 2 degranulation genes in this study had early childhood presentation of FHL. Those patients also displayed defective degranulation of lymphocytes as expected, unlike PRF1/Deg patients. This is the first retrospective study that suggests a potential digenic inheritance in FHL. However, it is worth noting several limitations of this study. Detailed clinical data were not available for several patients with suspected FLH. Data for analysis of NK function, CD107a degranulation, and perforin expression were also limited. In addition, parental and/or sibling samples were not available in most families for segregation studies. Ideally, future prospective studies should include a thorough clinical summary, outcome results, and data on genetic and functional studies performed on the proband and family members.

There is an Inside Blood Commentary on this article in this issue.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

The authors thank the staff of the Cincinnati Children’s Hospital Medical Center (a Federation of Clinical Immunology Societies Center of Excellence) and the diagnostic immunology and molecular genetics laboratories that performed the clinical tests for this study.

Contribution: K.Z. and A.H.F. designed the research, analyzed and interpreted data, and wrote the manuscript; S.C., C.A.V., and H.C. analyzed and interpreted data and wrote the manuscript; and D.K., J.A.J., and A.H. performed the research and collected and analyzed the data.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Kejian Zhang, Division of Human Genetics, Cincinnati Children’s Hospital Medical Center, 3333 Burnet Ave, ML 4006, Cincinnati, OH 45229; e-mail: kejian.zhang@cchmc.org.

1
Filipovich
 
AH
The expanding spectrum of hemophagocytic lymphohistiocytosis.
Curr Opin Allergy Clin Immunol
2011
, vol. 
11
 
6
(pg. 
512
-
516
)
2
Chandrakasan
 
S
Filipovich
 
AH
Hemophagocytic lymphohistiocytosis: advances in pathophysiology, diagnosis, and treatment.
J Pediatr
2013
, vol. 
163
 
5
(pg. 
1253
-
1259
)
3
Voskoboinik
 
I
Smyth
 
MJ
Trapani
 
JA
Perforin-mediated target-cell death and immune homeostasis.
Nat Rev Immunol
2006
, vol. 
6
 
12
(pg. 
940
-
952
)
4
Feldmann
 
J
Callebaut
 
I
Raposo
 
G
, et al. 
Munc13-4 is essential for cytolytic granules fusion and is mutated in a form of familial hemophagocytic lymphohistiocytosis (FHL3).
Cell
2003
, vol. 
115
 
4
(pg. 
461
-
473
)
5
zur Stadt
 
U
Rohr
 
J
Seifert
 
W
, et al. 
Familial hemophagocytic lymphohistiocytosis type 5 (FHL-5) is caused by mutations in Munc18-2 and impaired binding to syntaxin 11.
Am J Hum Genet
2009
, vol. 
85
 
4
(pg. 
482
-
492
)
6
zur Stadt
 
U
Schmidt
 
S
Kasper
 
B
, et al. 
Linkage of familial hemophagocytic lymphohistiocytosis (FHL) type-4 to chromosome 6q24 and identification of mutations in syntaxin 11.
Hum Mol Genet
2005
, vol. 
14
 
6
(pg. 
827
-
834
)
7
de Saint Basile
 
G
Ménasché
 
G
Fischer
 
A
Molecular mechanisms of biogenesis and exocytosis of cytotoxic granules.
Nat Rev Immunol
2010
, vol. 
10
 
8
(pg. 
568
-
579
)
8
Marcenaro
 
S
Gallo
 
F
Martini
 
S
, et al. 
Analysis of natural killer-cell function in familial hemophagocytic lymphohistiocytosis (FHL): defective CD107a surface expression heralds Munc13-4 defect and discriminates between genetic subtypes of the disease.
Blood
2006
, vol. 
108
 
7
(pg. 
2316
-
2323
)
9
Zhang
 
K
Johnson
 
JA
Biroschak
 
J
, et al. 
Familial haemophagocytic lymphohistiocytosis in patients who are heterozygous for the A91V perforin variation is often associated with other genetic defects.
Int J Immunogenet
2007
, vol. 
34
 
4
(pg. 
231
-
233
)
10
Zhang
 
K
Jordan
 
MB
Marsh
 
RA
, et al. 
Hypomorphic mutations in PRF1, MUNC13-4, and STXBP2 are associated with adult-onset familial HLH.
Blood
2011
, vol. 
118
 
22
(pg. 
5794
-
5798
)
11
Ming
 
JE
Muenke
 
M
Multiple hits during early embryonic development: digenic diseases and holoprosencephaly.
Am J Hum Genet
2002
, vol. 
71
 
5
(pg. 
1017
-
1032
)
12
Schäffer
 
AA
Digenic inheritance in medical genetics.
J Med Genet
2013
, vol. 
50
 
10
(pg. 
641
-
652
)
13
Egeler
 
RM
Shapiro
 
R
Loechelt
 
B
Filipovich
 
A
Characteristic immune abnormalities in hemophagocytic lymphohistiocytosis.
J Pediatr Hematol Oncol
1996
, vol. 
18
 
4
(pg. 
340
-
345
)
14
Kogawa
 
K
Lee
 
SM
Villanueva
 
J
Marmer
 
D
Sumegi
 
J
Filipovich
 
AH
Perforin expression in cytotoxic lymphocytes from patients with hemophagocytic lymphohistiocytosis and their family members.
Blood
2002
, vol. 
99
 
1
(pg. 
61
-
66
)
15
Molleran Lee
 
S
Villanueva
 
J
Sumegi
 
J
, et al. 
Characterisation of diverse PRF1 mutations leading to decreased natural killer cell activity in North American families with haemophagocytic lymphohistiocytosis.
J Med Genet
2004
, vol. 
41
 
2
(pg. 
137
-
144
)
16
Ng
 
PC
Henikoff
 
S
SIFT: Predicting amino acid changes that affect protein function.
Nucleic Acids Res
2003
, vol. 
31
 
13
(pg. 
3812
-
3814
)
17
Ramensky
 
V
Bork
 
P
Sunyaev
 
S
Human non-synonymous SNPs: server and survey.
Nucleic Acids Res
2002
, vol. 
30
 
17
(pg. 
3894
-
3900
)
18
Grantham
 
R
Amino acid difference formula to help explain protein evolution.
Science
1974
, vol. 
185
 
4154
(pg. 
862
-
864
)
19
Bryceson
 
YT
March
 
ME
Barber
 
DF
Ljunggren
 
HG
Long
 
EO
Cytolytic granule polarization and degranulation controlled by different receptors in resting NK cells.
J Exp Med
2005
, vol. 
202
 
7
(pg. 
1001
-
1012
)
Sign in via your Institution