DNA methylation levels vary markedly by cell-type makeup of a sample. Understanding these differences and estimating the cell-type makeup of a sample is an important aspect of studying DNA methylation. DNA from leukocytes in whole blood is simple to obtain and pervasive in research. However, leukocytes contain many distinct cell types and subtypes. We propose a two-stage model that estimates the proportions of six main cell types in whole blood (CD4+ T cells, CD8+ T cells, monocytes, B cells, granulocytes, and natural killer cells) as well as subtypes of T and B cells. Unlike previous methods that only estimate overall proportions of CD4+ T cell, CD8+ T cells, and B cells, our model is able to estimate proportions of naïve, memory, and regulatory CD4+ T cells as well as naïve and memory CD8+ T cells and naïve and memory B cells. Using real and simulated data, we are able to demonstrate that our model is able to reliably estimate proportions of these cell types and subtypes. In studies with DNA methylation data from Illumina's HumanMethylation450k arrays, our estimates will be useful both for testing for associations of cell type and subtype composition with phenotypes of interest as well as for adjustment purposes to prevent confounding in epigenetic association studies. Additionally, our method can be easily adapted for use with whole genome bisulfite sequencing (WGBS) data or any other genome-wide methylation data platform.
DNA methylation; whole blood; T cell subtypes; B cell subtypes; epigenetics; deconvolution; cell-type composition
Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association.
Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR.
The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10−4, OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10−7, OR 0.71; case-only pmeta=1.9×10−4, OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR.
These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.
Systemic Lupus Erythematosus; Autoantibodies; Gene Polymorphism; B cells
Multiple studies have demonstrated that single-nucleotide polymorphisms (SNPs) in the ITGAM locus (including the non-synonymous SNPs rs1143679, rs1143678, rs1143683) are associated with SLE. ITGAM encodes the protein CD11b, a subunit of the β2 integrin Mac-1. The purpose of this study was to determine the effects of ITGAM genetic variation on the biological functions of neutrophil Mac-1.
Neutrophils from ITGAM genotyped and sequenced healthy donors were isolated for functional studies. The phagocytic capacity of neutrophil ITGAM variants was probed with complement coated erythrocytes, serum treated zymosan, heat treated zymosan and IgG coated erythrocytes. The adhesion capacity of ITGAM variants, in adhering to either purified intercellular adhesion molecule 1 or tumor necrosis factor α-stimulated endothelial cells was assessed in a flow chamber. Expression levels of total CD11b and activation of CD11b were assessed by flow cytometry.
Mac-1–mediated neutrophil phagocytosis, determined in cultures with 2 different complement-coated particles, was significantly reduced in individuals with nonsynonymous variant alleles of ITGAM. This reduction in phagocytosis was related to variation at either rs1143679 (in the β-propeller region) or rs1143678/rs1143683 (highly linked SNPs in the cytoplasmic/calf-1 regions). Phagocytosis mediated by Fcγ receptors was also significantly reduced in donors with variant ITGAM alleles. Similarly, firm adhesion of neutrophils was significantly reduced in individuals with variant ITGAM alleles. These functional alterations were not attributable to differences in total receptor expression or activation.
The nonsynonymous ITGAM variants rs1143679 and rs1143678/rs113683 contribute to altered Mac-1 function on neutrophils. These results underscore the need to consider multiple nonsynonymous SNPs when assessing the functional consequences of ITGAM variation on immune cell processes and the risk of SLE.
We previously reported that a higher degree of methylation of CpG sites in the promoter (positions 31, 37, 43, 52 and 58) and enhancer site 7862 of HPV 16 was associated with a lower likelihood of being diagnosed with HPV 16-associated CIN 2+. The purpose of the current study was to replicate our previous findings and in addition to evaluate the influence of plasma concentrations of folate and vitamin B12 on the degree of HPV 16 methylation (HPV 16m). The study included 315 HPV 16 positive women diagnosed with either CIN 2+ or ≤ CIN 1. Pyrosequencing technology was used to quantify the degree of HPV 16m. We reproduced the previously reported inverse association between HPV 16m and risk of being diagnosed with CIN 2+. In addition, we observed that women with higher plasma folate and higher HPV 16m or those with higher plasma vitamin B12 and higher HPV 16m were 75% (P<0.01) and 60% (P=0.02) less likely to be diagnosed with CIN 2+, respectively. With a tertile increase in the plasma folate or vitamin B12, there was a 50% (P=0.03) and 40% (P=0.07) increase in the odds of having a higher degree of HPV 16m, respectively. The present study provides initial evidence that methyl donor micronutrients, folate and vitamin B12, may play an important role in maintaining a desirably high degree of methylation at specific CpG sites in the HPV E6 promoter and enhancer that are associated with the likelihood of being diagnosed with CIN 2+.
HPV; methylation; folate; vitamin B12; CIN
Neutrophil firm adhesion to endothelial cells plays a critical role in inflammation in both health and disease. The process of neutrophil firm adhesion involves many different adhesion molecules including members of the β2 integrin family and their counter-receptors of the ICAM family. Recently, naturally occurring genetic variants in both β2 integrins and ICAMs are reported to be associated with autoimmune disease. Thus, the quantitative adhesive capacity of neutrophils from individuals with varying allelic forms of these adhesion molecules is important to study in relation to mechanisms underlying development of autoimmunity. Adhesion studies in flow chamber systems can create an environment with fluid shear stress similar to that observed in the blood vessel environment in vivo. Here, we present a method using a flow chamber assay system to study the quantitative adhesive properties of human peripheral blood neutrophils to human umbilical vein endothelial cell (HUVEC) and to purified ligand substrates. With this method, the neutrophil adhesive capacities from donors with different allelic variants in adhesion receptors can be assessed and compared. This method can also be modified to assess adhesion of other primary cell types or cell lines.
Immunology; Issue 89; neutrophil adhesion; flow chamber; human umbilical vein endothelial cell (HUVEC); purified ligand
Kawasaki disease (KD), response to intravenous immunoglobulin (IVIG) therapy, and associated coronary artery disease progression have been associated with genetic polymorphisms in Fc gamma receptor (FcγR) genes. However, it is not known whether the existing gene copy number (GCN) variability relates to KD treatment response, susceptibility, or associated sequelae.
The copy number of individuals with KD (n = 510) and their family members (n = 808) for three variable FcγRs was assessed using pyrosequencing. We performed the transmission disequilibrium test to examine the association of GCN for FcγRs (FcγR2C, FcγR3A, and FcγR3B) with susceptibility and used logistic regression models to determine its association with IVIG treatment outcomes.
FcγR2C and FcγR3B GCN were significantly associated with KD susceptibility. IVIG response was associated with GCN variations of FcγR3B in Whites and FcγR2C in Hispanics, and gene risk score based on single nucleotide polymorphism and GCN in FcγRs were significantly different between IVIG responders and nonresponders among Whites. We found no significant associations between coronary artery disease and any of the FcγR copy numbers.
GCN of FcγR2C and FcγR3B influences IVIG treatment response and predisposes individuals to KD, providing potential insights into understanding the mechanism of the FcγR gene family in the IVIG pathway.
copy number variation; FcγR; genetic risk; intravenous immunoglobulin; Kawasaki disease
TNF ligand superfamily member 13B (B-lymphocyte stimulator (BLyS), B cell activating factor (BAFF)) promotes primary B cell proliferation and immunoglobulin production. While the soluble form of BLyS/BAFF is thought to be the primary biologically active form, little is known about the regulation of its cleavage and processing. We provide evidence that Fcγ receptor cross-linking triggers a rapid release of soluble, biologically active BLyS/BAFF from myeloid cells. Surprisingly, this function is primarily mediated by FcγRI, but not FcγRIIa as defined by specific mAb, and can be initiated by both IgG and C reactive protein (CRP) as ligands. The generation of a B cell proliferation and survival factor by both innate and adaptive immune opsonins through engagement of an Fcγ receptor, which can also enhance antigen uptake and presentation, provides a unique opportunity to facilitate antibody production. These results provide a mechanism by which Fcγ receptors can elevate circulating BLyS levels and promote autoantibody production in immune complex mediated autoimmune diseases.
Fc Receptors; Monocytes/Macrophages; Human; Autoimmunity
To identify genetic determinants of granulomatosis with polyangiitis (Wegener’s) (GPA).
We carried out a genome-wide association study (GWAS) of 492 GPA cases and 1,506 healthy controls (white subjects of European descent), followed by replication analysis of the most strongly associated signals in an independent cohort of 528 GPA cases and 1,228 controls.
Genome-wide significant associations were identified in 32 single-nucleotide polymorphic (SNP) markers across the HLA region, the majority of which were located in the HLA–DPB1 and HLA–DPA1 genes encoding the class II major histocompatibility complex (MHC) DPβ chain 1 and DPα chain 1 proteins, respectively. Peak association signals in these 2 genes, emanating from SNPs rs9277554 (for DPβ chain 1) and rs9277341 (DPα chain 1) were strongly replicated in an independent cohort (in the combined analysis of the initial cohort and the replication cohort, P = 1.92 × 10−50 and 2.18 × 10−39, respectively). Imputation of classic HLA alleles and conditional analyses revealed that the SNP association signal was fully accounted for by the classic HLA–DPB1*04 allele. An independent single SNP, rs26595, near SEMA6A (the gene for semaphorin 6A) on chromosome 5, was also associated with GPA, reaching genome-wide significance in a combined analysis of the GWAS and replication cohorts (P = 2.09 × 10−8).
We identified the SEMA6A and HLA–DP loci as significant contributors to risk for GPA, with the HLA–DPB1*04 allele almost completely accounting for the MHC association. These two associations confirm the critical role of immunogenetic factors in the development of GPA.
To investigate whether the FcγRIIIa-66R/H/L polymorphism influences net effective receptor function and to assess if the FCGR3A combined genotypes formed by FcγRIIIa-66R/H/L and FcγRIIIa-176F/V as well as copy number variation (CNV) confer risk for development of SLE and lupus nephritis.
FcγRIIIa variants, expressed on A20 IIA1.6 cells, were used in flow cytometry-based human IgG binding assays. FCGR3A SNP and CNV genotypes were determined by Pyrosequencing methodology in a cohort of 1728 SLE patients and 2404 healthy controls.
The FcγRIIIa-66L/H/R (rs10127939) polymorphism influences ligand binding capacity in the context of the FcγRIIIa-176V (rs396991) allele. The low binding FcγRIIIa-176F allele was associated with SLE nephritis (p = 0.0609) in African Americans but not in European Americans (p > 0.10). Nephritis among African American SLE subjects was associated with FcγRIIIa low binding haplotypes containing the 66R/H/L and 176F variants (p = 0.03) and with low binding genotype combinations (p = 0.002). No association was observed in European American SLE patients. The distribution of FCGR3A CNV was not significantly different between controls and SLE patients with or without nephritis.
FcγRIIIa-66R/H/L influences ligand binding. The low binding haplotypes formed by 66R/H/L and 176F confer enhanced risk for lupus nephritis in African Americans. FCGR3A CNVs are not associated with SLE or SLE nephritis in either African Americans or European Americans.
Thrombosis is a serious complication of systemic lupus erythematosus (SLE). Studies that have
investigated the genetics of thrombosis in SLE are limited. We undertook this study to assess the
association of previously implicated candidate genes, particularly Toll-like receptor (TLR) genes,
with pathogenesis of thrombosis.
We genotyped 3,587 SLE patients from 3 multiethnic populations for 77 single-nucleotide
polymorphisms (SNPs) in 10 genes, primarily in TLRs 2, 4, 7, and 9, and we also genotyped 64
ancestry-informative markers (AIMs). We first analyzed association with arterial and venous
thrombosis in the combined population via logistic regression, adjusting for top principal
components of the AIMs and other covariates. We also subjected an associated SNP, rs893629, to
meta-analysis (after stratification by ethnicity and study population) to confirm the association
and to test for study population or ethnicity effects.
In the combined analysis, the SNP rs893629 in the KIAA0922/TLR2 region was
significantly associated with arterial thrombosis (logistic P = 6.4 ×
10−5, false discovery rate P = 0.0044). Two additional SNPs in
TLR2 were also suggestive: rs1816702 (logistic P = 0.002) and
rs4235232 (logistic P = 0.009). In the meta-analysis by study population, the odds
ratio (OR) for arterial thrombosis with rs893629 was 2.44 (95% confidence interval
1.58–3.76), without evidence for heterogeneity (P = 0.78). By ethnicity, the
effect was most significant among African Americans (OR 2.42, P = 3.5 ×
10−4) and European Americans (OR 3.47, P = 0.024).
TLR2 gene variation is associated with thrombosis in SLE, particularly among
African Americans and European Americans. There was no evidence of association among Hispanics, and
results in Asian Americans were limited due to insufficient sample size. These results may help
elucidate the pathogenesis of this important clinical manifestation.
Over 30 years ago, receptors for the Fc region of IgG (FcγR) were implicated in the pathogenesis of SLE. Since those pioneering studies, our knowledge of the structure and function of these FcγRs has increased dramatically. We now know that FcγR contribute to regulation of acquired immunity and to regulation of innate immune responses where FcγRs act as specific receptors for innate opsonins (CRP and SAP). Our understanding of the genomic architecture of the genes encoding the FcγR has also witnessed remarkable advances. Numerous functionally relevant SNP variants and copy number (CN) variants have been characterized in the FcγR genes. Many of these variants have also been shown to associate with risk to development of SLE and some have been associated with disease progression. This review will provide an overview of the FcγR in relation to SLE including consideration of the role of genetic variants in FcγR in SLE pathogenesis. The difficulties in assessing genetic variation in these genes will be discussed. To enhance our understanding of the functional roles of these receptors in SLE, future research will need to integrate our knowledge of SNP variants, CN variants and the functional diversity of these receptors.
B cells are pivotal regulators of acquired immune responses and recent work in both experimental murine models and humans has demonstrated that subtle changes in the regulation of B cell function can significantly alter immunological responses. The balance of negative and positive signals in maintaining an appropriate B cell activation threshold is critical in B lymphocyte immune tolerance and autoreactivity. FcγRIIb (CD32B), the only recognized Fcγ receptor on B cells, provides IgG-mediated negative modulation through a tyrosine-based inhibition motif which down-regulates B cell receptor initiated signaling. These properties make FcγRIIb a promising target for antibody-based therapy. Here we report the discovery of allele-dependent expression of the activating FcγRIIc on B cells. Identical to FcγRIIb in the extracellular domain, FcγRIIc has a tyrosine-based activation motif in its cytoplasmic domain. In both human B cells and in B cells from mice transgenic for human FcγRIIc, FcγRIIc expression counterbalances the negative feedback of FcγRIIb and enhances humoral responses to immunization in mice and to BioThrax® vaccination in a human Anthrax vaccine trial. Moreover, the FCGR2C-ORF allele is associated with the risk of development of autoimmunity in humans. FcγRIIc expression on B cells challenges the prevailing paradigm of uni-directional negative feedback by IgG immune complexes via the inhibitory FcγRIIb, is a previously unrecognized determinant in human antibody/autoantibody responses, and opens the opportunity for more precise personalized use of B cell targeted antibody-based therapy.
Lupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE) that exhibits familial aggregation and may progress to end-stage renal disease (ESRD). LN is more prevalent among African Americans than among European Americans. This study was undertaken to investigate the hypothesis that the apolipoprotein L1 gene (APOL1) nephropathy risk alleles G1/G2, common in African Americans and rare in European Americans, contribute to the ethnic disparity in risk.
APOL1 G1 and G2 nephropathy alleles were genotyped in 855 African American SLE patients with LN-ESRD (cases) and 534 African American SLE patients without nephropathy (controls) and tested for association under a recessive genetic model, by logistic regression.
Ninety percent of the SLE patients were female. The mean ± SD age at SLE diagnosis was significantly lower in LN-ESRD cases than in SLE non-nephropathy controls (27.3 ± 10.9 years versus 39.5 ± 12.2 years). The mean ± SD time from SLE diagnosis to development of LN-ESRD in cases was 7.3 ± 7.2 years. The G1/G2 risk alleles were strongly associated with SLE-ESRD, with 25% of cases and 12% of controls having 2 nephropathy alleles (odds ratio [OR] 2.57, recessive model P = 1.49 × 10−9), and after adjustment for age, sex, and ancestry admixture (OR 2.72, P = 6.23 × 10−6). The age-, sex-, and admixture-adjusted population attributable risk for ESRD among patients with G1/G2 polymorphisms was 0.26, compared to 0.003 among European American patients. The mean time from SLE diagnosis to ESRD development was ~2 years earlier among individuals with APOL1 risk genotypes (P = 0.01).
APOL1 G1/G2 alleles strongly impact the risk of LN-ESRD in African Americans, as well as the time to progression to ESRD. The high frequency of these alleles in African Americans with near absence in European Americans explains an important proportion of the increased risk of LN-ESRD in African Americans.
African Americans (AA) disproportionately develop lupus nephritis (LN) relative to European Americans and familial clustering supports causative genes. Since MYH9 underlies approximately 40% of end-stage renal disease (ESRD) in AA, we tested for genetic association with LN.
Seven MYH9 single nucleotide polymorphisms (SNPs) and the E1 risk haplotype were tested for association with LN in three cohorts of AA.
A preliminary analysis revealed that the MYH9 E1 risk haplotype was associated with ESRD in 25 cases with presumed systemic lupus erythematosus (SLE)-associated ESRD, compared to 735 non-SLE controls (odds ratio 3.1; p = 0.010 recessive). Replication analyses were performed in 583 AA with SLE in the PROFILE cohort (318 with LN; 265 with SLE but without nephropathy) and 60 AA from the NIH (39 with LN; 21 with SLE but without nephropathy). Analysis of the NIH and larger PROFILE cohorts, as well as a combined analysis, did not support this association.
These results suggest that AA with ESRD and coincident SLE who were recruited from dialysis clinics more likely have kidney diseases in the MYH9-associated spectrum of focal segmental glomerulosclerosis. PROFILE and NIH participants, recruited from rheumatology practices, demonstrate that MYH9 does not contribute substantially to the development of LN in AA.
African Americans; Genetics; Lupus nephritis; Kidney; MYH9; Systemic lupus erythematosus
The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. In a case-control analysis, we investigated whether 33 established and novel single nucleotide polymorphisms (SNP) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis in the general population were risk factors for SLE development among Asians.
Patients in the discovery cohort were enrolled in one of two North American SLE cohorts. Patients in the replication cohort were enrolled in one of four Asian or two North American cohorts. SLE cases met American College of Rheumatology classification criteria. We first genotyped 263 Asian SLE and 357 healthy Asian control individuals for 33 SNPs using Luminex multiplex technology in the discovery phase, and then used Taqman and Immunochip assays to examine 5 SNPs in up to an additional 1496 cases and 993 controls in the Replication phase. SLE patients were compared to healthy controls for association with minor alleles in allelic models. Principal components analysis was used to control for intra-Asian ancestry in an analysis of the replication cohort.
Two genetic variants in the gene VKORC1, rs9934438 and rs9923231, were highly significant in both the discovery and replication cohorts: OR(disc) = 2.45 (p=2×10−9), OR(rep) = 1.53 (p=5×10−6) and OR(disc) = 2.40 (p=6×10−9), OR(rep) = 1.53 (p=5×10−6), respectively. These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: rs9934438 OR(adj) = 1.34 (p=0.0029) and rs9923231 OR(adj) = 1.34 (p=0.0032).
Genetic variants in VKORC1, involved in vitamin K reduction and associated with DVT, are associated with SLE development in Asians. These results suggest intersecting genetic pathways for the development of SLE and thrombosis.
systemic lupus erythematosus; single nucleotide polymorphisms; genetic risk factors
Ligand specificity characterizes receptors for antibody and many other immune receptors, but the common use of the FcR-γ-chain as their signaling subunit challenges the concept that these receptors are functionally distinct. We hypothesized that elements for specificity might be determined by the unique cytoplasmic domain (CY) sequences of the ligand-binding α-chains of γ-chain associated receptors. Among Fcγ receptors (FcRs), a protein kinase C (PKC) phosphorylation consensus motif, [RSSTR], identified within the FcγRIIIa (CD16A) CY by in silico analysis, is specifically phosphorylated by PKCs, unlike other FcRs. Phosphorylated CD16A mediates a more robust calcium flux, tyrosine phosphorylation of Syk and pro-inflammatory cytokine production while non-phosphorylatable CD16A is more effective at activation of the Gab2/PI3K pathway, leading to enhanced degranulation. S100A4, a specific protein binding partner for CD16A-CY newly identified by yeast two-hybrid analysis, inhibits phosphorylation of CD16A-CY by PKC in vitro, and reduction of S100A4 levels in vivo enhances receptor phosphorylation upon cross-linking. Taken together, PKC-mediated phosphorylation of CD16A modulates distinct signaling pathways engaged by the receptor. Calcium activated binding of S100A4 to CD16A, promoted by the initial calcium flux, attenuates the phosphorylation of CY, and acting as a molecular switch, may both serve as a negative feedback on cytokine production pathways during sustained receptor engagement and favor a shift to degranulation, consistent with the importance of granule release following conjugate formation between CD16A+ effector cells and target cells. This switch mechanism points to new therapeutic targets and provides a frame for understanding novel receptor polymorphisms.
Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified more than 30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting the NF-κB signaling. In order to better understand the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway, we analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog, TNIP2, in case-control sets of diverse ethnic origins.
We fine-mapped TNIP1, TNIP2, and TAX1BP1 in a total of 8372 SLE cases and 7492 healthy controls from European-ancestry, African-American, Hispanic, East Asian, and African-American Gullah populations. Levels of TNIP1 mRNA and ABIN1 protein were analyzed using quantitative RT-PCR and Western blotting, respectively, in EBV-transformed human B cell lines.
We found significant associations between genetic variants within TNIP1 and SLE but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified two independent risk haplotypes within TNIP1 in individuals of European-ancestry that were also present in African-American and Hispanic populations. These risk haplotypes produced lower levels of TNIP1 mRNA and ABIN1 protein suggesting they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.
Our results confirmed the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.
Amerindian-Europeans, Asians and African-Americans have an excess morbidity from SLE and higher prevalence of lupus nephritis than Caucasians. The aim of this study was to analyze the relationship between genetic ancestry and socio-demographic characteristics and clinical features in a large cohort of Amerindian-European SLE patients.
A total of 2116 SLE patients of Amerindian-European origin and 4001 SLE patients of European descent with clinical data were used in the study. Genotyping of 253 continental ancestry informative markers was performed on the Illumina platform. The STRUCTURE and ADMIXTURE software were used to determine genetic ancestry of each individual. Correlation between ancestry and socio-demographic and clinical data were analyzed using logistic regression.
The average Amerindian genetic ancestry of 2116 SLE patients was 40.7%. There was an increased risk of having renal involvement (P<0.0001, OR= 3.50 95%CI 2.63-4.63) and an early age of onset with the presence of Amerindian genetic ancestry (P<0.0001). Amerindian ancestry protected against photosensitivity (P<0.0001, OR= 0.58 95%CI 0.44-0.76), oral ulcers (P<0.0001, OR= 0.55 95%CI 0.42-0.72), and serositis (P<0.0001, OR= 0.56 95%CI 0.41-0.75) after adjustment by age, gender and age of onset. However, gender and age of onset had stronger effects on malar rash, discoid rash, arthritis and neurological involvement than genetic ancestry.
In general, genetic Amerindian ancestry correlates with lower socio-demographic status and increases the risk for developing renal involvement and SLE at an earlier age of onset.
Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10−8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.
Systemic lupus erythematosus (SLE), a debilitating autoimmune disease characterized by the production of pathogenic autoantibodies, has a strong genetic basis. Variants of the IL10 gene, which encodes cytokine interleukin-10 (IL-10) with known function of promoting B cell hyperactivity and autoantibody production, are associated with SLE and other autoimmune diseases, and serum IL-10 levels are elevated in SLE patients correlating with increased disease activity. In this study, to discover SLE-predisposing causal variant(s), we assessed variants within the genomic region containing IL10 and its gene family member IL19, IL20 and IL24 for association with SLE in case and control subjects from diverse ancestries. We identified SLE-associated SNP rs3122605 located at 9.2 kb upstream of IL10 as the most likely causal variant in subjects of European ancestry. The SLE-risk allele of rs3122605 was dose-dependently associated with elevated IL10 expression at both mRNA and protein levels in peripheral blood samples from SLE patients and controls, which could be explained, at least in part, by its preferential binding to Elk-1, a transcription factor activated in B cells during active disease of SLE patients. Elk-1-mediated IL-10 overexpression could be downregulated by inhibiting activation of mitogen-activated protein kinases, suggesting a potential therapeutic target for SLE.
To examine the association of previously identified autoimmune disease susceptibility loci with granulomatosis with polyangiitis (GPA, formerly known as Wegener’s granulomatosis), and determine whether genetic susceptibility profiles of other autoimmune diseases are associated with GPA
Genetic data from two cohorts were meta-analyzed. Genotypes for 168 previously identified single nucleotide polymorphisms (SNPs) associated with susceptibility to different autoimmune diseases were ascertained for a total of 880 GPA cases and 1969 controls of European descent. Single marker associations were identified using additive logistic regression models. Multi-SNP associations with GPA were assessed using genetic risk scores based on susceptibility loci for Crohn’s disease, type 1 diabetes, systemic lupus erythematosus, rheumatoid arthritis, celiac disease, and ulcerative colitis. Adjustment for population substructure was performed in all analyses using ancestry informative markers and principal components analysis.
Genetic polymorphisms in CTLA4 were significantly associated with GPA in the single-marker meta-analysis (OR 0.79. 95% CI 0.70–0.89, p=9.8×10−5). A genetic risk score based on rheumatoid arthritis susceptibility markers was significantly associated with GPA (OR 1.05 per 1-unit increase in genetic risk score, 95% CI 1.02–1.08, p=5.1×10−5).
Rheumatoid arthritis and GPA may arise from a similar genetic predisposition. Aside from CTLA4, other loci previously found to be associated with common autoimmune diseases were not statistically associated with GPA in this study.
genetics; vasculitis; granulomatosis with polyangiitis; rheumatoid arthritis; CTLA4
Protein tyrosine phosphatase non-receptor type 22 (PTPN22) is a negative regulator of T-cell activation associated with several autoimmune diseases, including systemic lupus erythematosus (SLE). Missense rs2476601 is associated with SLE in individuals with European ancestry. Since the rs2476601 risk allele frequency differs dramatically across ethnicities, we assessed robustness of PTPN22 association with SLE and its clinical sub-phenotypes across four ethnically diverse populations. Ten SNPs were genotyped in 8220 SLE cases and 7369 controls from in European-Americans (EA), African-Americans (AA), Asians (AS), and Hispanics (HS). We performed imputation-based association followed by conditional analysis to identify independent associations. Significantly associated SNPs were tested for association with SLE clinical sub-phenotypes, including autoantibody profiles. Multiple testing was accounted for by using false discovery rate. We successfully imputed and tested allelic association for 107 SNPs within the PTPN22 region and detected evidence of ethnic-specific associations from EA and HS. In EA, the strongest association was at rs2476601 (P = 4.7×10−9, OR = 1.40 (95% CI = 1.25–1.56)). Independent association with rs1217414 was also observed in EA, and both SNPs are correlated with increased European ancestry. For HS imputed intronic SNP, rs3765598, predicted to be a cis-eQTL, was associated (P = 0.007, OR = 0.79 and 95% CI = 0.67–0.94). No significant associations were observed in AA or AS. Case-only analysis using lupus-related clinical criteria revealed differences between EA SLE patients positive for moderate to high titers of IgG anti-cardiolipin (aCL IgG >20) versus negative aCL IgG at rs2476601 (P = 0.012, OR = 1.65). Association was reinforced when these cases were compared to controls (P = 2.7×10−5, OR = 2.11). Our results validate that rs2476601 is the most significantly associated SNP in individuals with European ancestry. Additionally, rs1217414 and rs3765598 may be associated with SLE. Further studies are required to confirm the involvement of rs2476601 with aCL IgG.
Little is known about the genetic etiology of systemic lupus erythematosus (SLE) in individuals of African ancestry, despite its higher prevalence and greater disease severity. Overproduction of nitric oxide (NO) and reactive oxygen species are implicated in the pathogenesis and severity of SLE, making NO synthases and other reactive intermediate related genes biological candidates for disease susceptibility. This study analyzed variation in reactive intermediate genes for association with SLE in two populations with African ancestry.
A total of 244 SNPs from 53 regions were analyzed in non-Gullah African Americans (AA; 1432 cases and 1687 controls) and the genetically more homogeneous Gullah of the Sea Islands of South Carolina (133 cases and 112 controls) and. Single-marker, haplotype, and two-locus interaction tests were computed for these populations.
The glutathione reductase gene GSR (rs2253409, P=0.0014, OR [95% CI]=1.26 [1.09–1.44]) was the most significant single-SNP association in AA. In the Gullah, the NADH dehydrogenase NDUFS4 (rs381575, P=0.0065, OR [95%CI]=2.10 [1.23–3.59]) and nitric oxide synthase gene NOS1 (rs561712, P=0.0072, OR [95%CI]=0.62 [0.44–0.88]) were most strongly associated with SLE. When both populations were analyzed together, GSR remained the most significant effect (rs2253409, P=0.00072, OR [95%CI]=1.26 [1.10–1.44]). Haplotype and two-locus interaction analyses also uncovered different loci in each population.
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (λS = ~30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 × 10−7 < Poverall < 1.6 × 10−23; odds ratio 0.82–1.62)in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 × 10−5) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at ≥9 other loci (P < 2 × 10−7). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.
Several lines of evidence have supported a host genetic contribution to vaccine response, but genome-wide assessments for specific determinants have been sparse. Here we describe a genome-wide association study (GWAS) of protective antigen-specific antibody (AbPA) responses among 726 European-Americans who received Anthrax Vaccine Adsorbed (AVA) as part of a clinical trial. After quality control, 736,996 SNPs were tested for association with the AbPA response to 3 or 4 AVA vaccinations given over a 6-month period. No SNP achieved the threshold of genome-wide significance (p=5x10−8), but suggestive associations (p<1x10−5) were observed for SNPs in or near the class II region of the major histocompatibility complex (MHC), in the promoter region of SPSB1, and adjacent to MEX3C. Multivariable regression modeling suggested that much of the association signal within the MHC corresponded to previously identified HLA DR-DQ haplotypes involving component HLA-DRB1 alleles of *15:01, *01:01, or *01:02. We estimated the proportion of additive genetic variance explained by common SNP variation for the AbPA response after the 6 month vaccination. This analysis indicated a significant, albeit imprecisely estimated, contribution of variation tagged by common polymorphisms (p=0.032). Future studies will be required to replicate these findings in European Americans and to further elucidate the host genetic factors underlying variable immune response to AVA.
Anthrax vaccines; Bacillus anthracis; bacterial vaccines; vaccination; Genome-wide association study
A functional polymorphism in the inhibitory IgG-Fc receptor FcγRIIB influences intravenous immunoglobulin (IVIG) response in Kawasaki Disease (KD) a vasculitis preferentially affecting the coronary arteries in children. We tested the hypothesis that the polymorphisms in the activating receptors (Fcγ RIIA, Fcγ RIIIA and Fcγ RIIIB) also influence susceptibility, IVIG treatment response, and coronary artery disease (CAD) in KD patients.
Methods and Results
We genotyped polymorphisms in the activating FcγRIIA, FcγRIIIA and FcγRIIIB genes using pyrosequencing in 443 KD patients, including 266 trios and 150 single parent-child pairs, in northwest US and genetically determined race with 155 ancestry information markers. We used the FBAT program to test for transmission disequilibrium and further generated pseudo-sibling controls for comparisons to the cases. The FcγRIIA-131H variant showed an association with KD (p = 0.001) with ORadditive = 1.51 [1.16–1.96], p = 0.002) for the primary combined population, which persisted in both Caucasian (p = .04) and Asian (p = .01) subgroups and is consistent with the recent genome-wide association study. We also identified over-transmission of FcγRIIIB-NA1 among IVIG non-responders (p = 0.0002), and specifically to Caucasian IVIG non-responders (p = 0.007). Odds ratios for overall and Caucasian non-responders were respectively 3.67 [1.75–7.66], p = 0.0006 and 3.60 [1.34–9.70], p = 0.01. Excess NA1 transmission also occurred to KD with CAD (ORadditive = 2.13 [1.11–4.0], p = 0.02).
A common variation in FcγRIIA is associated with increased KD susceptibility. The FcγRIIIB-NA1, which confers higher affinity for IgG compared to NA2, is a determining factor for treatment response. These activating FcγRs play an important role in KD pathogenesis and mechanism of IVIG anti-inflammatory.
coronary disease; pediatrics; Kawasaki disease; IVIG treatment response; FcγR