Rare,Low-Frequency,and Common Variants in the Protein-Coding Sequence of Biological Candidate Genes from GWASs Contribute to Risk of Rheumatoid Arthritis

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p_enrichment = 6.4 × 10⁻⁴). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10⁻⁶). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA.     

Low-Frequency Coding Variants at 6p21.33 and 20q11.21 Are Associated with Lung Cancer Risk in Chinese Populations

Genome-wide association studies have successfully identified a subset of common variants associated with lung cancer risk. However, these variants explain only a fraction of lung cancer heritability. It has been proposed that low-frequency or rare variants might have strong effects and contribute to the missing heritability. To assess the role of low-frequency or rare variants in lung cancer development, we analyzed exome chips representing 1,348 lung cancer subjects and 1,998 control subjects during the discovery stage and subsequently evaluated promising associations in an additional 4,699 affected subjects and 4,915 control subjects during the replication stages. Single-variant and gene-based analyses were carried out for coding variants with a minor allele frequency less than 0.05. We identified three low-frequency missense variants in BAT2 (rs9469031, c.1544C>T [p.Pro515Leu]; odds ratio [OR] = 0.55, p = 1.28 × 10⁻¹⁰), FKBPL (rs200847762, c.410C>T [p.Pro137Leu]; OR = 0.25, p = 9.79 × 10⁻¹²), and BPIFB1 (rs6141383, c.850G>A [p.Val284Met]; OR = 1.72, p = 1.79 × 10⁻⁷); these variants were associated with lung cancer risk. rs9469031 in BAT2 and rs6141383 in BPIFB1 were also associated with the age of onset of lung cancer (p = 0.001 and 0.006, respectively). BAT2 and FKBPL at 6p21.33 and BPIFB1 at 20q11.21 were differentially expressed in lung tumors and paired normal tissues. Gene-based analysis revealed that FKBPL, in which two independent variants were identified, might account for the association with lung cancer risk at 6p21.33. Our results highlight the important role low-frequency variants play in lung cancer susceptibility and indicate that candidate genes at 6p21.33 and 20q11.21 are potentially biologically relevant to lung carcinogenesis.     

Imputation of Exome Sequence Variants into Population- Based Samples and Blood-Cell-Trait-Associated Loci in African Americans: NHLBI GO Exome Sequencing Project

Researchers have successfully applied exome sequencing to discover causal variants in selected individuals with familial, highly penetrant disorders. We demonstrate the utility of exome sequencing followed by imputation for discovering low-frequency variants associated with complex quantitative traits. We performed exome sequencing in a reference panel of 761 African Americans and then imputed newly discovered variants into a larger sample of more than 13,000 African Americans for association testing with the blood cell traits hemoglobin, hematocrit, white blood count, and platelet count. First, we illustrate the feasibility of our approach by demonstrating genome-wide-significant associations for variants that are not covered by conventional genotyping arrays; for example, one such association is that between higher platelet count and an MPL c.117G>T (p.Lys39Asn) variant encoding a p.Lys39Asn amino acid substitution of the thrombpoietin receptor gene (p = 1.5 × 10⁻¹¹). Second, we identified an association between missense variants of LCT and higher white blood count (p = 4 × 10⁻¹³). Third, we identified low-frequency coding variants that might account for allelic heterogeneity at several known blood cell-associated loci: MPL c.754T>C (p.Tyr252His) was associated with higher platelet count; CD36 c.975T>G (p.Tyr325^∗) was associated with lower platelet count; and several missense variants at the α-globin gene locus were associated with lower hemoglobin. By identifying low-frequency missense variants associated with blood cell traits not previously reported by genome-wide association studies, we establish that exome sequencing followed by imputation is a powerful approach to dissecting complex, genetically heterogeneous traits in large population-based studies.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

Pooled Sequencing of 531 Genes in Inflammatory Bowel Disease Identifies an Associated Rare Variant in BTNL2 and Implicates Other Immune Related Genes

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn’s disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10⁻¹⁰, OR = 2.3[95% CI = 1.75–3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1–5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.     

WW-domain-containing oxidoreductase is associated with low plasma HDL-C levels

Low serum HDL-cholesterol (HDL-C) is a major risk factor for coronary artery disease. We performed targeted genotyping of a 12.4 Mb linked region on 16q to test for association with low HDL-C by using a regional-tag SNP strategy. We identified one SNP, rs2548861, in the WW-domain-containing oxidoreductase (WWOX) gene with region-wide significance for low HDL-C in dyslipidemic families of Mexican and European descent and in low-HDL-C cases and controls of European descent (p = 6.9 × 10⁻⁷). We extended our investigation to the population level by using two independent unascertained population-based Finnish cohorts, the cross-sectional METSIM cohort of 4,463 males and the prospective Young Finns cohort of 2,265 subjects. The combined analysis provided p = 4 × 10⁻⁴ to 2 × 10⁻⁵. Importantly, in the prospective cohort, we observed a significant longitudinal association of rs2548861 with HDL-C levels obtained at four different time points over 21 years (p = 0.003), and the T risk allele explained 1.5% of the variance in HDL-C levels. The rs2548861 resides in a highly conserved region in intron 8 of WWOX. Results from our in vitro reporter assay and electrophoretic mobility-shift assay demonstrate that this region functions as a cis-regulatory element whose associated rs2548861 SNP has a specific allelic effect and that the region forms an allele-specific DNA-nuclear-factor complex. In conclusion, analyses of 9,798 subjects show significant association between HDL-C and a WWOX variant with an allele-specific cis-regulatory function.     

Meta-analysis of 65,734 Individuals Identifies TSPAN15 and SLC44A2 as Two Susceptibility Loci for Venous Thromboembolism

Venous thromboembolism (VTE), the third leading cause of cardiovascular mortality, is a complex thrombotic disorder with environmental and genetic determinants. Although several genetic variants have been found associated with VTE, they explain a minor proportion of VTE risk in cases. We undertook a meta-analysis of genome-wide association studies (GWASs) to identify additional VTE susceptibility genes. Twelve GWASs totaling 7,507 VTE case subjects and 52,632 control subjects formed our discovery stage where 6,751,884 SNPs were tested for association with VTE. Nine loci reached the genome-wide significance level of 5 × 10⁻⁸ including six already known to associate with VTE (ABO, F2, F5, F11, FGG, and PROCR) and three unsuspected loci. SNPs mapping to these latter were selected for replication in three independent case-control studies totaling 3,009 VTE-affected individuals and 2,586 control subjects. This strategy led to the identification and replication of two VTE-associated loci, TSPAN15 and SLC44A2, with lead risk alleles associated with odds ratio for disease of 1.31 (p = 1.67 × 10⁻¹⁶) and 1.21 (p = 2.75 × 10⁻¹⁵), respectively. The lead SNP at the TSPAN15 locus is the intronic rs78707713 and the lead SLC44A2 SNP is the non-synonymous rs2288904 previously shown to associate with transfusion-related acute lung injury. We further showed that these two variants did not associate with known hemostatic plasma markers. TSPAN15 and SLC44A2 do not belong to conventional pathways for thrombosis and have not been associated to other cardiovascular diseases nor related quantitative biomarkers. Our findings uncovered unexpected actors of VTE etiology and pave the way for novel mechanistic concepts of VTE pathophysiology.     

A Large-Scale Genetic Analysis Reveals a Strong Contribution of the HLA Class II Region to Giant Cell Arteritis Susceptibility

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10⁻⁴⁰, OR = 1.73). A multivariate model including class II amino acids of HLA-DRβ1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1^∗04. An omnibus test on polymorphic amino acid positions highlighted DRβ1 13 (p = 4.08 × 10⁻⁴³) and HLA-DQα1 47 (p = 4.02 × 10⁻⁴⁶), 56, and 76 (both p = 1.84 × 10⁻⁴⁵) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10⁻⁶, OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10⁻⁶, OR = 1.20), and REL (rs115674477, p = 1.10 × 10⁻⁵, OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.     

Identification of Allelic Heterogeneity at Type-2 Diabetes Loci and Impact on Prediction

Although over 60 single nucleotide polymorphisms (SNPs) have been identified by meta-analysis of genome-wide association studies for type-2 diabetes (T2D) among individuals of European descent, much of the genetic variation remains unexplained. There are likely many more SNPs that contribute to variation in T2D risk, some of which may lie in the regions surrounding established SNPs - a phenomenon often referred to as allelic heterogeneity. Here, we use the summary statistics from the DIAGRAM consortium meta-analysis of T2D genome-wide association studies along with linkage disequilibrium patterns inferred from a large reference sample to identify novel SNPs associated with T2D surrounding each of the previously established risk loci. We then examine the extent to which the use of these additional SNPs improves prediction of T2D risk in an independent validation dataset. Our results suggest that multiple SNPs at each of 3 loci contribute to T2D susceptibility (TCF7L2, CDKN2A/B, and KCNQ1; p<5×10⁻⁸). Using a less stringent threshold (p<5×10⁻⁴), we identify 34 additional loci with multiple associated SNPs. The addition of these SNPs slightly improves T2D prediction compared to the use of only the respective lead SNPs, when assessed using an independent validation cohort. Our findings suggest that some currently established T2D risk loci likely harbor multiple polymorphisms which contribute independently and collectively to T2D risk. This opens a promising avenue for improving prediction of T2D, and for a better understanding of the genetic architecture of T2D.     

A Genome-Wide Association Study Reveals Variants in ARL15 that Influence Adiponectin Levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10⁻⁸). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10⁻¹⁹ for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10⁻⁸, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10⁻⁶, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10⁻³, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.     

ANGPTL8 protein-truncating variant associated with lower serum triglycerides and risk of coronary disease

Protein-truncating variants (PTVs) affecting dyslipidemia risk may point to therapeutic targets for cardiometabolic disease. Our objective was to identify PTVs that were associated with both lipid levels and the risk of coronary artery disease (CAD) or type 2 diabetes (T2D) and assess their possible associations with risks of other diseases. To achieve this aim, we leveraged the enrichment of PTVs in the Finnish population and tested the association of low-frequency PTVs in 1,209 genes with serum lipid levels in the Finrisk Study (n = 23,435). We then tested which of the lipid-associated PTVs were also associated with the risks of T2D or CAD, as well as 2,683 disease endpoints curated in the FinnGen Study (n = 218,792). Two PTVs were associated with both lipid levels and the risk of CAD or T2D: triglyceride-lowering variants in ANGPTL8 (-24.0[-30.4 to -16.9] mg/dL per rs760351239-T allele, P = 3.4 × 10⁻⁹) and ANGPTL4 (-14.4[-18.6 to -9.8] mg/dL per rs746226153-G allele, P = 4.3 × 10⁻⁹). The risk of T2D was lower in carriers of the ANGPTL4 PTV (OR = 0.70[0.60–0.81], P = 2.2 × 10⁻⁶) than noncarriers. The odds of CAD were 47% lower in carriers of a PTV in ANGPTL8 (OR = 0.53[0.37–0.76], P = 4.5 × 10⁻⁴) than noncarriers. Finally, the phenome-wide scan of the ANGPTL8 PTV showed that the ANGPTL8 PTV carriers were less likely to use statin therapy (68,782 cases, OR = 0.52[0.40–0.68], P = 1.7 × 10⁻⁶) compared to noncarriers. Our findings provide genetic evidence of potential long-term efficacy and safety of therapeutic targeting of dyslipidemias.     

Identification of New Genetic Risk Variants for Type 2 Diabetes

Although more than 20 genetic susceptibility loci have been reported for type 2 diabetes (T2D), most reported variants have small to moderate effects and account for only a small proportion of the heritability of T2D, suggesting that the majority of inter-person genetic variation in this disease remains to be determined. We conducted a multistage, genome-wide association study (GWAS) within the Asian Consortium of Diabetes to search for T2D susceptibility markers. From 590,887 SNPs genotyped in 1,019 T2D cases and 1,710 controls selected from Chinese women in Shanghai, we selected the top 2,100 SNPs that were not in linkage disequilibrium (r²<0.2) with known T2D loci for in silico replication in three T2D GWAS conducted among European Americans, Koreans, and Singapore Chinese. The 5 most promising SNPs were genotyped in an independent set of 1,645 cases and 1,649 controls from Shanghai, and 4 of them were further genotyped in 1,487 cases and 3,316 controls from 2 additional Chinese studies. Consistent associations across all studies were found for rs1359790 (13q31.1), rs10906115 (10p13), and rs1436955 (15q22.2) with P-values (per allele OR, 95%CI) of 6.49×10⁻⁹ (1.15, 1.10–1.20), 1.45×10⁻⁸ (1.13, 1.08–1.18), and 7.14×10⁻⁷ (1.13, 1.08–1.19), respectively, in combined analyses of 9,794 cases and 14,615 controls. Our study provides strong evidence for a novel T2D susceptibility locus at 13q31.1 and the presence of new independent risk variants near regions (10p13 and 15q22.2) reported by previous GWAS.     

Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Meta-Analysis of Genome-Wide Association Studies in African Americans Provides Insights into the Genetic Architecture of Type 2 Diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.     

Novel association in chromosome 4q27 region with rheumatoid arthritis and confirmation of type 1 diabetes point to a general risk locus for autoimmune diseases

Recently, association of celiac disease with common single-nucleotide polymorphism (SNP) variants in an extensive linkage-disequilibrium block of 480 kb containing the KIAA1109, Tenr, IL2, and IL21 genes has been demonstrated in three independent populations (rs6822844 P_combined=1.3×10^-14). The KIAA1109/Tenr/IL2/IL21 block corresponds to the Idd3 locus in the nonobese diabetic mouse model of type 1 diabetes (T1D). This block was recently found to be associated with T1D in a genomewide association study, although this finding lacks unequivocal confirmation. We therefore aimed to investigate whether the KIAA1109/Tenr/IL2/IL21 region is involved in susceptibility to multiple autoimmune diseases. We tested SNP rs6822844 for association with disease in 350 T1D-affected and 1,047 rheumatoid arthritis (RA)–affected Dutch patients and in 929 controls. We replicated the association with T1D (P=.0006; OR 0.64 [95% CI 0.50–0.83]), and revealed a similar novel association with RA (P=.0002; OR 0.72 [95% CI 0.61–0.86]). Our results replicate and extend the association found in the KIAA1109/Tenr/IL2/IL21 gene region with autoimmune diseases, implying that this locus is a general risk factor for multiple autoimmune diseases.     

Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones—Role of Age and Comorbid Diseases

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P = 4.1×10⁻¹⁰). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P = 1.3×10⁻²³) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P = 3.0×10⁻¹⁷) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P = 1.0×10⁻⁶), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P = 5.7×10⁻⁵).     

Fine Mapping Major Histocompatibility Complex Associations in Psoriasis and Its Clinical Subtypes

Psoriasis vulgaris (PsV) risk is strongly associated with variation within the major histocompatibility complex (MHC) region, but its genetic architecture has yet to be fully elucidated. Here, we conducted a large-scale fine-mapping study of PsV risk in the MHC region in 9,247 PsV-affected individuals and 13,589 controls of European descent by imputing class I and II human leukocyte antigen (HLA) genes from SNP genotype data. In addition, we imputed sequence variants for MICA, an MHC HLA-like gene that has been associated with PsV, to evaluate association at that locus as well. We observed that HLA-C^∗06:02 demonstrated the lowest p value for overall PsV risk (p = 1.7 × 10⁻³⁶⁴). Stepwise analysis revealed multiple HLA-C^∗06:02-independent risk variants in both class I and class II HLA genes for PsV susceptibility (HLA-C^∗12:03, HLA-B amino acid positions 67 and 9, HLA-A amino acid position 95, and HLA-DQα1 amino acid position 53; p < 5.0 × 10⁻⁸), but no apparent risk conferred by MICA. We further evaluated risk of two major clinical subtypes of PsV, psoriatic arthritis (PsA; n = 3,038) and cutaneous psoriasis (PsC; n = 3,098). We found that risk heterogeneity between PsA and PsC might be driven by HLA-B amino acid position 45 (p_omnibus = 2.2 × 10⁻¹¹), indicating that different genetic factors underlie the overall risk of PsV and the risk of specific PsV subphenotypes. Our study illustrates the value of high-resolution HLA and MICA imputation for fine mapping causal variants in the MHC.