Phenome-wide association study maps new diseases to the human major histocompatibility complex region

Lichen Planus: What is New in Diagnosis and Treatment?

Lichen planus (LP), an idiopathic, multifaceted chronic inflammatory disease with a heterogeneous clinical presentation, affects approximately 0.5-1% of the population. The various clinical manifestations of LP fall into three broad categories, namely cutaneous, appendageal, and mucosal, with further subclassification depending on the morphology and distribution patterns of individual lesions. There is mounting evidence that LP has systemic associations, including autoimmune conditions, glucose intolerance, dyslipidemia, and cardiovascular disorders. Cutaneous hypertrophic and mucosal forms of LP are at a heightened risk for malignant transformation. Familiarity with these potential associations in conjunction with long-term follow-up and regular screening could lead to a timely diagnosis and management of concomitant conditions. In addition, the frequent quality of life (QoL) impairment in LP underscores the need for a comprehensive approach including psychological evaluation and support. Several treatment strategies have been attempted, though most of them have not been adopted in clinical practice because of suboptimal benefit-to-risk ratios or lack of evidence. More recent studies toward pathogenesis-driven treatments have identified Janus kinase inhibitors such as tofacitinib, phosphodiesterase-4 inhibitors such as apremilast, and biologics targeting the interleukin-23/interleukin-17 pathway as novel therapeutic options, resulting in a dramatic change of the treatment landscape of LP. This contemporary review focuses on the diagnosis and management of LP, and places emphasis on more recently described targeted treatment options.

Significant Correlation Between Cutaneous Abundance of Streptococcus and Psoriasis Severity in Patients with FBXL19 Gene Variants

Psoriasis results from both genetic predisposition and environmental triggers, such as Streptococcal infections. This study aimed to explore the correlation between the abundance of the Streptococcus genus on the skin and psoriasis severity in individuals carrying specific psoriasis-associated genetic variants. Studying 39 chronic plaque psoriasis patients, the elbow skin microbiome and 49 psoriasis-related single nucleotide polymorphisms (SNPs) were analysed using a MiSeq instrument for 16S rDNA sequencing, and CLC Genomic Workbench for processing and analysis. Through multivariate linear regression analysis, a positive correlation was found between Streptococcus genus abundance and psoriasis severity in patients with certain FBXL19 gene-related heterozygous SNPs (rs12924903, rs10782001, rs12445568). Conversely, a negative association was observed in patients with homozygous genotypes. Moreover, we identified an association between Streptococcus abundance and psoriasis severity in patients with genetic variants related to IL-22, ERAP1, NOS2, and ILF3. This is the first study highlighting a positive association between Streptococcus skin colonization and psoriasis severity in patients with heterozygous genotypes within the FBXL19 gene region. FXBL19 targets the IL-33/IL1RL1 axis, crucial in infectious diseases and innate immunity promotion. These novel results suggests an intricate interaction among host genetics, Streptococcus skin colonization, and psoriasis inflammation, offering potential avenues for novel treatment approaches.

Major histocompatibility complex (MHC) gene frequency in acquired dermal macular hyperpigmentation: a case control study

BACKGROUND

Human leukocyte antigen (HLA) allele frequencies have a known association with the pathogenesis of various autoimmune diseases.

METHODS

We recruited 31 Indian patients of acquired dermal macular hyperpigmentation (ADMH) and 60 unrelated, age-and-gender-matched healthy controls. After history and clinical examination, 5 ml of blood in EDTA vials was collected. These samples were subjected to DNA extraction and the expression of HLA A, B, C, DR, DQ-A, and DQ-B was studied.

RESULTS

There was a predominance of females with a gender ratio of 23 : 8 and the most common phototype was Fitzpatrick type IV (83.9%). There was a significant association of HLA A*03:01 (OR: 5.8, CI: 1.7-17.0, P = 0.005), HLA B*07:02 (OR: 5.3, CI: 1.9-14.6, P = 0.003), HLA C*07:02 (OR: 4.3, CI: 1.8-9.6, P = 0.001), HLA DRB1*10:01 (OR: 7.6, CI: 1.7-38.00, P = 0.022), and HLA DRB1*15:02 (OR: 31.0, CI: 4.4-341.8, P < 0.001) with patients compared to controls, whereas HLA DQB*03:01 was less associated with patients compared to controls (OR: 0.2, CI: 0.0-0.6, P = 0.009).

CONCLUSION

Patients with ADMH are more likely to have the HLA A*03:01, HLA B 07*02, HLA C*07:02, HLA DRB1*10:01, HLA DRB1*15:02 and less likely to have the HLA DQB*03:01 allele. Larger cohort studies may thus be conducted studying these specific alleles.

Oral Lichen Planus: An Updated Review of Etiopathogenesis, Clinical Presentation, and Management

Lichen planus (LP) is a chronic idiopathic immune-mediated inflammatory condition. LP is a heterogeneous disease with varied clinical presentations having different natural history, prognosis, sequelae, and outcomes. It can affect skin, hair, nails, and mucosae. Mucosal LP (including oral LP) tends to be persistent and resistant to treatment, compared to cutaneous LP. Oral LP (OLP) is broadly divided into two main categories: hyperkeratotic (usually asymptomatic) and erosive (commonly symptomatic). It can present with symptoms including odynophagia, dysphagia, dysgeusia, and sensitivity to hot spicy foods. Apart from the superficial epidermal changes, which vary with the type of clinical presentation, histopathologically oral LP shows a unifying similar and consistent feature of a lichenoid interface dermatitis. Recently, researchers have highlighted the critical role played by IL-17 in the pathogenesis of OLP. World Health Organization has categorized oral LP as one of the oral potentially malignant disorders (OPMD), albeit with a low risk of malignant transformation. Also, in the last couple of years there have been various reports on the usage of newer drugs like anti-IL17, anti-IL12/23, anti-IL 23, PDE4 inhibitors, and JAK inhibitors in the management of refractory OLP. The principal aim of treatment still remains to resolve the symptoms, prolong the symptoms free period, and reduce the risk of potential malignant transformation. We have described many new revelations made in recent times regarding the etiopathogenesis, associated conditions as well as management of OLP. Thus, the objective of this review is to present a comprehensive up-to-date knowledge including the recent advances made regarding OLP.

Phenome-wide association study on miRNA-related sequence variants: the UK Biobank

BACKGROUND

Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown.

METHODS

Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis.

RESULTS

We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia.

CONCLUSIONS

Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Estimating the efficacy of pharmacogenomics over a lifetime

It is well known that common variants in specific genes influence drug metabolism and response, but it is currently unknown what fraction of patients are given prescriptions over a lifetime that could be contraindicated by their pharmacogenomic profiles. To determine the clinical utility of pharmacogenomics over a lifetime in a general patient population, we sequenced the genomes of 300 deceased Marshfield Clinic patients linked to lifelong medical records. Genetic variants in 33 pharmacogenes were evaluated for their lifetime impact on drug prescribing using extensive electronic health records. Results show that 93% of the 300 deceased patients carried clinically relevant variants. Nearly 80% were prescribed approximately three medications on average that may have been impacted by these variants. Longitudinal data suggested that the optimal age for pharmacogenomic testing was prior to age 50, but the optimal age is greatly influenced by the stability of the population in the healthcare system. This study emphasizes the broad clinical impact of pharmacogenomic testing over a lifetime and demonstrates the potential application of genomic medicine in a general patient population for the advancement of precision medicine.

Comprehensive Insight into Lichen Planus Immunopathogenesis

Lichen planus is a chronic disease affecting the skin, appendages, and mucous membranes. A cutaneous lichen planus is a rare disease occurring in less than 1% of the general population, while oral illness is up to five times more prevalent; still, both forms equally impair the patient's quality of life. The etiology of lichen planus is not entirely understood. Yet, immune-mediated mechanisms have been recognized since environmental factors such as hepatitis virus infection, mechanical trauma, psychological stress, or microbiome changes can trigger the disease in genetically susceptible individuals. According to current understanding, lichen planus immunopathogenesis is caused by cell-mediated cytotoxicity, particularly cytotoxic T lymphocytes, whose activity is further influenced by Th1 and IL-23/Th-17 axis. However, other immunocytes and inflammatory pathways complement these mechanisms. This paper presents a comprehensive insight into the actual knowledge about lichen planus, with the causal genetic and environmental factors being discussed, the immunopathogenesis described, and the principal effectors of its inflammatory circuits identified.

Phenome-wide association study of the major histocompatibility complex region in the Korean population identifies novel association signals

Human leukocyte antigen (HLA) gene variants in the major histocompatibility complex (MHC) region are associated with numerous complex human diseases and quantitative traits. Previous phenome-wide association studies (PheWAS) for this region demonstrated that HLA association patterns to the phenome have both population-specific and population-shared components. We performed MHC PheWAS in the Korean population by analyzing associations between phenotypes and genetic variants in the MHC region using the Korea Biobank Array project data samples from the Korean Genome and Epidemiology Study (KoGES) cohorts. Using this single-population dataset, we curated and analyzed 82 phenotypes for 125 673 Korean individuals after imputing HLA using CookHLA, a recently developed imputation framework. More than one-third of these phenotypes showed significant associations, confirming 56 known associations and discovering 13 novel association signals that were not reported previously. In addition, we analyzed heritability explained by the variants in the MHC region and genetic correlations among phenotypes based on the MHC variants.

pyPheWAS: A Phenome-Disease Association Tool for Electronic Medical Record Analysis

Along with the increasing availability of electronic medical record (EMR) data, phenome-wide association studies (PheWAS) and phenome-disease association studies (PheDAS) have become a prominent, first-line method of analysis for uncovering the secrets of EMR. Despite this recent growth, there is a lack of approachable software tools for conducting these analyses on large-scale EMR cohorts. In this article, we introduce pyPheWAS, an open-source python package for conducting PheDAS and related analyses. This toolkit includes 1) data preparation, such as cohort censoring and age-matching; 2) traditional PheDAS analysis of ICD-9 and ICD-10 billing codes; 3) PheDAS analysis applied to a novel EMR phenotype mapping: current procedural terminology (CPT) codes; and 4) novelty analysis of significant disease-phenotype associations found through PheDAS. The pyPheWAS toolkit is approachable and comprehensive, encapsulating data prep through result visualization all within a simple command-line interface. The toolkit is designed for the ever-growing scale of available EMR data, with the ability to analyze cohorts of 100,000 + patients in less than 2 h. Through a case study of Down Syndrome and other intellectual developmental disabilities, we demonstrate the ability of pyPheWAS to discover both known and potentially novel disease-phenotype associations across different experiment designs and disease groups. The software and user documentation are available in open source at https://github.com/MASILab/pyPheWAS .

Lichen Planus

Lichen planus (LP) is a T cell-mediated disease affecting the stratified squamous epithelia of the skin and/or mucus membrane. Histologically, the disease is characterized by a lichenoid inflammatory infiltrate and vacuolar degeneration of the basal layer of the epidermis. LP has three major subtypes: Cutaneous, mucosal and appendageal LP. Rarely, it may affect the nails in the absence of skin and/or mucosal changes. LP may also be induced by several drugs, typically anti-hypertensive medication or be associated with infections, particularly viral hepatitis. The diagnosis is based on the clinical presentation and characteristic histological findings. Although the disease is often self-limiting, the intractable pruritus and painful mucosal erosions result in significant morbidity. The current first-line treatment are topical and/or systemic corticosteroids. In addition, immunosuppressants may be used as corticosteroid-sparing agents. These, however are often not sufficient to control disease. Janus kinase inhibitors and biologics (anti-IL-12/23, anti-IL17) have emerged as novel future treatment options. Thus, one may expect a dramatic change of the treatment landscape of LP in the near future.

Targeted RNA-Based Oxford Nanopore Sequencing for Typing 12 Classical HLA Genes

Identification of human leukocyte antigen (HLA) alleles from next-generation sequencing (NGS) data is challenging because of the high polymorphism and mosaic nature of HLA genes. Owing to the complex nature of HLA genes and consequent challenges in allele assignment, Oxford Nanopore Technologies' (ONT) single-molecule sequencing technology has been of great interest due to its fitness for sequencing long reads. In addition to the read length, ONT's advantages are its portability and possibility for a rapid real-time sequencing, which enables a simultaneous data analysis. Here, we describe a targeted RNA-based method for HLA typing using ONT sequencing and SeqNext-HLA SeqPilot software (JSI Medical Systems GmbH). Twelve classical HLA genes were enriched from cDNA of 50 individuals, barcoded, pooled, and sequenced in 10 MinION R9.4 SpotON flow cell runs producing over 30,000 reads per sample. Using barcoded 2D reads, SeqPilot assigned HLA alleles to two-field typing resolution or higher with the average read depth of 1750x. Sequence analysis resulted in 99-100% accuracy at low-resolution level (one-field) and in 74-100% accuracy at high-resolution level (two-field) with the expected alleles. There are still some limitations with ONT RNA sequencing, such as noisy reads, homopolymer errors, and the lack of robust algorithms, which interfere with confident allele assignment. These issues need to be inspected carefully in the future to improve the allele call rates. Nevertheless, here we show that sequencing of multiplexed cDNA amplicon libraries on ONT MinION can produce accurate high-resolution typing results of 12 classical HLA loci. For HLA research, ONT RNA sequencing is a promising method due to its capability to sequence full-length HLA transcripts. In addition to HLA genotyping, the technique could also be applied for simultaneous expression analysis.

PheProb: probabilistic phenotyping using diagnosis codes to improve power for genetic association studies

Objective

Standard approaches for large scale phenotypic screens using electronic health record (EHR) data apply thresholds, such as ≥2 diagnosis codes, to define subjects as having a phenotype. However, the variation in the accuracy of diagnosis codes can impair the power of such screens. Our objective was to develop and evaluate an approach which converts diagnosis codes into a probability of a phenotype (PheProb). We hypothesized that this alternate approach for defining phenotypes would improve power for genetic association studies.

Methods

The PheProb approach employs unsupervised clustering to separate patients into 2 groups based on diagnosis codes. Subjects are assigned a probability of having the phenotype based on the number of diagnosis codes. This approach was developed using simulated EHR data and tested in a real world EHR cohort. In the latter, we tested the association between low density lipoprotein cholesterol (LDL-C) genetic risk alleles known for association with hyperlipidemia and hyperlipidemia codes (ICD-9 272.x). PheProb and thresholding approaches were compared.

Results

Among n = 1462 subjects in the real world EHR cohort, the threshold-based p-values for association between the genetic risk score (GRS) and hyperlipidemia were 0.126 (≥1 code), 0.123 (≥2 codes), and 0.142 (≥3 codes). The PheProb approach produced the expected significant association between the GRS and hyperlipidemia: p = .001.

Conclusions

PheProb improves statistical power for association studies relative to standard thresholding approaches by leveraging information about the phenotype in the billing code counts. The PheProb approach has direct applications where efficient approaches are required, such as in Phenome-Wide Association Studies.

Mapping ICD-10 and ICD-10-CM Codes to Phecodes: Workflow Development and Initial Evaluation

BACKGROUND

The phecode system was built upon the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) for phenome-wide association studies (PheWAS) using the electronic health record (EHR).

OBJECTIVE

The goal of this paper was to develop and perform an initial evaluation of maps from the International Classification of Diseases, 10th Revision (ICD-10) and the International Classification of Diseases, 10th Revision, Clinical Modification (ICD-10-CM) codes to phecodes.

METHODS

We mapped ICD-10 and ICD-10-CM codes to phecodes using a number of methods and resources, such as concept relationships and explicit mappings from the Centers for Medicare & Medicaid Services, the Unified Medical Language System, Observational Health Data Sciences and Informatics, Systematized Nomenclature of Medicine-Clinical Terms, and the National Library of Medicine. We assessed the coverage of the maps in two databases: Vanderbilt University Medical Center (VUMC) using ICD-10-CM and the UK Biobank (UKBB) using ICD-10. We assessed the fidelity of the ICD-10-CM map in comparison to the gold-standard ICD-9-CM phecode map by investigating phenotype reproducibility and conducting a PheWAS.

RESULTS

We mapped >75% of ICD-10 and ICD-10-CM codes to phecodes. Of the unique codes observed in the UKBB (ICD-10) and VUMC (ICD-10-CM) cohorts, >90% were mapped to phecodes. We observed 70-75% reproducibility for chronic diseases and <10% for an acute disease for phenotypes sourced from the ICD-10-CM phecode map. Using the ICD-9-CM and ICD-10-CM maps, we conducted a PheWAS with a Lipoprotein(a) genetic variant, rs10455872, which replicated two known genotype-phenotype associations with similar effect sizes: coronary atherosclerosis (ICD-9-CM: P<.001; odds ratio (OR) 1.60 [95% CI 1.43-1.80] vs ICD-10-CM: P<.001; OR 1.60 [95% CI 1.43-1.80]) and chronic ischemic heart disease (ICD-9-CM: P<.001; OR 1.56 [95% CI 1.35-1.79] vs ICD-10-CM: P<.001; OR 1.47 [95% CI 1.22-1.77]).

CONCLUSIONS

This study introduces the beta versions of ICD-10 and ICD-10-CM to phecode maps that enable researchers to leverage accumulated ICD-10 and ICD-10-CM data for PheWAS in the EHR.

Variation in Cilia Protein Genes and Progression of Lung Disease in Cystic Fibrosis

RATIONALE

Cystic fibrosis, like primary ciliary dyskinesia, is an autosomal recessive disorder characterized by abnormal mucociliary clearance and obstructive lung disease. We hypothesized that genes underlying the development or function of cilia may modify lung disease severity in persons with cystic fibrosis.

OBJECTIVES

To test this hypothesis, we compared variants in 93 candidate genes in both upper and lower tertiles of lung function in a large cohort of children and adults with cystic fibrosis with those of a population control dataset.

METHODS

Variants within candidate genes were tested for association using the SKAT-O test, comparing cystic fibrosis cases defined by poor (n = 127) or preserved (n = 127) lung function with population controls (n = 3,269 or 3,148, respectively). Associated variants were then tested for association with related phenotypes in independent datasets.

RESULTS

Variants in DNAH14 and DNAAF3 were associated with poor lung function in cystic fibrosis, whereas variants in DNAH14 and DNAH6 were associated with preserved lung function in cystic fibrosis. Associations between DNAH14 and lung function were replicated in disease-related phenotypes characterized by obstructive lung disease in adults.

CONCLUSIONS

Genetic variants within DNAH6, DNAH14, and DNAAF3 are associated with variation in lung function among persons with cystic fibrosis.

Lichen planus and lichenoid dermatoses: Clinical overview and molecular basis

Deriving from the Greek word λειχήν for "tree moss" and the Latin word planus for "planar," lichen planus is a relatively uncommon and heterogeneous cutaneous disorder that typically develops in middle-aged adults. Despite the significant clinical burden associated with the disorder, little well-conducted molecular research has been undertaken, possibly because of heterogeneity impeding consistent and confident phenotyping. The multiple variants of lichenoid disease bear overlapping clinical and pathologic features despite manifesting as distinct clinical disorders. The first article in this 2-part continuing medical education series provides a comprehensive overview of the clinical and pathologic characteristics of cutaneous lichenoid dermatoses and links these manifestations to recent advances in our understanding of the underlying pathobiology of such diseases.

Genomic and Phenomic Research in the 21st Century

The field of human genomics has changed dramatically over time. Initial genomic studies were predominantly restricted to rare disorders in small families. Over the past decade, researchers changed course from family-based studies and instead focused on common diseases and traits in populations of unrelated individuals. With further advancements in biobanking, computer science, electronic health record (EHR) data, and more affordable high-throughput genomics, we are experiencing a new paradigm in human genomic research. Rapidly changing technologies and resources now make it possible to study thousands of diseases simultaneously at the genomic level. This review will focus on these advancements as scientists begin to incorporate phenome-wide strategies in human genomic research to understand the etiology of human diseases and develop new drugs to treat them.

Vulvo-vaginal lichen planus: A focussed review for the clinician

Vulvo-vaginal lichen planus is a condition which imposes a significant burden of symptoms and sequale. There is a paucity of knowledge and evidenced-based management with significant diagnostic delay prior to appropriate treatment being common. There remains great variability in clinician practice in the context of limited knowledge. This clinical review presents current evidence on the clinicopathological features, practical assessment and management options. Learning points include characterization of this chronic, burdensome clinical entity that has no standardized diagnostic or management protocols. We provide practical clinical conclusions based on current knowledge and identify areas for future research to improve patient outcomes.

Statistical methods to detect pleiotropy in human complex traits

In recent years pleiotropy, the phenomenon of one genetic locus influencing several traits, has become a widely researched field in human genetics. With the increasing availability of genome-wide association study summary statistics, as well as the establishment of deeply phenotyped sample collections, it is now possible to systematically assess the genetic overlap between multiple traits and diseases. In addition to increasing power to detect associated variants, multi-trait methods can also aid our understanding of how different disorders are aetiologically linked by highlighting relevant biological pathways. A plethora of available tools to perform such analyses exists, each with their own advantages and limitations. In this review, we outline some of the currently available methods to conduct multi-trait analyses. First, we briefly introduce the concept of pleiotropy and outline the current landscape of pleiotropy research in human genetics; second, we describe analytical considerations and analysis methods; finally, we discuss future directions for the field.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

Current Scope and Challenges in Phenome-Wide Association Studies

PURPOSE OF REVIEW

Over many decades, researchers have been designing studies to investigate the relationship between genotypes and phenotypes to gain an understanding about the effect of genetics on disease. Recently, a high-throughput approach called phenome-wide associations studies (PheWAS) have been extensively used to identify associations between genetic variants and many diseases and traits simultaneously. In this review, we describe the value of PheWAS along with methodological issues and challenges in interpretation for current applications of PheWAS.

RECENT FINDINGS

PheWAS have uncovered a paradigm to identify new associations for genetic loci across many diseases. The application of PheWAS have been effective with phenotype data from electronic health records, epidemiological studies, and clinical trials data.

SUMMARY

The key strength of a PheWAS is to identify the association of one or more genetic variants with multiple phenotypes, which can showcase interconnections among the phenotypes due to shared genetic associations. While the PheWAS approach appears promising, there are a number of challenges that need to be addressed to provide additional robustness to PheWAS findings.

Phenome-wide scanning identifies multiple diseases and disease severity phenotypes associated with HLA variants

Although many phenotypes have been associated with variants in human leukocyte antigen (HLA) genes, the full phenotypic impact of HLA variants across all diseases is unknown. We imputed HLA genomic variation from two populations of 28,839 and 8431 European ancestry individuals and tested association of HLA variation with 1368 phenotypes. A total of 104 four-digit and 92 two-digit HLA allele phenotype associations were significant in both discovery and replication cohorts, the strongest being HLA-DQB1*03:02 and type 1 diabetes. Four previously unidentified associations were identified across the spectrum of disease with two- and four-digit HLA alleles and 10 with nonsynonymous variants. Some conditions associated with multiple HLA variants and stronger associations with more severe disease manifestations were identified. A comprehensive, publicly available catalog of clinical phenotypes associated with HLA variation is provided. Examining HLA variant disease associations in this large data set allows comprehensive definition of disease associations to drive further mechanistic insights.