Genome-wide association analysis reveals the associations of NPHP4, TYW1-AUTS2 and SEMA6D for Behçet's disease and HLA-B*46:01 for its intestinal involvement.

BACKGROUND: Intestinal involvement in Behçet's disease (BD) is associated with poor prognosis and is more prevalent in East Asian than in Mediterranean populations. Identifying the genetic causes of intestinal BD is important for understanding the pathogenesis and for appropriate treatment of BD patients. METHODS: We performed genome-wide association studies (GWAS) and imputation/replication genotyping of human leukocyte antigen (HLA) alleles for 1,689 Korean and Turkish patients with BD (including 379 patients with intestinal BD) and 2,327 healthy controls, followed by replication using 593 Japanese patients with BD (101 patients with intestinal BD) and 737 healthy controls. Stratified cross-phenotype analyses were performed for 1) overall BD, 2) intestinal BD, and 3) intestinal BD without association of overall BD. RESULTS: We identified three novel genome-wide significant susceptibility loci including NPHP4 (rs74566205; P=1.36 × 10-8), TYW1-AUTS2 (rs60021986; P=1.14 × 10-9), and SEMA6D (rs4143322; P=5.54 × 10-9) for overall BD, and a new association with HLA-B*46:01 for intestinal BD (P=1.67 × 10-8) but not for BD without intestinal involvement. HLA peptide binding analysis revealed that Mycobacterial peptides, have a stronger binding affinity to HLA-B*46:01 compared to the known risk allele HLA-B*51:01. CONCLUSIONS: HLA-B*46:01 is associated with the development of intestinal BD; NPHP4, TYW1-AUTS2, and SEMA6D are susceptibility loci for overall BD.

The interleukin-11 receptor variant p.W307R results in craniosynostosis in humans.

Craniosynostosis is characterized by the premature fusion and ossification of one or more of the sutures of the calvaria, often resulting in abnormal features of the face and the skull. In cases in which growth of the brain supersedes available space within the skull, developmental delay or cognitive impairment can occur. A complex interplay of different cell types and multiple signaling pathways are required for correct craniofacial development. In this study, we report on two siblings with craniosynostosis and a homozygous missense pathogenic variant within the IL11RA gene (c.919 T > C; p.W307R). The patients present with craniosynostosis, exophthalmos, delayed tooth eruption, mild platybasia, and a basilar invagination. The p.W307R variant is located within the arginine-tryptophan-zipper within the D3 domain of the IL-11R, a structural element known to be important for the stability of the cytokine receptor. Expression of IL-11R-W307R in cells shows impaired maturation of the IL-11R, no transport to the cell surface and intracellular retention. Accordingly, cells stably expressing IL-11R-W307R do not respond when stimulated with IL-11, arguing for a loss-of-function mutation. In summary, the IL-11R-W307R variant, reported here for the first time to our knowledge, is most likely the causative variant underlying craniosynostosis in these patients.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants.

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

Shared Genetics of Multiple System Atrophy and Inflammatory Bowel Disease.

BACKGROUND: Multiple system atrophy (MSA) is a rare neurodegenerative disease characterized by intracellular accumulations of α-synuclein and nerve cell loss in striatonigral and olivopontocerebellar structures. Epidemiological and clinical studies have reported potential involvement of autoimmune mechanisms in MSA pathogenesis. However, genetic etiology of this interaction remains unknown. We aimed to investigate genetic overlap between MSA and 7 autoimmune diseases and to identify shared genetic loci. METHODS: Genome-wide association study summary statistics of MSA and 7 autoimmune diseases were combined in cross-trait conjunctional false discovery rate analysis to explore overlapping genetic background. Expression of selected candidate genes was compared in transgenic MSA mice and wild-type mice. Genetic variability of candidate genes was further investigated using independent whole-exome genotyping data from large cohorts of MSA and autoimmune disease patients and healthy controls. RESULTS: We observed substantial polygenic overlap between MSA and inflammatory bowel disease and identified 3 shared genetic loci with leading variants upstream of the DENND1B and RSP04 genes, and in intron of the C7 gene. Further, the C7 gene showed significantly dysregulated expression in the degenerating midbrain of transgenic MSA mice compared with wild-type mice and had elevated burden of protein-coding variants in independent MSA and inflammatory bowel disease cohorts. CONCLUSION: Our study provides evidence of shared genetic etiology between MSA and inflammatory bowel disease with an important role of the C7 gene in both phenotypes, with the implication of immune and gut dysfunction in MSA pathophysiology. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The Macrophage Migration Inhibitory Factor (MIF) Promoter Polymorphisms (rs3063368, rs755622) Predict Acute Kidney Injury and Death after Cardiac Surgery.

BACKGROUND: Macrophage Migration Inhibitory Factor (MIF) is highly elevated after cardiac surgery and impacts the postoperative inflammation. The aim of this study was to analyze whether the polymorphisms CATT5-7 (rs5844572/rs3063368,"-794") and G>C single-nucleotide polymorphism (rs755622,-173) in the MIF gene promoter are related to postoperative outcome. METHODS: In 1116 patients undergoing cardiac surgery, the MIF gene polymorphisms were analyzed and serum MIF was measured by ELISA in 100 patients. RESULTS: Patients with at least one extended repeat allele (CATT7) had a significantly higher risk of acute kidney injury (AKI) compared to others (23% vs. 13%; OR 2.01 (1.40-2.88), p = 0.0001). Carriers of CATT7 were also at higher risk of death (1.8% vs. 0.4%; OR 5.12 (0.99-33.14), p = 0.026). The GC genotype was associated with AKI (20% vs. GG/CC:13%, OR 1.71 (1.20-2.43), p = 0.003). Multivariate analyses identified CATT7 predictive for AKI (OR 2.13 (1.46-3.09), p < 0.001) and death (OR 5.58 (1.29-24.04), p = 0.021). CATT7 was associated with higher serum MIF before surgery (79.2 vs. 50.4 ng/mL, p = 0.008). CONCLUSION: The CATT7 allele associates with a higher risk of AKI and death after cardiac surgery, which might be related to chronically elevated serum MIF. Polymorphisms in the MIF gene may constitute a predisposition for postoperative complications and the assessment may improve risk stratification and therapeutic guidance.

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression.

BACKGROUND: Fifteen percent of atopic dermatitis (AD) liability-scale heritability could be attributed to 31 susceptibility loci identified by using genome-wide association studies, with only 3 of them (IL13, IL-6 receptor [IL6R], and filaggrin [FLG]) resolved to protein-coding variants. OBJECTIVE: We examined whether a significant portion of unexplained AD heritability is further explained by low-frequency and rare variants in the gene-coding sequence. METHODS: We evaluated common, low-frequency, and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 control subjects combined with whole-transcriptome data on lesional, nonlesional, and healthy skin samples of 27 patients and 38 control subjects. RESULTS: An additional 12.56% (SE, 0.74%) of AD heritability is explained by rare protein-coding variation. We identified docking protein 2 (DOK2) and CD200 receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with AD are further enriched in 5 genes (IL-4 receptor [IL4R], IL13, Janus kinase 1 [JAK1], JAK2, and tyrosine kinase 2 [TYK2]) of the IL13 pathway, all of which are targets for novel systemic AD therapeutics. Multiomics-based network and RNA sequencing analysis revealed DOK2 as a central hub interacting with, among others, CD200R1, IL6R, and signal transducer and activator of transcription 3 (STAT3). Multitissue gene expression profile analysis for 53 tissue types from the Genotype-Tissue Expression project showed that disease-associated protein-coding variants exert their greatest effect in skin tissues. CONCLUSION: Our discoveries highlight a major role of rare coding variants in AD acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of the novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility.

ZNF133 is associated with infliximab responsiveness in patients with inflammatory bowel diseases.

BACKGROUND AND AIM: Infliximab has been widely prescribed for treating inflammatory bowel disease (IBD). However, the response rates to infliximab differ among patients. Therefore, we aimed to identify the genetic and clinical markers that predict infliximab response. METHODS: A total of 139 Korean patients with IBD who received infliximab were classified according to infliximab response as follows: (i) primary response vs nonresponse and (ii) sustained response vs loss of response. We performed an association study using whole-exome sequencing data to identify genetic variants associated with infliximab response. Candidate variants were validated in 77 German patients with IBD. Stepwise multivariate logistic regression was performed to identify predictors. RESULTS: We found five candidate variants that were associated with primary nonresponse to infliximab (P < 5 × 10-6 ). Of the five variants, rs2228273 in ZNF133 was validated in German (combined P = 6.49 × 10-7 ). We also identified the best genetic variant (rs9144, P = 4.60 × 10-6 ) associated with the loss of infliximab response. In multivariate regression analysis, rs2228273 (P = 2.10 × 10-5 ), concurrent azathioprine/6-mercaptopurine use, and bodyweight at the first infliximab use (< 50 kg) were associated with primary nonresponse. In addition, the Crohn's disease activity index at the first infliximab use and rs9144 (P = 0.001) were independently associated with the loss of response in patients with Crohn's disease. CONCLUSIONS: We identified clinical and genetic markers associated with infliximab response in IBD patients. Our findings could provide insights to maximize the efficacy of infliximab therapy in IBD patients.

Genetic markers associated with long-term cardiovascular outcome in kidney transplant recipients.

There is a clear genetic contribution to the risk of cardiovascular diseases, and a composite genetic risk score (GRS) based on 27 single nucleotide polymorphisms (SNPs) was reported to predict risk of cardiovascular events in the general population. We aimed to evaluate this risk score in renal transplant recipients, a population with heightened cardiovascular risk, with a yet unknown genetic contribution. A total of 1640 participants from the ALERT trial (Assessment of Lescol in Renal Transplantation), a study comparing fluvastatin with placebo in stable renal transplant recipients, were genotyped for all SNPs making up the GRS. Risk alleles were weighted by the log of odds ratios reported in genome wide association studies and summed. Associations between GRS and time from study inclusion to first major cardiovascular event (MACE) were analyzed by Cox regression. In analyses adjusted for cardiovascular risk factors, GRS was significantly associated with MACE (hazard ratio [HR] 1.81, P = .006) when comparing genetic high-risk patients (quartile 4) with genetic low-risk participants (quartile 1). A 27-SNP GRS, which predicted cardiovascular events in the nontransplant population, appears to have predictive value also in kidney allograft recipients. Refining the score to better fit the transplant population seems feasible.

Exome-wide association study reveals novel psoriasis susceptibility locus at TNFSF15 and rare protective alleles in genes contributing to type I IFN signalling.

Psoriasis is a common inflammatory skin disorder for which multiple genetic susceptibility loci have been identified, but few resolved to specific functional variants. In this study, we sought to identify common and rare psoriasis-associated gene-centric variation. Using exome arrays we genotyped four independent cohorts, totalling 11 861 psoriasis cases and 28 610 controls, aggregating the dataset through statistical meta-analysis. Single variant analysis detected a previously unreported risk locus at TNFSF15 (rs6478108; P = 1.50 × 10-8, OR = 1.10), and association of common protein-altering variants at 11 loci previously implicated in psoriasis susceptibility. We validate previous reports of protective low-frequency protein-altering variants within IFIH1 (encoding an innate antiviral receptor) and TYK2 (encoding a Janus kinase), in each case establishing a further series of protective rare variants (minor allele frequency < 0.01) via gene-wide aggregation testing (IFIH1: pburden = 2.53 × 10-7, OR = 0.707; TYK2: pburden = 6.17 × 10-4, OR = 0.744). Both genes play significant roles in type I interferon (IFN) production and signalling. Several of the protective rare and low-frequency variants in IFIH1 and TYK2 disrupt conserved protein domains, highlighting potential mechanisms through which their effect may be exerted.

Large scale meta-analysis characterizes genetic architecture for common psoriasis associated variants.

Psoriasis is a complex disease of skin with a prevalence of about 2%. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for psoriasis to date, including data from eight different Caucasian cohorts, with a combined effective sample size >39,000 individuals. We identified 16 additional psoriasis susceptibility loci achieving genome-wide significance, increasing the number of identified loci to 63 for European-origin individuals. Functional analysis highlighted the roles of interferon signalling and the NFκB cascade, and we showed that the psoriasis signals are enriched in regulatory elements from different T cells (CD8+ T-cells and CD4+ T-cells including TH0, TH1 and TH17). The identified loci explain ∼28% of the genetic heritability and generate a discriminatory genetic risk score (AUC=0.76 in our sample) that is significantly correlated with age at onset (p=2 × 10-89). This study provides a comprehensive layout for the genetic architecture of common variants for psoriasis.

Genome-wide association study of primary sclerosing cholangitis identifies new risk loci and quantifies the genetic relationship with inflammatory bowel disease.

Primary sclerosing cholangitis (PSC) is a rare progressive disorder leading to bile duct destruction; ∼75% of patients have comorbid inflammatory bowel disease (IBD). We undertook the largest genome-wide association study of PSC (4,796 cases and 19,955 population controls) and identified four new genome-wide significant loci. The most associated SNP at one locus affects splicing and expression of UBASH3A, with the protective allele (C) predicted to cause nonstop-mediated mRNA decay and lower expression of UBASH3A. Further analyses based on common variants suggested that the genome-wide genetic correlation (rG) between PSC and ulcerative colitis (UC) (rG = 0.29) was significantly greater than that between PSC and Crohn's disease (CD) (rG = 0.04) (P = 2.55 × 10-15). UC and CD were genetically more similar to each other (rG = 0.56) than either was to PSC (P < 1.0 × 10-15). Our study represents a substantial advance in understanding of the genetics of PSC.

CD4+ T cells from patients with primary sclerosing cholangitis exhibit reduced apoptosis and down-regulation of proapoptotic Bim in peripheral blood.

The pathogenesis of the progressive liver disease, primary sclerosing cholangitis (PSC), remains largely elusive. The strong genetic association with HLA loci suggests that T cell-dependent, adaptive immune reactions could contribute to disease pathogenesis. Recent studies have indicated that PSC is also associated with polymorphisms in the locus encoding for proapoptotic Bim (BCL2L11). Bim is crucial for the maintenance of immunologic tolerance through induction of apoptosis in activated T cells. Of interest with regard to PSC is the finding that BCL2L11-deficient mice develop periductular infiltrates. We, therefore, investigated, whether defective apoptosis of T cells might contribute to the phenotype of PSC. Thus, we induced apoptosis of T cells from patients with PSC and controls by repeated T cell receptor (TCR) stimulation or cytokine withdrawal. We found that CD4+ T cells, but not CD8+ T cells, from patients with PSC exhibited significantly reduced apoptosis in response to both, TCR restimulation or cytokine withdrawal. This increased apoptosis resistance was associated with significantly reduced up-regulation of proapoptotic Bim in T cells from patients with PSC. However, T cell apoptosis did not seem to be influenced by the previously described BCL2L11 polymorphisms. Reduced CD4+ T cell apoptosis in patients with PSC was not due to reduced cell activation, as indicated by a similar surface expression of the activation markers CD69, CD25, and CD28 in T cells from patients and controls. Thus, decreased apoptosis of activated CD4+ T cells may be part of the immune dysregulation observed in patients with PSC.

A protein-truncating R179X variant in RNF186 confers protection against ulcerative colitis.

Protein-truncating variants protective against human disease provide in vivo validation of therapeutic targets. Here we used targeted sequencing to conduct a search for protein-truncating variants conferring protection against inflammatory bowel disease exploiting knowledge of common variants associated with the same disease. Through replication genotyping and imputation we found that a predicted protein-truncating variant (rs36095412, p.R179X, genotyped in 11,148 ulcerative colitis patients and 295,446 controls, MAF=up to 0.78%) in RNF186, a single-exon ring finger E3 ligase with strong colonic expression, protects against ulcerative colitis (overall P=6.89 × 10(-7), odds ratio=0.30). We further demonstrate that the truncated protein exhibits reduced expression and altered subcellular localization, suggesting the protective mechanism may reside in the loss of an interaction or function via mislocalization and/or loss of an essential transmembrane domain.

Pooled Resequencing of 122 Ulcerative Colitis Genes in a Large Dutch Cohort Suggests Population-Specific Associations of Rare Variants in MUC2.

Genome-wide association studies have revealed several common genetic risk variants for ulcerative colitis (UC). However, little is known about the contribution of rare, large effect genetic variants to UC susceptibility. In this study, we performed a deep targeted re-sequencing of 122 genes in Dutch UC patients in order to investigate the contribution of rare variants to the genetic susceptibility to UC. The selection of genes consists of 111 established human UC susceptibility genes and 11 genes that lead to spontaneous colitis when knocked-out in mice. In addition, we sequenced the promoter regions of 45 genes where known variants exert cis-eQTL-effects. Targeted pooled re-sequencing was performed on DNA of 790 Dutch UC cases. The Genome of the Netherlands project provided sequence data of 500 healthy controls. After quality control and prioritization based on allele frequency and pathogenicity probability, follow-up genotyping of 171 rare variants was performed on 1021 Dutch UC cases and 1166 Dutch controls. Single-variant association and gene-based analyses identified an association of rare variants in the MUC2 gene with UC. The associated variants in the Dutch population could not be replicated in a German replication cohort (1026 UC cases, 3532 controls). In conclusion, this study has identified a putative role for MUC2 on UC susceptibility in the Dutch population and suggests a population-specific contribution of rare variants to UC.

Vy-PER: eliminating false positive detection of virus integration events in next generation sequencing data.

Several pathogenic viruses such as hepatitis B and human immunodeficiency viruses may integrate into the host genome. These virus/host integrations are detectable using paired-end next generation sequencing. However, the low number of expected true virus integrations may be difficult to distinguish from the noise of many false positive candidates. Here, we propose a novel filtering approach that increases specificity without compromising sensitivity for virus/host chimera detection. Our detection pipeline termed Vy-PER (Virus integration detection bY Paired End Reads) outperforms existing similar tools in speed and accuracy. We analysed whole genome data from childhood acute lymphoblastic leukemia (ALL), which is characterised by genomic rearrangements and usually associated with radiation exposure. This analysis was motivated by the recently reported virus integrations at genomic rearrangement sites and association with chromosomal instability in liver cancer. However, as expected, our analysis of 20 tumour and matched germline genomes from ALL patients finds no significant evidence for integrations by known viruses. Nevertheless, our method eliminates 12,800 false positives per genome (80× coverage) and only our method detects singleton human-phiX174-chimeras caused by optical errors of the Illumina HiSeq platform. This high accuracy is useful for detecting low virus integration levels as well as non-integrated viruses.

Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility.

Psoriasis is a common inflammatory skin disease with complex genetics and different degrees of prevalence across ethnic populations. Here we present the largest trans-ethnic genome-wide meta-analysis (GWMA) of psoriasis in 15,369 cases and 19,517 controls of Caucasian and Chinese ancestries. We identify four novel associations at LOC144817, COG6, RUNX1 and TP63, as well as three novel secondary associations within IFIH1 and IL12B. Fine-mapping analysis of MHC region demonstrates an important role for all three HLA class I genes and a complex and heterogeneous pattern of HLA associations between Caucasian and Chinese populations. Further, trans-ethnic comparison suggests population-specific effect or allelic heterogeneity for 11 loci. These population-specific effects contribute significantly to the ethnic diversity of psoriasis prevalence. This study not only provides novel biological insights into the involvement of immune and keratinocyte development mechanism, but also demonstrates a complex and heterogeneous genetic architecture of psoriasis susceptibility across ethnic populations.

Psoriasis and cardiometabolic traits: modest association but distinct genetic architectures.

Psoriasis has been linked to cardiometabolic diseases, but epidemiological findings are inconsistent. We investigated the association between psoriasis and cardiometabolic outcomes in a German cross-sectional study (n=4,185) and a prospective cohort of German Health Insurance beneficiaries (n=1,811,098). A potential genetic overlap was explored using genome-wide data from >22,000 coronary artery disease and >4,000 psoriasis cases, and with a dense genotyping study of cardiometabolic risk loci on 927 psoriasis cases and 3,717 controls. After controlling for major confounders, in the cross-sectional analysis psoriasis was significantly associated with type 2 diabetes (T2D, adjusted odds ratio (OR)=2.36; 95% confidence interval CI=1.26-4.41) and myocardial infarction (MI, OR=2.26; 95% CI=1.03-4.96). In the longitudinal study, psoriasis slightly increased the risk for incident T2D (adjusted relative risk (RR)=1.11; 95% CI=1.08-1.14) and MI (RR=1.14; 95% CI=1.06-1.22), with highest risk increments in systemically treated psoriasis, which accounted for 11 and 17 excess cases of T2D and MI per 10,000 person-years. Except for weak signals from within the major histocompatibility complex, there was no evidence of genetic risk loci shared between psoriasis and cardiometabolic traits. Our findings suggest that psoriasis, in particular severe psoriasis, increases the risk for T2D and MI, and that the genetic architecture of psoriasis and cardiometabolic traits is largely distinct.