Association of JAG1 with Bone Mineral Density and Osteoporotic Fractures: A Genome-wide Association Study and Follow-up Replication Studies

Bone mineral density (BMD), a diagnostic parameter for osteoporosis and a clinical predictor of fracture, is a polygenic trait with high heritability. To identify genetic variants that influence BMD in different ethnic groups, we performed a genome-wide association study (GWAS) on 800 unrelated Southern Chinese women with extreme BMD and carried out follow-up replication studies in six independent study populations of European descent and Asian populations including 18,098 subjects. In the meta-analysis, rs2273061 of the Jagged1 (JAG1) gene was associated with high BMD (p = 5.27 × 10⁻⁸ for lumbar spine [LS] and p = 4.15 × 10⁻⁵ for femoral neck [FN], n = 18,898). This SNP was further found to be associated with the low risk of osteoporotic fracture (p = 0.009, OR = 0.7, 95% CI 0.57–0.93, n = 1881). Region-wide and haplotype analysis showed that the strongest association evidence was from the linkage disequilibrium block 5, which included rs2273061 of the JAG1 gene (p = 8.52 × 10⁻⁹ for LS and 3.47 × 10⁻⁵ at FN). To assess the function of identified variants, an electrophoretic mobility shift assay demonstrated the binding of c-Myc to the “G” but not “A” allele of rs2273061. A mRNA expression study in both human bone-derived cells and peripheral blood mononuclear cells confirmed association of the high BMD-related allele G of rs2273061 with higher JAG1 expression. Our results identify the JAG1 gene as a candidate for BMD regulation in different ethnic groups, and it is a potential key factor for fracture pathogenesis.     

Sequence Variants in Three Loci Influence Monocyte Counts and Erythrocyte Volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p = 8.9 × 10⁻¹⁴). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p = 1.2 × 10⁻⁹, 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p = 5.3 × 10⁻⁹, 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values = 0.001, 9.9 × 10⁻⁵, and 7 × 10⁻⁵, respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.     

Genetic variants at 13q12.12 are associated with high myopia in the Han Chinese population

High myopia, which is extremely prevalent in the Chinese population, is one of the leading causes of blindness in the world. Genetic factors play a critical role in the development of the condition. To identify the genetic variants associated with high myopia in the Han Chinese, we conducted a genome-wide association study (GWAS) of 493,947 SNPs in 1088 individuals (419 cases and 669 controls) from a Han Chinese cohort and followed up on signals that were associated with p < 1.0 × 10⁻⁴ in three independent cohorts (combined, 2803 cases and 5642 controls). We identified a significant association between high myopia and a variant at 13q12.12 (rs9318086, combined p = 1.91 × 10⁻¹⁶, heterozygous odds ratio = 1.32, and homozygous odds ratio = 1.64). Furthermore, five additional SNPs (rs9510902, rs3794338, rs1886970, rs7325450, and rs7331047) in the same linkage disequilibrium (LD) block with rs9318086 also proved to be significantly associated with high myopia in the Han Chinese population; p values ranged from 5.46 × 10⁻¹¹ to 6.16 × 10⁻¹⁶. This associated locus contains three genes—MIPEP, C1QTNF9B-AS1, and C1QTNF9B. MIPEP and C1QTNF9B were found to be expressed in the retina and retinal pigment epithelium (RPE) and are more likely than C1QTNF9B-AS1 to be associated with high myopia given the evidence of retinal signaling that controls eye growth. Our results suggest that the variants at 13q12.12 are associated with high myopia.     

A Weighted False Discovery Rate Control Procedure Reveals Alleles at FOXA2 that Influence Fasting Glucose Levels

Association signals in GWAS are usually prioritized solely by p values. Here, we attempt to improve the power of GWAS by using a weighted false discovery rate control procedure to detect associations of low-frequency variants with effect sizes similar to or even larger than those of common variants. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to test for association with fasting glucose levels in the Atherosclerosis Risk in Communities Study (ARIC) population. In addition to finding several previously identified sequence variations, we identified a low-frequency variant (rs1209523; minor allele frequency = 0.043) near FOXA2 that was associated with fasting glucose levels in European Americans (EAs) (n = 7428, p value = 1.3 × 10⁻⁵). The association between rs1209523 and glucose levels was also significant in African Americans (AAs) (n = 2029, p value = 6.7 × 10⁻³) of the ARIC and was confirmed by replication in both EAs and AAs of the Dallas Heart Study (n = 963 and 1571, respectively; p values = 5.3 × 10⁻³ and 5.8 × 10⁻⁴, respectively) and in EAs of the Cooper Center Longitudinal Study (n = 2862; p value = 1.6 × 10⁻²). A meta-analysis of these five populations yielded an estimated effect size of −1.31 mg/dl per minor allele (p value = 2.2 × 10⁻¹¹). This study reveals that there is a cache of less-frequent variants in GWAS arrays that can be identified via analytical approaches accounting for allele frequencies.     

A Genome-wide Association Study of Lung Cancer Identifies a Region of Chromosome 5p15 Associated with Risk for Adenocarcinoma

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR] = 1.23, 95% confidence interval [CI] = 1.13–1.33, p = 3.02 × 10⁻⁷), but not with other histologic types (OR = 1.01, p = 0.84 and OR = 1.00, p = 0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR = 1.24, 95% CI = 1.17–1.31, p = 3.74 × 10⁻¹⁴ for AD; OR = 0.99, p = 0.69 and OR = 0.97, p = 0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.     

The IL-33-ST2L Pathway Is Associated with Coronary Artery Disease in a Chinese Han Population

The effects of interleukin-33 (IL-33) on the immune system have been clearly demonstrated; however, in cardiovascular diseases, especially in coronary artery disease (CAD), these effects have not yet been clarified. In this study, we investigate the genetic role of the IL-33-ST2L pathway in CAD. We performed three-stage case-control association analyses on a total of 4,521 individuals with CAD and 4,809 controls via tag SNPs in the genes encoding IL-33 and ST2L-IL-1RL1. One tag SNP in each gene was significantly associated with CAD (rs7025417^T in IL33, p_adj = 1.19 × 10⁻²⁸, OR = 1.39, 95% CI: 1.31–1.47; rs11685424^G in IL1RL1, p_adj = 6.93 × 10⁻³⁰, OR = 1.40, 95% CI: 1.32–1.48). Combining significant variants in two genes, the risk for CAD increased nearly 5-fold (p_adj = 8.90 × 10⁻²¹, OR = 4.98, 95% CI: 3.56–6.97). Traditional risk factors for CAD were adjusted for the association studies by SPSS with logistic regression analysis. With the two variants above, both located within the gene promoter regions, reporter gene analysis indicated that the rs7025417 C>T and rs11685424 A>G changes resulted in altered regulation of IL33 and IL1RL1 gene expression, respectively (p < 0.005). Further studies revealed that the rs7025417 genotype was significantly associated with plasma IL-33 levels in the detectable subjects (n = 227, R² = 0.276, p = 1.77 × 10⁻¹⁷): the level of IL-33 protein increased with the number of rs7025417 risk (T) alleles. Based on genetic evidence in humans, the IL-33-ST2L pathway appears to have a causal role in the development of CAD, highlighting this pathway as a valuable target for the prevention and treatment of CAD.     

Meta-analysis of Dense Genecentric Association Studies Reveals Common and Uncommon Variants Associated with Height

Height is a classic complex trait with common variants in a growing list of genes known to contribute to the phenotype. Using a genecentric genotyping array targeted toward cardiovascular-related loci, comprising 49,320 SNPs across approximately 2000 loci, we evaluated the association of common and uncommon SNPs with adult height in 114,223 individuals from 47 studies and six ethnicities. A total of 64 loci contained a SNP associated with height at array-wide significance (p < 2.4 × 10⁻⁶), with 42 loci surpassing the conventional genome-wide significance threshold (p < 5 × 10⁻⁸). Common variants with minor allele frequencies greater than 5% were observed to be associated with height in 37 previously reported loci. In individuals of European ancestry, uncommon SNPs in IL11 and SMAD3, which would not be genotyped with the use of standard genome-wide genotyping arrays, were strongly associated with height (p < 3 × 10⁻¹¹). Conditional analysis within associated regions revealed five additional variants associated with height independent of lead SNPs within the locus, suggesting allelic heterogeneity. Although underpowered to replicate findings from individuals of European ancestry, the direction of effect of associated variants was largely consistent in African American, South Asian, and Hispanic populations. Overall, we show that dense coverage of genes for uncommon SNPs, coupled with large-scale meta-analysis, can successfully identify additional variants associated with a common complex trait.     

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (n_rs495828 = 53 and n_rs505922 = 59); strongest association with venous embolism, odds ratio (OR_rs495828 1.33 (p = 1.32 x 10⁻¹⁹⁹), and thrombosis OR_rs505922 1.33, p = 2.2 x10^-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10⁻¹⁹¹; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10⁻¹². The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10^-23, lupus OR 0.84, p = 3.97 x 10⁻⁰⁶. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10⁻¹³ among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.     

Integrating autoimmune risk loci with gene-expression data identifies specific pathogenic immune cell subsets

Although genome-wide association studies have implicated many individual loci in complex diseases, identifying the exact causal alleles and the cell types within which they act remains greatly challenging. To ultimately understand disease mechanism, researchers must carefully conceive functional studies in relevant pathogenic cell types to demonstrate the cellular impact of disease-associated genetic variants. This challenge is highlighted in autoimmune diseases, such as rheumatoid arthritis, where any of a broad range of immunological cell types might potentially be impacted by genetic variation to cause disease. To this end, we developed a statistical approach to identify potentially pathogenic cell types in autoimmune diseases by using a gene-expression data set of 223 murine-sorted immune cells from the Immunological Genome Consortium. We found enrichment of transitional B cell genes in systemic lupus erythematosus (p = 5.9 × 10⁻⁶) and epithelial-associated stimulated dendritic cell genes in Crohn disease (p = 1.6 × 10⁻⁵). Finally, we demonstrated enrichment of CD4+ effector memory T cell genes within rheumatoid arthritis loci (p < 10⁻⁶). To further validate the role of CD4+ effector memory T cells within rheumatoid arthritis, we identified 436 loci that were not yet known to be associated with the disease but that had a statistically suggestive association in a recent genome-wide association study (GWAS) meta-analysis (p_GWAS < 0.001). Even among these putative loci, we noted a significant enrichment for genes specifically expressed in CD4+ effector memory T cells (p = 1.25 × 10⁻⁴). These cell types are primary candidates for future functional studies to reveal the role of risk alleles in autoimmunity. Our approach has application in other phenotypes, outside of autoimmunity, where many loci have been discovered and high-quality cell-type-specific gene expression is available.     

Abdominal aortic aneurysm is associated with a variant in low-density lipoprotein receptor-related protein 1

Abdominal aortic aneurysm (AAA) is a common cause of morbidity and mortality and has a significant heritability. We carried out a genome-wide association discovery study of 1866 patients with AAA and 5435 controls and replication of promising signals (lead SNP with a p value < 1 × 10⁻⁵) in 2871 additional cases and 32,687 controls and performed further follow-up in 1491 AAA and 11,060 controls. In the discovery study, nine loci demonstrated association with AAA (p < 1 × 10⁻⁵). In the replication sample, the lead SNP at one of these loci, rs1466535, located within intron 1 of low-density-lipoprotein receptor-related protein 1 (LRP1) demonstrated significant association (p = 0.0042). We confirmed the association of rs1466535 and AAA in our follow-up study (p = 0.035). In a combined analysis (6228 AAA and 49182 controls), rs1466535 had a consistent effect size and direction in all sample sets (combined p = 4.52 × 10⁻¹⁰, odds ratio 1.15 [1.10–1.21]). No associations were seen for either rs1466535 or the 12q13.3 locus in independent association studies of coronary artery disease, blood pressure, diabetes, or hyperlipidaemia, suggesting that this locus is specific to AAA. Gene-expression studies demonstrated a trend toward increased LRP1 expression for the rs1466535 CC genotype in arterial tissues; there was a significant (p = 0.029) 1.19-fold (1.04–1.36) increase in LRP1 expression in CC homozygotes compared to TT homozygotes in aortic adventitia. Functional studies demonstrated that rs1466535 might alter a SREBP-1 binding site and influence enhancer activity at the locus. In conclusion, this study has identified a biologically plausible genetic variant associated specifically with AAA, and we suggest that this variant has a possible functional role in LRP1 expression.     

Identification of KCNJ15 as a Susceptibility Gene in Asian Patients with Type 2 Diabetes Mellitus

Recent advances in genome research have enabled the identification of new genomic variations that are associated with type 2 diabetes mellitus (T2DM). Via fine mapping of SNPs in a candidate region of chromosome 21q, the current study identifies potassium inwardly-rectifying channel, subfamily J, member 15 (KCNJ15) as a new T2DM susceptibility gene. KCNJ15 is expressed in the β cell of the pancreas, and a synonymous SNP, rs3746876, in exon 4 (C566T) of this gene, with T allele frequency among control subjects of 3.1%, showed a significant association with T2DM affecting lean individuals in three independent Japanese sample sets (p = 2.5 × 10⁻⁷, odds ratio [OR] = 2.54, 95% confidence interval [CI] = 1.76–3.67) and with unstratified T2DM (p = 6.7 × 10⁻⁶, OR = 1.76, 95% CI = 1.37–2.25). The diabetes risk allele frequency was, however, very low among Europeans in whom no association between this variant and T2DM could be shown. Functional analysis in human embryonic kidney 293 cells demonstrated that the risk allele of the synonymous SNP in exon 4 increased KCNJ15 expression via increased mRNA stability, which resulted in the higher expression of protein as compared to that of the nonrisk allele. We also showed that KCNJ15 is expressed in human pancreatic β cells. In conclusion, we demonstrated a significant association between a synonymous variant in KCNJ15 and T2DM in lean Japanese patients with T2DM, suggesting that KCNJ15 is a previously unreported susceptibility gene for T2DM among Asians.     

Genome-wide screen identifies rs646776 near sortilin as a regulator of progranulin levels in human plasma

Recent studies suggest progranulin (GRN) is a neurotrophic factor. Loss-of-function mutations in the progranulin gene (GRN) cause frontotemporal lobar degeneration (FTLD), a progressive neurodegenerative disease affecting ∼10% of early-onset dementia patients. Using an enzyme-linked immunosorbent assay, we previously showed that GRN is detectable in human plasma and can be used to predict GRN mutation status. This study also showed a wide range in plasma GRN levels in non-GRN mutation carriers, including controls. We have now performed a genome-wide association study of 313,504 single-nucleotide polymorphisms (SNPs) in 533 control samples and identified on chromosome 1p13.3 two SNPs with genome-wide significant association with plasma GRN levels (top SNP rs646776; p = 1.7 × 10⁻³⁰). The association of rs646776 with plasma GRN levels was replicated in two independent series of 508 controls (p = 1.9 × 10⁻¹⁹) and 197 FTLD patients (p = 6.4 × 10⁻¹²). Overall, each copy of the minor C allele decreased GRN levels by ∼15%. SNP rs646776 is located near sortilin (SORT1), and the minor C allele of rs646776 was previously associated with increased SORT1 mRNA levels. Supporting these findings, overexpression of SORT1 in cultured HeLa cells dramatically reduced GRN levels in the conditioned media, whereas knockdown of SORT1 increased extracellular GRN levels. In summary, we identified significant association of a locus on chromosome 1p13.3 with plasma GRN levels through an unbiased genome-wide screening approach and implicated SORT1 as an important regulator of GRN levels. This finding opens avenues for future research into GRN biology and the pathophysiology of neurodegenerative diseases.     

Quantitative Trait Loci for CD4:CD8 Lymphocyte Ratio Are Associated with Risk of Type 1 Diabetes and HIV-1 Immune Control

Abnormal expansion or depletion of particular lymphocyte subsets is associated with clinical manifestations such as HIV progression to AIDS and autoimmune disease. We sought to identify genetic predictors of lymphocyte levels and reasoned that these may play a role in immune-related diseases. We tested 2.3 million variants for association with five lymphocyte subsets, measured in 2538 individuals from the general population, including CD4+ T cells, CD8+ T cells, CD56+ natural killer (NK) cells, and the derived measure CD4:CD8 ratio. We identified two regions of strong association. The first was located in the major histocompatibility complex (MHC), with multiple SNPs strongly associated with CD4:CD8 ratio (rs2524054, p = 2.1 × 10⁻²⁸). The second region was centered within a cluster of genes from the Schlafen family and was associated with NK cell levels (rs1838149, p = 6.1 × 10⁻¹⁴). The MHC association with CD4:CD8 replicated convincingly (p = 1.4 × 10⁻⁹) in an independent panel of 988 individuals. Conditional analyses indicate that there are two major independent quantitative trait loci (QTL) in the MHC region that regulate CD4:CD8 ratio: one is located in the class I cluster and influences CD8 levels, whereas the second is located in the class II cluster and regulates CD4 levels. Jointly, both QTL explained 8% of the variance in CD4:CD8 ratio. The class I variants are also strongly associated with durable host control of HIV, and class II variants are associated with type-1 diabetes, suggesting that genetic variation at the MHC may predispose one to immune-related diseases partly through disregulation of T cell homeostasis.     

Genome-wide Association Study Reveals Multiple Nasopharyngeal Carcinoma-Associated Loci within the HLA Region at Chromosome 6p21.3

Nasopharyngeal carcinoma (NPC) is a multifactorial malignancy closely associated with genetic factors and Epstein-Barr virus infection. To identify the common genetic variants linked to NPC susceptibility, we conducted a genome-wide association study (GWAS) in 277 NPC patients and 285 healthy controls within the Taiwanese population, analyzing 480,365 single-nucleotide polymorphisms (SNPs). Twelve statistically significant SNPs were identified and mapped to chromosome 6p21.3. Associations were replicated in two independent sets of case-control samples. Two of the most significant SNPs (rs2517713 and rs2975042; p_combined = 3.9 × 10⁻²⁰ and 1.6 × 10⁻¹⁹, respectively) were located in the HLA-A gene. Moreover, we detected significant associations between NPC and two genes: specifically, gamma aminobutyric acid b receptor 1 (GABBR1) (rs29232; p_combined = 8.97 × 10⁻¹⁷) and HLA-F (rs3129055 and rs9258122; p_combined = 7.36 × 10⁻¹¹ and 3.33 × 10⁻¹⁰, respectively). Notably, the association of rs29232 remained significant (residual p < 5 × 10⁻⁴) after adjustment for age, gender, and HLA-related SNPs. Furthermore, higher GABA_B receptor 1 expression levels can be found in the tumor cells in comparison to the adjacent epithelial cells (p < 0.001) in NPC biopsies, implying a biological role of GABBR1 in NPC carcinogenesis. To our knowledge, it is the first GWAS report of NPC showing that multiple loci (HLA-A, HLA-F, and GABBR1) within chromosome 6p21.3 are associated with NPC. Although some of these relationships may be attributed to linkage disequilibrium between the loci, the findings clearly provide a fresh direction for the study of NPC development.     

PRKCA: A Positional Candidate Gene for Body Mass Index and Asthma

Asthma incidence and prevalence are higher in obese individuals. A potential mechanistic basis for this relationship is pleiotropy. We hypothesized that significant linkage and candidate-gene association would be found for body mass index (BMI) in a population ascertained on asthma affection status. Linkage analysis for BMI was performed on 657 subjects in eight Costa Rican families enrolled in a study of asthma. Family-based association studies were conducted for BMI with SNPs within a positional candidate gene, PRKCA. SNPs within PRKCA were also tested for association with asthma. Association studies were conducted in 415 Costa Rican parent-child trios and 493 trios participating in the Childhood Asthma Management Program (CAMP). Although only modest evidence of linkage for BMI was obtained for the whole cohort, significant linkage was noted for BMI in females on chromosome 17q (peak LOD = 3.39). Four SNPs in a candidate gene in this region (PRKCA) had unadjusted association p values < 0.05 for BMI in both cohorts, with the joint p value for two SNPs remaining significant after adjustment for multiple comparisons (rs228883 and rs1005651, joint p values = 9.5 × 10⁻⁵ and 5.6 × 10⁻⁵). Similarly, eight SNPs had unadjusted association p values < 0.05 for asthma in both populations, with one SNP remaining significant after adjustment for multiple comparisons (rs11079657, joint p value = 2.6 × 10⁻⁵). PRKCA is a pleiotropic locus that is associated with both BMI and asthma and that has been identified via linkage analysis of BMI in a population ascertained on asthma.     

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p_meta-Euro = 2.08 × 10⁻¹⁰), transmembrane protein 39A (TMEM39A; rs1132200; p_meta-all = 8.62 × 10⁻⁹), and 17q21 (rs1453560; p_meta-all = 3.48 × 10⁻¹⁰) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10⁻⁸ < p_meta-Euro < 9.99 × 10⁻⁵) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.     

Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10⁻⁵ were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10⁻⁸, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10⁻⁷ and rs9275245, P = 1.39×10⁻⁷. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10⁻⁵, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10⁻⁵) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10⁻¹⁰, OR:1.25), TNIP1 (P = 4.68×10⁻⁹, OR:1.31), and RHOB loci (P = 3.17×10⁻⁶, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.     

A Locus on Chromosome 1p36 Is Associated with Thyrotropin and Thyroid Function as Identified by Genome-wide Association Study

Thyroid hormones are key regulators of cellular growth, development, and metabolism, and thyroid disorders are a common cause of ill health in the community. Circulating concentrations of thyrotropin (TSH), thyroxine (T4) and triiodothyronine (T3) have a strong heritable component and are thought to be under polygenic control, but the genes responsible are mostly unknown. In order to identify genetic loci associated with these metabolic phenotypes, we performed a genome-wide association study of 2,120,505 SNPs in 2014 female twins from the TwinsUK study and found a significant association between rs10917469 on chromosome 1p36.13 and serum TSH (p = 3.2 × 10⁻⁸). The association of rs10917469 with serum TSH was replicated (p = 2.0 × 10⁻⁴) in an independent community-based sample of 1154 participants in the Busselton Health Study. This SNP is located near CAPZB, which might be a regulator of TSH secretion and thus of pituitary-thyroid axis function. Twenty-nine percent of white individuals carry the variant, and the difference in mean TSH concentrations between wild-type individuals and those homozygous for the minor G allele was 0.5 mU/l, which is likely to be clinically relevant. We also provide evidence of suggestive association (p < 5.0 × 10⁻⁶) of other SNPs with serum TSH, free T4, and free T3 concentrations, and these SNPs might be good targets for further studies. These results advance understanding of the genetic basis of pituitary-thyroid axis function and metabolic regulation.     

Combined analysis of genome-wide association studies for Crohn disease and psoriasis identifies seven shared susceptibility loci

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 × 10⁻⁸) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p_rs1250544 = 3.53 × 10⁻⁸, 11q13 near PRDX5, p_rs694739 = 3.71 × 10⁻⁰⁹, 22q11 at YDJC, p_rs181359 = 8.02 × 10⁻¹⁰). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p_rs4780355 = 4.99 × 10⁻⁸). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.     

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⁻⁸).