Key Points
Whole-exome sequencing may identify specific therapeutic opportunities for patients with HLH.
HLH should be conceptualized as a critical illness phenotype driven by toxic activation of immune cells from different underlying mechanisms.
Abstract
The HLH-2004 criteria are used to diagnose hemophagocytic lymphohistiocytosis (HLH), yet concern exists for their misapplication, resulting in suboptimal treatment of some patients. We sought to define the genomic spectrum and associated outcomes of a diverse cohort of children who met the HLH-2004 criteria. Genetic testing was performed clinically or through research-based whole-exome sequencing. Clinical metrics were analyzed with respect to genomic results. Of 122 subjects enrolled over the course of 17 years, 101 subjects received genetic testing. Biallelic familial HLH (fHLH) gene defects were identified in only 19 (19%) and correlated with presentation at younger than 1 year of age (P < .0001). Digenic fHLH variants were observed but lacked statistical support for disease association. In 28 (58%) of 48 subjects, research whole-exome sequencing analyses successfully identified likely molecular explanations, including underlying primary immunodeficiency diseases, dysregulated immune activation and proliferation disorders, and potentially novel genetic conditions. Two-thirds of patients identified by the HLH-2004 criteria had underlying etiologies for HLH, including genetic defects, autoimmunity, and malignancy. Overall survival was 45%, and increased mortality correlated with HLH triggered by infection or malignancy (P < .05). Differences in survival did not correlate with genetic profile or extent of therapy. HLH should be conceptualized as a phenotype of critical illness characterized by toxic activation of immune cells from different underlying mechanisms. In most patients with HLH, targeted sequencing of fHLH genes remains insufficient for identifying pathogenic mechanisms. Whole-exome sequencing, however, may identify specific therapeutic opportunities and affect hematopoietic stem cell transplantation options for these patients.
Introduction
Isolated cases of children with fever, fatal systemic inflammation, and bone marrow or tissue specimens demonstrating hemophagocytosis were first described in 1939,1 followed by familial reports that suggested a potential genetic predisposition toward hemophagocytic lymphohistiocytosis (HLH).2 Consensus criteria were developed to standardize clinical trial enrollment (supplemental Table 1, available on the Blood Web site)3 and now serve as de facto diagnostic criteria for HLH. HLH, however, might be better conceptualized as a convergent entity or common final pathway/phenotype characterized by toxic immune activation secondary to a range of diverse underlying molecular diagnoses rather than as a distinct disease.
Historically, patients have long been categorized as having either “primary” HLH from an inherited defect in cytolytic immune function or “secondary” HLH resulting from immune activation incited by various antigenic challenges, including autoimmune disease, persistent infection, or malignancy. For primary or familial HLH (fHLH), mutations in genes needed for cytolytic immune function have been identified and serve as part of the HLH-2004 diagnostic criteria.3,4 The clinical distinctions between primary and secondary HLH may not be this straightforward; in HLH-94, no difference in survival was observed for patients with presumed primary vs secondary HLH.5 Thus, it becomes apparent that clinical fulfillment of the HLH-2004 criteria does not necessarily imply defective cytolytic activity for all patients. Similarly, the spectrum of inherited genetic defects that inform a need for immune suppression or hematopoietic stem cell transplantation (HSCT) are inadequately captured by primary or secondary HLH categories.
The uncertain performance of the HLH-2004 criteria and dire consequences of missed or incorrect HLH diagnosis have led to widely variable approaches to diagnosis and treatment of patients.6,7 If untreated, fHLH is fatal in 90% of children, and prompt initiation of immunochemotherapy in patients with active HLH remains critical for survival.8-10 HSCT is necessary for long-term survival in patients at risk for recurrent episodes of HLH.7 Alternately, it has been proposed that clinical diagnosis of HLH using the HLH-2004 criteria with no confirmed gene mutations may prove dangerous and that these criteria fail to distinguish HLH from sepsis or systemic inflammatory response syndrome, for which immunosuppression or HSCT may be contraindicated.6 Targeted sequencing of fHLH genes in ethnically enriched cohorts has identified disease-causing gene mutations in 34% to 100% of affected individuals with greater enrichment when genetic testing is restricted to patients with abnormal expression of HLH-related proteins or defective cytolytic function.11-14 In this study, we sought to evaluate the genomic spectrum and associated outcomes of a prospectively enrolled, multiethnic cohort of children clinically diagnosed with HLH, using the HLH-2004 criteria at a large tertiary care center.
Methods
Subject enrollment
Patients at Texas Children’s Hospital or collaborating referral centers who met HLH-2004 criteria for the diagnosis of HLH were offered participation in this research study, approved by the Baylor College of Medicine Institutional Review Board. All subjects and family members provided written informed consent to have their clinical and genetic information published in medical or scientific journals. All procedures involving human participants were performed in accordance with institutional and international ethical standards.
Ascertainment of clinical data
Fulfillment of the HLH-2004 diagnostic criteria was confirmed for all subjects. Comorbidities, clinical outcomes, and clinical genetic test results, including targeted individual gene sequencing, HLH gene panel testing, and clinical whole-exome sequencing (WES), were derived from the medical record. Clinical WES was conducted by Baylor Genetics Laboratories (Houston, TX). None of the subjects had known features of or mutations in genes associated with primary immunodeficiency diseases or dysregulated immune activation or proliferation before HLH onset and diagnosis.
Research-based WES and analysis
WES of probands who lacked clinical molecular diagnoses and their family members was conducted within the Baylor-Hopkins Center for Mendelian Genomics, as previously described (supplemental Methods).15,16 For all subjects who underwent HSCT, pretransplantation samples were sequenced.
Variant selection methods have been previously reported.15,16 Briefly, variants were prioritized using 7 criteria: existence in known HLH-causing genes; appropriate cosegregation with phenotype among family members; allelic frequencies in internal and public (eg, Exome Aggregation Consortium) exome databases less than 0.005 for autosomal and X-linked recessive inheritance models and less than 0.001 for autosomal dominant inheritance models; in silico conservation and damage prediction (eg, Combined Annotation Dependent Depletion score17 greater than the Mutational Significance Cutoff18 ); expression in immune cells, especially T cells, natural killer (NK) cells, monocytes, and macrophages (eg, http://biogps.org); relevance of the affected gene to known diseases that predispose toward infection, autoimmunity, or malignancy (eg, in the Online Mendelian Inheritance for Man catalog); and known function of the affected protein in immune function and inflammation (eg, in PubMed). Stratified disease-associated variants were confirmed by Sanger sequencing. Variants in common fHLH genes (PRF1, STX11, STXBP2, and UNC13D plus LYST and RAB27A) were scrutinized for possible contributions to disease.
Categorization of variants in primary immunodeficiency disease genes identified by WES
For variants in primary immunodeficiency disease (PIDD)-associated genes, molecular diagnoses were designated as “likely” (supplemental Table 4) if HLH had been previously demonstrated in other individuals with the associated PIDD, the associated PIDD was known to impair cytolytic T or NK cell function, or the molecular diagnosis was consistent with relevant features of the clinical phenotype of the subject. If the molecular diagnosis was otherwise solely supported by the exome data analyses, it was considered a “possible” genetic explanation for HLH in the subject.
Statistical analysis
Demographic and clinical information were abstracted from medical records. The χ-squared test was used if counts exceeded 5; otherwise Fisher’s Exact test was implemented. Logistic regression was used to obtain odds ratios and 95% confidence intervals for comparisons with counts exceeding 5. Kaplan-Meier survival curves were generated to estimate survival from time of diagnosis to end of follow-up, and a log-rank test estimated differences across strata of interest. Subjects whose genetic profiles aligned with multiple categories (n = 3) were excluded from analyses stratified by genetic profile (n = 119). Otherwise, all 122 subjects were assessed. All statistical analyses were conducted in STATA 13.v1.
Results
Demographics
Overall, 122 subjects clinically diagnosed with HLH at Texas Children’s Hospital and referring centers from 1999 to 2016 were enrolled (Table 1). All subjects met HLH-2004 diagnostic criteria for HLH. Median age at diagnosis was 6.1 years (range, 2.1 weeks-18.5 years) with equal proportions of boys (n = 62) and girls (n = 60). Eight subjects had family histories of consanguinity. Six others denied consanguinity, but distant parental kinship could not be excluded after an increased number of rare homozygous variants were observed in the WES data.
Genetic testing workflow
Of the 122 subjects, a total of 101 received genetic testing. Clinical genetic testing was performed for 87 subjects, including 60 who were tested only by targeted Sanger sequencing, 15 who were assessed only by HLH gene panel testing, 11 who underwent clinical WES, and 1 who received only chromosomal microarray testing (Figures 1A; supplemental Figure 1). Research-based WES analyses were performed for 48 subjects. Genetic testing was not performed in 21 subjects for various reasons, including early death and clinical testing availability at the time of evaluation.
Genetic testing confirmed the presence of disease-causing fHLH gene mutations in a limited number of subjects
Nineteen of 101 subjects tested (19%) had biallelic mutations in 6 designated fHLH genes (Figure 1B). Clinical sequencing of individual fHLH genes identified 13 cases. Research-based WES of 20 of the remaining 48 subjects within this group revealed a rare, computationally predicted damaging homozygous RAB27A variant (supplemental Table 2) in a boy who presented with presumed malignancy-associated HLH and lacked the hair hypopigmentation associated with Griscelli syndrome.19,20 Clinical HLH gene panel testing identified 2 cases of Munc13-4 deficiency, and clinical WES successfully diagnosed 2 cases of perforin deficiency and 1 case of Munc13-4 deficiency.
Variants of unknown significance in fHLH genes were observed in 26 subjects who underwent research-based WES (supplemental Table 2). Most of these variants were either observed frequently in public and internal exome databases or not located in protein-encoding regions. The known pathogenic intronic UNC13D c.118-308C>T mutation11,21 was excluded by Sanger sequencing in the 3 subjects who had single heterozygous UNC13D variants (supplemental Table 2). Four subjects had clinical microarray testing performed to exclude copy number variations. Of the 26 subjects, 9 (HLH_014, HLH_016, HLH_017, HLH_063, HLH_087, HLH_097, HLH_109, HLH_117, and HLH_122) had rare and potentially protein sequence-altering single allelic variants of unknown significance in fHLH genes. Similarly, 5 subjects who received only clinical genetic testing carried solitary single heterozygous fHLH gene variants. The possibility that some of these monogenic variants may cause dominant-negative loss of function of the corresponding protein22 cannot be fully excluded.
Seven subjects who underwent research-based WES (HLH_014, HLH_016, HLH_017, HLH_063, HLH_087, HLH_097, and HLH_117) had single allelic variants in 2 fHLH genes that were potentially protein sequence-altering (supplemental Table 2). Computational analyses (supplemental Methods) found that heterozygous variants in the various “digenic” dual-gene combinations, except for one (LYST/RAB27A in subject HLH_117, who carries a disease-associated homozygous RAG1 mutation), were not uncommon (supplemental Table 3). Another subject received only clinical individual gene sequencing and was found to have single allelic variants in STX11 and UNC13D. This variant combination was associated with the HLH phenotype (P = .03). These analyses suggest that select combined variants23 (ie, mutational burden) in fHLH genes may be clinically relevant, but the majority of such variants are present in the general population. Thus, the potential digenic model for variant contribution to HLH lacked overall statistical support.
WES reveals the presence of other molecular diagnoses in subjects who meet HLH-2004 criteria
Excluding the subject discovered to have RAB27A deficiency, 47 subjects underwent research WES analyses and lacked biallelic mutations in fHLH genes, and potential underlying PIDDs were identified in 14 (30%; supplemental Table 4). Four of the 14 subjects had variants in genes relevant for cytolytic function in T and NK cells: WAS,24 TTC7A25 /LRBA,26 DOCK8,27 and CARMIL2.28-30 For the subject with compound heterozygous variants in both TTC7A and LRBA, it remains unclear whether defects in both of these genes are necessary for development of HLH. For DOCK8, deletion of the 5′ portion of the gene was discovered through bioinformatic detection of loss of WES reads and confirmed by array comparative genomic hybridization testing (see Stray-Pedersen et al.,15 Figure E4, Proband 30.1). We also identified biallelic mutations in RAG1 and RAG2 in 2 subjects who in hindsight had Omenn syndrome, which shares similar features with HLH.31 Two other subjects had underlying chronic granulomatous disease.31 Both subjects carried variants in second genes relevant toward inflammation or cytolytic immune function that may have further increased the risk for HLH. Four subjects had defects in key T or NK cell-signaling molecules: STAT1,32 STAT2,33 STAT3,34 and PIK3CD.35 WES identified another subject with biallelic variants in MCM3AP and MCM9. MCM3-associated protein deficiency has been reported to cause PIDD associated with a DNA repair defect.36 The precise mechanism for HLH in this subject remains unclear. Finally, 1 subject was discovered to have autoimmune lymphoproliferative syndrome type 2 because of a mutation in CASP10. Patients with autoimmune lymphoproliferative syndrome may present with features that overlap with HLH.37
Here, we introduce a new term, DIAP, or dysregulated immune activation or proliferation, to indicate a group of genes with putative or proven roles in dysregulated immune activation or proliferation. Of the 47 subjects who did not have biallelic mutations in fHLH genes and underwent research WES analyses, 11 (23%) had alterations in DIAP genes (supplemental Table 4). Three of the subjects (HLH_120, HLH_121, and HLH_122) also carried pathogenic or contributing variants in PIDD-associated genes. We bioinformatically tested the significance of differences in allelic counts between cases and controls for variants in these DIAP genes (supplemental Methods). This analysis demonstrated that the single allelic variants in NLRC438 and NLRP1239 and biallelic variants in NLRP4,40 NLRC3,41,42 and NLRP13 were significantly associated with HLH (supplemental Table 4; Figure 1C). The correlation between the NLRC4 variants and HLH corroborates the known link between NLRC4 mutations and HLH.38,43 These novel findings underscore a role for mutations in inflammation-associated genes in producing a phenotype that meets criteria for HLH. Finally, 1 subject was found to have a mutation in NRAS that is known to produce RAS-associated autoimmune leukoproliferative disorder.44 Because the variant exists de novo in this subject and all reported mutations in NRAS that cause RAS-associated autoimmune leukoproliferative disorder are somatic, this case is included in the DIAP and not the PIDD category.
Other potentially novel genetic explanations were identified in 5 subjects (10%; supplemental Table 4). Because none of these candidates have yet been fully biologically validated, assignment of causality remains uncertain. Nonetheless, they remain leading hypotheses because of their relevance in pathways known to be defective in HLH, including cytotoxicity, regulation of immune activation, and control of proliferation of hematopoietic cells.
Genetic findings are not enriched in subjects with abnormal CD107a mobilization or decreased expression of HLH-related proteins
HLH immunologic characterization was performed in 76 subjects: 33 were tested both for mobilization of CD107a and expression of HLH-related proteins (perforin, signaling lymphocytic activation molecule-associated protein, and/or baculoviral inhibitor of apoptosis repeat-containing protein 4), 37 were assessed for protein expression alone, and 6 received isolated CD107a functional testing. Genetic testing of 39 of the 43 individuals with abnormal results identified likely explanations in only 20 (51%). Although 10 fHLH probands were captured using this approach, subjects with PIDD (n = 3), DIAP (n = 4), and combined PIDD/DIAP (n = 3) were also discovered among the cases.
Correlations identified between genetic findings and clinical presentations
Genetic findings correlated with aspects of clinical presentation (Table 2 and Figure 2). Among the 19 subjects with biallelic fHLH genetic mutations, 58% had no identified trigger for HLH and 37% had associated infections; no fHLH subjects had rheumatologic triggers (Figure 2A). For the 11 subjects who had PIDD-associated gene variants and no variants in DIAP-associated genes, 36% had no identified trigger for HLH and 46% had associated infections (Figure 2A). Among the 8 subjects with DIAP-associated variants and no underlying PIDD, the majority (63%) had rheumatologic triggers (Figure 2A). The incidence of HLH-associated triggers did not vary significantly with age aside from perhaps rheumatologic triggers, which were suggested as increased in subjects older than 12 years compared with subjects younger than 1 year (36% vs 11%; P = .06; Figure 2B). Interestingly, however, subjects with biallelic fHLH mutations were significantly younger (P = .02) than subjects who lacked genetic findings, unlike subjects in the other genetically defined categories (Figure 2C; Table 2). In fact, subjects with fHLH were more likely than other subjects (P < .0001) to present younger than 1 year (Figure 2D). In contrast, all subjects who had underlying variants in DIAP-associated genes presented at older than 1 year, whereas subjects with underlying PIDD-associated gene defects presented with similar frequency at all ages.
Clinical outcomes are associated with genetic findings
Overall survival for the entire cohort was 45% with a median 18.6 months of follow-up (range, 0.04-201.4; Figure 3A). Overall survival estimates were not suggested to differ by HLH genetic profile (Figure 3B; supplemental Table 5). Individuals with fHLH were more likely to receive HSCT than other therapies compared with subjects with no genetic explanation (odds ratio, 7.33; 95% confidence interval, 2.35-22.89). Disease trigger was significantly associated with survival: subjects with rheumatologic disease or no obvious trigger exhibited better survival than subjects with HLH triggered by infection or malignancy (P < .05; Figure 3C; supplemental Table 6). All 9 subjects whose milder clinical course required only observation survived. No significant difference in survival was observed between subjects treated with biologics/steroids, immunochemotherapy, and HSCT, and no deaths occurred after 1 year in subjects treated with biologics/steroids (Table 3).
Discussion
Diagnostic considerations
These results demonstrate that patients clinically diagnosed with HLH using the HLH-2004 criteria represent a distinct category of patients with persistent, toxic immune activation at high risk for death. We found that 80 (66%) of 122 subjects with HLH had likely or potential disease-causing genetic, autoimmune, or oncologic etiologies underlying the phenotype. Importantly, 55 subjects received genetic testing and lacked biallelic fHLH gene mutations or other genetic explanations, yet 35 (64%) were not evaluated by WES. For 48 subjects who underwent research WES analyses, potential disease-causing genetic defects were identified in 28 (58%).
In this ethnically diverse cohort, although the HLH-2004 criteria did not demonstrate the high prevalence of fHLH exhibited in the European studies,11-14 they identified patients with HLH resulting from other genetic conditions. Overall, genetic testing of 101 subjects with HLH secured likely molecular diagnoses in 46 (46%; or 70% excluding the 35 subjects not tested by WES). Biallelic fHLH gene mutations were identified in only 19% of 101 subjects. This result may reflect the fact that genetic testing was performed on the basis of overall fulfillment of HLH-2004 criteria and not simply focused toward subjects with abnormal expression of HLH-related proteins or defective cytolytic function.11-14 Furthermore, modeling of mutational burden in fHLH genes found that most potential digenic combinations were unlikely to arise with greater frequency in HLH than in the general population. Clinical laboratory results reporting digenic variants should therefore be interpreted with caution. Meanwhile, variants associated with PIDDs and DIAP were found in 22% of 101 cases (33% excluding the 35 subjects not tested by WES). Thus, in this cohort, HLH-2004 criteria identified more cases of potential PIDD or DIAP than fHLH. This observation emphasizes the need to conceptualize HLH as a phenotype rather than a single disease. For example, in patients with PIDD, impaired ability to eliminate pathogens can drive formation of the HLH phenotype. This phenomenon has been clearly reported by others.31 In contrast, the HLH phenotype may arise in patients with DIAP because of exaggerated or poorly controlled immunologic responses. Neither of these conditions necessitates defective cytotoxicity. Thus, targeted sequencing or testing of cytolytic function may not be sufficient for identifying inherited causes of HLH.
An emerging concept creates a spectrum of genotype and phenotype rather than dichotomous primary and secondary categories.4 Our results argue that even patients with presumed secondary HLH may benefit from investigations for underlying genetic conditions. Among our subjects with presumed secondary HLH, disease-associated variants in fHLH, PIDD, and DIAP genes were identified in 22% (7 of 32) of subjects with autoimmune disease and 21% (3 of 14) of subjects with malignancy.
Our results also suggest that addition of an age parameter to the HLH-2004 criteria may increase the specificity for fHLH. We observed that fHLH was statistically enriched in subjects who were diagnosed at younger than 1 year of age. This phenomenon is consistent with findings reported by others.12,13,45 Meanwhile, the majority of subjects who had underlying variants in DIAP-associated genes were older. Nevertheless, age ranges for these groups overlap.
Therapeutic implications
In this cohort, overall survival was 45%, which is comparable to other reports.5 The management of HLH at this time is dictated by the clinical presentation and course of each patient. Patients with minimal or short-lived symptoms typically receive only observation or minimally toxic therapies. The current approach to patients who require immunochemotherapy involves initiation of etoposide and dexamethasone (for patients lacking co-incident autoimmune disease).7 These agents, however, are not feasible for long-term immunosuppression because of cumulative toxicity.46 In general, HSCT is reserved for patients presenting with HLH because of fHLH gene defects, exhibiting persistent or progressive pathologic inflammation after HLH-94 induction, demonstrating persistent cytolytic immune function defects, or having central nervous system inflammation.7 Thus, patients who exhibit greater severity of disease and likelihood for nonsurvival receive more intense therapies that convey increased risk for poor outcomes (Table 3).
The current practice of relying on failure of therapy to dictate escalation of treatment places patients at risk for increased morbidity and mortality. Although our results suggest that WES may be deferred in patients who are managed by observation alone, this study also supports the use of WES to guide long-term clinical decision-making for patients with clinically diagnosed HLH on initiation of therapy. First, these data revealed multiple potential molecular pathways leading to HLH. Early identification of these altered pathways may offer alternative therapeutic opportunities, such as targeted biologic therapy. Second, WES can confirm the presence of fHLH or an underlying PIDD for which HSCT is indicated. Finally, WES can identify unanticipated inherited gene defects that may be shared by potential allograft donors.
Genetic testing in children with HLH
Our genetic findings reflect a very diverse population at a large referral center and differ from results elsewhere that suggest focused use of genetic testing toward subjects who demonstrate decreased mobilization of CD107a or expression of HLH-related proteins.14 Although such an approach led to genetic diagnoses in 85% of subjects in 1 cohort, we found that only 51% of individuals with such laboratory abnormalities in our population had likely genetic explanations. The discrepancy may be partially attributed to differences in study populations and design, as 46 of the subjects in this cohort, including 6 subjects with fHLH, did not receive testing of CD107a mobilization or expression of HLH-related proteins. Nonetheless, the defective genes in our cohort extended well beyond PRF1, UNC13D, STX11, STXBP2, LYST, RAB27A, AP3B1, SH2D1A, and XIAP. Furthermore, such a focused approach would have failed to capture the subject with CASP10 deficiency and another subject with a potential novel genetic diagnosis, as both probands had normal CD107a function and HLH-related protein expression.
We propose a broader approach for use of WES, if rapidly available, in patients with HLH (Figure 4). This algorithm combines our findings with experiences reported by others.12,14,47 First, HLH must be considered in ill children, and the HLH-2004 diagnostic criteria should be interrogated accordingly. If sufficient criteria are not fulfilled, clinical judgment must be exercised to determine which individuals should be observed closely for potential progression such that the criteria eventually become met. Children who satisfy 5 or more of the 8 criteria should receive immunologic laboratory testing for cytolytic function and expression of HLH-related proteins. If these tests demonstrate abnormalities, targeted gene sequencing should be performed, especially if the patient is younger than 1 year. We argue that WES should follow if no genetic diagnosis has been secured. Children with HLH who have normal or unhelpful immunologic laboratory testing results and require minimal therapy may be observed without genetic testing, although targeted gene sequencing or WES should be considered. We did not find any cases of familial HLH, PIDD, or DIAP in this small subgroup of individuals, but further studies will be needed to determine whether WES should be used to search for other underlying genetic conditions. In contrast, patients with persistent or recurrent HLH or who require initiation of steroids, biologics, or more should receive WES.
It remains important to note that although WES currently serves as the most expedient broad genetic test in terms of cost, analytic difficulty, turnaround time, and availability, other techniques may eventually supplant this method because of several of its limitations. WES is known to miss structural variations, and despite enhanced bioinformatic tools,15,48 identification of gross copy number variations remains difficult. In addition, WES is not designed to interrogate intronic variants that may affect either promoter or enhancer elements, resulting in altered protein expression, or introduce cryptic start codons or splicing sites. These abilities are better captured by whole-genome sequencing and RNA sequencing, which may ultimately become preferred diagnostic tests in the future.
In conclusion, in this study the HLH-2004 criteria identified a cohort of critically ill children with a broad spectrum of genetic conditions, autoimmune disease, and malignancy that result in the HLH phenotype. Clinical suspicion for HLH and application of the HLH-2004 criteria remain the most critical steps in identifying patients with pathologic immune activation. After clinical identification, subsequent evaluations should include studies to identify inherited molecular defects that underlie inability to regulate immunologic homeostasis. Additional prospective trials enrolling at time of suspicion of HLH will prove helpful for determining therapeutic strategies based on pathogenic molecular mechanisms.
The online version of this article contains a data supplement.
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.
Acknowledgments
This work was funded in part by grants from the HistioCure Foundation (K.L.M. and C.E.A.), St. Baldrick’s Foundation (to NACHO Consortium, K.L.M., and C.E.A.), National Institutes of Health, National Cancer Institute (R25CA160078: Training Program in Pediatric Cancer Epidemiology and Control Grant to E.C.P.-G.), American Society of Hematology (Scholar Award in Clinical Research to E.C.P.-G.), the Thrasher Research Fund (Early Career Award to E.C.P.-G.), and the Feldman Family. This work was also supported by the following: National Institutes of Health, National Institute of Allergy and Infectious Diseases grants R01AI067946 and R01AI120989 (J.S.O.), National Institutes of Health, National Human Genome Research Institute/National Heart, Lung, and Blood Institute grant UM1HG006542 (to the Baylor-Hopkins Center for Mendelian Genomics), National Institutes of Health, National Institute of Neurological Disorders and Stroke grant R01NS058529 (J.R.L.), and the Jeffrey Modell Foundation Translational Research Program Grant (I.K.C.).
Authorship
Contribution: I.K.C. analyzed data and wrote the manuscript; O.S.E. and B.R.G. analyzed clinical data; E.C.P.-G. and Z.H.C.-A. performed statistical and bioinformatic analyses, respectively, and wrote portions of the manuscript; L.R.F., S.K.N., H.A.A., M.I.D., H.E.H., R.A.K., C.A.M., T.C.N., D.A.B., J.R.G., L.A.P., M.C.P., J.C.A.-B., S.A.M., W.A.-H., A.C., T.N.C., and D.N.H. collected samples and clinical or research data; E.M.M., T.P.V., and A.S.-P. analyzed genetic data; S.N.J., D.M.M., and R.A.G. performed WES; J.R.L., J.S.O., and K.L.M. provided expertise and supervision of the research; and C.E.A. supervised the research study and wrote portions of the manuscript.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Carl E. Allen, Feigin Center, Suite 730.06, 1102 Bates St, Houston, TX 77030; e-mail: ceallen@txch.org; and Ivan K. Chinn, Feigin Center, Suite 330, 1102 Bates St, Houston, TX 77030; e-mail: chinn@bcm.edu.