NFKB1 and MANBA Confer Disease Susceptibility to Primary Biliary Cholangitis via Independent Putative Primary Functional Variants

A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis

BACKGROUND AND AIMS

Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-nonspecific) and nonshared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC.

APPROACH AND RESULTS

Protein tyrosine phosphatase nonreceptor type 2 ( PTPN2) was identified as a novel PBC susceptibility gene locus through GWAS and subsequent genome-wide meta-analysis involving 2181 cases and 2699 controls from the Japanese population (GWAS-lead variant: rs8098858, p = 2.6 × 10 -8 ). In silico and in vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells and plasmacytoid dendritic cells. Infiltration of PTPN2-positive T-cells and plasmacytoid dendritic cells was confirmed in the portal area of the PBC liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758.

CONCLUSIONS

PTPN2 , a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC through an insufficient negative feedback loop caused by the risk allele of rs2292758 in IFN-γ signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

Noninvasive assessment of organ-specific and shared pathways in multi-organ fibrosis using T1 mapping

Fibrotic diseases affect multiple organs and are associated with morbidity and mortality. To examine organ-specific and shared biologic mechanisms that underlie fibrosis in different organs, we developed machine learning models to quantify T1 time, a marker of interstitial fibrosis, in the liver, pancreas, heart and kidney among 43,881 UK Biobank participants who underwent magnetic resonance imaging. In phenome-wide association analyses, we demonstrate the association of increased organ-specific T1 time, reflecting increased interstitial fibrosis, with prevalent diseases across multiple organ systems. In genome-wide association analyses, we identified 27, 18, 11 and 10 independent genetic loci associated with liver, pancreas, myocardial and renal cortex T1 time, respectively. There was a modest genetic correlation between the examined organs. Several loci overlapped across the examined organs implicating genes involved in a myriad of biologic pathways including metal ion transport (SLC39A8, HFE and TMPRSS6), glucose metabolism (PCK2), blood group antigens (ABO and FUT2), immune function (BANK1 and PPP3CA), inflammation (NFKB1) and mitosis (CENPE). Finally, we found that an increasing number of organs with T1 time falling in the top quintile was associated with increased mortality in the population. Individuals with a high burden of fibrosis in ≥3 organs had a 3-fold increase in mortality compared to those with a low burden of fibrosis across all examined organs in multivariable-adjusted analysis (hazard ratio = 3.31, 95% confidence interval 1.77-6.19; P = 1.78 × 10-4). By leveraging machine learning to quantify T1 time across multiple organs at scale, we uncovered new organ-specific and shared biologic pathways underlying fibrosis that may provide therapeutic targets.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Mapping of the gene network that regulates glycan clock of ageing

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

A Systematic Review of Host Genomic Variation and Neuropsychological Outcomes for Pediatric Cancer Survivors

Pediatric survivors of brain tumors and acute lymphoblastic leukemia (ALL) are at risk for long-term deficits in their neuropsychological functioning. Researchers have begun examining associations between germline single nucleotide polymorphisms (SNPs), which interact with cancer treatment, and neuropsychological outcomes. This review synthesizes the impact of treatment-related toxicity from germline SNPs by neuropsychological domain (i.e., working memory, processing speed, psychological functioning) in pediatric survivors. By focusing on specific neuropsychological domains, this review will examine outcome measurement and critique methodology. Fourteen studies were identified and included in this review using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). All studies were published in peer-reviewed journals in English by November 24^th, 2021. Reviewed studies were not of sufficient quality for a meta-analysis due to varying measurement strategies, gaps in reported descriptive variables, and low power. All neuropsychological domains evaluated in this review had associations with SNPs, except fine motor and visual integration abilities. Only five SNPs had consistent neuropsychological findings in more than one study or cohort. Future research and replication studies should use validated measures of discrete skills that are central to empirically validated models of survivors' long-term outcomes (i.e., attention, working memory, processing speed). Researchers should examine SNPs across pathophysiological pathways to investigate additive genetic risk in pediatric cancer survivors. Two SNPs were identified that confer resiliency in neuropsychological functioning, and future work should investigate resiliency genotypes and their underlying biological mechanisms.

The Genetics of Primary Biliary Cholangitis: A GWAS and Post-GWAS Update

Primary biliary cholangitis (PBC) is a chronic, progressive cholestatic liver disease in which the small intrahepatic bile ducts are destroyed by autoimmune reactions. Among autoimmune diseases, which are polygenic complex traits caused by the combined contribution of genetic and environmental factors, PBC exhibits the strongest involvement of genetic heritability in disease development. As at December 2022, genome-wide association studies (GWASs) and associated meta-analyses identified approximately 70 PBC susceptibility gene loci in various populations, including those of European and East Asian descent. However, the molecular mechanisms through which these susceptibility loci affect the pathogenesis of PBC are not fully understood. This study provides an overview of current data regarding the genetic factors of PBC as well as post-GWAS approaches to identifying primary functional variants and effector genes in disease-susceptibility loci. Possible mechanisms of these genetic factors in the development of PBC are also discussed, focusing on four major disease pathways identified by in silico gene set analyses, namely, (1) antigen presentation by human leukocyte antigens, (2) interleukin-12-related pathways, (3) cellular responses to tumor necrosis factor, and (4) B cell activation, maturation, and differentiation pathways.

The Clostridium Metabolite P-Cresol Sulfate Relieves Inflammation of Primary Biliary Cholangitis by Regulating Kupffer Cells

OBJECTIVE

To study the effect and mechanism of the Clostridium metabolite p-Cresol sulfate (PCS) in primary biliary cholangitis (PBC).

METHODS

Gas chromatography-mass spectrometry (GC-MS) was used to detect differences in tyrosine, phenylalanine, tryptophan, PCS, and p-Cresyl glucuronide (PCG) between the serum of PBC patients and healthy controls. In vivo experiments, mice were divided into the normal control, PBC group, and PBC tyrosine group. GC-MS was used to detect PCS and PCG. Serum and liver inflammatory factors were compared between groups along with the polarization of liver Kupffer cells. Additionally, PCS was cultured with normal bile duct epithelial cells and Kupffer cells, respectively. PCS-stimulated Kupffer cells were co-cultured with lipopolysaccharide-injured bile duct epithelial cells to detect changes in inflammatory factors.

RESULTS

Levels of tyrosine and phenylalanine were increased, but PCS level was reduced in PBC patients, with PCG showing a lower concentration distribution in both groups. PCS in PBC mice was also lower than those in normal control mice. After oral administration of tyrosine feed to PBC mice, PCS increased, liver inflammatory factors were decreased, and anti-inflammatory factors were increased. Furthermore, Kupffer cells in the liver polarized form M1 transitioned to M2. PCS can damage normal bile duct epithelial cells and suppress the immune response of Kupffer cells. But PCS protects bile duct epithelial cells damaged by LPS through Kupffer cells.

CONCLUSIONS

PCS produced by Clostridium-metabolized tyrosine reduced PBC inflammation, suggesting that intervention by food, or supplementation with PCS might represent an effective clinical strategy for treating PBC.

Endothelin receptor B enhances liver injury and pro-inflammatory responses by increasing G-protein-coupled receptor kinase-2 expression in primary biliary cholangitis

Severe diseases like cirrhosis and liver failure can be developed from primary biliary cholangitis (PBC). Endothelin-2 (EDN2) and endothelin receptor B (EDNRB) are related to the pathogenesis of PBC. However, the roles of EDN2 and EDNRB in PBC-related liver injury and inflammation along with molecular mechanisms are poorly defined. In this study, histopathologic alterations of liver tissues were assessed through hematoxylin-eosin staining. Alanine transaminase (ALT), alkaline phosphatase (ALP), aspartate transaminase (AST), and γ-Glutamyltranspetidase (GGT) (4 liver function indexes) serum levels were detected with corresponding activity assay kits. Also, we determined the levels of M2 subtype anti-mitochondrial antibody (AMA-M2), interferon-gamma (IFN-γ), and tumor-necrosis factor alpha (TNFα) in serum with ELISA assay. Later, RT-qPCR assay was used to measure the expression of genes at mRNA levels, while western blotting and immunohistochemical techniques were used to detect protein levels of genes. Our results showed that the liver tissues of PBC patients and mice presented with severe hepatocyte injury and inflammatory cell infiltration as well as destruction of intrahepatic small bile ducts. ALP, AST, ALT, GGT, AMA-M2, IFN-γ, and TNF-α serum levels were higher in PBC patients and mice. Besides, EDN2 and EDNRB were highly expressed in serums and livers of PBC patients and mice. EDNRB potentiated PBC-related liver injury and pro-inflammatory responses, as evidenced by observation of serious liver pathologic injury and increased serum levels of ALP, AST, ALT, AMA-M2, IFN-γ, and TNF-α in PBC mice following EDNRB overexpression. EDNRB overexpression or activation via its agonist IRL-1620 TFA triggered liver injury and pro-inflammatory responses, increased GRK2 expression and induced NF-κB expression and activation in wild-type mice. EDNRB knockdown or inhibition by Bosentan alleviated liver damage and inflammation, reduced GRK2 expression, and inhibited NF-κB in PBC mice. These findings suggested EDNRB loss or inhibition weakened liver injury and pro-inflammatory responses by down-regulating GRK2 and inhibiting the NF-κB pathway in PBC mice.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms

Background

Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule.

Results

Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms.

Conclusion

The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00419-7.

NFκB1 Polymorphisms Are Associated with Severe Influenza A (H1N1) Virus Infection in a Canadian Population

Background: We examined associations between NFκB1 polymorphisms and influenza A (H1N1) clinical outcomes in Canadian. Methods: A total of thirty-six Caucasian patients admitted to the intensive care unit (ICU) in hospitals in Canada were recruited during the 2009 H1N1 pandemic. Genomic DNA was extracted from the whole blood samples. The NFkB1 gene was targeted for genotyping using next-generation sequencing technology—Roche 454. Results: A total of 136 single nucleotide polymorphisms (SNPs) were discovered within the NFκB1 gene. Among them, 63 SNPs were significantly enriched in patients admitted in the ICU (p < 0.05) compared with the British Caucasian population in the 1000 Genomes study. These enriched SNPs are mainly intron variants, and only two are exon SNPs from the non-transcribing portion of the NFκB1 gene. Conclusions: Genetic variations in the NFκB1 gene could influence clinical outcomes of pandemic H1N1 infections. Our findings showed that sequence variations of the NFκB1 gene might influence patient response to influenza infection.

Impact of Multiple Sclerosis Risk Polymorphism rs7665090 on MANBA Activity, Lysosomal Endocytosis, and Lymphocyte Activation

Deficiencies in Mannosidase β (MANBA) are associated with neurological abnormalities and recurrent infections. The single nucleotide polymorphism located in the 3′UTR of MANBA, rs7665090, was found to be associated with multiple sclerosis (MS) susceptibility. We aimed to study the functional impact of this polymorphism in lymphocytes isolated from MS patients and healthy controls. A total of 152 MS patients and 112 controls were genotyped for rs7665090. MANBA mRNA expression was quantified through qPCR and MANBA enzymatic activity was analyzed. Upon phytohemagglutinin stimulation, immune activation was evaluated by flow cytometry detection of CD69, endocytic function, and metabolic rates with Seahorse XFp Analyzer, and results were stratified by variation in rs7665090. A significantly reduced gene expression (p < 0.0001) and enzymatic activity (p = 0.018) of MANBA were found in lymphocytes of MS patients compared to those of controls. The rs7665090*GG genotype led to a significant β-mannosidase enzymatic deficiency correlated with lysosomal dysfunction, as well as decreased metabolic activation in lymphocytes of MS patients compared to those of rs7665090*GG controls. In contrast, lymphocytes of MS patients and controls carrying the homozygous AA genotype behaved similarly. Our work provides new evidence highlighting the impact of the MS-risk variant, rs7665090, and the role of MANBA in the immunopathology of MS.

Differentiation and comorbidity of bipolar disorder and attention deficit and hyperactivity disorder in children, adolescents, and adults: A clinical and nosological perspective

Bipolar Disorder (BD) and Attention Deficit and Hyperactivity Disorder (ADHD) are mental disorders with high degree of lifetime comorbidity. Both BD and ADHD are disorders with onset in childhood and early adolescence. Both disorders are often undiagnosed, misdiagnosed, and sometimes overdiagnosed, leading to high rates of morbidity and disability. The psychiatric and behavioral symptoms associated with ADHD and BD have significant overlap. Albeit the existence of a large body of literature, it is far from being clear whether comorbidity can be explained by the confounding overlap of operationally defined criteria or whether it reflects a genuine comorbidity of two biologically distinct disorders. The aim of this paper is to recognize and/or differentiate the pattern of ADHD across the course of BD from a nosological point of view, focusing on specific clinical and neurobiological dimensions. We found that some critical issues may help to fulfill the purpose of our perspective. We suggest that the relationship between ADHD and BD, based on clinical, developmental, and epidemiological commonalities, can be better clarified using four different scenarios.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

Emerging Insights Into Chronic Renal Disease Pathogenesis in Hypertension From Human and Animal Genomic Studies

The pathogenic links between elevated blood pressure and chronic kidney disease remain obscure. This article examines progress in population genetics and in animal models of hypertension and chronic kidney disease. It also provides a critique of the application of genome-wide association studies to understanding the heritability of renal function. Emerging themes identified indicate that heritable risk of chronic kidney disease in hypertension can arise from genetic variation in (1) glomerular and tubular protein handling mechanisms; (2) autoregulatory capacity of the renal vasculature; and (3) innate and adaptive immune mechanisms. Increased prevalence of hypertension-associated chronic kidney disease that occurs with aging may reflect amplification of heritable risks by normal aging processes affecting immunity and autoregulation.

Regional heritability mapping identifies several novel loci (STAT4, ULK4, and KCNH5) for primary biliary cholangitis in the Japanese population

While the advent of GWAS more than a decade ago has ushered in remarkable advances in our understanding of complex traits, the limitations of single-SNP analysis have also led to the development of several other approaches. Simulation studies have shown that the regional heritability mapping (RHM) method, which makes use of multiple adjacent SNPs jointly to estimate the genetic effect of a given region of the genome, generally has higher detection power than single-SNP GWAS. However, thus far its use has been mostly limited to agricultural settings, and its potential for the discovery of new genes in human diseases is yet to be fully exploited. In this study, by applying the RHM method to primary biliary cholangitis (PBC) in the Japanese population, we identified three novel loci (STAT4, ULK4, and KCNH5) at the genome-wide significance level, two of which (ULK4 and KCNH5) have not been found associated with PBC in any population previously. Notably, these genes could not be detected by using conventional single-SNP GWAS, highlighting the potential of the RHM method for the detection of new susceptibility loci in human diseases. These findings thereby provide strong empirical evidence that RHM is an effective and practical complementary approach to GWAS in this context. Also, liver tissue mRNA microarray analysis revealed higher gene expression levels in ULK4 in PBC patients (P < 0.01). Lastly, we estimated the common SNP heritability of PBC in the Japanese population (0.210 ± 0.026).

Identification of genetic overlap and novel risk loci for attention-deficit/hyperactivity disorder and bipolar disorder

Differential diagnosis between childhood onset attention-deficit/hyperactivity disorder (ADHD) and bipolar disorder (BD) remains a challenge, mainly due to overlapping symptoms and high rates of comorbidity. Despite this, genetic correlation reported for these disorders is low and non-significant. Here we aimed to better characterize the genetic architecture of these disorders utilizing recent large genome-wide association studies (GWAS). We analyzed independent GWAS summary statistics for ADHD (19,099 cases and 34,194 controls) and BD (20,352 cases and 31,358 controls) applying the conditional/conjunctional false discovery rate (condFDR/conjFDR) statistical framework that increases the power to detect novel phenotype-specific and shared loci by leveraging the combined power of two GWAS. We observed cross-trait polygenic enrichment for ADHD conditioned on associations with BD, and vice versa. Leveraging this enrichment, we identified 19 novel ADHD risk loci and 40 novel BD risk loci at condFDR <0.05. Further, we identified five loci jointly associated with ADHD and BD (conjFDR < 0.05). Interestingly, these five loci show concordant directions of effect for ADHD and BD. These results highlight a shared underlying genetic risk for ADHD and BD which may help to explain the high comorbidity rates and difficulties in differentiating between ADHD and BD in the clinic. Improving our understanding of the underlying genetic architecture of these disorders may aid in the development of novel stratification tools to help reduce these diagnostic difficulties.

Moderate Exercise Inhibits Age-Related Inflammation, Liver Steatosis, Senescence, and Tumorigenesis

Age-related chronic inflammation promotes cellular senescence, chronic disease, cancer, and reduced lifespan. In this study, we wanted to explore the effects of a moderate exercise regimen on inflammatory liver disease and tumorigenesis. We used an established model of spontaneous inflammaging, steatosis, and cancer (nfkb1-/- mouse) to demonstrate whether 3 mo of moderate aerobic exercise was sufficient to suppress liver disease and cancer development. Interventional exercise when applied at a relatively late disease stage was effective at reducing tissue inflammation (liver, lung, and stomach), oxidative damage, and cellular senescence, and it reversed hepatic steatosis and prevented tumor development. Underlying these benefits were transcriptional changes in enzymes driving the conversion of tryptophan to NAD+, this leading to increased hepatic NAD+ and elevated activity of the NAD+-dependent deacetylase sirtuin. Increased SIRT activity was correlated with enhanced deacetylation of key transcriptional regulators of inflammation and metabolism, NF-κB (p65), and PGC-1α. We propose that moderate exercise can effectively reprogram pre-established inflammatory and metabolic pathologies in aging with the benefit of prevention of disease.

Mapping of susceptible variants for cold medicine-related Stevens-Johnson syndrome by whole-genome resequencing

Stevens–Johnson syndrome (SJS) and its severe condition with extensive skin detachment and a poor prognosis, toxic epidermal necrolysis (TEN), are immunologically mediated severe cutaneous reactions of the skin and mucous membranes such as the ocular surface. Genetic variations on the HLA-A and other autosomal genes have been identified as risk factors for cold medicine-related SJS/TEN with severe ocular complications (CM-SJS/TEN with SOC). Using a whole-genome sequencing (WGS) approach, we explored other susceptible variants of CM-SJS/TEN with SOC, especially among rare variants and structural variants (SVs). WGS was performed on samples from 133 patients with CM-SJS/TEN with SOC and 418 healthy controls to obtain single nucleotide polymorphisms (SNPs) and SVs. Genome-wide association tests were performed with these variants. Our genome-wide association test reproduced the associations of the common variants of HLA-A and loci on chromosome 16q12.1. We also identified novel associations of SVs on these loci and an aggregation of rare coding variants on the TPRM8 gene. In silico gene expression analysis on the HLA-A locus revealed that the SNP (rs12202296), which was significantly associated with susceptibility to CM-SJS/TEN with SOC, was correlated to an increase in HLA-A expression levels in the whole blood (P = 2.9 × 10⁻¹⁷), from the GTEx database. The majority of variants that were significantly associated with CM-SJS/TEN with SOC were found in non-coding regions, indicating the regulatory role of genetic variations in the pathogenesis of CM-SJS/TEN with SOC.