Exploring the genetic diversity of the Japanese population: Insights from a large-scale whole genome sequencing analysis.

The Japanese archipelago is a terminal location for human migration, and the contemporary Japanese people represent a unique population whose genomic diversity has been shaped by multiple migrations from Eurasia. We analyzed the genomic characteristics that define the genetic makeup of the modern Japanese population from a population genetics perspective from the genomic data of 9,287 samples obtained by high-coverage whole-genome sequencing (WGS) by the National Center Biobank Network. The dataset comprised populations from the Ryukyu Islands and other parts of the Japanese archipelago (Hondo). The Hondo population underwent two episodes of population decline during the Jomon period, corresponding to the Late Neolithic, and the Edo period, corresponding to the Early Modern era, while the Ryukyu population experienced a population decline during the shell midden period of the Late Neolithic in this region. Haplotype analysis suggested increased allele frequencies for genes related to alcohol and fatty acid metabolism, which were reported as loci that had experienced positive natural selection. Two genes related to alcohol metabolism were found to be 12,500 years out of phase with the time when they began to increase in the allele frequency; this finding indicates that the genomic diversity of Japanese people has been shaped by events closely related to agriculture and food production.

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

BACKGROUND AIMS: Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-non-specific) and non-shared (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 RESULTS: Protein tyrosine phosphatase non-receptor type 2 (PTPN2) was identified as a novel PBC susceptibility gene locus through a GWAS and subsequent genome-wide meta-analysis involving 2,181 cases and 2,699 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 (Tfh) and plasmacytoid dendritic cells (pDCs). Infiltration of PTPN2-positive T-cells and pDCs were 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 via an insufficient negative feedback loop caused by the PTPN2 risk allele of rs2292758 in IFN signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

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.

HLA Class I Analysis Provides Insight Into the Genetic and Epigenetic Background of Immune Evasion in Colorectal Cancer With High Microsatellite Instability.

BACKGROUND & AIMS: A detailed understanding of antitumor immunity is essential for optimal cancer immune therapy. Although defective mutations in the B2M and HLA-ABC genes, which encode molecules essential for antigen presentation, have been reported in several studies, the effects of these defects on tumor immunity have not been quantitatively evaluated. METHODS: Mutations in HLA-ABC genes were analyzed in 114 microsatellite instability-high colorectal cancers using a long-read sequencer. The data were further analyzed in combination with whole-exome sequencing, transcriptome sequencing, DNA methylation array, and immunohistochemistry data. RESULTS: We detected 101 truncating mutations in 57 tumors (50%) and loss of 61 alleles in 21 tumors (18%). Based on the integrated analysis that enabled the immunologic subclassification of microsatellite instability-high colorectal cancers, we identified a subtype of tumors in which lymphocyte infiltration was reduced, partly due to reduced expression of HLA-ABC genes in the absence of apparent genetic alterations. Survival time of patients with such tumors was shorter than in patients with other tumor types. Paradoxically, tumor mutation burden was highest in the subtype, suggesting that the immunogenic effect of accumulating mutations was counterbalanced by mutations that weakened immunoreactivity. Various genetic and epigenetic alterations, including frameshift mutations in RFX5 and promoter methylation of PSMB8 and HLA-A, converged on reduced expression of HLA-ABC genes. CONCLUSIONS: Our detailed immunogenomic analysis provides information that will facilitate the improvement and development of cancer immunotherapy.

Response to COVID-19 vaccine is reduced in patients with inflammatory bowel disease, but improved with additional dose.

BACKGROUND: Coronavirus disease 2019 (COVID-19) vaccination is recommended for patients with inflammatory bowel disease (IBD); however, suppressed immune responses have been reported for fully vaccinated patients under immunosuppressive therapy, mainly from Western countries. We prospectively analyzed antibody titers of IBD patients in Asia induced by two-dose and additional dose of messengerRNA COVID-19 vaccine. METHODS: After measuring high-affinity antibody titers, factors associated with antibody titers were identified by multiple regression analyses using the following covariates: sex, age (≥60 or <60 years), disease type (Crohn's disease or ulcerative colitis), vaccine type (BNT162b2 or mRNA-1273), time from second/third vaccination, molecular-targeted agent (anti-tumor necrosis factor [TNF] agents, ustekinumab, vedolizumab, tofacitinib, or no molecular-targeted agents), thiopurine, steroid, and 5-aminosalicylic acid. RESULTS: Among 409 patients analyzed, mean titer was 1316.7 U/mL (SD, 1799.3); 403 (98.5%) were judged to be seropositive (≥0.8 U/mL), and 389 (95.1%) had neutralizing antibodies (≥15 U/mL). After the third vaccination, mean titer raised up to 21 123.8 U/mL (SD, 23 474.5); all 179 were seropositive, and 178 (99.4%) had neutralizing antibodies. In 248 patients with genetic data, there was no difference in mean titer after two/third doses between carriers and non-carriers of HLA-A24 associated with severe disease during COVID-19 infection. A multiple regression analyses using covariates revealed that older age, vaccine type (BNT162b2), time from second/third dose, anti-TNF agent, tofacitinib, and thiopurine were independently associated with lower antibody titers. CONCLUSIONS: Our findings further support the recommendation for COVID-19 vaccination in patients under immunosuppressive therapy, especially additional third dose for patients receiving anti-TNF agents and/or thiopurine or tofacitinib.

Prediction Model with HLA-A*33:03 Reveals Number of Days to Develop Liver Cancer from Blood Test.

The development of liver cancer in patients with hepatitis B is a major problem, and several models have been reported to predict the development of liver cancer. However, no predictive model involving human genetic factors has been reported to date. For the items incorporated in the prediction model reported so far, we selected items that were significant in predicting liver carcinogenesis in Japanese patients with hepatitis B and constructed a prediction model of liver carcinogenesis by the Cox proportional hazard model with the addition of Human Leukocyte Antigen (HLA) genotypes. The model, which included four items-sex, age at the time of examination, alpha-fetoprotein level (log10AFP) and presence or absence of HLA-A*33:03-revealed an area under the receiver operating characteristic curve (AUROC) of 0.862 for HCC prediction within 1 year and an AUROC of 0.863 within 3 years. A 1000 repeated validation test resulted in a C-index of 0.75 or higher, or sensitivity of 0.70 or higher, indicating that this predictive model can distinguish those at high risk of developing liver cancer within a few years with high accuracy. The prediction model constructed in this study, which can distinguish between chronic hepatitis B patients who develop hepatocellular carcinoma (HCC) early and those who develop HCC late or not, is clinically meaningful.

Protective association of HLA-DPB1*04:01:01 with acute encephalopathy with biphasic seizures and late reduced diffusion identified by HLA imputation.

Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a severe syndrome of acute encephalopathy that affects infants and young children. AESD is a polygenic disorder preceded by common viral infections with high fever. We conducted an association study of human leukocyte antigen (HLA) regions with AESD using HLA imputation. SNP genotyping was performed on 254 Japanese patients with AESD and 799 healthy controls. We conducted 3-field HLA imputation for 14 HLA genes based on Japanese-specific references using data from our previous genome-wide association study. After quality control, 208 patients and 737 controls were included in the analysis of HLA alleles. We then compared the carrier frequencies of HLA alleles and haplotypes between the patients and controls. HLA-DPB1*04:01:01 showed a significant association with AESD, exerting a protective effect against the disease (p = 0.0053, pcorrected = 0.042, odds ratio = 0.43, 95% confidence interval = 0.21-0.80). The allele frequency of HLA-DPB1*04:01:01 was lower in East Asians than in Caucasians, which may partially account for the higher incidence of AESD in the Japanese population. The present results demonstrate the importance of fine-mapping of the HLA region to investigate disease susceptibilities and elucidate the pathogenesis of AESD.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs.

BACKGROUNDS & AIMS: Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening. METHODS: We combined new and existing genotype data for 10,516 cases and 20,772 controls from 5 European and 2 East Asian cohorts. RESULTS: We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are well-established in the treatment of other autoimmune disorders. CONCLUSIONS: This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders. LAY SUMMARY: Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10,516 people with PBC and 20,772 healthy individuals recruited in Canada, China, Italy, Japan, the UK, or the USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these 'candidate genes' to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC.

Common risk variants in NPHS1 and TNFSF15 are associated with childhood steroid-sensitive nephrotic syndrome.

To understand the genetics of steroid-sensitive nephrotic syndrome (SSNS), we conducted a genome-wide association study in 987 childhood SSNS patients and 3,206 healthy controls with Japanese ancestry. Beyond known associations in the HLA-DR/DQ region, common variants in NPHS1-KIRREL2 (rs56117924, P=4.94E-20, odds ratio (OR) =1.90) and TNFSF15 (rs6478109, P=2.54E-8, OR=0.72) regions achieved genome-wide significance and were replicated in Korean, South Asian and African populations. Trans-ethnic meta-analyses including Japanese, Korean, South Asian, African, European, Hispanic and Maghrebian populations confirmed the significant associations of variants in NPHS1-KIRREL2 (Pmeta=6.71E-28, OR=1.88) and TNFSF15 (Pmeta=5.40E-11, OR=1.33) loci. Analysis of the NPHS1 risk alleles with glomerular NPHS1 mRNA expression from the same person revealed allele specific expression with significantly lower expression of the transcript derived from the risk haplotype (Wilcox test p=9.3E-4). Because rare pathogenic variants in NPHS1 cause congenital nephrotic syndrome of the Finnish type (CNSF), the present study provides further evidence that variation along the allele frequency spectrum in the same gene can cause or contribute to both a rare monogenic disease (CNSF) and a more complex, polygenic disease (SSNS).

Amino acid position 37 of HLA-DRß1 affects susceptibility to Crohn's disease in Asians.

Crohn's disease (CD) and ulcerative colitis (UC) are the major types of chronic inflammatory bowel disease (IBD) characterized by recurring episodes of inflammation of the gastrointestinal tract. Although it is well established that human leukocyte antigen (HLA) is a major risk factor for IBD, it is yet to be determined which HLA alleles or amino acids drive the risks of CD and UC in Asians. To define the roles of HLA for IBD in Asians, we fine-mapped HLA in 12 568 individuals from Korea and Japan (3294 patients with CD, 1522 patients with UC and 7752 controls). We identified that the amino acid position 37 of HLA-DRß1 plays a key role in the susceptibility to CD (presence of serine being protective, P = 3.6 × 10-67, OR = 0.48 [0.45-0.52]). For UC, we confirmed the known association of the haplotype spanning HLA-C*12:02, HLA-B*52:01 and HLA-DRB1*1502 (P = 1.2 × 10-28, OR = 4.01 [3.14-5.12]).

A high-resolution HLA imputation system for the Taiwanese population: a study of the Taiwan Biobank.

An imputation algorithm for human leukocyte antigen (HLA) is helpful for exploring novel disease associations. However, population-specific HLA imputation references are essential for achieving high imputation accuracy. In this study, a subset of 1012 individuals from the Taiwan Biobank (TWB) who underwent both whole-genome SNP array and NGS-based HLA typing were used to establish Taiwanese HLA imputation references. The HIBAG package was used to generate the imputation references for eight HLA loci at a two- and three-field resolution. Internal validation was carried out to evaluate the call threshold and accuracy for each HLA gene. HLA class II genes found to be associated with rheumatoid arthritis (RA) were validated in this study by the imputed HLA alleles. Our Taiwanese population-specific references achieved average HLA imputation accuracies of 98.11% for two-field and 98.08% for three-field resolution. The frequency distribution of imputed HLA alleles among 23,972 TWB subjects were comparable with PCR-based HLA alleles in general Taiwanese reported in the allele frequency net database. We replicated four common HLA alleles (HLA-DRB1*03:01, DRB1*04:05, DQA1*03:03, and DQB1*04:01) significantly associated with RA. The population-specific references provide an informative tool to investigate the associations of HLA variants and human diseases in large-scale population-based studies.

Human leukocyte antigen variants and risk of hepatocellular carcinoma modified by hepatitis C virus genotypes: A genome-wide association study.

We conducted a genome-wide association study to discover genetic variants associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC). We genotyped 502 HCC cases and 749 non-HCC controls using the Axiom-CHB genome-wide array. After identifying single-nucleotide polymorphism clusters located in the human leukocyte antigen (HLA) region which were potentially associated with HCC, HLA-DQB1 genotyping was performed to analyze 994 anti-HCV seropositives collected in the period 1991-2013 in a community-based cohort for evaluating long-term predictability of HLA variants for identifying the risk of HCC. Cox proportional hazards models were used to estimate the hazard ratios and 95% confidence intervals of HLA genotypes for determining the aforementioned HCC risk. Eight single-nucleotide polymorphisms in the proximity of HLA-DQB1 were associated with HCC (P < 8.7 × 10-8 ) in the genome-wide association study. Long-term follow-up showed a significant association with HLA-DQB1*03:01 and DQB1*06:02 (P < 0.05). The adjusted hazard ratios associated with HCC were 0.45 (0.30-0.68) and 2.11 (1.34-3.34) for DQB1*03:01 and DQB1*06:02, respectively. After stratification by HCV genotypes, DQB1*03:01 showed protective effects only in patients with HCV genotype 1, whereas DQB1*06:02 conferred risk of HCC only in patients with HCV non-1 genotypes. HLA imputation analyses revealed that HLA-DRB1*15:01, which is in linkage disequilibrium with DQB1*06:02, also increased the risk of HCC (odds ratio, 1.96; 95% confidence interval, 1.31-2.93). Haplotype analysis supported that DQB1*03:01 and DQB1*06:02 are primarily protective and susceptible variants, respectively. Conclusion: HLA-DQB1 was independently associated with HCC; HCV genotypes modified the effects of HLA-DQB1 on the risk of HCC. (Hepatology 2018;67:651-661).

Key HLA-DRB1-DQB1 haplotypes and role of the BTNL2 gene for response to a hepatitis B vaccine.

Approximately 5-10% of individuals who are vaccinated with a hepatitis B (HB) vaccine designed based on the hepatitis B virus (HBV) genotype C fail to acquire protective levels of antibodies. Here, host genetic factors behind low immune response to this HB vaccine were investigated by a genome-wide association study (GWAS) and Human Leukocyte Antigen (HLA) association tests. The GWAS and HLA association tests were carried out using a total of 1,193 Japanese individuals including 107 low responders, 351 intermediate responders, and 735 high responders. Classical HLA class II alleles were statistically imputed using the genome-wide SNP typing data. The GWAS identified independent associations of HLA-DRB1-DQB1, HLA-DPB1 and BTNL2 genes with immune response to a HB vaccine designed based on the HBV genotype C. Five HLA-DRB1-DQB1 haplotypes and two DPB1 alleles showed significant associations with response to the HB vaccine in a comparison of three groups of 1,193 HB vaccinated individuals. When frequencies of DRB1-DQB1 haplotypes and DPB1 alleles were compared between low immune responders and HBV patients, significant associations were identified for three DRB1-DQB1 haplotypes, and no association was identified for any of the DPB1 alleles. In contrast, no association was identified for DRB1-DQB1 haplotypes and DPB1 alleles in a comparison between high immune responders and healthy individuals. Conclusion: The findings in this study clearly show the importance of HLA-DR-DQ (i.e., recognition of a vaccine related HB surface antigen (HBsAg) by specific DR-DQ haplotypes) and BTNL2 molecules (i.e., high immune response to HB vaccine) for response to a HB vaccine designed based on the HBV genotype C. (Hepatology 2018).

Strong Association of the HLA-DR/DQ Locus with Childhood Steroid-Sensitive Nephrotic Syndrome in the Japanese Population.

Background Nephrotic syndrome is the most common cause of chronic glomerular disease in children. Most of these patients develop steroid-sensitive nephrotic syndrome (SSNS), but the loci conferring susceptibility to childhood SSNS are mainly unknown.Methods We conducted a genome-wide association study (GWAS) in the Japanese population; 224 patients with childhood SSNS and 419 adult healthy controls were genotyped using the Affymetrix Japonica Array in the discovery stage. Imputation for six HLA genes (HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1) was conducted on the basis of Japanese-specific references. We performed genotyping for HLA-DRB1/-DQB1 using a sequence-specific oligonucleotide-probing method on a Luminex platform. Whole-genome imputation was conducted using a phased reference panel of 2049 healthy Japanese individuals. Replication was performed in an independent Japanese sample set including 216 patients and 719 healthy controls. We genotyped candidate single-nucleotide polymorphisms using the DigiTag2 assay.Results The most significant association was detected in the HLA-DR/DQ region and replicated (rs4642516 [minor allele G], combined Pallelic=7.84×10-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined Pallelic=1.72×10-25; OR, 0.29; 95% CI, 0.23 to 0.37). HLA-DRB1*08:02 (Pc=1.82×10-9; OR, 2.62; 95% CI, 1.94 to 3.54) and HLA-DQB1*06:04 (Pc=2.09×10-12; OR, 0.10; 95% CI, 0.05 to 0.21) were considered primary HLA alleles associated with childhood SSNS. HLA-DRB1*08:02-DQB1*03:02 (Pc=7.01×10-11; OR, 3.60; 95% CI, 2.46 to 5.29) was identified as the most significant genetic susceptibility factor.Conclusions The most significant association with childhood SSNS was detected in the HLA-DR/DQ region. Further HLA allele/haplotype analyses should enhance our understanding of molecular mechanisms underlying SSNS.

A variant at 9q34.11 is associated with HLA-DQB1*06:02 negative essential hypersomnia.

Essential hypersomnia (EHS) is a lifelong disorder characterized by excessive daytime sleepiness without cataplexy. EHS is associated with human leukocyte antigen (HLA)-DQB1*06:02, similar to narcolepsy with cataplexy (narcolepsy). Previous studies suggest that DQB1*06:02-positive and -negative EHS are different in terms of their clinical features and follow different pathological pathways. DQB1*06:02-positive EHS and narcolepsy share the same susceptibility genes. In the present study, we report a genome-wide association study with replication for DQB1*06:02-negative EHS (408 patients and 2247 healthy controls, all Japanese). One single-nucleotide polymorphism, rs10988217, which is located 15-kb upstream of carnitine O-acetyltransferase (CRAT), was significantly associated with DQB1*06:02-negative EHS (P = 7.5 × 10-9, odds ratio = 2.63). The risk allele of the disease-associated SNP was correlated with higher expression levels of CRAT in various tissues and cell types, including brain tissue. In addition, the risk allele was associated with levels of succinylcarnitine (P = 1.4 × 10-18) in human blood. The leading SNP in this region was the same in associations with both DQB1*06:02-negative EHS and succinylcarnitine levels. The results suggest that DQB1*06:02-negative EHS may be associated with an underlying dysfunction in energy metabolic pathways.

New susceptibility variants to narcolepsy identified in HLA class II region.

Narcolepsy, a sleep disorder characterized by excessive daytime sleepiness, cataplexy and rapid eye movement sleep abnormalities, is tightly associated with human leukocyte antigen HLA-DQB1*06:02. DQB1*06:02 is common in the general population (10-30%); therefore, additional genetic factors are needed for the development of narcolepsy. In the present study, HLA-DQB1 in 664 Japanese narcoleptic subjects and 3131 Japanese control subjects was examined to determine whether HLA-DQB1 alleles located in trans of DQB1*06:02 are associated with narcolepsy. The strongest association was with DQB1*06:01 (P = 1.4 × 10(-10), odds ratio, OR = 0.39), as reported in previous studies. Additional predisposing effects of DQB1*03:02 were also found (P = 2.5 × 10(-9), OR = 1.97). A comparison between DQB1*06:02 heterozygous cases and controls revealed dominant protective effects of DQB1*06:01 and DQB1*05:01. In addition, a single-nucleotide polymorphism-based conditional analysis controlling for the effect of HLA-DQB1 was performed to determine whether there were other independent HLA associations outside of HLA-DQB1. This analysis revealed associations at HLA-DPB1 in the HLA class II region (rs3117242, P = 4.1 × 10(-5), OR = 2.45; DPB1*05:01, P = 8.1 × 10(-3), OR = 1.39). These results indicate that complex HLA class II associations contribute to the genetic predisposition to narcolepsy.

The first genome-wide association study identifying new susceptibility loci for obstetric antiphospholipid syndrome.

Antiphospholipid syndrome (APS) is the most important treatable cause of recurrent pregnancy loss. The live birth rate is limited to only 70-80% in patients with APS undergoing established anticoagulant therapy. Lupus anticoagulant (LA), but not anticardiolipin antibody (aCL), was found to predict adverse pregnancy outcome. Recent genome-wide association studies (GWAS) of APS focusing on aCL have shown that several molecules may be involved. This is the first GWAS for obstetric APS focusing on LA. A GWAS was performed to compare 115 Japanese patients with obstetric APS, diagnosed according to criteria of the International Congress on APS, and 419 healthy individuals. Allele or genotype frequencies were compared in a total of 426 344 single-nucleotide polymorphisms (SNPs). Imputation analyses were also performed for the candidate regions detected by the GWAS. One SNP (rs2288493) located on the 3'-UTR of TSHR showed an experiment-wide significant APS association (P=7.85E-08, OR=6.18) under a recessive model after Bonferroni correction considering the number of analyzed SNPs. Another SNP (rs79154414) located around the C1D showed a genome-wide significant APS association (P=4.84E-08, OR=6.20) under an allelic model after applying the SNP imputation. Our findings demonstrate that a specific genotype of TSHR and C1D genes can be a risk factor for obstetric APS.

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

BACKGROUND: Over 160 disease phenotypes have been mapped to the major histocompatibility complex (MHC) region on chromosome 6 by genome-wide association study (GWAS), suggesting that the MHC region as a whole may be involved in the aetiology of many phenotypes, including unstudied diseases. The phenome-wide association study (PheWAS), a powerful and complementary approach to GWAS, has demonstrated its ability to discover and rediscover genetic associations. The objective of this study is to comprehensively investigate the MHC region by PheWAS to identify new phenotypes mapped to this genetically important region. METHODS: In the current study, we systematically explored the MHC region using PheWAS to associate 2692 MHC-linked variants (minor allele frequency ≥0.01) with 6221 phenotypes in a cohort of 7481 subjects from the Marshfield Clinic Personalized Medicine Research Project. RESULTS: Findings showed that expected associations previously identified by GWAS could be identified by PheWAS (eg, psoriasis, ankylosing spondylitis, type I diabetes and coeliac disease) with some having strong cross-phenotype associations potentially driven by pleiotropic effects. Importantly, novel associations with eight diseases not previously assessed by GWAS (eg, lichen planus) were also identified and replicated in an independent population. Many of these associated diseases appear to be immune-related disorders. Further assessment of these diseases in 16 484 Marshfield Clinic twins suggests that some of these diseases, including lichen planus, may have genetic aetiologies. CONCLUSIONS: These results demonstrate that the PheWAS approach is a powerful and novel method to discover SNP-disease associations, and is ideal when characterising cross-phenotype associations, and further emphasise the importance of the MHC region in human health and disease.