Inborn Errors of Immunity Contribute to the Burden of Skin Disease and Create Opportunities for Improving the Practice of Dermatology.

Opportunities to improve the clinical management of skin disease are being created by advances in genomic medicine. Large-scale sequencing increasingly challenges notions about single-gene disorders. It is now apparent that monogenic etiologies make appreciable contributions to the population burden of disease and that they are underrecognized in clinical practice. A genetic diagnosis informs on molecular pathology and may direct targeted treatments and tailored prevention strategies for patients and family members. It also generates knowledge about disease pathogenesis and management that is relevant to patients without rare pathogenic variants. Inborn errors of immunity are a large class of monogenic etiologies that have been well-studied and contribute to the population burden of inflammatory diseases. To further delineate the contributions of inborn errors of immunity to the pathogenesis of skin disease, we performed a set of analyses that identified 316 inborn errors of immunity associated with skin pathologies, including common skin diseases. These data suggest that clinical sequencing is underutilized in dermatology. We next use these data to derive a network that illuminates the molecular relationships of these disorders and suggests an underlying etiological organization to immune-mediated skin disease. Our results motivate the further development of a molecularly derived and data-driven reorganization of clinical diagnoses of skin disease.

Participant-guided development of bilingual genomic educational infographics for Electronic Medical Records and Genomics Phase IV study.

OBJECTIVE: Developing targeted, culturally competent educational materials is critical for participant understanding of engagement in a large genomic study that uses computational pipelines to produce genome-informed risk assessments. MATERIALS AND METHODS: Guided by the Smerecnik framework that theorizes understanding of multifactorial genetic disease through 3 knowledge types, we developed English and Spanish infographics for individuals enrolled in the Electronic Medical Records and Genomics Network. Infographics were developed to explain concepts in lay language and visualizations. We conducted iterative sessions using a modified "think-aloud" process with 10 participants (6 English, 4 Spanish-speaking) to explore comprehension of and attitudes towards the infographics. RESULTS: We found that all but one participant had "awareness knowledge" of genetic disease risk factors upon viewing the infographics. Many participants had difficulty with "how-to" knowledge of applying genetic risk factors to specific monogenic and polygenic risks. Participant attitudes towards the iteratively-refined infographics indicated that design saturation was reached. DISCUSSION: There were several elements that contributed to the participants' comprehension (or misunderstanding) of the infographics. Visualization and iconography techniques best resonated with those who could draw on prior experiences or knowledge and were absent in those without. Limited graphicacy interfered with the understanding of absolute and relative risks when presented in graph format. Notably, narrative and storytelling theory that informed the creation of a vignette infographic was most accessible to all participants. CONCLUSION: Engagement with the intended audience who can identify strengths and points for improvement of the intervention is necessary to the development of effective infographics.

Inborn Errors of Immunity in Hidradenitis Suppurativa Pathogenesis and Disease Burden.

Hidradenitis suppurativa (HS), also known as Verneuil's disease and acne inversa, is a prevalent, debilitating, and understudied inflammatory skin disease. It is marked by repeated bouts of pathological inflammation causing pain, hyperplasia, aberrant healing, and fibrosis. HS is difficult to manage and has many unmet medical needs. There is clinical and pharmacological evidence for extensive etiological heterogeneity with HS, suggesting that this clinical diagnosis is capturing a spectrum of disease entities. Human genetic studies provide robust insight into disease pathogenesis. They also can be used to resolve etiological heterogeneity and to identify drug targets. However, HS has not been extensively investigated with well-powered genetic studies. Here, we review what is known about its genetic architecture. We identify overlap in molecular, cellular, and clinical features between HS and inborn errors of immunity (IEI). This evidence indicates that HS may be an underrecognized component of IEI and suggests that undiagnosed IEI are present in HS cohorts. Inborn errors of immunity represent a salient opportunity for rapidly resolving the immunological landscape of HS pathogenesis, for prioritizing drug repurposing studies, and for improving the clinical management of HS.

Quantifying the phenome-wide disease burden of obesity using electronic health records and genomics.

OBJECTIVE: High BMI is associated with many comorbidities and mortality. This study aimed to elucidate the overall clinical risk of obesity using a genome- and phenome-wide approach. METHODS: This study performed a phenome-wide association study of BMI using a clinical cohort of 736,726 adults. This was followed by genetic association studies using two separate cohorts: one consisting of 65,174 adults in the Electronic Medical Records and Genomics (eMERGE) Network and another with 405,432 participants in the UK Biobank. RESULTS: Class 3 obesity was associated with 433 phenotypes, representing 59.3% of all billing codes in individuals with severe obesity. A genome-wide polygenic risk score for BMI, accounting for 7.5% of variance in BMI, was associated with 296 clinical diseases, including strong associations with type 2 diabetes, sleep apnea, hypertension, and chronic liver disease. In all three cohorts, 199 phenotypes were associated with class 3 obesity and polygenic risk for obesity, including novel associations such as increased risk of renal failure, venous insufficiency, and gastroesophageal reflux. CONCLUSIONS: This combined genomic and phenomic systematic approach demonstrated that obesity has a strong genetic predisposition and is associated with a considerable burden of disease across all disease classes.

Whole exome sequencing in Alopecia Areata identifies rare variants in KRT82.

Alopecia areata is a complex genetic disease that results in hair loss due to the autoimmune-mediated attack of the hair follicle. We previously defined a role for both rare and common variants in our earlier GWAS and linkage studies. Here, we identify rare variants contributing to Alopecia Areata using a whole exome sequencing and gene-level burden analyses approach on 849 Alopecia Areata patients compared to 15,640 controls. KRT82 is identified as an Alopecia Areata risk gene with rare damaging variants in 51 heterozygous Alopecia Areata individuals (6.01%), achieving genome-wide significance (p = 2.18E-07). KRT82 encodes a hair-specific type II keratin that is exclusively expressed in the hair shaft cuticle during anagen phase, and its expression is decreased in Alopecia Areata patient skin and hair follicles. Finally, we find that cases with an identified damaging KRT82 variant and reduced KRT82 expression have elevated perifollicular CD8 infiltrates. In this work, we utilize whole exome sequencing to successfully identify a significant Alopecia Areata disease-relevant gene, KRT82, and reveal a proposed mechanism for rare variant predisposition leading to disrupted hair shaft integrity.

Medical Records-Based Genetic Studies of the Complement System.

BACKGROUND: Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. METHODS: We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. RESULTS: In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; ß=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; ß=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; ß=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (ß=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (ß=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. CONCLUSIONS: We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

Blockade of IL-7 signaling suppresses inflammatory responses and reverses alopecia areata in C3H/HeJ mice.

The interleukin-7 (IL-7) signaling pathway plays an important role in regulation of T cell function and survival. We detected overexpression of IL-7 in lesional skin from both humans and C3H/HeJ mice with alopecia areata (AA), a T cell-mediated autoimmune disease of the hair follicle. We found that exogenous IL-7 accelerated the onset of AA by augmenting the expansion of alopecic T cells. Conversely, blockade of IL-7 stopped the progression of AA and reversed early AA in C3H/HeJ mice. Mechanistically, we observed that IL-7Rα blockade substantially reduced the total number of most T cell subsets, but relative sparing of regulatory T cells (Tregs). We postulated that short-term anti-IL-7Rα treatment in combination with a low dose of Treg-tropic cytokines might improve therapeutic efficacy in AA. We demonstrated that short-term IL-7Rα blockade in combination with low doses of Treg-tropic cytokines enhanced therapeutic effects in the treatment of AA, and invite further clinical investigation.

A semi-supervised model to predict regulatory effects of genetic variants at single nucleotide resolution using massively parallel reporter assays.

MOTIVATION: Predicting regulatory effects of genetic variants is a challenging but important problem in functional genomics. Given the relatively low sensitivity of functional assays, and the pervasiveness of class imbalance in functional genomic data, popular statistical prediction models can sharply underestimate the probability of a regulatory effect. We describe here the presence-only model (PO-EN), a type of semi-supervised model, to predict regulatory effects of genetic variants at sequence-level resolution in a context of interest by integrating a large number of epigenetic features and massively parallel reporter assays (MPRAs). RESULTS: Using experimental data from a variety of MPRAs we show that the presence-only model produces better calibrated predicted probabilities and has increased accuracy relative to state-of-the-art prediction models. Furthermore, we show that the predictions based on pre-trained PO-EN models are useful for prioritizing functional variants among candidate eQTLs and significant SNPs at GWAS loci. In particular, for the costimulatory locus, associated with multiple autoimmune diseases, we show evidence of a regulatory variant residing in an enhancer 24.4 kb downstream of CTLA4, with evidence from capture Hi-C of interaction with CTLA4. Furthermore, the risk allele of the regulatory variant is on the same risk increasing haplotype as a functional coding variant in exon 1 of CTLA4, suggesting that the regulatory variant acts jointly with the coding variant leading to increased risk to disease. AVAILABILITY: The presence-only model is implemented in the R package 'PO.EN', freely available on CRAN. A vignette describing a detailed demonstration of using the proposed PO-EN model can be found on github at https://github.com/Iuliana-Ionita-Laza/PO.EN/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An Imperative Need for Further Genetic Studies of Alopecia Areata.

Human genetic studies of diseases that are multifactorial and prevalent have generated a wealth of knowledge about the genetic architecture of chronic diseases. Generalizable attributes are shaping the development of models to explain how the human genome influences our health and can be leveraged to improve it. Importantly, both rare and common genetic variants contribute to disease risk and provide complementary information. Although initial genetic studies of alopecia areata have yielded insight with high clinical impact, there remains a number of important unanswered questions pertaining to disease biology and patient care that could be addressed by further genetic investigations.

Integrative analysis of rare copy number variants and gene expression data in alopecia areata implicates an aetiological role for autophagy.

Alopecia areata (AA) is a highly prevalent autoimmune disease that attacks the hair follicle and leads to hair loss that can range from small patches to complete loss of scalp and body hair. Our previous linkage and genome-wide association studies (GWAS) generated strong evidence for aetiological contributions from inherited genetic variants at different population frequencies, including both rare mutations and common polymorphisms. Additionally, we conducted gene expression (GE) studies on scalp biopsies of 96 patients and controls to establish signatures of active disease. In this study, we performed an integrative analysis on these two datasets to test the hypothesis that rare CNVs in patients with AA could be leveraged to identify drivers of disease in our AA GE signatures. We analysed copy number variants (CNVs) in a case-control cohort of 673 patients with AA and 16 311 controls independent of the case-control cohort of 96 research participants used in our GE study. Using an integrative computational analysis, we identified 14 genes whose expression levels were altered by CNVs in a consistent direction of effect, corresponding to gene expression changes in lesional skin of patients. Four of these genes were affected by CNVs in three or more unrelated patients with AA, including ATG4B and SMARCA2, which are involved in autophagy and chromatin remodelling, respectively. Our findings identified new classes of genes with potential contributions to AA pathogenesis.

Association of Genetic Risk of Obesity with Postoperative Complications Using Mendelian Randomization.

BACKGROUND: The extent to which obesity and genetics determine postoperative complications is incompletely understood. METHODS: We performed a retrospective study using two population cohorts with electronic health record (EHR) data. The first included 736,726 adults with body mass index (BMI) recorded between 1990 and 2017 at Vanderbilt University Medical Center. The second cohort consisted of 65,174 individuals from 12 institutions contributing EHR and genome-wide genotyping data to the Electronic Medical Records and Genomics (eMERGE) Network. Pairwise logistic regression analyses were used to measure the association of BMI categories with postoperative complications derived from International Classification of Disease-9 codes, including postoperative infection, incisional hernia, and intestinal obstruction. A genetic risk score was constructed from 97 obesity-risk single-nucleotide polymorphisms for a Mendelian randomization study to determine the association of genetic risk of obesity on postoperative complications. Logistic regression analyses were adjusted for sex, age, site, and race/principal components. RESULTS: Individuals with overweight or obese BMI (≥25 kg/m2) had increased risk of incisional hernia (odds ratio [OR] 1.7-5.5, p < 3.1 × 10-20), and people with obesity (BMI ≥ 30 kg/m2) had increased risk of postoperative infection (OR 1.2-2.3, p < 2.5 × 10-5). In the eMERGE cohort, genetically predicted BMI was associated with incisional hernia (OR 2.1 [95% CI 1.8-2.5], p = 1.4 × 10-6) and postoperative infection (OR 1.6 [95% CI 1.4-1.9], p = 3.1 × 10-6). Association findings were similar after limitation of the cohorts to those who underwent abdominal procedures. CONCLUSIONS: Clinical and Mendelian randomization studies suggest that obesity, as measured by BMI, is associated with the development of postoperative incisional hernia and infection.

Big Data Reveal Insights into Alopecia Areata Comorbidities.

Autoimmune diseases create a substantial burden of disease, and alopecia areata is among the more prevalent forms. Comorbidities are medical conditions that tend to occur together and may provide etiologic insights, suggest novel therapeutic strategies, and help patients and family members understand the risk of other health conditions. It is well established that having one autoimmune disease increases risk for others because of an underlying shared biology. Precision medicine initiatives are creating vast amounts of data that allow us to efficiently identify comorbidities. A survey across various datasets suggests that patients with autoimmune disease, and patients with alopecia areata in particular, may have comorbid neuropsychiatric and metabolic conditions.

Genome-Wide MicroRNA Analysis Implicates miR-30b/d in the Etiology of Alopecia Areata.

Alopecia areata (AA) is one of the most common forms of human hair loss. Although genetic studies have implicated autoimmune processes in AA etiology, understanding of the etiopathogenesis is incomplete. Recent research has implicated microRNAs, a class of small noncoding RNAs, in diverse autoimmune diseases. To our knowledge, no study has investigated the role of microRNAs in AA. In this study, gene-based analyses were performed for microRNAs using data of the largest genome-wide association meta-analysis of AA to date. Nominally, significant P-values were obtained for 78 of the 617 investigated microRNAs. After correction for multiple testing, three of the 78 microRNAs remained significant. Of these, miR-30b/d was the most significant microRNA for the follow-up analyses, which also showed lower expression in the hair follicle of AA patients. Target gene analyses for the three microRNAs showed 42 significantly associated target genes. These included IL2RA, TNXB, and ERBB3, which had been identified as susceptibility loci in previous genome-wide association studies. Using luciferase assay, site-specific miR-30b regulation of the AA risk genes IL2RA, STX17, and TNXB was validated. This study implicates microRNAs in the pathogenesis of AA. This finding may facilitate the development of future treatment strategies.

Computational derivation of a molecular framework for hair follicle biology from disease genes.

Knowledge about genetic drivers of disease increases the efficiency of interpreting patient DNA sequence and helps to identify and prioritize biological points of intervention. Discoveries of genes with single mutations exerting substantial phenotypic impact reliably provide new biological insight, although such approaches tend to generate knowledge that is disjointed from the complexity of biological systems governed by elaborate networks. Here we sought to facilitate diagnostic sequencing for hair disorders and assess the underlying biology by compiling an archive of 684 genes discovered in studies of monogenic disorders and identifying molecular annotations enriched by them. To demonstrate utility for this dataset, we performed two data driven analyses. First, we extracted and analyzed data implicating enriched signaling pathways and identified previously unrecognized contributions from Hippo signaling. Second, we performed hierarchical clustering on the entire dataset to investigate the underlying causal structure of hair disorders. We identified 35 gene clusters representing genetically derived biological modules that provide a foundation for the development of a new disease taxonomy grounded in biology, rather than clinical presentations alone. This Resource will be useful for diagnostic sequencing in patients with diseases affecting the hair follicle, improved characterization of hair follicle biology, and methods development in precision medicine.