Search tips
Search criteria

Results 1-25 (76)

Clipboard (0)

Select a Filter Below

more »
Year of Publication
more »
1.  Estimation of Cell-Type Composition Including T and B Cell Subtypes for Whole Blood Methylation Microarray Data 
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.
PMCID: PMC4757643  PMID: 26925097
DNA methylation; whole blood; T cell subtypes; B cell subtypes; epigenetics; deconvolution; cell-type composition
2.  Preferential association of a functional variant in complement receptor 2 with antibodies to double-stranded DNA 
Annals of the Rheumatic Diseases  2014;75(1):242-252.
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.
PMCID: PMC4717392  PMID: 25180293
Systemic Lupus Erythematosus; Autoantibodies; Gene Polymorphism; B cells
3.  The IRF5–TNPO3 association with systemic lupus erythematosus has two components that other autoimmune disorders variably share 
Kottyan, Leah C. | Zoller, Erin E. | Bene, Jessica | Lu, Xiaoming | Kelly, Jennifer A. | Rupert, Andrew M. | Lessard, Christopher J. | Vaughn, Samuel E. | Marion, Miranda | Weirauch, Matthew T. | Namjou, Bahram | Adler, Adam | Rasmussen, Astrid | Glenn, Stuart | Montgomery, Courtney G. | Hirschfield, Gideon M. | Xie, Gang | Coltescu, Catalina | Amos, Chris | Li, He | Ice, John A. | Nath, Swapan K. | Mariette, Xavier | Bowman, Simon | Rischmueller, Maureen | Lester, Sue | Brun, Johan G. | Gøransson, Lasse G. | Harboe, Erna | Omdal, Roald | Cunninghame-Graham, Deborah S. | Vyse, Tim | Miceli-Richard, Corinne | Brennan, Michael T. | Lessard, James A. | Wahren-Herlenius, Marie | Kvarnström, Marika | Illei, Gabor G. | Witte, Torsten | Jonsson, Roland | Eriksson, Per | Nordmark, Gunnel | Ng, Wan-Fai | Anaya, Juan-Manuel | Rhodus, Nelson L. | Segal, Barbara M. | Merrill, Joan T. | James, Judith A. | Guthridge, Joel M. | Hal Scofield, R. | Alarcon-Riquelme, Marta | Bae, Sang-Cheol | Boackle, Susan A. | Criswell, Lindsey A. | Gilkeson, Gary | Kamen, Diane L. | Jacob, Chaim O. | Kimberly, Robert | Brown, Elizabeth | Edberg, Jeffrey | Alarcón, Graciela S. | Reveille, John D. | Vilá, Luis M. | Petri, Michelle | Ramsey-Goldman, Rosalind | Freedman, Barry I. | Niewold, Timothy | Stevens, Anne M. | Tsao, Betty P. | Ying, Jun | Mayes, Maureen D. | Gorlova, Olga Y. | Wakeland, Ward | Radstake, Timothy | Martin, Ezequiel | Martin, Javier | Siminovitch, Katherine | Moser Sivils, Kathy L. | Gaffney, Patrick M. | Langefeld, Carl D. | Harley, John B. | Kaufman, Kenneth M.
Human Molecular Genetics  2014;24(2):582-596.
Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5–TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5–TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10−49; OR = 1.38–1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10−27–10−32, OR = 1.7–1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5–TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5–TNPO3.
PMCID: PMC4275071  PMID: 25205108
4.  Folate and vitamin B12 may play a critical role in lowering the HPV 16 methylation associated risk of developing higher grades of CIN 
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+.
PMCID: PMC4711818  PMID: 25145486
HPV; methylation; folate; vitamin B12; CIN
5.  Association of Granulomatosis With Polyangiitis (Wegener’s) With HLA–DPB1*04 and SEMA6A Gene Variants Evidence From Genome-Wide Analysis 
Arthritis and rheumatism  2013;65(9):2457-2468.
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.
PMCID: PMC4471994  PMID: 23740775
6.  FcγRIIIa SNPs and haplotypes affect human IgG binding and association with lupus nephritis in African Americans 
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.
PMCID: PMC4069204  PMID: 24782186
7.  FcγR gene copy number in Kawasaki disease and intravenous immunoglobulin treatment response 
Pharmacogenetics and genomics  2013;23(9):455-462.
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.
PMCID: PMC4400828  PMID: 23778324
copy number variation; FcγR; genetic risk; intravenous immunoglobulin; Kawasaki disease
8.  A Polymorphism in TLR2 Is Associated With Arterial Thrombosis in a Multiethnic Population of Patients With Systemic Lupus Erythematosus 
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.
PMCID: PMC4269184  PMID: 24578102
9.  Multiple Lupus Associated ITGAM Variants Alter Mac-1 Function on Neutrophils 
Arthritis and rheumatism  2013;65(11):2907-2916.
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.
PMCID: PMC3969028  PMID: 23918739
10.  Human Neutrophil Flow Chamber Adhesion Assay 
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.
PMCID: PMC4166108  PMID: 25045887
Immunology; Issue 89; neutrophil adhesion; flow chamber; human umbilical vein endothelial cell (HUVEC); purified ligand
11.  Allelic Dependent Expression of an Activating Fc receptor on B cells Enhances Humoral Immune Responses 
Science translational medicine  2013;5(216):216ra175.
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.
PMCID: PMC3982386  PMID: 24353158
12.  End-Stage Renal Disease in African Americans With Lupus Nephritis Is Associated With APOL1 
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.
PMCID: PMC4002759  PMID: 24504811
13.  SNPs in VKORC1 are Risk Factors for Systemic Lupus Erythematosus in Asians 
Arthritis and rheumatism  2013;65(1):211-215.
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.
PMCID: PMC3670944  PMID: 23124848
systemic lupus erythematosus; single nucleotide polymorphisms; genetic risk factors
14.  The unique cytoplasmic domain of human FcγRIIIA regulates receptor mediated function 
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.
PMCID: PMC3478424  PMID: 23024279
15.  Two Independent Functional Risk Haplotypes in TNIP1 are Associated with Systemic Lupus Erythematosus 
Arthritis and rheumatism  2012;64(11):3695-3705.
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.
PMCID: PMC3485412  PMID: 22833143
16.  Impact of Genetic Ancestry and Socio-Demographic Status on the Clinical Expression of Systemic Lupus Erythematosus in Amerindian-European Populations 
Arthritis and rheumatism  2012;64(11):3687-3694.
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.
PMCID: PMC3485439  PMID: 22886787
17.  Preferential Binding to Elk-1 by SLE-Associated IL10 Risk Allele Upregulates IL10 Expression 
PLoS Genetics  2013;9(10):e1003870.
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.
Author Summary
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.
PMCID: PMC3794920  PMID: 24130510
18.  Meta-analysis in granulomatosis with polyangiitis reveals shared susceptibility loci with rheumatoid arthritis 
Arthritis and rheumatism  2012;64(10):3463-3471.
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.
PMCID: PMC3425721  PMID: 22508400
genetics; vasculitis; granulomatosis with polyangiitis; rheumatoid arthritis; CTLA4
19.  Genome-Wide DNA Methylation Analysis of Systemic Lupus Erythematosus Reveals Persistent Hypomethylation of Interferon Genes and Compositional Changes to CD4+ T-cell Populations 
PLoS Genetics  2013;9(8):e1003678.
Systemic lupus erythematosus (SLE) is an autoimmune disease with known genetic, epigenetic, and environmental risk factors. To assess the role of DNA methylation in SLE, we collected CD4+ T-cells, CD19+ B-cells, and CD14+ monocytes from 49 SLE patients and 58 controls, and performed genome-wide DNA methylation analysis with Illumina Methylation450 microarrays. We identified 166 CpGs in B-cells, 97 CpGs in monocytes, and 1,033 CpGs in T-cells with highly significant changes in DNA methylation levels (p<1×10−8) among SLE patients. Common to all three cell-types were widespread and severe hypomethylation events near genes involved in interferon signaling (type I). These interferon-related changes were apparent in patients collected during active and quiescent stages of the disease, suggesting that epigenetically-mediated hypersensitivity to interferon persists beyond acute stages of the disease and is independent of circulating interferon levels. This interferon hypersensitivity was apparent in memory, naïve and regulatory T-cells, suggesting that this epigenetic state in lupus patients is established in progenitor cell populations. We also identified a widespread, but lower amplitude shift in methylation in CD4+ T-cells (>16,000 CpGs at FDR<1%) near genes involved in cell division and MAPK signaling. These cell type-specific effects are consistent with disease-specific changes in the composition of the CD4+ population and suggest that shifts in the proportion of CD4+ subtypes can be monitored at CpGs with subtype-specific DNA methylation patterns.
Author Summary
We have analyzed DNA methylation, an epigenetic modification that influences gene expression, in lupus patients and control subjects. Our analysis was run in three different immune cell types, T-cells, B-cells, and monocytes, to discern common epigenetic effects in lupus from cell type-specific effects. We have identified a lupus-related reduction in methylation around genes that respond to interferon, a cytokine that induces inflammation in response to pathogens. This hypomethylation suggests that lupus patients are hypersensitive to interferon, as DNA methylation is typically an inhibitor of gene expression. We also find that this hypersensitivity is preserved in lupus patients beyond active stages of the disease, and this may help explain the chronic, recurrent nature of the disease. In addition, we have identified DNA methylation changes in T-cells that suggest an alteration in the proportions of these cells in lupus patients, which may help explain the disease process.
PMCID: PMC3738443  PMID: 23950730
20.  PTPN22 Association in Systemic Lupus Erythematosus (SLE) with Respect to Individual Ancestry and Clinical Sub-Phenotypes 
PLoS ONE  2013;8(8):e69404.
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.
PMCID: PMC3737240  PMID: 23950893
21.  Variable association of reactive intermediate genes with systemic lupus erythematosus (SLE) in populations with different African ancestry 
The Journal of rheumatology  2013;40(6):842-849.
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.
PMCID: PMC3735344  PMID: 23637325
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
22.  Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci 
Nature genetics  2008;40(2):204-210.
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.
PMCID: PMC3712260  PMID: 18204446
23.  A Genome-wide Association Study of Host Genetic Determinants of the Antibody Response to Anthrax Vaccine Adsorbed 
Vaccine  2012;30(32):4778-4784.
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.
PMCID: PMC3387748  PMID: 22658931
Anthrax vaccines; Bacillus anthracis; bacterial vaccines; vaccination; Genome-wide association study
24.  Immune opsonins modulate BLyS/BAFF release in a receptor-specific fashion* 
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.
PMCID: PMC3684394  PMID: 18606652
Fc Receptors; Monocytes/Macrophages; Human; Autoimmunity
25.  Evaluation of TRAF6 in a Large Multi-Ancestral Lupus Cohort 
Arthritis and Rheumatism  2012;64(6):1960-1969.
Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with significant immune system aberrations resulting from complex heritable genetics as well as environmental factors. TRAF6 is a candidate gene for SLE, which has a major role in several signaling pathways that are important for immunity and organ development.
Fifteen single-nucleotide polymorphisms (SNPs), across TRAF6 were evaluated in 7,490 SLE and 6,780 control subjects from different ancestries. Population-based case-control association analyses and meta-analyses were performed. P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated.
Evidence of associations in multiple SNPs was detected. The best overall p values were obtained for SNPs rs5030437 and rs4755453 (p=7.85×10−5 and p=4.73×10−5, respectively) without significant heterogeneity among populations (p=0.67 and p=0.50 in Q-statistic). In addition, rs540386 previously reported to be associated with RA was found to be in LD with these two SNPs (r2= 0.95) and demonstrated evidence of association with SLE in the same direction (meta-analysis p=9.15×10−4, OR=0.89, 95%CI=0.83–0.95). Thrombocytopenia improved the overall results in different populations (meta-analysis p=1.99×10−6, OR=0.57, 95%CI=0.45–0.72, for rs5030470). Finally evidence of family based association in 34 African-American pedigrees with the presence of thrombocytopenia were detected in one available SNP rs5030437 with Z score magnitude of 2.28 (p=0.02) under a dominant model.
Our data indicate the presence of association of TRAF6 with SLE in agreement with the previous report of association with RA. These data provide further support for the involvement of TRAF6 in the pathogenesis of autoimmunity.
PMCID: PMC3380425  PMID: 22231568
TRAF6; polymorphism; systemic lupus erythematosus

Results 1-25 (76)