Local Ancestry at the Major Histocompatibility Complex Region is Not a Major Contributor to Disease Heterogeneity in a Multiethnic Lupus Cohort

OBJECTIVE

Systemic lupus erythematosus (SLE) is an autoimmune disease resulting in debilitating clinical manifestations that vary in severity by race and ethnicity with a disproportionate burden in African American, Mestizo, and Asian populations compared with populations of European descent. Differences in global and local genetic ancestry may shed light on the underlying mechanisms contributing to these disparities, including increased prevalence of lupus nephritis, younger age of symptom onset, and presence of autoantibodies.

METHODS

A total of 1,139 European, African American, and Mestizos patients with SLE were genotyped using the Affymetrix LAT1 World array. Global ancestry proportions were estimated using ADMIXTURE, and local ancestry was estimated using RFMIXv2.0. We investigated associations between lupus nephritis, age at onset, and autoantibody status with both global and local ancestry proportions within the Major Histocompatibility Complex region.

RESULTS

Our results showed small effect sizes that did not meet the threshold for statistical significance for global or local ancestry proportions in either African American or Mestizo patients with SLE who presented with the clinical manifestations of interest compared with those who did not.

CONCLUSION

These findings suggest that local genetic ancestry within the Major Histocompatibility Complex region is not a major contributor to these SLE manifestations among patients with SLE from admixed populations.

Screening of Environmental Chemicals to Characterize Exposures in Participants With Systemic Lupus Erythematosus

OBJECTIVE

There is a need to characterize exposures associated with the pathogenesis of systemic lupus erythematosus (SLE). In this pilot study, we explore a hypothesis-free approach that can measure thousands of exogenous chemicals in blood ("exposome") in patients with SLE and unaffected controls.

METHODS

This cross-sectional study analyzed a cohort of patients with prevalent SLE (n = 285) and controls (n = 106). Plasma was analyzed by liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF/MS). Mass spectrometry features present in at least 25% of all samples were selected for association analysis (n = 2,737). Features were matched to potential chemicals using available databases. Association analysis of abundances of features with SLE status was performed, adjusting for age and sex. We also explored features associated with SLE phenotypes, sociodemographic factors, and current medication use.

RESULTS

We found 30 features significantly associated with SLE status (Bonferroni P < 0.05). Of these, seven matched chemical names based on databases. These seven features included phthalate metabolites, a formetanate metabolite, and eugenol. The abundance of acid pesticides differed between patients with SLE and controls (Bonferroni P < 0.05). Two unmatched features were associated with a history of lupus nephritis, and one with anti-double-stranded DNA antibody production (Bonferroni P < 0.05). Seventeen features varied by self-reported race and ethnicity, including a polyfluoroalkyl substance (analysis of variance P < 1.69 × 10-5 ). Eleven features correlated with antimalarials, 6 with mycophenolate mofetil, and 29 with prednisone use.

CONCLUSION

This proof-of-concept study demonstrates that LC-QTOF/MS is a powerful tool that agnostically detects circulating exogenous compounds. These analyses can generate hypotheses of disease-related exposures for future prospective, longitudinal studies.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis

OBJECTIVE

We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment.

METHODS

Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software.

RESULTS

We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples.

CONCLUSION

We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

Mer-tyrosine kinase: a novel susceptibility gene for SLE related end-stage renal disease

OBJECTIVE

Lupus nephritis (LN) is a common and severe manifestation of SLE. The genetic risk for nephritis and progression to end-stage renal disease (ESRD) in patients with LN remains unclear. Herein, we aimed to identify novel genetic associations with LN, focusing on subphenotypes and ESRD.

METHODS

We analysed genomic data on 958 patients with SLE (discovery cohort: LN=338) with targeted sequencing data from 1832 immunological pathway genes. We used an independent multiethnic cohort comprising 1226 patients with SLE (LN=603) as a replication dataset. Detailed functional annotation and functional epigenomic enrichment analyses were applied to predict functional effects of the candidate variants.

RESULTS

A genetic variant (rs56097910) within the MERTK gene was associated with ESRD in both cohorts, meta-analysis OR=5.4 (2.8 to 10.6); p=1.0×10-6. We observed decreased methylation levels in peripheral blood cells from SLE patients with ESRD, compared with patients without renal SLE (p=2.7×10-4), at one CpG site (cg16333401) in close vicinity to the transcription start site of MERTK and located in a DNAse hypersensitivity region in T and B cells. Rs56097910 is linked to altered MERTK expression in kidney tissue in public eQTL databases. Two loci were replicated for association with proliferative LN: PRDM1 (rs6924535, pmeta=1.6×10-5, OR=0.58) and APOA1BP (NAXE) (rs942960, pmeta=1.2×10-5, OR=2.64).

CONCLUSION

We identified a novel genetic risk locus, MERTK, associated with SLE-ESRD using the data from two large SLE cohorts. Through DNA methylation analysis and functional annotation, we showed that the risk could be mediated through regulation of gene expression. Our results suggest that variants in the MERTK gene are important for the risk of developing SLE-ESRD and suggest a role for PRDM1 and APOA1BP in proliferative LN.

Dynamics of Methylation of CpG Sites Associated With Systemic Lupus Erythematosus Subtypes in a Longitudinal Cohort

OBJECTIVE

Findings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points.

METHODS

A total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components.

RESULTS

The majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P < 2.70 × 10^-5 ), and 15 CpGs had at least 1 significant association with immunosuppressant use.

CONCLUSION

In this longitudinal SLE cohort, we identified a subset of SLE subtype-associated CpGs that remained stable over time and may be useful as biomarkers of disease subtypes. Another subset of SLE subtype-associated CpGs changed at a higher proportion compared to the genome-wide methylome. Additional studies are needed to understand the etiology and impact of these changes on methylation of SLE-associated CpGs.

A Human MSH6 Germline Variant Associated With Systemic Lupus Erythematosus Induces Lupus-like Disease in Mice

OBJECTIVE

To determine if single-nucleotide polymorphisms (SNPs) in DNA repair genes are enriched in individuals with systemic lupus erythematosus (SLE) and if they are sufficient to confer a disease phenotype in a mouse model.

METHODS

Human exome chip data of 2499 patients with SLE and 1230 healthy controls were analyzed to determine if variants in 10 different mismatch repair genes (MSH4, EXO1, MSH2, MSH6, MLH1, MSH3, POLH, PMS2, ML3, and APEX2) were enriched in individuals with SLE. A mouse model of the MSH6 SNP, which was found to be enriched in individuals with SLE, was created using CRISPR/Cas9 gene targeting. Wildtype mice and mice heterozygous and homozygous for the MSH6 variant were then monitored for 2 years for the development of autoimmune phenotypes, including the presence of high levels of antinuclear antibodies (ANA). Additionally, somatic hypermutation frequencies and spectra of the intronic region downstream of the VH J558-rearranged JH4 immunoglobulin gene was characterized from Peyer's patches.

RESULTS

Based on the human exome chip data, the MSH6 variant (rs63750897, p.Ser503Cys) is enriched among patients with SLE versus controls after we corrected for ancestry (odds ratio = 8.39, P = 0.0398). Mice homozygous for the MSH6 variant (Msh6S502C/S502C ) harbor significantly increased levels of ANA. Additionally, the Msh6S502C/S502C mice display a significant increase in the infiltration of CD68+ cells (a marker for monocytes and macrophages) into the lung alveolar space as well as apoptotic cells. Furthermore, characterization of somatic hypermutation in these mice reveals an increase in the DNA polymerase η mutational signature.

CONCLUSION

An MSH6 mutation that is enriched in humans diagnosed with lupus was identified. Mice harboring this Msh6 mutation develop increased autoantibodies and an inflammatory lung disease. These results suggest that the human MSH6 variant is linked to the development of SLE.

P2RY8 variants in lupus patients uncover a role for the receptor in immunological tolerance

B cell self-tolerance is maintained through multiple checkpoints, including restraints on intracellular signaling and cell trafficking. P2RY8 is a receptor with established roles in germinal center (GC) B cell migration inhibition and growth regulation. Somatic P2RY8 variants are common in GC-derived B cell lymphomas. Here, we identify germline novel or rare P2RY8 missense variants in lupus kindreds or the related antiphospholipid syndrome, including a "de novo" variant in a child with severe nephritis. All variants decreased protein expression, F-actin abundance, and GPCR-RhoA signaling, and those with stronger effects increased AKT and ERK activity and cell migration. Remarkably, P2RY8 was reduced in B cell subsets from some SLE patients lacking P2RY8 gene variants. Low P2RY8 correlated with lupus nephritis and increased age-associated B cells and plasma cells. By contrast, P2RY8 overexpression in cells and mice restrained plasma cell development and reinforced negative selection of DNA-reactive developing B cells. These findings uncover a role of P2RY8 in immunological tolerance and lupus pathogenesis.

Variants in BANK1 are associated with lupus nephritis of European ancestry

The genetic background of lupus nephritis (LN) has not been completely elucidated. We performed a case-only study of 2886 SLE patients, including 947 (33%) with LN. Renal biopsies were available from 396 patients. The discovery cohort (Sweden, n = 1091) and replication cohort 1 (US, n = 962) were genotyped on the Immunochip and replication cohort 2 (Denmark/Norway, n = 833) on a custom array. Patients with LN, proliferative nephritis, or LN with end-stage renal disease were compared with SLE without nephritis. Six loci were associated with LN (p < 1 × 10⁻⁴, NFKBIA, CACNA1S, ITGA1, BANK1, OR2Y, and ACER3) in the discovery cohort. Variants in BANK1 showed the strongest association with LN in replication cohort 1 (p = 9.5 × 10⁻⁴) and proliferative nephritis in a meta-analysis of discovery and replication cohort 1. There was a weak association between BANK1 and LN in replication cohort 2 (p = 0.052), and in the meta-analysis of all three cohorts the association was strengthened (p = 2.2 × 10⁻⁷). DNA methylation data in 180 LN patients demonstrated methylation quantitative trait loci (meQTL) effects between a CpG site and BANK1 variants. To conclude, we describe genetic variations in BANK1 associated with LN and evidence for genetic regulation of DNA methylation within the BANK1 locus. This indicates a role for BANK1 in LN pathogenesis.

Large-Scale Imputation of KIR Copy Number and HLA Alleles in North American and European Psoriasis Case-Control Cohorts Reveals Association of Inhibitory KIR2DL2 With Psoriasis

Killer cell immunoglobulin-like receptors (KIR) regulate immune responses in NK and CD8+ T cells via interaction with HLA ligands. KIR genes, including KIR2DS1, KIR3DL1, and KIR3DS1 have previously been implicated in psoriasis susceptibility. However, these previous studies were constrained to small sample sizes, in part due to the time and expense required for direct genotyping of KIR genes. Here, we implemented KIR*IMP to impute KIR copy number from single-nucleotide polymorphisms (SNPs) on chromosome 19 in the discovery cohort (n=11,912) from the PAGE consortium, University of California San Francisco, and the University of Dundee, and in a replication cohort (n=66,357) from Kaiser Permanente Northern California. Stratified multivariate logistic regression that accounted for patient ancestry and high-risk HLA alleles revealed that KIR2DL2 copy number was significantly associated with psoriasis in the discovery cohort (p ≤ 0.05). The KIR2DL2 copy number association was replicated in the Kaiser Permanente replication cohort. This is the first reported association of KIR2DL2 copy number with psoriasis and highlights the importance of KIR genetics in the pathogenesis of psoriasis.

Cross-Tissue Transcriptomic Analysis Leveraging Machine Learning Approaches Identifies New Biomarkers for Rheumatoid Arthritis

There is an urgent need to identify biomarkers for diagnosis and disease activity monitoring in rheumatoid arthritis (RA). We leveraged publicly available microarray gene expression data in the NCBI GEO database for whole blood (N=1,885) and synovial (N=284) tissues from RA patients and healthy controls. We developed a robust machine learning feature selection pipeline with validation on five independent datasets culminating in 13 genes: TNFAIP6, S100A8, TNFSF10, DRAM1, LY96, QPCT, KYNU, ENTPD1, CLIC1, ATP6V0E1, HSP90AB1, NCL and CIRBP which define the RA score and demonstrate its clinical utility: the score tracks the disease activity DAS28 (p = 7e-9), distinguishes osteoarthritis (OA) from RA (OR 0.57, p = 8e-10) and polyJIA from healthy controls (OR 1.15, p = 2e-4) and monitors treatment effect in RA (p = 2e-4). Finally, the immunoblotting analysis of six proteins on an independent cohort confirmed two proteins, TNFAIP6/TSG6 and HSP90AB1/HSP90.

Transcriptomic analysis of immune cells in a multi-ethnic cohort of systemic lupus erythematosus patients identifies ethnicity- and disease-specific expression signatures

Systemic lupus erythematosus (SLE) is an autoimmune disease in which outcomes vary among different racial groups. We leverage cell-sorted RNA-seq data (CD14+ monocytes, B cells, CD4+ T cells, and NK cells) from 120 SLE patients (63 Asian and 57 White individuals) and apply a four-tier approach including unsupervised clustering, differential expression analyses, gene co-expression analyses, and machine learning to identify SLE subgroups within this multiethnic cohort. K-means clustering on each cell-type resulted in three clusters for CD4 and CD14, and two for B and NK cells. To understand the identified clusters, correlation analysis revealed significant positive associations between the clusters and clinical parameters including disease activity as well as ethnicity. We then explored differentially expressed genes between Asian and White groups for each cell-type. The shared differentially expressed genes across cells were involved in SLE or other autoimmune-related pathways. Co-expression analysis identified similarly regulated genes across samples and grouped these genes into modules. Finally, random forest classification of disease activity in the White and Asian cohorts showed the best classification in CD4+ T cells in White individuals. The results from these analyses will help stratify patients based on their gene expression signatures to enable SLE precision medicine.

Evaluation of a Genetic Risk Score for Diagnosis of Psoriatic Arthritis

Background

Diagnosis of psoriatic arthritis (PsA) can be challenging, resulting in delays that contribute to irreversible joint damage, reduced quality of life, and increased mortality.

Objective

Use genetic markers to develop and evaluate a PsA genetic risk score (GRS) for its ability to discriminate between psoriasis (PsO) only and PsO with PsA among a psoriatic cohort with full genome-wide genotype data.

Methods

Genome-wide single-nucleotide polymorphism genotyping was performed on 724 psoriatic patients. A set of 11 candidate risk genes previously shown to be preferentially associated with PsO or PsA were selected. To evaluate the cumulative effects of these risk loci, a PsA GRS was developed using an unweighted risk allele count (cGRS) and a weighted (wGRS) approach. Additional analyses included only human leukocyte antigen (HLA) risk alleles.

Results

The discriminative power attributable to each GRS was evaluated by calculating the areas under the receiver operator characteristic curve (AUROC). The AUROC for the wGRS is 56.2% versus 54.1% for the cGRS, and the AUROC for the HLA-only wGRS model was 56.9% versus 55.7% for the HLA-only cGRS.

Conclusion

The AUROC of 56.9% for HLA-only wGRS indicates that this approach has the greatest power in discriminating PsA from PsO among these models. Given that an AUROC of 56.9% is quite modest, this study suggests that using a small number of well-validated genetic loci provides limited predictive power for PsA, and that future approaches may benefit from using a larger number of genetic loci.

A phenotypic and genomics approach in a multi-ethnic cohort to subtype systemic lupus erythematosus

Systemic lupus erythematous (SLE) is a heterogeneous autoimmune disease in which outcomes vary among different racial groups. Here, we aim to identify SLE subgroups within a multiethnic cohort using an unsupervised clustering approach based on the American College of Rheumatology (ACR) classification criteria. We identify three patient clusters that vary according to disease severity. Methylation association analysis identifies a set of 256 differentially methylated CpGs across clusters, including 101 CpGs in genes in the Type I Interferon pathway, and we validate these associations in an external cohort. A cis-methylation quantitative trait loci analysis identifies 744 significant CpG-SNP pairs. The methylation signature is enriched for ethnic-associated CpGs suggesting that genetic and non-genetic factors may drive outcomes and ethnic-associated methylation differences. Our computational approach highlights molecular differences associated with clusters rather than single outcome measures. This work demonstrates the utility of applying integrative methods to address clinical heterogeneity in multifactorial multi-ethnic disease settings.

Increased alloreactive and autoreactive antihuman leucocyte antigen antibodies associated with systemic lupus erythematosus and rheumatoid arthritis

Objectives

Rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) disproportionately affect women during and following childbearing years. Antihuman leucocyte antigen (HLA) alloantibody responses are common in healthy parous women, and as these diseases are both linked with HLA and immune dysregulation, we sought to evaluate anti-HLA antibodies in RA and SLE.

Methods

Anti-HLA antibodies were measured among parous SLE cases (n=224), parous RA cases (n=202) and healthy parous controls (n=239) and compared with each other as well as with nulliparous female and male controls. Antibody specificities were identified and compared against subject HLA types to determine autoreactivity versus alloreactivity. The association of anti-HLA antibodies with clinical outcomes was evaluated.

Results

Levels and frequencies of anti-HLA antibodies were significantly higher among parous females with SLE (52%) or RA (46%) compared with controls (26%), and anti-HLA antibodies were also found among nulliparous females and males with SLE and RA. Autoreactive anti-HLA antibodies were observed among SLE and RA antibody-positive subjects, but not healthy controls, with the highest frequency of autoreactive anti-HLA antibodies found in the SLE subjects. Higher levels of anti-HLA antibodies were associated with nephritis among the nulliparous SLE cases (p<0.01). The presence of anti-class I HLA antibodies was associated with younger age at diagnosis among both the RA and SLE nulliparous cases.

Conclusions

Both autoreactive and alloreactive anti-HLA antibodies were found at high levels in RA and SLE subjects. These occurred even in the absence of alloexposure, particularly among SLE subjects and may be linked with disease severity.

Transcriptional landscape of epithelial and immune cell populations revealed through FACS-seq of healthy human skin

Human skin consists of multiple cell types, including epithelial, immune, and stromal cells. Transcriptomic analyses have previously been performed from bulk skin samples or from epithelial and immune cells expanded in cell culture. However, transcriptomic analysis of bulk skin tends to drown out expression signals from relatively rare cells while cell culture methods may significantly alter cellular phenotypes and gene expression profiles. To identify distinct transcriptomic profiles of multiple cell populations without substantially altering cell phenotypes, we employed a fluorescence activated cell sorting method to isolate keratinocytes, dendritic cells, CD4+ T effector cells, and CD8+ T effector cells from healthy skin samples, followed by RNA-seq of each cell population. Principal components analysis revealed distinct clustering of cell types across samples, while differential expression and coexpression network analyses revealed transcriptional profiles of individual cell populations distinct from bulk skin, most strikingly in the least abundant CD8+ T effector population. Our work provides a high resolution view of cutaneous cellular gene expression and suggests that transcriptomic profiling of bulk skin may inadequately capture the contribution of less abundant cell types.

Genome-wide profiling identifies associations between lupus nephritis and differential methylation of genes regulating tissue hypoxia and type 1 interferon responses

OBJECTIVE

Previous studies have shown that differential DNA methylation is associated with SLE susceptibility. How DNA methylation affects the clinical heterogeneity of SLE has not been fully defined. We conducted this study to identify differentially methylated CpG sites associated with nephritis among women with SLE.

METHODS

The methylation status of 428 229 CpG sites across the genome was characterised for peripheral blood cells from 322 women of European descent with SLE, 80 of whom had lupus nephritis, using the Illumina HumanMethylation450 BeadChip. Multivariable linear regression adjusting for population substructure and leucocyte cell proportions was used to identify differentially methylated sites associated with lupus nephritis. The influence of genetic variation on methylation status was investigated using data from a genome-wide association study of SLE. Pathway analyses were used to identify biological processes associated with lupus nephritis.

RESULTS

We identified differential methylation of 19 sites in 18 genomic regions that was associated with nephritis among patients with SLE (false discovery rate q<0.05). Associations for four sites in HIF3A, IFI44 and PRR4 were replicated when examining methylation data derived from CD4+ T cells collected from an independent set of patients with SLE. These associations were not driven by genetic variation within or around the genomic regions. In addition, genes associated with lupus nephritis in a prior genome-wide association study were not differentially methylated in this epigenome-wide study. Pathway analysis indicated that biological processes involving type 1 interferon responses and the development of the immune system were associated with nephritis in patients with SLE.

CONCLUSIONS

Differential methylation of genes regulating the response to tissue hypoxia and interferon-mediated signalling is associated with lupus nephritis among women with SLE. These findings have not been identified in genetic studies of lupus nephritis, suggesting that epigenome-wide association studies can help identify the genomic differences that underlie the clinical heterogeneity of SLE.

Psoriasis risk SNPs and their association with HIV-1 control

Human evolution has resulted in selection for genetic polymorphisms beneficial in the defense against pathogens. However, such polymorphisms may have the potential to heighten the risk of autoimmune disease. Here, we investigated whether psoriasis-associated single nucleotide polymorphisms influence host control of HIV-1 infection. We studied psoriasis and viral immune response variants in three HIV-positive cohorts: (1) HIV-1 controllers and non-controllers in the Study of the Consequences of the Protease Inhibitor Era (SCOPE) cohort (n=366), (2) Individuals with primary HIV infection in the Options cohort (n=675), and (3) HIV-positive injection drug users from the Urban Health Study (UHS) (n=987). We found a strong association of two psoriasis MHC variants, rs9264942 and rs3021366, with both HIV-1 controller status and viral load, and identified another Class III MHC variant rs9368699 to be strongly associated with viral load. A number of genetic variants outside the MHC (SOX5, TLR9, SDC4, PROX1, IL12B, TLR4, MBL-2, TYK2, IFIH1) demonstrated nominal significance. Overall, our data suggest that several psoriasis variants within the MHC have a robust impact on HIV-1 control, while variants outside the MHC require further investigation.

Genome-Wide Assessment of Differential DNA Methylation Associated with Autoantibody Production in Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is characterized by the development of autoantibodies associated with specific clinical manifestations. Previous studies have shown an association between differential DNA methylation and SLE susceptibility, but have not investigated SLE-related autoantibodies. Our goal was to determine whether DNA methylation is associated with production of clinically relevant SLE-related autoantibodies, with an emphasis on the anti-dsDNA autoantibody. In this study, we characterized the methylation status of 467,314 CpG sites in 326 women with SLE. Using a discovery and replication study design, we identified and replicated significant associations between anti-dsDNA autoantibody production and the methylation status of 16 CpG sites (p_discovery<1.07E-07 and p_replication<0.0029) in 11 genes. Associations were further investigated using multivariable regression to adjust for estimated leukocyte cell proportions and population substructure. The adjusted mean DNA methylation difference between anti-dsDNA positive and negative cases ranged from 1.2% to 19%, and the adjusted odds ratio for anti-dsDNA autoantibody production comparing the lowest and highest methylation tertiles ranged from 6.8 to 18.2. Differential methylation for these CpG sites was also associated with anti-SSA, anti-Sm, and anti-RNP autoantibody production. Overall, associated CpG sites were hypomethylated in autoantibody positive compared to autoantibody negative cases. Differential methylation of CpG sites within the major histocompatibility region was not strongly associated with autoantibody production. Genes with differentially methylated CpG sites represent multiple biologic pathways, and have not been associated with autoantibody production in genetic association studies. In conclusion, hypomethylation of CpG sites within genes from different pathways is associated with anti-dsDNA, anti-SSA, anti-Sm, and anti-RNP production in SLE, and these associations are not explained by genetic variation. Thus, studies of epigenetic mechanisms such as DNA methylation represent a complementary method to genetic association studies to identify biologic pathways that may contribute to the clinical heterogeneity of autoimmune diseases.

A polymorphism in TLR2 is associated with arterial thrombosis in a multiethnic population of patients with systemic lupus erythematosus

OBJECTIVE

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.

METHODS

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.

RESULTS

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

CONCLUSION

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.

NCR3/NKp30 contributes to pathogenesis in primary Sjogren's syndrome

Primary Sjögren's syndrome (pSS) is a chronic autoimmune disease characterized by a lymphocytic exocrinopathy. However, patients often have evidence of systemic autoimmunity, and they are at markedly increased risk for the development of non- Hodgkin's lymphoma. Similar to other autoimmune disorders, a strong interferon (IFN) signature is present among subsets of pSS patients, although the precise etiology remains uncertain. NCR3/NKp30 is a natural killer (NK)-specific activating receptor regulating the cross talk between NK and dendritic cells and type II IFN secretion. We performed a case-control study of genetic polymorphisms of the NCR3/NKp30 gene and found that rs11575837 (G>A) residing in the promoter was associated with reduced gene transcription and function as well as protection to pSS. We also demonstrated that circulating levels of NCR3/NKp30 were significantly increased among pSS patients compared with controls and correlated with higher NCR3/NKp30 but not CD16-dependent IFN-γ secretion by NK cells. Excess accumulation of NK cells in minor salivary glands correlated with the severity of the exocrinopathy. B7H6, the ligand of NKp30, was expressed by salivary epithelial cells. These findings suggest that NK cells may promote an NKp30-dependent inflammatory state in salivary glands and that blockade of the B7H6/NKp30 axis could be clinically relevant in pSS.

Differentiation of follicular helper T cells by salivary gland epithelial cells in primary Sjögren's syndrome

Follicular helper T cells (Tfh), which play a pivotal role in B cell activation and differentiation in lymphoid structures, secrete IL-21 whose augmented secretion is a hallmark of several autoimmune diseases. To decipher the cellular and molecular interactions occurring in salivary glands of patients suffering from primary Sjögren's syndrome (pSS), we investigated whether salivary gland epithelial cells (SGECs) were capable to induce Tfh differentiation. Co-cultures of naïve CD4(+) T cells and SGECs from both patients with pSS and controls were performed. Here, we report that IL-6 and ICOSL expression by SGECs contributes to naïve CD4(+) T differentiation into Tfh cells, as evidenced by their acquisition of a specific phenotype, characterized by Bcl-6, ICOS and CXCR5 expression and IL-21 secretion, but also but by their main functional feature: the capacity to enhance B lymphocytes survival. We demonstrated an increase of serum IL-21 with systemic activity. Finally, we analyzed the potential occurrence of a genetic association between IL-21 or IL-21R gene polymorphisms and pSS or elevated IL-21 secretion. This study, which demonstrates a direct induction of Tfh differentiation by SGECs, emphasizes a yet unknown pathogenic role of SGECs and suggests that Tfh and IL-21 might be relevant biomarkers and/or therapeutic targets in primary Sjögren's syndrome.

Assessing breast cancer risk models in Marin County, a population with high rates of delayed childbirth

INTRODUCTION

This study was designed to compare the Breast Cancer Risk Assessment Tool (BCRAT; Gail), International Breast Intervention Study (IBIS; Tyrer-Cuzick), and BRCAPRO breast cancer risk assessment models using data from the Marin Women's Study, a cohort of women within Marin County, California, with high rates of breast cancer, nulliparity, and delayed childbirth. Existing models have not been well-validated in these high-risk populations.

METHODS

Discrimination was assessed using the area under the receiver operating characteristic curve (AUC) and calibration by estimating the ratio of expected-to-observed (E/O) cases. The models were assessed using data from 12,843 participants, of whom 203 had developed cancer during a 5-year period. All tests of statistical significance were 2-sided.

RESULTS

The IBIS model achieved an AUC of 0.65 (95% confidence interval [CI], 0.61-0.68) compared with 0.62 (95% CI, 0.59-0.66) for BCRAT and 0.60 (95% CI, 0.56-0.63) for BRCAPRO. The corresponding estimated E/O ratios for the models were 1.08 (95% CI, 0.95-1.25), 0.81 (95% CI, 0.71-0.93), and 0.59 (95% CI, 0.52-0.68). In women with age at first birth > 30 years, the AUC for the IBIS, BCRAT, and BRCAPRO models was 0.69 (95% CI, 0.62-0.75), 0.63 (95% CI, 0.56-0.70), and 0.62 (95% CI, 0.56-0.68) and the E/O ratio was 1.15 (95% CI, 0.89-1.47), 0.81 (95% CI, 0.63-1.05), and 0.53 (95% CI, 0.41-0.68), respectively.

CONCLUSIONS

The IBIS model was well calibrated for the high-risk Marin mammography population and demonstrated the best calibration of the 3 models in nulliparous women. The IBIS model also achieved the greatest overall discrimination and displayed superior discrimination for women with age at first birth > 30 years.

Brief Report: Single-nucleotide polymorphisms in VKORC1 are risk factors for systemic lupus erythematosus in Asians

OBJECTIVE

The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. The present study was undertaken to investigate whether 33 established and novel single-nucleotide polymorphisms (SNPs) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis (DVT) in the general population are risk factors for SLE among Asian subjects.

METHODS

Patients in the discovery cohort were enrolled in 1 of 2 North American SLE cohorts. Patients in the replication cohort were enrolled in 1 of 4 Asian or 2 North American cohorts. We first genotyped 263 Asian patients with SLE and 357 healthy Asian control subjects for 33 SNPs in the discovery phase, and then genotyped 5 SNPs in up to an additional 1,496 patients and 993 controls in the replication phase. Patients were compared to controls for bivariate association with minor alleles. Principal components analysis was used to control for intra-Asian ancestry in the replication cohort.

RESULTS

Two genetic variants in the gene VKORC1 were highly significant in both the discovery and replication cohorts: rs9934438 (in the discovery cohort, odds ratio [OR] 2.45, P=2×10(-9); in the replication cohort, OR 1.54, P=4×10(-6)) and rs9923231 (in the discovery cohort, OR 2.40, P=6×10(-9); in the replication cohort, OR 1.53, P=5×10(-6)). These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: for rs9934438, OR 1.34, P=0.0029; for rs9923231, OR 1.34, P=0.0032.

CONCLUSION

Genetic variants in VKORC1, which are involved in vitamin K reduction and associated with DVT, correlate with SLE development in Asian subjects. These results suggest that there may be intersecting genetic pathways for the development of SLE and thrombosis.

Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti-dsDNA + SLE and were not associated with anti-dsDNA - SLE. In conclusion, we identified differential genetic associations with SLE based on anti-dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti-dsDNA - SLE.

Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor alpha therapy

OBJECTIVE

Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy.

METHODS

A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models.

RESULTS

Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99).

CONCLUSION

Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10⁻¹³. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and involvement of multiple organ systems. Although the cause of SLE remains unknown, several lines of evidence underscore the importance of genetic factors. As is true for most autoimmune diseases, a substantial genetic contribution to disease risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6. This region of the genome contains a large number of genes that participate in the immune response. However, the full contribution of this genomic region to SLE risk has not yet been defined. In the current study we characterize a large number of SLE patients and family members for approximately 2,000 MHC region variants to identify the specific genes that influence disease risk. Our results, for the first time, implicate four different MHC regions in SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.

European population substructure correlates with systemic lupus erythematosus endophenotypes in North Americans of European descent

Previous work has demonstrated that Northern and Southern European ancestries are associated with specific systemic lupus erythematosus (SLE) manifestations. In this study, 1855 SLE cases of European descent were genotyped for 4965 single-nucleotide polymorphisms and principal components analysis of genotype information was used to define population substructure. The first principal component (PC1) distinguished Northern from Southern European ancestry, PC2 differentiated Eastern from Western European ancestry and PC3 delineated Ashkenazi Jewish ancestry. Compared with Northern European ancestry, Southern European ancestry was associated with autoantibody production (odds ratio (OR)=1.40, 95% confidence interval (CI) 1.07-1.83) and renal involvement (OR 1.41, 95% CI 1.06-1.87), and was protective for discoid rash (OR=0.51, 95% CI 0.32-0.82) and photosensitivity (OR=0.74, 95% CI 0.56-0.97). Both serositis (OR=1.46, 95% CI 1.12-1.89) and autoantibody production (OR=1.38, 95% CI 1.06-1.80) were associated with Western compared to Eastern European ancestry. Ashkenazi Jewish ancestry was protective against neurologic manifestations of SLE (OR=0.62, 95% CI 0.40-0.94). Homogeneous clusters of cases defined by multiple PCs demonstrated stronger phenotypic associations. Genetic ancestry may contribute to the development of SLE endophenotypes and should be accounted for in genetic studies of disease characteristics.

European population substructure is associated with mucocutaneous manifestations and autoantibody production in systemic lupus erythematosus

OBJECTIVE

To determine whether genetic substructure in European-derived populations is associated with specific manifestations of systemic lupus erythematosus (SLE), including mucocutaneous phenotypes, autoantibody production, and renal disease.

METHODS

SLE patients of European descent (n=1,754) from 8 case collections were genotyped for >1,400 ancestry informative markers that define a north-south gradient of European substructure. Using the Structure program, each SLE patient was characterized in terms of percent Northern (versus percent Southern) European ancestry based on these genetic markers. Nonparametric methods, including tests for trend, were used to identify associations between Northern European ancestry and specific SLE manifestations.

RESULTS

In multivariate analyses, increasing levels of Northern European ancestry were significantly associated with photosensitivity (Ptrend=0.0021, odds ratio for highest quartile of Northern European ancestry versus lowest quartile [ORhigh-low] 1.64, 95% confidence interval [95% CI] 1.13-2.35) and discoid rash (Ptrend=0.014, ORhigh-low 1.93, 95% CI 0.98-3.83). In contrast, increasing levels of Northern European ancestry had a protective effect against the production of anticardiolipin autoantibodies (Ptrend=1.6x10(-4), ORhigh-low 0.46, 95% CI 0.30-0.69) and anti-double-stranded DNA autoantibodies (Ptrend=0.017, ORhigh-low 0.67, 95% CI 0.46-0.96).

CONCLUSION

This study demonstrates that specific SLE manifestations vary according to Northern versus Southern European ancestry. Thus, genetic ancestry may contribute to the clinical heterogeneity and variation in disease outcomes among SLE patients of European descent. Moreover, these results suggest that genetic studies of SLE subphenotypes will need to carefully address issues of population substructure based on genetic ancestry.

Clinical and genetic risk factors for pneumonia in systemic lupus erythematosus

OBJECTIVE

To define the contribution of polymorphisms in genes encoding tumor necrosis factor (TNF), mannose-binding lectin (MBL), and Fcgamma receptor IIa (FCGR2A) as well as clinical factors, to the development of pneumonia in patients with systemic lupus erythematosus (SLE).

METHODS

We studied 282 SLE patients from a multiethnic cohort. Pneumonia events and clinical risk factors for pneumonia were identified through medical record review. Genotyping was performed for MBL (+223, +230, and +239), TNF (-308, -238, and +488), and FCGR2A (-131H/R) polymorphisms. Univariate analyses were performed to identify clinical and genetic risk factors for pneumonia. Covariates for multivariate analysis included sex, ethnicity, treatment with immunomodulators, and leukopenia.

RESULTS

Forty-two patients (15%) had at least 1 episode of pneumonia. Polymorphism of the TNF gene, particularly the -238A allele and a related haplotype, revealed the most striking and consistent association with pneumonia in univariate analyses. Results of multivariate analyses indicated an odds ratio (OR) for the TNF -238A allele of 3.5 (P = 0.007) and an OR for the related haplotype of 5.4 (P = 0.001). Male sex, treatment with immunomodulators, and leukopenia also influenced the risk of pneumonia.

CONCLUSION

These findings suggest that specific TNF variants may identify SLE patients who are at particularly high risk of developing pneumonia. Given the prevalence and excessive morbidity associated with pneumonia in SLE, these findings have clinical relevance and provide insight into the pathogenesis.