Proteogenomic analysis integrated with electronic health records data reveals disease-associated variants in Black Americans.

BACKGROUND: Most genome wide association studies (GWAS) of plasma proteomics have focused on White individuals of European ancestry, limiting biological insight from other ancestry enriched protein quantitative loci (pQTL). METHODS: We conducted a discovery GWAS of ~3,000 plasma proteins measured by the antibody based Olink platform in 1,054 Black adults from the Jackson Heart Study (JHS), and validated our findings in the Multi-Ethnic Study of Atherosclerosis (MESA). The genetic architecture of identified pQTLs were further explored through fine mapping and admixture association analysis. Finally, using our pQTL findings, we performed a phenome wide association study (PheWAS) across two large multi-ethnic electronic health record (EHR) systems in All of Us and BioMe. RESULTS: We identified 1002 pQTLs for 925 proteins. Fine mapping and admixture analyses suggested allelic heterogeneity of the plasma proteome across diverse populations. We identified associations for variants enriched in African ancestry, many in diseases that lack precise biomarkers, including cis-pQTLs for Cathepsin L (CTSL) and Siglec-9 that were linked with sarcoidosis and non-Hodgkin's lymphoma, respectively. We found concordant associations across clinical diagnoses and laboratory measurements, elucidating disease pathways, including a cis-pQTL associated with circulating CD58, white blood cell count, and multiple sclerosis. CONCLUSIONS: Our findings emphasize the value of leveraging diverse populations to enhance biological insights from proteomics GWAS, and we have made this resource readily available as an interactive web portal.

Inherited human RelB deficiency impairs innate and adaptive immunity to infection.

We report two unrelated adults with homozygous (P1) or compound heterozygous (P2) private loss-of-function variants of V-Rel Reticuloendotheliosis Viral Oncogene Homolog B (RELB). The resulting deficiency of functional RelB impairs the induction of NFKB2 mRNA and NF-κB2 (p100/p52) protein by lymphotoxin in the fibroblasts of the patients. These defects are rescued by transduction with wild-type RELB complementary DNA (cDNA). By contrast, the response of RelB-deficient fibroblasts to Tumor Necrosis Factor (TNF) or IL-1ß via the canonical NF-κB pathway remains intact. P1 and P2 have low proportions of naïve CD4+ and CD8+ T cells and of memory B cells. Moreover, their naïve B cells cannot differentiate into immunoglobulin G (IgG)- or immunoglobulin A (IgA)-secreting cells in response to CD40L/IL-21, and the development of IL-17A/F-producing T cells is strongly impaired in vitro. Finally, the patients produce neutralizing autoantibodies against type I interferons (IFNs), even after hematopoietic stem cell transplantation, attesting to a persistent dysfunction of thymic epithelial cells in T cell selection and central tolerance to some autoantigens. Thus, inherited human RelB deficiency disrupts the alternative NF-κB pathway, underlying a T- and B cell immunodeficiency, which, together with neutralizing autoantibodies against type I IFNs, confers a predisposition to viral, bacterial, and fungal infections.

IRS2 Signaling Protects Against Stress-Induced Arrhythmia by Maintaining Ca2+ Homeostasis.

BACKGROUND: The docking protein IRS2 (insulin receptor substrate protein-2) is an important mediator of insulin signaling and may also regulate other signaling pathways. Murine hearts with cardiomyocyte-restricted deletion of IRS2 (cIRS2-KO) are more susceptible to pressure overload-induced cardiac dysfunction, implying a critical protective role of IRS2 in cardiac adaptation to stress through mechanisms that are not fully understood. There is limited evidence regarding the function of IRS2 beyond metabolic homeostasis regulation, particularly in the context of cardiac disease. METHODS: A retrospective analysis of an electronic medical record database was conducted to identify patients with IRS2 variants and assess their risk of cardiac arrhythmias. Arrhythmia susceptibility was examined in cIRS2-KO mice. The underlying mechanisms were investigated using confocal calcium imaging of ex vivo whole hearts and isolated cardiomyocytes to assess calcium handling, Western blotting to analyze the involved signaling pathways, and pharmacological and genetic interventions to rescue arrhythmias in cIRS2-KO mice. RESULTS: The retrospective analysis identified patients with IRS2 variants of uncertain significance with a potential association to an increased risk of cardiac arrhythmias compared with matched controls. cIRS2-KO hearts were found to be prone to catecholamine-sensitive ventricular tachycardia and reperfusion ventricular tachycardia. Confocal calcium imaging of ex vivo whole hearts and single isolated cardiomyocytes from cIRS2-KO hearts revealed decreased Ca²+ transient amplitudes, increased spontaneous Ca²+ sparks, and reduced sarcoplasmic reticulum Ca²+ content during sympathetic stress, indicating sarcoplasmic reticulum dysfunction. We identified that overactivation of the AKT1/NOS3 (nitric oxide synthase 3)/CaMKII (Ca2+/calmodulin-dependent protein kinase II)/RyR2 (type 2 ryanodine receptor) signaling pathway led to calcium mishandling and catecholamine-sensitive ventricular tachycardia in cIRS2-KO hearts. Pharmacological AKT inhibition or genetic stabilization of RyR2 rescued catecholamine-sensitive ventricular tachycardia in cIRS2-KO mice. CONCLUSIONS: Cardiac IRS2 inhibits sympathetic stress-induced AKT/NOS3/CaMKII/RyR2 overactivation and calcium-dependent arrhythmogenesis. This novel IRS2 signaling axis, essential for maintaining cardiac calcium homeostasis under stress, presents a promising target for developing new antiarrhythmic therapies.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency.

G6PC3 deficiency is a monogenic immunometabolic disorder that causes syndromic congenital neutropenia. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Here, we investigated the origin and functional consequence of the G6PC3 c.210delC variant found in patients of Mexican origin. Based on the shared haplotypes amongst carriers of the c.210delC mutation, we estimated that this variant originated from a founder effect in a common ancestor. Furthermore, by ancestry analysis, we concluded that it originated in the indigenous Mexican population. At the protein level, we showed that this frameshift mutation leads to an aberrant protein expression in overexpression and patient-derived cells. G6PC3 pathology is driven by the intracellular accumulation of the metabolite 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. We characterized how the variant c.210delC impacts glycolysis by performing extracellular flux assays on patient-derived cells. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Finally, we compared the clinical presentation of patients with the mutation c.210delC and all other G6PC3 deficient patients reported in the literature to date, and we found that c.210delC carriers display all prominent clinical features observed in prior G6PC3 deficient patients. In conclusion, G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

The landscape of rare genetic variation associated with inflammatory bowel disease and Parkinson's disease comorbidity.

BACKGROUND: Inflammatory bowel disease (IBD) and Parkinson's disease (PD) are chronic disorders that have been suggested to share common pathophysiological processes. LRRK2 has been implicated as playing a role in both diseases. Exploring the genetic basis of the IBD-PD comorbidity through studying high-impact rare genetic variants can facilitate the identification of the novel shared genetic factors underlying this comorbidity. METHODS: We analyzed whole exomes from the BioMe BioBank and UK Biobank, and whole genomes from a cohort of 67 European patients diagnosed with both IBD and PD to examine the effects of LRRK2 missense variants on IBD, PD and their co-occurrence (IBD-PD). We performed optimized sequence kernel association test (SKAT-O) and network-based heterogeneity clustering (NHC) analyses using high-impact rare variants in the IBD-PD cohort to identify novel candidate genes, which we further prioritized by biological relatedness approaches. We conducted phenome-wide association studies (PheWAS) employing BioMe BioBank and UK Biobank whole exomes to estimate the genetic relevance of the 14 prioritized genes to IBD-PD. RESULTS: The analysis of LRRK2 missense variants revealed significant associations of the G2019S and N2081D variants with IBD-PD in addition to several other variants as potential contributors to increased or decreased IBD-PD risk. SKAT-O identified two significant genes, LRRK2 and IL10RA, and NHC identified 6 significant gene clusters that are biologically relevant to IBD-PD. We observed prominent overlaps between the enriched pathways in the known IBD, PD, and candidate IBD-PD gene sets. Additionally, we detected significantly enriched pathways unique to the IBD-PD, including MAPK signaling, LPS/IL-1 mediated inhibition of RXR function, and NAD signaling. Fourteen final candidate IBD-PD genes were prioritized by biological relatedness methods. The biological importance scores estimated by protein-protein interaction networks and pathway and ontology enrichment analyses indicated the involvement of genes related to immunity, inflammation, and autophagy in IBD-PD. Additionally, PheWAS provided support for the associations of candidate genes with IBD and PD. CONCLUSIONS: Our study confirms and uncovers new LRRK2 associations in IBD-PD. The identification of novel inflammation and autophagy-related genes supports and expands previous findings related to IBD-PD pathogenesis, and underscores the significance of therapeutic interventions for reducing systemic inflammation.

Molecular and clinical characterization of a founder mutation causing G6PC3 deficiency.

Background: G6PC3 deficiency is a rare genetic disorder that causes syndromic congenital neutropenia. It is driven by the intracellular accumulation of a metabolite named 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Objective: The G6PC3 c.210delC variant has been identified in patients of Mexican origin. We set out to study the origin and functional consequence of this mutation. Furthermore, we sought to characterize the clinical phenotypes caused by it. Methods: Using whole-genome sequencing data, we conducted haplotype analysis to estimate the age of this allele and traced its ancestral origin. We examined how this mutation affected G6PC3 protein expression and performed extracellular flux assays on patient-derived cells to characterize how this mutation impacts glycolysis. Finally, we compared the clinical presentations of patients with the c.210delC mutation relative to other G6PC3 deficient patients published to date. Results: Based on the length of haplotypes shared amongst ten carriers of the G6PC3 c.210delC mutation, we estimated that this variant originated in a common ancestor of indigenous American origin. The mutation causes a frameshift that introduces a premature stop codon, leading to a complete loss of G6PC3 protein expression. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient-derived cells exhibited markedly reduced engagement of glycolysis. Clinically, c.210delC carriers display all the clinical features of syndromic severe congenital neutropenia type 4 observed in prior reports of G6PC3 deficiency. Conclusion: The G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Individuals with JAK1 variants are affected by syndromic features encompassing autoimmunity, atopy, colitis, and dermatitis.

Inborn errors of immunity lead to autoimmunity, inflammation, allergy, infection, and/or malignancy. Disease-causing JAK1 gain-of-function (GoF) mutations are considered exceedingly rare and have been identified in only four families. Here, we use forward and reverse genetics to identify 59 individuals harboring one of four heterozygous JAK1 variants. In vitro and ex vivo analysis of these variants revealed hyperactive baseline and cytokine-induced STAT phosphorylation and interferon-stimulated gene (ISG) levels compared with wild-type JAK1. A systematic review of electronic health records from the BioME Biobank revealed increased likelihood of clinical presentation with autoimmunity, atopy, colitis, and/or dermatitis in JAK1 variant-positive individuals. Finally, treatment of one affected patient with severe atopic dermatitis using the JAK1/JAK2-selective inhibitor, baricitinib, resulted in clinically significant improvement. These findings suggest that individually rare JAK1 GoF variants may underlie an emerging syndrome with more common presentations of autoimmune and inflammatory disease (JAACD syndrome). More broadly, individuals who present with such conditions may benefit from genetic testing for the presence of JAK1 GoF variants.

Development of a human genetics-guided priority score for 19,365 genes and 399 drug indications.

Studies have shown that drug targets with human genetic support are more likely to succeed in clinical trials. Hence, a tool integrating genetic evidence to prioritize drug target genes is beneficial for drug discovery. We built a genetic priority score (GPS) by integrating eight genetic features with drug indications from the Open Targets and SIDER databases. The top 0.83%, 0.28% and 0.19% of the GPS conferred a 5.3-, 9.9- and 11.0-fold increased effect of having an indication, respectively. In addition, we observed that targets in the top 0.28% of the score were 1.7-, 3.7- and 8.8-fold more likely to advance from phase I to phases II, III and IV, respectively. Complementary to the GPS, we incorporated the direction of genetic effect and drug mechanism into a directional version of the score called the GPS with direction of effect. We applied our method to 19,365 protein-coding genes and 399 drug indications and made all results available through a web portal.

Genome-wide prediction of pathogenic gain- and loss-of-function variants from ensemble learning of a diverse feature set.

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available at https://itanlab.shinyapps.io/goflof/ .

Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients.

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

Genetic susceptibility to diabetic kidney disease is linked to promoter variants of XOR.

The lifetime risk of kidney disease in people with diabetes is 10-30%, implicating genetic predisposition in the cause of diabetic kidney disease (DKD). Here we identify an expression quantitative trait loci (QTLs) in the cis-acting regulatory region of the xanthine dehydrogenase, or xanthine oxidoreductase (Xor), a binding site for C/EBPß, to be associated with diabetes-induced podocyte loss in DKD in male mice. We examine mouse inbred strains that are susceptible (DBA/2J) and resistant (C57BL/6J) to DKD, as well as a panel of recombinant inbred BXD mice, to map QTLs. We also uncover promoter XOR orthologue variants in humans associated with high risk of DKD. We introduced the risk variant into the 5'-regulatory region of XOR in DKD-resistant mice, which resulted in increased Xor activity associated with podocyte depletion, albuminuria, oxidative stress and damage restricted to the glomerular endothelium, which increase further with type 1 diabetes, high-fat diet and ageing. Therefore, differential regulation of Xor contributes to phenotypic consequences with diabetes and ageing.

Identifying shared genetic factors underlying epilepsy and congenital heart disease in Europeans.

Epilepsy (EP) and congenital heart disease (CHD) are two apparently unrelated diseases that nevertheless display substantial mutual comorbidity. Thus, while congenital heart defects are associated with an elevated risk of developing epilepsy, the incidence of epilepsy in CHD patients correlates with CHD severity. Although genetic determinants have been postulated to underlie the comorbidity of EP and CHD, the precise genetic etiology is unknown. We performed variant and gene association analyses on EP and CHD patients separately, using whole exomes of genetically identified Europeans from the UK Biobank and Mount Sinai BioMe Biobank. We prioritized biologically plausible candidate genes and investigated the enriched pathways and other identified comorbidities by biological proximity calculation, pathway analyses, and gene-level phenome-wide association studies. Our variant- and gene-level results point to the Voltage-Gated Calcium Channels (VGCC) pathway as being a unifying framework for EP and CHD comorbidity. Additionally, pathway-level analyses indicated that the functions of disease-associated genes partially overlap between the two disease entities. Finally, phenome-wide association analyses of prioritized candidate genes revealed that cerebral blood flow and ulcerative colitis constitute the two main traits associated with both EP and CHD.

Novel brown adipose tissue candidate genes predicted by the human gene connectome.

Brown adipose tissue (BAT) is a promising therapeutic target against obesity. Therefore, research on the genetic architecture of BAT could be key for the development of successful therapies against this complex phenotype. Hypothesis-driven candidate gene association studies are useful for studying genetic determinants of complex traits, but they are dependent upon the previous knowledge to select candidate genes. Here, we predicted 107 novel-BAT candidate genes in silico using the uncoupling protein one (UCP1) as the hallmark of BAT activity. We first identified the top 1% of human genes predicted by the human gene connectome to be biologically closest to the UCP1, estimating 167 additional pathway genes (BAT connectome). We validated this prediction by showing that 60 genes already associated with BAT were included in the connectome and they were biologically closer to each other than expected by chance (p < 2.2 × 10-16). The rest of genes (107) are potential candidates for BAT, being also closer to known BAT genes and more expressed in BAT biopsies than expected by chance (p < 2.2 × 10-16; p = 4.39 × 10-02). The resulting new list of predicted human BAT genes should be useful for the discovery of novel BAT genes and metabolic pathways.

VIPPID: a gene-specific single nucleotide variant pathogenicity prediction tool for primary immunodeficiency diseases.

Distinguishing pathogenic variants from non-pathogenic ones remains a major challenge in clinical genetic testing of primary immunodeficiency (PID) patients. Most of the existing mutation pathogenicity prediction tools treat all mutations as homogeneous entities, ignoring the differences in characteristics of different genes, and use the same model for genes in different diseases. In this study, we developed a single nucleotide variant (SNV) pathogenicity prediction tool, Variant Impact Predictor for PIDs (VIPPID; https://mylab.shinyapps.io/VIPPID/), which was tailored for PIDs genes and used a specific model for each of the most prevalent PID known genes. It employed a Conditional Inference Forest model and utilized information of 85 features of SNVs and scores from 20 existing prediction tools. Evaluation of VIPPID showed that it had superior performance (area under the curve = 0.91) over non-specific conventional tools. In addition, we also showed that the gene-specific model outperformed the non-gene-specific models. Our study demonstrated that disease-specific and gene-specific models can improve SNV pathogenicity prediction performance. This observation supports the notion that each feature of mutations in the model can be potentially used, in a new algorithm, to investigate the characteristics and function of the encoded proteins.

Identification of discriminative gene-level and protein-level features associated with pathogenic gain-of-function and loss-of-function variants.

Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

The genetic structure of the Turkish population reveals high levels of variation and admixture.

The construction of population-based variomes has contributed substantially to our understanding of the genetic basis of human inherited disease. Here, we investigated the genetic structure of Turkey from 3,362 unrelated subjects whose whole exomes (n = 2,589) or whole genomes (n = 773) were sequenced to generate a Turkish (TR) Variome that should serve to facilitate disease gene discovery in Turkey. Consistent with the history of present-day Turkey as a crossroads between Europe and Asia, we found extensive admixture between Balkan, Caucasus, Middle Eastern, and European populations with a closer genetic relationship of the TR population to Europeans than hitherto appreciated. We determined that 50% of TR individuals had high inbreeding coefficients (≥0.0156) with runs of homozygosity longer than 4 Mb being found exclusively in the TR population when compared to 1000 Genomes Project populations. We also found that 28% of exome and 49% of genome variants in the very rare range (allele frequency < 0.005) are unique to the modern TR population. We annotated these variants based on their functional consequences to establish a TR Variome containing alleles of potential medical relevance, a repository of homozygous loss-of-function variants and a TR reference panel for genotype imputation using high-quality haplotypes, to facilitate genome-wide association studies. In addition to providing information on the genetic structure of the modern TR population, these data provide an invaluable resource for future studies to identify variants that are associated with specific phenotypes as well as establishing the phenotypic consequences of mutations in specific genes.

Sudden Cardiac Arrest in a Patient With Mitral Valve Prolapse and LMNA and SCN5A Mutations.

Bileaflet mitral valve prolapse (Bi-MVP) is associated with increased risk for cardiac arrest. We describe a patient who presented after a cardiac arrest with Bi-MVP and variants in Lamin A/C (LMNA) and the sodium channel alpha-subunit 5a (SCN5A). Genetic variants may be the culprit for arrhythmogenesis in Bi-MVP patients. (Level of Difficulty: Intermediate.).

Burden of Cardiomyopathic Genetic Variation in Lethal Pediatric Myocarditis.

BACKGROUND: Acute myocarditis (AM) is a well-known cause of sudden death and heart failure, often caused by prevalent viruses. We previously showed that some pediatric AM correlates with putatively damaging variants in genes related to cardiomyocyte structure and function. We sought to evaluate whether deleterious cardiomyopathic variants were enriched among fatal pediatric AM cases in New York City compared with ancestry-matched controls. METHODS: Twenty-four children (aged 3 weeks to 20 years) with death due to AM were identified through autopsy records; histologies were reviewed to confirm that all cases met Dallas criteria for AM and targeted panel sequencing of 57 cardiomyopathic genes was performed. Controls without cardiovascular disease were identified from a pediatric database and matched by genetic ancestry to cases using principal components from exome sequencing. Rates of putative deleterious variations (DV) were compared between cases and controls. Where available, AM tissues underwent viral analysis by polymerase chain reaction. RESULTS: DV were identified in 4 of 24 AM cases (16.7%), compared with 2 of 96 age and ancestry-matched controls (2.1%, P=0.014). Viral causes were proven for 6 of 8 AM cases (75%), including the one DV+ case where tissue was available for testing. DV+ cases were more likely to be female, have no evidence of chronic inflammation, and associate with sudden cardiac death than DV- cases. CONCLUSIONS: Deleterious variants in genes related to cardiomyocyte integrity are more common in children with fatal AM than controls, likely conferring susceptibility. Additionally, genetically mediated AM may progress more rapidly and be more severe.

A computational approach for detecting physiological homogeneity in the midst of genetic heterogeneity.

The human genetic dissection of clinical phenotypes is complicated by genetic heterogeneity. Gene burden approaches that detect genetic signals in case-control studies are underpowered in genetically heterogeneous cohorts. We therefore developed a genome-wide computational method, network-based heterogeneity clustering (NHC), to detect physiological homogeneity in the midst of genetic heterogeneity. Simulation studies showed our method to be capable of systematically converging genes in biological proximity on the background biological interaction network, and capturing gene clusters harboring presumably deleterious variants, in an efficient and unbiased manner. We applied NHC to whole-exome sequencing data from a cohort of 122 individuals with herpes simplex encephalitis (HSE), including 13 individuals with previously published monogenic inborn errors of TLR3-dependent IFN-α/ß immunity. The top gene cluster identified by our approach successfully detected and prioritized all causal variants of five TLR3 pathway genes in the 13 previously reported individuals. This approach also suggested candidate variants of three reported genes and four candidate genes from the same pathway in another ten previously unstudied individuals. TLR3 responsiveness was impaired in dermal fibroblasts from four of the five individuals tested, suggesting that the variants detected were causal for HSE. NHC is, therefore, an effective and unbiased approach for unraveling genetic heterogeneity by detecting physiological homogeneity.