PMCC PMCC

Search tips
Search criteria

Advanced
Results 1-25 (116)
 

Clipboard (0)
None

Select a Filter Below

Journals
more »
Authors
more »
Year of Publication
more »
1.  The Effect of Inversion at 8p23 on BLK Association with Lupus in Caucasian Population 
PLoS ONE  2014;9(12):e115614.
To explore the potential influence of the polymorphic 8p23.1 inversion on known autoimmune susceptibility risk at or near BLK locus, we validated a new bioinformatics method that utilizes SNP data to enable accurate, high-throughput genotyping of the 8p23.1 inversion in a Caucasian population. Methods: Principal components analysis (PCA) was performed using markers inside the inversion territory followed by k-means cluster analyses on 7416 European derived and 267 HapMaP CEU and TSI samples. A logistic regression conditional analysis was performed. Results: Three subgroups have been identified; inversion homozygous, heterozygous and non-inversion homozygous. The status of inversion was further validated using HapMap samples that had previously undergone Fluorescence in situ hybridization (FISH) assays with a concordance rate of above 98%. Conditional analyses based on the status of inversion were performed. We found that overall association signals in the BLK region remain significant after controlling for inversion status. The proportion of lupus cases and controls (cases/controls) in each subgroup was determined to be 0.97 for the inverted homozygous group (1067 cases and 1095 controls), 1.12 for the inverted heterozygous group (1935 cases 1717 controls) and 1.36 for non-inverted subgroups (924 cases and 678 controls). After calculating the linkage disequilibrium between inversion status and lupus risk haplotype we found that the lupus risk haplotype tends to reside on non-inversion background. As a result, a new association effect between non-inversion status and lupus phenotype has been identified ((p = 8.18×10−7, OR = 1.18, 95%CI = 1.10–1.26). Conclusion: Our results demonstrate that both known lupus risk haplotype and inversion status act additively in the pathogenesis of lupus. Since inversion regulates expression of many genes in its territory, altered expression of other genes might also be involved in the development of lupus.
doi:10.1371/journal.pone.0115614
PMCID: PMC4278715  PMID: 25545785
2.  Whole Exome Sequencing Reveals Overlap Between Macrophage Activation Syndrome in Systemic Juvenile Idiopathic Arthritis and Familial Hemophagocytic Lymphohistiocytosis 
Objective
Macrophage activation syndrome (MAS), a life-threatening complication of systemic Juvenile Idiopathic Arthritis (SJIA), resembles Familial Hemophagocytic Lymphohistiocytosis (FHLH), a constellation of autosomal recessive immune disorders resulting from deficiency in cytolytic pathway proteins. We hypothesized that MAS predisposition in SJIA could be attributed to rare gene sequence variants affecting the cytotolytic pathway.
Methods
Whole exome sequencing (WES) was used in 14 SJIA/MAS patients and their parents to identify protein altering SNPs/indels in the known HLH-associated genes. To discover new candidate genes, the entire WES data were filtered to identify protein altering, rare recessive homozygous, compound heterozygous, and de novo variants with the potential to affect the cytolytic pathway.
Results
Heterozygous protein-altering rare variants in the known genes (LYST, MUNC13-4, and STXBP2) were found in 5 of 14 SJIA/MAS patients (35.7%). This was in contrast to only 4 variants in 4 of 29 (13,7%) SJIA patients without MAS. Homozygosity and compound heterozygosity analysis applied to the entire WES data in SJIAMAS, revealed 3 recessive pairs in 3 genes, and 76 compound heterozygotes in 75 genes. We also identified 22 heterozygous rare protein altering variants that occurred in at least two patients. Many of the identified genes encode proteins with a role in actin and microtubule reorganization and vesicle-mediated transport. “Cellular assembly and organization” was the top cellular function category based on Ingenuity Pathways Analysis (p<3.10E-05).
Conclusion
WES performed in SJIA/MAS patients identified rare protein altering variants in the known HLH associated genes as well as new candidate genes.
doi:10.1002/art.38793
PMCID: PMC4321811  PMID: 25047945
3.  Genome-wide association analysis of eosinophilic esophagitis provides insight into the tissue specificity of this allergic disease 
Nature genetics  2014;46(8):895-900.
Eosinophilic esophagitis (EoE) is a chronic inflammatory disorder associated with allergic hypersensitivity to food. We interrogated >1.5 million genetic variants in European EoE cases and subsequently in a multi-site cohort with local and out-of-study control subjects. In addition to replication of the 5q22 locus (meta-analysis p = 1.9×10−16), we identified association at 2p23 (encoding CAPN14, p = 2.5×10−10). CAPN14 was specifically expressed in the esophagus, dynamically upregulated as a function of disease activity and genetic haplotype and after exposure of epithelial cells to IL-13, and located in an epigenetic hotspot modified by IL-13. There was enriched esophageal expression for the genes neighboring the top 208 EoE sequence variants. Multiple allergic sensitization loci were associated with EoE susceptibility (4.8×10−2 < p < 5.1×10−11). We propose a model that elucidates the tissue specific nature of EoE that involves the interplay of allergic sensitization with an EoE-specific, IL-13–inducible esophageal response involving CAPN14.
doi:10.1038/ng.3033
PMCID: PMC4121957  PMID: 25017104
4.  Lupus risk variants in the PXK locus alter B-cell receptor internalization 
Frontiers in Genetics  2015;5:450.
Genome wide association studies have identified variants in PXK that confer risk for humoral autoimmune diseases, including systemic lupus erythematosus (SLE or lupus), rheumatoid arthritis and more recently systemic sclerosis. While PXK is involved in trafficking of epidermal growth factor Receptor (EGFR) in COS-7 cells, mechanisms linking PXK to lupus pathophysiology have remained undefined. In an effort to uncover the mechanism at this locus that increases lupus-risk, we undertook a fine-mapping analysis in a large multi-ancestral study of lupus patients and controls. We define a large (257kb) common haplotype marking a single causal variant that confers lupus risk detected only in European ancestral populations and spans the promoter through the 3′ UTR of PXK. The strongest association was found at rs6445972 with P < 4.62 × 10−10, OR 0.81 (0.75–0.86). Using stepwise logistic regression analysis, we demonstrate that one signal drives the genetic association in the region. Bayesian analysis confirms our results, identifying a 95% credible set consisting of 172 variants spanning 202 kb. Functionally, we found that PXK operates on the B-cell antigen receptor (BCR); we confirmed that PXK influenced the rate of BCR internalization. Furthermore, we demonstrate that individuals carrying the risk haplotype exhibited a decreased rate of BCR internalization, a process known to impact B cell survival and cell fate. Taken together, these data define a new candidate mechanism for the genetic association of variants around PXK with lupus risk and highlight the regulation of intracellular trafficking as a genetically regulated pathway mediating human autoimmunity.
doi:10.3389/fgene.2014.00450
PMCID: PMC4288052  PMID: 25620976
lupus; PXK; fine-mapping; B cells; BCR
5.  Phenome-wide association study (PheWAS) in EMR-linked pediatric cohorts, genetically links PLCL1 to speech language development and IL5-IL13 to Eosinophilic Esophagitis 
Frontiers in Genetics  2014;5:401.
Objective: We report the first pediatric specific Phenome-Wide Association Study (PheWAS) using electronic medical records (EMRs). Given the early success of PheWAS in adult populations, we investigated the feasibility of this approach in pediatric cohorts in which associations between a previously known genetic variant and a wide range of clinical or physiological traits were evaluated. Although computationally intensive, this approach has potential to reveal disease mechanistic relationships between a variant and a network of phenotypes.
Method: Data on 5049 samples of European ancestry were obtained from the EMRs of two large academic centers in five different genotyped cohorts. Recently, these samples have undergone whole genome imputation. After standard quality controls, removing missing data and outliers based on principal components analyses (PCA), 4268 samples were used for the PheWAS study. We scanned for associations between 2476 single-nucleotide polymorphisms (SNP) with available genotyping data from previously published GWAS studies and 539 EMR-derived phenotypes. The false discovery rate was calculated and, for any new PheWAS findings, a permutation approach (with up to 1,000,000 trials) was implemented.
Results: This PheWAS found a variety of common variants (MAF > 10%) with prior GWAS associations in our pediatric cohorts including Juvenile Rheumatoid Arthritis (JRA), Asthma, Autism and Pervasive Developmental Disorder (PDD) and Type 1 Diabetes with a false discovery rate < 0.05 and power of study above 80%. In addition, several new PheWAS findings were identified including a cluster of association near the NDFIP1 gene for mental retardation (best SNP rs10057309, p = 4.33 × 10−7, OR = 1.70, 95%CI = 1.38 − 2.09); association near PLCL1 gene for developmental delays and speech disorder [best SNP rs1595825, p = 1.13 × 10−8, OR = 0.65(0.57 − 0.76)]; a cluster of associations in the IL5-IL13 region with Eosinophilic Esophagitis (EoE) [best at rs12653750, p = 3.03 × 10−9, OR = 1.73 95%CI = (1.44 − 2.07)], previously implicated in asthma, allergy, and eosinophilia; and association of variants in GCKR and JAZF1 with allergic rhinitis in our pediatric cohorts [best SNP rs780093, p = 2.18 × 10−5, OR = 1.39, 95%CI = (1.19 − 1.61)], previously demonstrated in metabolic disease and diabetes in adults.
Conclusion: The PheWAS approach with re-mapping ICD-9 structured codes for our European-origin pediatric cohorts, as with the previous adult studies, finds many previously reported associations as well as presents the discovery of associations with potentially important clinical implications.
doi:10.3389/fgene.2014.00401
PMCID: PMC4235428  PMID: 25477900
PheWAS; ICD-9 code; genetic polymorphism
8.  Variants at multiple loci implicated in both innate and adaptive immune responses are associated with Sjögren’s syndrome 
Nature genetics  2013;45(11):10.1038/ng.2792.
Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (Pmeta=7.65×10−114), we establish associations with IRF5-TNPO3 (Pmeta=2.73×10−19), STAT4 (Pmeta=6.80×10−15), IL12A (Pmeta =1.17×10−10), FAM167A-BLK (Pmeta=4.97×10−10), DDX6-CXCR5 (Pmeta=1.10×10−8), and TNIP1 (Pmeta=3.30×10−8). Suggestive associations with Pmeta<5×10−5 were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome.
doi:10.1038/ng.2792
PMCID: PMC3867192  PMID: 24097067
9.  End-Stage Renal Disease in African Americans With Lupus Nephritis Is Associated With APOL1 
Objective
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.
Methods
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.
Results
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).
Conclusion
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.
doi:10.1002/art.38220
PMCID: PMC4002759  PMID: 24504811
10.  Excess female siblings and male fetal loss in families with systemic lupus erythematosus 
The Journal of rheumatology  2013;40(4):430-434.
Background
Systemic lupus erythematosus (SLE) occurs more frequently among woman than men. We undertook the present study to determine whether the male-female ratio in SLE families is different than expected by chance, and whether excess male fetal loss is found.
Methods
All SLE patients met the revised American College of Rheumatology Classification criteria, while SLE-unaffected subjects were shown not to satisfy these same criteria. Putative family relationships were confirmed by genetic testing. Pregnancy history was obtained from all subjects, including unrelated control woman. Adjusted Wald binomial confidence intervals (CI) were calculated for ratio of boys to girls in families and compared to the expected ratio of 1.06
Results
There were 2578 subjects with SLE with 6056 siblings. Considering all subjects, we found 3201 boys and 5434 girls (ratio=0.59, of 95% CI 0.576–0.602). When considering only the SLE-unaffected siblings, there were 2919 boys and 3137 girls (ratio=0.93, 95%CI 0.92–0.94). In both cases, the ratio of males-to-females is statistically different than the known birth rate. Among SLE patients as well as among their sisters and mothers there was an excess of male fetal loss compared to the controls.
Conclusion
Siblings of SLE patients are more likely to be girls than expected. This finding may be in part explained by excess male fetal loss, which is found among SLE patients and their first degree relatives.
doi:10.3899/jrheum.120643
PMCID: PMC3693848  PMID: 23378464
Systemic lupus erythematosus; sex ratio; fetal loss; pregnancy
11.  The struggle to find reliable results in exome sequencing data: filtering out Mendelian errors 
Next Generation Sequencing studies generate a large quantity of genetic data in a relatively cost and time efficient manner and provide an unprecedented opportunity to identify candidate causative variants that lead to disease phenotypes. A challenge to these studies is the generation of sequencing artifacts by current technologies. To identify and characterize the properties that distinguish false positive variants from true variants, we sequenced a child and both parents (one trio) using DNA isolated from three sources (blood, buccal cells, and saliva). The trio strategy allowed us to identify variants in the proband that could not have been inherited from the parents (Mendelian errors) and would most likely indicate sequencing artifacts. Quality control measurements were examined and three measurements were found to identify the greatest number of Mendelian errors. These included read depth, genotype quality score, and alternate allele ratio. Filtering the variants on these measurements removed ~95% of the Mendelian errors while retaining 80% of the called variants. These filters were applied independently. After filtering, the concordance between identical samples isolated from different sources was 99.99% as compared to 87% before filtering. This high concordance suggests that different sources of DNA can be used in trio studies without affecting the ability to identify causative polymorphisms. To facilitate analysis of next generation sequencing data, we developed the Cincinnati Analytical Suite for Sequencing Informatics (CASSI) to store sequencing files, metadata (eg. relatedness information), file versioning, data filtering, variant annotation, and identify candidate causative polymorphisms that follow either de novo, rare recessive homozygous or compound heterozygous inheritance models. We conclude the data cleaning process improves the signal to noise ratio in terms of variants and facilitates the identification of candidate disease causative polymorphisms.
doi:10.3389/fgene.2014.00016
PMCID: PMC3921572  PMID: 24575121
whole exome sequencing; variant filtering; next-generation sequencing; disease causative polymorphisms; Mendelian errors; Mendel errors; CASSI
12.  Comparison of the American-European Consensus Group Sjögren's syndrome classification criteria to newly proposed American College of Rheumatology criteria in a large, carefully characterized sicca cohort 
Annals of the rheumatic diseases  2013;73(1):10.1136/annrheumdis-2013-203845.
Objective
To compare the performance of the American-European Consensus Group (AECG) and the newly proposed American College of Rheumatology (ACR) classification criteria for Sjögren's syndrome in a well-characterized sicca cohort, given ongoing efforts to resolve discrepancies and weaknesses in the systems.
Methods
In a multidisciplinary clinic for the evaluation of sicca, we assessed features of salivary and lacrimal gland dysfunction and autoimmunity as defined by tests of both AECG and ACR criteria in 646 participants. Global gene expression profiles were compared in a subset of 180 participants.
Results
Application of the AECG and ACR criteria resulted in classification of 279 and 268 participants with SS, respectively. Both criteria were met by 244 participants (81%). In 26 of the 35 AECG+/ACR- participants, the minor salivary gland biopsy focal score was ≥1 (74%), while 9 had positive anti-Ro/La (26%). There were 24 AECG-/ACR+ who met ACR criteria mainly due to differences in the scoring of corneal staining. All patients with SS, regardless of classification, had similar gene expression profiles, which were distinct from the healthy controls.
Conclusion
The two sets of classification criteria yield concordant results in the majority of cases and gene expression profiling suggests that patients meeting either set of criteria are more similar to other SS participants than to healthy controls. Thus, there is no clear evidence for increased value of the new ACR criteria over the old AECG criteria from the clinical or biological perspective. It is our contention, supported by this report, that improvements in diagnostic acumen will require a more fundamental understanding of the pathogenic mechanisms than is at present available.
doi:10.1136/annrheumdis-2013-203845
PMCID: PMC3855629  PMID: 23968620
Sjögren's syndrome; Classification; Diagnosis
13.  EMR-linked GWAS study: investigation of variation landscape of loci for body mass index in children 
Frontiers in Genetics  2013;4:268.
Common variations at the loci harboring the fat mass and obesity gene (FTO), MC4R, and TMEM18 are consistently reported as being associated with obesity and body mass index (BMI) especially in adult population. In order to confirm this effect in pediatric population five European ancestry cohorts from pediatric eMERGE-II network (CCHMC-BCH) were evaluated.
Method: Data on 5049 samples of European ancestry were obtained from the Electronic Medical Records (EMRs) of two large academic centers in five different genotyped cohorts. For all available samples, gender, age, height, and weight were collected and BMI was calculated. To account for age and sex differences in BMI, BMI z-scores were generated using 2000 Centers of Disease Control and Prevention (CDC) growth charts. A Genome-wide association study (GWAS) was performed with BMI z-score. After removing missing data and outliers based on principal components (PC) analyses, 2860 samples were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and BMI was tested using linear regression adjusting for age, gender, and PC by cohort. The effects of SNPs were modeled assuming additive, recessive, and dominant effects of the minor allele. Meta-analysis was conducted using a weighted z-score approach.
Results: The mean age of subjects was 9.8 years (range 2–19). The proportion of male subjects was 56%. In these cohorts, 14% of samples had a BMI ≥95 and 28 ≥ 85%. Meta analyses produced a signal at 16q12 genomic region with the best result of p = 1.43 × 10-7 [p(rec) = 7.34 × 10-8) for the SNP rs8050136 at the first intron of FTO gene (z = 5.26) and with no heterogeneity between cohorts (p = 0.77). Under a recessive model, another published SNP at this locus, rs1421085, generates the best result [z = 5.782, p(rec) = 8.21 × 10-9]. Imputation in this region using dense 1000-Genome and Hapmap CEU samples revealed 71 SNPs with p < 10-6, all at the first intron of FTO locus. When hetero-geneity was permitted between cohorts, signals were also obtained in other previously identified loci, including MC4R (rs12964056, p = 6.87 × 10-7, z = -4.98), cholecystokinin CCK (rs8192472, p = 1.33 × 10-6, z = -4.85), Interleukin 15 (rs2099884, p = 1.27 × 10-5, z = 4.34), low density lipoprotein receptor-related protein 1B [LRP1B (rs7583748, p = 0.00013, z = -3.81)] and near transmembrane protein 18 (TMEM18) (rs7561317, p = 0.001, z = -3.17). We also detected a novel locus at chromosome 3 at COL6A5 [best SNP = rs1542829, minor allele frequency (MAF) of 5% p = 4.35 × 10-9, z = 5.89].
Conclusion: An EMR linked cohort study demonstrates that the BMI-Z measurements can be successfully extracted and linked to genomic data with meaningful confirmatory results. We verified the high prevalence of childhood rate of overweight and obesity in our cohort (28%). In addition, our data indicate that genetic variants in the first intron of FTO, a known adult genetic risk factor for BMI, are also robustly associated with BMI in pediatric population.
doi:10.3389/fgene.2013.00268
PMCID: PMC3847941  PMID: 24348519
BMI; obesity; polymorphism; GWAS
14.  Two Independent Functional Risk Haplotypes in TNIP1 are Associated with Systemic Lupus Erythematosus 
Arthritis and rheumatism  2012;64(11):3695-3705.
Objective
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.
Methods
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.
Results
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.
Conclusion
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.
doi:10.1002/art.34642
PMCID: PMC3485412  PMID: 22833143
15.  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.
Objective
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.
Methods
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.
Results
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.
Conclusion
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.
doi:10.1002/art.34650
PMCID: PMC3485439  PMID: 22886787
16.  Comparison of autoantibody specificities between traditional and bead-based assays in a large, diverse collection of SLE patients and family members 
Arthritis and rheumatism  2012;64(11):3677-3686.
Objective
The replacement of standard immunofluorescence anti-nuclear antibody (ANA) methods with bead-based assays is a new clinical option. A large, multi-racial cohort of SLE patients, blood relatives and unaffected control individuals was evaluated for familial aggregation and subset clustering of autoantibodies by high-throughput serum screening technology and traditional methods.
Methods
Serum samples (1,540 SLE patients, 1,127 unaffected relatives, and 906 healthy, population-based controls) were analyzed for SLE autoantibodies using a bead-based assay, immunofluorescence, and immunodiffusion. Autoantibody prevalence, disease sensitivity, clustering, and association with standard immunodiffusion results were evaluated.
Results
ANA frequency in SLE patient sera were 89%, 73%, and 67% by BioPlex 2200 and 94%, 84%, and 86% by immunofluorescence in African-American, Hispanic, and European-American patients respectively. 60kD Ro, La, Sm, nRNP A, and ribosomal P prevalence were compared across assays, with sensitivities ranging from 0.92 to 0.83 and specificities ranging from 0.90 to 0.79. Cluster autoantibody analysis showed association of three subsets: 1) 60kD Ro, 52kD Ro and La, 2) spliceosomal proteins, and 3) dsDNA, chromatin, and ribosomal P. Familial aggregation of Sm/RNP, ribosomal P, and 60kD Ro in SLE patient sibling pairs was observed (p ≤ 0.004). Simplex pedigree patients had a greater prevalence for dsDNA (p=0.0003) and chromatin (p=0.005) autoantibodies than multiplex patients.
Conclusion
ANA frequencies detected by a bead-based assay are lower in European-American SLE patients compared to immunofluorescence. These assays have strong positive predictive values across racial groups, provide useful information for clinical care, and provide unique insights to familial aggregation and autoantibody clustering.
doi:10.1002/art.34651
PMCID: PMC3490432  PMID: 23112091
systemic lupus erythematosus; autoantibodies; ancestry
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.
doi:10.1371/journal.pgen.1003870
PMCID: PMC3794920  PMID: 24130510
18.  Familial Aggregation of High Tumor Necrosis Factor Alpha Levels in Systemic Lupus Erythematosus 
Systemic lupus erythematosus (SLE) patients frequently have high circulating tumor necrosis factor alpha (TNF-α) levels. We explored circulating TNF-α levels in SLE families to determine whether high levels of TNF-α were clustered in a heritable pattern. We measured TNF-α in 242 SLE patients, 361 unaffected family members, 23 unaffected spouses of SLE patients, and 62 unrelated healthy controls. Familial correlations and relative recurrence risk rates for the high TNF-α trait were assessed. SLE-affected individuals had the highest TNF-α levels, and TNF-α was significantly higher in unaffected first degree relatives than healthy unrelated subjects (P = 0.0025). No Mendelian patterns were observed, but 28.4% of unaffected first degree relatives of SLE patients had high TNF-α levels, resulting in a first degree relative recurrence risk of 4.48 (P = 2.9 × 10−5). Interestingly, the median TNF-α value in spouses was similar to that of the first degree relatives. Concordance of the TNF-α trait (high versus low) in SLE patients and their spouses was strikingly high at 78.2%. These data support a role for TNF-α in SLE pathogenesis, and TNF-α levels may relate with heritable factors. The high degree of concordance in SLE patients and their spouses suggests that environmental factors may also play a role in the observed familial aggregation.
doi:10.1155/2013/267430
PMCID: PMC3800640  PMID: 24187561
20.  Large Scale Analysis of Tumor Necrosis Factor Alpha Levels in Systemic Lupus Erythematosus 
Arthritis and rheumatism  2012;64(9):2947-2952.
Background
SLE disease manifestations are highly variable between patients, and the prevalence of individual clinical features differs significantly by ancestry. Serum tumor necrosis factor alpha (TNF-α) is elevated in some SLE patients, and may play a role in disease pathogenesis. We detected associations between serum TNF-α, clinical manifestations, autoantibodies, and serum IFN-α in a large multi-ancestral SLE cohort.
Methods
We studied serum TNF-α in 653 SLE patients, including 214 African-American, 298 European-Americans and 141 Hispanic-American subjects. TNF-α was measured using ELISA, and IFN-α was measured with a functional reporter cell assay. Stratified and multivariate analyses were used to detect associations in each ancestral background separately, with meta-analysis when appropriate.
Results
Serum TNF-α levels were significantly higher in SLE patients than in nonautoimmune controls (p<5.0×10−3 for each ancestral background). High serum TNF-α was positively correlated with high serum IFN-α when tested in the same sample across all ancestral backgrounds (meta-analysis OR=1.8, p=1.2×10−3). While serum TNF-α levels alone did not differ significantly between SLE patients of different ancestral backgrounds, the proportion of patients with concurrently high TNF-α and high IFN-α was highest in African-Americans and lowest in European-Americans (p=5.0×10−3). Serum TNF-α was not associated with autoantibodies, clinical criteria for the diagnosis of SLE, or age at time of sample.
Conclusions
Serum TNF-α levels are high in many SLE patients, and we observed a positive correlation between serum TNF-α and IFN-α. These data support a role for TNF-α in SLE pathogenesis across all ancestral backgrounds, and suggest important cytokine subgroups within the disease.
doi:10.1002/art.34483
PMCID: PMC3396783  PMID: 22488302
systemic lupus erythematosus; tumor necrosis factor alpha; autoantibodies, ancestry
21.  Evidence of Dynamically Dysregulated Gene Expression Pathways in Hyperresponsive B Cells from African American Lupus Patients 
PLoS ONE  2013;8(8):e71397.
Recent application of gene expression profiling to the immune system has shown a great potential for characterization of complex regulatory processes. It is becoming increasingly important to characterize functional systems through multigene interactions to provide valuable insights into differences between healthy controls and autoimmune patients. Here we apply an original systematic approach to the analysis of changes in regulatory gene interconnections between in Epstein-Barr virus transformed hyperresponsive B cells from SLE patients and normal control B cells. Both traditional analysis of differential gene expression and analysis of the dynamics of gene expression variations were performed in combination to establish model networks of functional gene expression. This Pathway Dysregulation Analysis identified known transcription factors and transcriptional regulators activated uniquely in stimulated B cells from SLE patients.
doi:10.1371/journal.pone.0071397
PMCID: PMC3744560  PMID: 23977035
22.  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.
doi:10.1371/journal.pone.0069404
PMCID: PMC3737240  PMID: 23950893
23.  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.
Objective
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.
Methods
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.
Results
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.
Conclusion
These results suggest distinct patterns of association with SLE in African-derived populations; specific loci may be more strongly associated within select population groups.
doi:10.3899/jrheum.120989
PMCID: PMC3735344  PMID: 23637325
systemic lupus erythematosus; African Americans; genetic association studies; oxygen compounds; single nucleotide polymorphism
24.  Trans-Ancestral Studies Fine Map the SLE-Susceptibility Locus TNFSF4 
PLoS Genetics  2013;9(7):e1003554.
We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10−34, OR = 1.43[1.26–1.60]) and rs1234317-T (P = 1.16×10−28, OR = 1.38[1.24–1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.
Author Summary
Systemic lupus erythematosus (SLE/lupus) is a complex disease in which the body's immune cells cause inflammation in one or more systems to cause the associated morbidity. Hormones, the environment and genes are all causal contributors to SLE and over the past several years the genetic component of SLE has been firmly established. Several genes which are regulators of the immune system are associated with disease risk. We have established one of these, the tumour-necrosis family superfamily member 4 (TNFSF4) gene, as a lupus susceptibility gene in Northern Europeans. A major obstacle in pinpointing the marker(s) at TNFSF4 which best explain the risk of SLE has been the strong correlation (linkage disequilibrium, LD) between adjacent markers across the TNFSF4 region in this population. To address this, we have typed polymorphisms in several populations in addition to the European groups. The mixed ancestry of these populations gives a different LD pattern than that found in Europeans, presenting a method of pinpointing the section of the TNFSF4 region which results in SLE susceptibility. The Non-European populations have allowed identification of a polymorphism likely to regulate expression of TNFSF4 to increase susceptibility to SLE.
doi:10.1371/journal.pgen.1003554
PMCID: PMC3715547  PMID: 23874208
25.  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.
doi:10.1038/ng.81
PMCID: PMC3712260  PMID: 18204446

Results 1-25 (116)