Less than 7 hours of sleep per night is associated with transitioning to systemic lupus erythematosus

Background The role of sleep in the etiology of systemic lupus erythematosus (SLE) has not been well studied. We examined whether sleep duration was associated with subsequent transitioning to SLE in individuals at risk for SLE. Methods Four hundred and thirty-six relatives of SLE patients who did not have SLE themselves at baseline were evaluated again an average of 6.3 (± 3.9) years later. Fifty-six individuals transitioned to SLE (≥ 4 cumulative American College of Rheumatology (ACR) criteria). Sleep duration, medication use and medical history were assessed by questionnaire; ACR criteria were confirmed by medical record review. Vitamin D was measured by ELISA. Generalized estimating equations, accounting for correlation within families, assessed associations between baseline sleep and the outcome of transitioning to SLE. Results Reporting sleeping less than 7 hours per night at baseline was more common in those who subsequently transitioned than those who did not transition to SLE (55% versus 32%, p = 0.0005; OR: 2.8, 95% CI 1.6-4.9). Those who transitioned to SLE were more likely to sleep less than 7 hours per night than those who did not transition to SLE adjusting for age, sex and race (OR: 2.8, 95% CI 1.6-5.1). This association remained after individual adjustment for conditions and early symptoms that could affect sleep, including prednisone use, vitamin D deficiency and number of ACR criteria (OR: 2.0, 95% CI 1.1-4.2). Conclusion Lack of sleep may be associated with transitioning to SLE, independent of early clinical manifestations of SLE that may influence sleep duration. Further evaluation of sleeping patterns and biomarkers in at-risk individuals is warranted.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. We show that transcription factors (TFs) occupy multiple loci of individual complex genetic disorders using novel computational methods. Application to 213 phenotypes and 1,544 TF binding datasets identifies 2,264 relationships between hundreds of TFs and 94 phenotypes, including AR in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of the systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many co-clustering human TFs, revealing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis, and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent upon EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new paradigms for disease origins.

Unique clinical characteristics, autoantibodies and medication use in Native American patients with systemic lupus erythematosus

Objective

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with varied morbidity and mortality. We assessed clinical presentations, autoantibody specificities and therapeutic interventions in Native American (NA) patients with SLE.

Methods

Patients with SLE meeting 1997 American College of Rheumatology classification criteria (n=3148) were enrolled between 1992 and 2010 in the multiethnic, cross-sectional Lupus Family Registry and Repository. Clinical, demographic and therapeutic information were extracted from medical records using a standardised form and formalised training. Autoantibodies were assessed by indirect immunofluorescence (antinuclear antibodies (ANA) and antidouble-stranded DNA), precipitin (ENA) and ELISA (IgG and IgM anticardiolipins).

Results

NA patients met SLE classification at a younger age (29.89±12.3 years) than European Americans (EA; 32.02±12.87, P=0.0157) and a similar age to African-Americans (AAs) and Hispanics (HIS). More NA patients had concurrent rheumatic diseases or symptoms, such as Raynaud’s phenomenon, interstitial lung disease, Sjӧgren’s syndrome and systemic sclerosis. Compared with EAs, NAs were more likely to have high-titre ANA (≥1:3240; P<0.0001) and had more SLE-associated autoantibodies. Autoantibodies with unknown specificities were more common in NAs (41%) compared with other racial/ethnic groups in this collection (AA: 24%, P=0.0006; EA: 17%, P<0.0001; HIS: 23%, P=0.0050). Fewer NA patients used hydroxychloroquine (68%) compared with others (AA: 74%, P=0.0308; EA: 79%, P=0.0001, HIS: 77%, P=0.0173); this was influenced by lower hydroxychloroquine use in NA patients from Latin America (32%). NA patients had higher rates of methotrexate use (28%) compared with AA (18%, P=0.0006) and HIS patients (14%, P=0.0003), higher azathioprine use (38%) compared with EA patients (30%, P=0.0105) and higher mycophenolate mofetil use (26%) compared with EA (17%, P=0.0012) and HIS patients (11%, P<0.0001).

Conclusions

NA patients are diagnosed with SLE earlier in life and present worse concurrent rheumatic disease symptoms than EA patients. NA patients also are more likely to have expanded autoantibody profiles and precipitins of unknown specificities.

Characterization and classification of lupus patients based on plasma thermograms

Objective

Plasma thermograms (thermal stability profiles of blood plasma) are being utilized as a new diagnostic approach for clinical assessment. In this study, we investigated the ability of plasma thermograms to classify systemic lupus erythematosus (SLE) patients versus non SLE controls using a sample of 300 SLE and 300 control subjects from the Lupus Family Registry and Repository. Additionally, we evaluated the heterogeneity of thermograms along age, sex, ethnicity, concurrent health conditions and SLE diagnostic criteria.

Methods

Thermograms were visualized graphically for important differences between covariates and summarized using various measures. A modified linear discriminant analysis was used to segregate SLE versus control subjects on the basis of the thermograms. Classification accuracy was measured based on multiple training/test splits of the data and compared to classification based on SLE serological markers.

Results

Median sensitivity, specificity, and overall accuracy based on classification using plasma thermograms was 86%, 83%, and 84% compared to 78%, 95%, and 86% based on a combination of five antibody tests. Combining thermogram and serology information together improved sensitivity from 78% to 86% and overall accuracy from 86% to 89% relative to serology alone. Predictive accuracy of thermograms for distinguishing SLE and osteoarthritis / rheumatoid arthritis patients was comparable. Both gender and anemia significantly interacted with disease status for plasma thermograms (p<0.001), with greater separation between SLE and control thermograms for females relative to males and for patients with anemia relative to patients without anemia.

Conclusion

Plasma thermograms constitute an additional biomarker which may help improve diagnosis of SLE patients, particularly when coupled with standard diagnostic testing. Differences in thermograms according to patient sex, ethnicity, clinical and environmental factors are important considerations for application of thermograms in a clinical setting.

Clinical and Serologic Features in Patients With Incomplete Lupus Classification Versus Systemic Lupus Erythematosus Patients and Controls

OBJECTIVE

Incomplete lupus erythematosus (ILE) involves clinical and/or serologic manifestations consistent with but insufficient for systemic lupus erythematosus (SLE) classification. Because the nature of ILE is poorly understood and no treatment recommendations exist, we examined the clinical manifestations, medication history, and immunologic features in a diverse collection of ILE and SLE patients.

METHODS

Medical records of subjects enrolled in the Lupus Family Registry and Repository were reviewed for medication history and American College of Rheumatology (ACR) classification criteria to identify ILE patients (3 ACR criteria; n = 440) and SLE patients (≥4 ACR criteria; n = 3,397). Participants completed the Connective Tissue Disease Screening Questionnaire. Anticardiolipin and plasma B lymphocyte stimulator (BLyS) were measured by enzyme-linked immunosorbent assay, antinuclear antibodies (ANAs) by indirect immunofluorescence, and 13 autoantibodies by bead-based assays.

RESULTS

On average, ILE patients were older than SLE patients (46.2 years versus 42.0 years; P < 0.0001), and fewer ILE patients were African American (23.9% versus 32.2%; P < 0.001). ILE patients exhibited fewer autoantibody specificities than SLE patients (1.3 versus 2.6; P < 0.0001) and were less likely to have ANA titers ≥1:1,080 (10.5% versus 19.5%; P < 0.0001). BLyS levels were intermediate in ILE patients (controls < ILE; P = 0.016; ILE < SLE; P = 0.008). Pericarditis, renal, or neurologic manifestations occurred in 12.5% of ILE patients and were associated with non-European American race/ethnicity (P = 0.012). Hydroxychloroquine use increased over time, but was less frequent in ILE than SLE patients (65.2% versus 83.1%; P < 0.0001).

CONCLUSION

Although usually characterized by milder symptoms, ILE manifestations may require immunomodulatory treatments. Longitudinal studies are necessary to understand how ILE affects organ damage and future SLE risk, and to delineate molecular pathways unique to ILE.

Rare X Chromosome Abnormalities in Systemic Lupus Erythematosus and Sjögren's Syndrome

OBJECTIVE

Sjögren's syndrome (SS) and systemic lupus erythematosus (SLE) are related by clinical and serologic manifestations as well as genetic risks. Both diseases are more commonly found in women than in men, at a ratio of ~10 to 1. Common X chromosome aneuploidies, 47,XXY and 47,XXX, are enriched among men and women, respectively, in either disease, suggesting a dose effect on the X chromosome.

METHODS

We examined cohorts of SS and SLE patients by constructing intensity plots of X chromosome single-nucleotide polymorphism alleles, along with determining the karyotype of selected patients.

RESULTS

Among ~2,500 women with SLE, we found 3 patients with a triple mosaic, consisting of 45,X/46,XX/47,XXX. Among ~2,100 women with SS, 1 patient had 45,X/46,XX/47,XXX, with a triplication of the distal p arm of the X chromosome in the 47,XXX cells. Neither the triple mosaic nor the partial triplication was found among the controls. In another SS cohort, we found a mother/daughter pair with partial triplication of this same region of the X chromosome. The triple mosaic occurs in ~1 in 25,000-50,000 live female births, while partial triplications are even rarer.

CONCLUSION

Very rare X chromosome abnormalities are present among patients with either SS or SLE and may inform the location of a gene(s) that mediates an X dose effect, as well as critical cell types in which such an effect is operative.

Identification of Four Novel Loci in Asthma in European American and African American Populations

RATIONALE

Despite significant advances in knowledge of the genetic architecture of asthma, specific contributors to the variability in the burden between populations remain uncovered.

OBJECTIVES

To identify additional genetic susceptibility factors of asthma in European American and African American populations.

METHODS

A phenotyping algorithm mining electronic medical records was developed and validated to recruit cases with asthma and control subjects from the Electronic Medical Records and Genomics network. Genome-wide association analyses were performed in pediatric and adult asthma cases and control subjects with European American and African American ancestry followed by metaanalysis. Nominally significant results were reanalyzed conditioning on allergy status.

MEASUREMENTS AND MAIN RESULTS

The validation of the algorithm yielded an average of 95.8% positive predictive values for both cases and control subjects. The algorithm accrued 21,644 subjects (65.83% European American and 34.17% African American). We identified four novel population-specific associations with asthma after metaanalyses: loci 6p21.31, 9p21.2, and 10q21.3 in the European American population, and the PTGES gene in African Americans. TEK at 9p21.2, which encodes TIE2, has been shown to be involved in remodeling the airway wall in asthma, and the association remained significant after conditioning by allergy. PTGES, which encodes the prostaglandin E synthase, has also been linked to asthma, where deficient prostaglandin E₂ synthesis has been associated with airway remodeling.

CONCLUSIONS

This study adds to understanding of the genetic architecture of asthma in European Americans and African Americans and reinforces the need to study populations of diverse ethnic backgrounds to identify shared and unique genetic predictors of asthma.

Strong viral associations with SLE among Filipinos

Objectives

Epstein-Barr virus (EBV) is considered an important environmental factor in SLE aetiology, but the relationship between SLE and EBV in the Filipino population is unknown. We tested associations between SLE, lupus-associated autoantibodies and seropositivity for EBV and other herpes viruses in the Filipino population.

Methods

Sera from Filipino patients with SLE (n=233), unaffected first-degree relatives (FDRs, n=543) and unrelated controls (n=221) were tested for antibodies against EBV, cytomegalovirus (CMV) and herpes simplex viruses (HSV-1 and HSV-2) by standardised ELISAs. Humoral specificity against EBV nuclear antigen (EBNA)-1 was compared by solid-phase epitope mapping. Autoantibodies were detected by a bead-based multiplex assay. Results were analysed by Fisher's exact test, Student's t-test, χ² test and one-way analysis of variance, as appropriate for the question.

Results

Filipino patients with SLE had increased seroprevalence and elevated antibody concentrations against EBV viral capsid antigen (EBV-VCA), CMV, HSV-1 and HSV-2 compared with unrelated controls (p<0.05). Seropositivity for anti-EBV early antigen (EA), a marker of EBV reactivation, was dramatically increased in patients with SLE compared with unrelated controls (92.3% vs 40.4%; OR 17.15(95% CI 10.10, 30.66), p<0.0001) or unaffected FDRs (49.4%; OR 12.04(7.42, 20.74), p<0.0001), despite similar seroprevalence of EBV-VCA in patients and FDRs. In patients with SLE, EBV-EA seropositivity correlated with lupus-associated autoantibodies (p<0.001), most notably with autoantibodies against dsDNA, chromatin, Sm, SmRNP and RNP A (p<0.01). Patient and unrelated control sera reacted to the highly repetitive glycine and alanine domain of EBNA-1. An epitope spanning EBNA-1_410-420 was identified in sera of patients with SLE and showed limited binding by FDR and control sera.

Conclusions

Filipino patients with SLE have elevated prevalence and concentrations of antibodies against EBV, CMV, HSV-1 and HSV-2 antigens, along with altered anti-EBNA-1 specificities. EBV reactivation is more common among Filipino patients with SLE compared with healthy Filipinos and may contribute to SLE pathogenesis in this population.

Monocyte MicroRNA Expression in Active Systemic Juvenile Idiopathic Arthritis Implicates MicroRNA-125a-5p in Polarized Monocyte Phenotypes

OBJECTIVE

Systemic juvenile idiopathic arthritis (JIA) is an inflammatory disease of childhood in which cells of the monomyelocytoid lineage are thought to be key effector cells. Monocytes from patients with systemic JIA have a distinct phenotype, with features of both M1 and M2 alternative activation. MicroRNAs are critical regulators of monocyte polarization and function, but cellular microRNAs in systemic JIA have not been examined systematically.

METHODS

MicroRNA TaqMan arrays were used to determine the expression profiles of monocytes from children with systemic JIA. Expression of microRNA-125a-5p (miR-125a-5p) and its contribution to monocyte polarization were examined using in vitro-polarized THP-1 cells and primary human monocytes.

RESULTS

A total of 110 microRNAs were found to be differentially expressed in monocytes from patients with active systemic JIA, including molecules implicated in rheumatoid arthritis pathogenesis, cytokine production, and monocyte polarization. MicroRNA-125a-5p was identified as being highly up-regulated in monocytes from children with active systemic JIA, as compared to those from children with clinically inactive JIA or those with active polyarticular JIA, and correlated with systemic features of the disease. In vitro, monocyte miR125a-5p expression was increased after polarization under M2b or M2c conditions. Inhibition of miR-125a-5p showed that this microRNA contributed to full polarization of M2b regulatory macrophages. In contrast, miR-125a-5p overexpression enhanced M2b polarization and altered other polarized populations, including increasing the production of M2 markers. Indeed, in vitro overexpression of this microRNA altered the macrophage phenotype toward that observed in systemic JIA.

CONCLUSION

Children with active systemic JIA have profound alterations in the expression of microRNAs that are implicated in monocyte function and polarization. One of these microRNAs, miR-125a-5p, is also a regulator of immunoregulatory M2b macrophages.

X Chromosome Dose and Sex Bias in Autoimmune Diseases: Increased Prevalence of 47,XXX in Systemic Lupus Erythematosus and Sjögren's Syndrome

OBJECTIVE

More than 80% of autoimmune disease predominantly affects females, but the mechanism for this female bias is poorly understood. We suspected that an X chromosome dose effect accounts for this, and we undertook this study to test our hypothesis that trisomy X (47,XXX; occurring in ∼1 in 1,000 live female births) would be increased in patients with female-predominant diseases (systemic lupus erythematosus [SLE], primary Sjögren's syndrome [SS], primary biliary cirrhosis, and rheumatoid arthritis [RA]) compared to patients with diseases without female predominance (sarcoidosis) and compared to controls.

METHODS

All subjects in this study were female. We identified subjects with 47,XXX using aggregate data from single-nucleotide polymorphism arrays, and, when possible, we confirmed the presence of 47,XXX using fluorescence in situ hybridization or quantitative polymerase chain reaction.

RESULTS

We found 47,XXX in 7 of 2,826 SLE patients and in 3 of 1,033 SS patients, but in only 2 of 7,074 controls (odds ratio in the SLE and primary SS groups 8.78 [95% confidence interval 1.67-86.79], P = 0.003 and odds ratio 10.29 [95% confidence interval 1.18-123.47], P = 0.02, respectively). One in 404 women with SLE and 1 in 344 women with SS had 47,XXX. There was an excess of 47,XXX among SLE and SS patients.

CONCLUSION

The estimated prevalence of SLE and SS in women with 47,XXX was ∼2.5 and ∼2.9 times higher, respectively, than that in women with 46,XX and ∼25 and ∼41 times higher, respectively, than that in men with 46,XY. No statistically significant increase of 47,XXX was observed in other female-biased diseases (primary biliary cirrhosis or RA), supporting the idea of multiple pathways to sex bias in autoimmunity.

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Use of SLICC criteria in a large, diverse lupus registry enables SLE classification of a subset of ACR-designated subjects with incomplete lupus

Objective

SLE is traditionally classified using the American College of Rheumatology (ACR) criteria. The Systemic Lupus International Collaborating Clinics (SLICC) recently validated an alternative system. This study examined large cohorts of subjects with SLE and incomplete lupus erythematosus (ILE) to compare the impact of ACR and SLICC criteria.

Methods

Medical records of subjects in the Lupus Family Registry and Repository were reviewed for documentation of 1997 ACR classification criteria, SLICC classification criteria and medication usage. Autoantibodies were assessed by indirect immunofluorescence (ANA, antidouble-stranded DNA), precipitin (Sm) and ELISA (anticardiolipin). Other relevant autoantibodies were detected by precipitin and with a bead-based multiplex assay.

Results

Of 3575 subjects classified with SLE under at least one system, 3312 (92.6%) were classified as SLE by both systems (SLE^both), 85 only by ACR criteria (SLE^ACR-only) and 178 only by SLICC criteria (SLE^SLICC-only). Of 440 subjects meeting 3 ACR criteria, 33.9% (149/440) were SLE^SLICC-only, while 66.1% (n=291, designated ILE) did not meet the SLICC classification criteria. Under the SLICC system, the complement criterion and the individual autoantibody criteria enabled SLE classification of SLE^SLICC-only subjects, while SLE^ACR-only subjects failed to meet SLICC classification due to the combined acute/subacute cutaneous criterion. The SLICC criteria classified more African-American subjects by the leucopenia/lymphopenia criterion than did ACR criteria. Compared with SLE^ACR-only subjects, SLE^SLICC-only subjects exhibited similar numbers of affected organ systems, rates of major organ system involvement (∼30%: pulmonary, cardiovascular, renal, neurological) and medication history.

Conclusions

The SLICC criteria classify more subjects with SLE than ACR criteria; however, individuals with incomplete lupus still exist under SLICC criteria. Subjects who gain SLE classification through SLICC criteria exhibit heterogeneous disease, including potential major organ involvement. These results provide supportive evidence that SLICC criteria may be more inclusive of SLE subjects for clinical studies.

Discerning Risk of Disease Transition in Relatives of Systemic Lupus Erythematosus Patients Utilizing Soluble Mediators and Clinical Features

OBJECTIVE

Systemic lupus erythematosus (SLE) and other autoimmune diseases cause significant morbidity. Identifying populations at risk of developing SLE is essential for curtailing irreversible inflammatory damage. The aim of this study was to identify factors associated with transition to classified disease that would inform our understanding of the risk of SLE.

METHODS

Previously identified blood relatives of patients with SLE, who had <4 American College of Rheumatology (ACR) classification criteria for SLE at baseline, were enrolled in this follow-up study (n = 409 unaffected relatives). Participants provided detailed family, demographic, and clinical information, including the SLE-specific portion of the Connective Tissue Disease Screening Questionnaire (SLE-CSQ). Serum and plasma samples were tested for the presence of lupus-associated autoantibodies and 52 soluble mediators. Generalized estimating equations (GEEs) were applied to identify factors predictive of disease transition.

RESULTS

Of the 409 unaffected relatives of SLE patients, 45 (11%) had transitioned to classified SLE at follow-up (mean time to follow-up 6.4 years). Relatives who transitioned to SLE displayed more lupus-associated autoantibody specificities and higher SLE-CSQ scores (P < 0.0001) at baseline than did relatives who did not transition. Importantly, those who had developed SLE during the follow-up period also had elevated baseline plasma levels of inflammatory mediators, including B lymphocyte stimulator, stem cell factor (SCF), and interferon-associated chemokines (P ≤ 0.02), with concurrent decreases in the levels of regulatory mediators, transforming growth factor β (TGFβ), and interleukin-10 (P ≤ 0.03). GEE analyses revealed that baseline SLE-CSQ scores or ACR scores (number of ACR criteria satisfied) and plasma levels of SCF and TGFβ, but not autoantibodies, were significant and independent predictors of SLE transition (P ≤ 0.03).

CONCLUSION

Preclinical alterations in levels of soluble mediators may predict transition to classified disease in relatives of SLE patients. Thus, immune perturbations precede SLE classification and can help identify high-risk relatives for rheumatology referral and potential enrollment in prevention trials.

Altered type II interferon precedes autoantibody accrual and elevated type I interferon activity prior to systemic lupus erythematosus classification

OBJECTIVES

The relationship of immune dysregulation and autoantibody production that may contribute to systemic lupus erythematosus (SLE) pathogenesis is unknown. This study evaluates the individual and combined contributions of autoantibodies, type I interferon (IFN-α) activity, and IFN-associated soluble mediators to disease development leading to SLE.

METHODS

Serial serum specimens from 55 individuals collected prior to SLE classification (average timespan=4.3 years) and unaffected healthy controls matched by age (±5 years), gender, race and time of sample procurement were obtained from the Department of Defense Serum Repository. Levels of serum IFN-α activity, IFN-associated mediators and autoantibodies were evaluated and temporal relationships assessed by growth curve modelling, path analysis, analysis of covariance and random forest models.

RESULTS

In cases, but not matched controls, autoantibody specificities and IFN-associated mediators accumulated over a period of years, plateauing near the time of disease classification (p<0.001). Autoantibody positivity coincided with or followed type II IFN dysregulation, preceding IFN-α activity in growth curve models, with elevated IFN-α activity and B-lymphocyte stimulator levels occurring shortly before SLE classification (p≤0.005). Cases were distinguished by multivariate random forest models incorporating IFN-γ, macrophage chemoattractant protein (MCP)-3, anti-chromatin and anti-spliceosome antibodies (accuracy 93% >4 years pre-classification; 97% within 2 years of SLE classification).

CONCLUSIONS

Years before SLE classification, enhancement of the type II IFN pathway allows for accumulation of autoantibodies and subsequent elevations in IFN-α activity immediately preceding SLE classification. Perturbations in select immunological processes may help identify at-risk individuals for further clinical evaluation or participation in prospective intervention trials.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage.

METHODS

We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated.

RESULTS

The IRF5-TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control-adjusted P [Pgcadj ] = 2.61 × 10(-29), OR 2.12 [95% CI 1.88-2.39]), followed by HLA class II on the DQA2-DQB1 loci (rs9275572: Pgcadj  = 1.11 × 10(-16), OR 1.62 [95% CI 1.46-1.80] and rs9271366: Pgcadj  = 6.46 × 10(-12), OR 2.06 [95% CI 1.71-2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj  = 1.39 × 10(-8)) with an expression quantitative trait locus (eQTL) effect (Peqtl  = 8.0 × 10(-37) at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE.

CONCLUSION

Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

In Vivo Therapeutic Success of MicroRNA-155 Antagomir in a Mouse Model of Lupus Alveolar Hemorrhage

OBJECTIVE

Diffuse alveolar hemorrhage (DAH) is a rare but life-threatening complication of systemic lupus erythematosus (SLE). Pristane-treated B6 mice develop severe DAH within 2 weeks of treatment. MicroRNA-155 (miR-155) is a pleiotropic microRNA that plays a crucial role in the regulation of immune responses. Recent studies have revealed a pathogenic role of miR-155 in various autoimmune disorders. The purpose of this study was to examine the role of miR-155 in the development of DAH in pristane-induced lupus using miR-155-knockout (miR-155(-/-)) mice and miR-155 antagomir to silence miR-155.

METHODS

DAH was induced by an intraperitoneal injection of 0.5 ml of pristane. MicroRNA-155 antagomir was administered intravenously to silence miR-155 expression. Lung tissues were collected for RNA extraction and were embedded in paraffin for sectioning. Gene expression profiling data were analyzed using Ingenuity Pathway Analysis. Real-time quantitative polymerase chain reaction analysis was used for single-gene validation. Luciferase reporter assay and argonaute 2 immunoprecipitation were performed for target validation.

RESULTS

MicroRNA-155 expression was significantly increased during the development of DAH. Disease progression was reduced in miR-155(-/-) mice as well as by in vivo silencing of miR-155 using a miR-155 antagomir. MicroRNA-155 silencing dampened pristane-induced ectopic activation of multiple inflammatory pathways and reduced the expression of proinflammatory cytokines. Several negative regulators of NF-κB signaling were inhibited by pristane and were reactivated in miR-155(-/-) mice. In particular, the antiinflammatory factor peroxisome proliferator-activated receptor α was identified as a direct target of miR-155.

CONCLUSION

MicroRNA-155 promotes pristane-induced lung inflammation. It contributes to ectopic activation of NF-κB signaling pathways by targeting multiple negative regulators. MicroRNA-155 antagomir may be a promising therapeutic strategy for treating acute lung inflammation in lupus.

The curiously suspicious: a role for Epstein-Barr virus in lupus

While the events initiating the development of autoantibodies in systemic lupus erythematosus (SLE) have not yet been convincingly established, newly developed tools for molecular investigation make such an undertaking increasingly practical. Applied to the earliest events in the sequence culminating in lupus autoimmunity, we present a critical potential role for Epstein-Barr virus (EBV) in the development and perhaps perpetuation of SLE. The expected properties for an environmental risk factor for SLE are found in this virus and the human host response against it. Existing data show the molecular progression to autoimmunity observed in SLE patient sera, the discovery of the first autoimmune epitopes in the Sm and Ro autoantigen systems, and the possible emergence of these autoantibodies from the heterologous antibodies against Epstein-Barr nuclear antigen-1 (EBNA-1). Further, existing data demonstrate association of SLE with EBV infection, even preceding the development of autoimmunity. Finally, the data are consistent with a proposed model of lupus pathogenesis that begins with antibodies to EBNA-1, predisposing to immune responses that develop crossreactive autoantibodies that culminate in the development of SLE autoimmunity.

Dysregulation of innate and adaptive serum mediators precedes systemic lupus erythematosus classification and improves prognostic accuracy of autoantibodies

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a poorly understood preclinical stage of immune dysregulation and symptom accrual. Accumulation of antinuclear autoantibody (ANA) specificities is a hallmark of impending clinical disease. Yet, many ANA-positive individuals remain healthy, suggesting that additional immune dysregulation underlies SLE pathogenesis. Indeed, we have recently demonstrated that interferon (IFN) pathways are dysregulated in preclinical SLE. To determine if other forms of immune dysregulation contribute to preclinical SLE pathogenesis, we measured SLE-associated autoantibodies and soluble mediators in samples from 84 individuals collected prior to SLE classification (average timespan = 5.98 years), compared to unaffected, healthy control samples matched by race, gender, age (±5 years), and time of sample procurement. We found that multiple soluble mediators, including interleukin (IL)-5, IL-6, and IFN-γ, were significantly elevated in cases compared to controls more than 3.5 years pre-classification, prior to or concurrent with autoantibody positivity. Additional mediators, including innate cytokines, IFN-associated chemokines, and soluble tumor necrosis factor (TNF) superfamily mediators increased longitudinally in cases approaching SLE classification, but not in controls. In particular, levels of B lymphocyte stimulator (BLyS) and a proliferation-inducing ligand (APRIL) were comparable in cases and controls until less than 10 months pre-classification. Over the entire pre-classification period, random forest models incorporating ANA and anti-Ro/SSA positivity with levels of IL-5, IL-6, and the IFN-γ-induced chemokine, MIG, distinguished future SLE patients with 92% (±1.8%) accuracy, compared to 78% accuracy utilizing ANA positivity alone. These data suggest that immune dysregulation involving multiple pathways contributes to SLE pathogenesis. Importantly, distinct immunological profiles are predictive for individuals who will develop clinical SLE and may be useful for delineating early pathogenesis, discovering therapeutic targets, and designing prevention trials.

Combined role of vitamin D status and CYP24A1 in the transition to systemic lupus erythematosus

OBJECTIVE

We examined whether measures of vitamin D were associated with transitioning to systemic lupus erythematosus (SLE) in individuals at risk for SLE.

METHODS

436 individuals who reported having a relative with SLE but who did not have SLE themselves were evaluated at baseline and again an average of 6.3 (±3.9) years later. Fifty-six individuals transitioned to SLE (≥4 cumulative American College of Rheumatology criteria). 25-Hydroxyvitamin D (25[OH]D) levels were measured by ELISA. Six single-nucleotide polymorphisms in four vitamin D genes were genotyped. Generalised estimating equations, adjusting for correlation within families, were used to test associations between the vitamin D variables and the outcome of transitioning to SLE.

RESULTS

Mean baseline 25[OH]D levels (p=0.42) and vitamin D supplementation (p=0.65) were not different between those who did and did not transition to SLE. Vitamin D deficiency (25[OH]D <20 ng/mL) was greater in those who transitioned compared with those who did not transition to SLE (46% vs 33%, p=0.05). The association between 25[OH]D and SLE was modified by CYP24A1 rs4809959, where for each additional minor allele increased 25[OH]D was associated with decreased SLE risk: zero minor alleles (adjusted OR: 1.03, CI 0.98 to 1.09), one minor allele (adjusted OR: 1.01, CI 0.97 to 1.05) and two minor alleles (adjusted OR: 0.91, CI 0.84 to 0.98). Similarly, vitamin D deficiency significantly increased the risk of transitioning to SLE in those with two minor alleles at rs4809959 (adjusted OR: 4.90, CI 1.33 to 18.04).

CONCLUSIONS

Vitamin D status and CYP24A1 may have a combined role in the transition to SLE in individuals at increased genetic risk for SLE.

Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network

Genetic variation can affect drug response in multiple ways, although it remains unclear how rare genetic variants affect drug response. The electronic Medical Records and Genomics (eMERGE) Network, collaborating with the Pharmacogenomics Research Network, began eMERGE‐PGx, a targeted sequencing study to assess genetic variation in 82 pharmacogenes critical for implementation of “precision medicine.” The February 2015 eMERGE‐PGx data release includes sequence‐derived data from ∼5,000 clinical subjects. We present the variant frequency spectrum categorized by variant type, ancestry, and predicted function. We found 95.12% of genes have variants with a scaled Combined Annotation‐Dependent Depletion score above 20, and 96.19% of all samples had one or more Clinical Pharmacogenetics Implementation Consortium Level A actionable variants. These data highlight the distribution and scope of genetic variation in relevant pharmacogenes, identifying challenges associated with implementing clinical sequencing for drug treatment at a broader level, underscoring the importance for multifaceted research in the execution of precision medicine.

Dysregulation of T helper-type cytokines appears prior to and concurrent with appearance of autoantibody specificities during early SLE pathogenesis and contributes to prognosis of transition to classified disease

Autoantibodies and immune dysregulation are associated with systemic lupus erythematosus (SLE), but how these factors influence disease development is unknown. This study evaluates timing and changes in SLE-linked autoantibodies and soluble mediator pathways in the years preceding SLE classification. Serial sera from 84 SLE cases spanning pre- and post-SLE classification (average time = 5.98 years) and matched healthy controls (HC) were obtained from the Department of Defense Serum Repository. Adjusting for multiple comparison, a number of soluble mediators, including IL-5 (q=4.35×10−6), IL-6 (q=8.26×10−6), and IFN-γ (q=0.037), were significantly elevated in cases vs. HC &gt;3.5 years pre-classification, prior to (IL-5 and IL-6) and concurrent with (IFN-γ) the earliest SLE-specific autoantibody specificity, anti-Ro/SSA. Th1-type, Th17-type, and TNF superfamily soluble mediators increased longitudinally in cases approaching SLE classification, but not in HC (q&lt;0.05). In particular, serum levels of BLyS and APRIL were comparable in cases and HC until &lt;10 months pre-classification (q=0.003 and q=0.019, respectively). Random forest models incorporating IL-5, IL-6, and IFN-γ levels (79–82% accuracy) distinguished future SLE patients better than models with ANA alone (58% accuracy) &gt;3.5 years before classification. These data indicate that pre-clinical immune system perturbations allow for the accumulation of autoantibodies preceding clinical symptoms and SLE classification.

Heightened preclinical dysregulation of distinct adaptive and renal-associated mediators in patients who develop nephritis as they transition to systemic lupus erythematosus classification

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease marked by immune dysregulation. Why some patients have only moderate symptoms and others develop organ-threatening manifestations is unclear. This study evaluates the temporal expression of autoantibodies and cytokines in sera from the Department of Defense Serum Repository during transition from preclinical lupus to SLE in patients (n=83) who do or do not present with nephritis. Patients who met renal criteria (n=30; 36%) did so within 5.2 (±5.5) months of SLE classification (range −5.2 to +12 months). Renal cases experienced earlier onset of autoantibody positivity vs. non-renal cases (mean −3.8 vs −2.6 years relative to SLE classification, p=0.0442). Prior to developing nephritis, renal cases exhibited elevated levels of distinct soluble mediators compared to matched non-renal cases and healthy controls, including adaptive mediators IL-4, IL-5, IL-12, and IFN-γ, chemokines IP-10 and MCP-3, and nephritis-associated mediators SCF and shed TNFRII (all p&lt;0.05), increasing again at the time of nephritis (all p&lt;0.02 compared to pre-nephritis levels). These same mediators were significantly increased at SLE classification in patients with nephritis vs. non-renal patients (all p&lt;0.0001). These data indicate that perturbations in distinct immune mediated inflammatory processes may help identify individuals at high risk of renal involvement for early and continued monitoring and intervention.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population.

METHODS

A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China.

RESULTS

Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified.

CONCLUSION

Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

Blood relatives with incomplete lupus are distinguished from unaffected relatives and relatives with classified disease by altered autoantibodies, soluble inflammatory and regulatory mediators, and clinical features

Blood relatives (Rel) of lupus patients have increased risk of systemic lupus erythematosus (SLE). Some have autoantibodies and clinical features, but do not meet ≥ 4 ACR criteria needed to reach SLE classification (incomplete lupus, ILE). This study examined individuals enrolled in the Lupus Family Registry and Repository who only met 3 ACR classification criteria during medical record review. ILE patients (n=77) were matched to unaffected Rel and unrelated controls (Ctls), with a subset of ILE patients (n=55) matched to medical record-confirmed SLE patients. ILE patients had significantly higher SLE-specific screening questionnaire (CSQ) scores than unaffected Rel and Ctls (p≤0.0001). A number of soluble mediators that positively correlated with CSQ scores were highest in ILE patients, including SCF (p=0.0001), BLyS (p=0.0018), MCP-3 (p=0.0167), and TNFRI (p=0.0196), as well as ANA titer (p&lt;0.01) and number of lupus-associated autoantibodies (p&lt;0.01). Comparing ILE vs. SLE patients, SLE patients had higher rates of arthritis, serositis, and renal disease (p≤0.004), number of autoantibody specificities (p&lt;0.004), and levels of BLyS (p=0.0138), IL-2Rα (p=0.0201), IP-10 (p=0.0269), and TNFRII (p=0.0309). Yet, ILE patients had higher levels of the regulatory mediator TGF-β (p=0.0454). Identification of factors which discern relatives at increased risk of transitioning to classified SLE may be beneficial to curtail inflammatory damage.

X chromosome gene Ddx3x is involved in B-cell development

Ddx3x, located on the X chromosome, belongs to the DEAD box RNA helicase family. Ddx3x is previously implicated in RNA metabolism, cell cycle control, apoptosis, tumorigeneses and viral infections.

In the process of exploring the role of ddx3x in the female predominance of systemic lupus erythematosus, we found abnormalities in B-cell development within conditional ddx3x-deficient mice.

Breeding of ddx3x-floxed mice with Vav1-Cre mice led to deletion of ddx3x in hematopoietic cells (Vav1ddx3x). Ddx3x expression was reduced 80% in the bone marrow (BM) of Vav1ddx3x mice. Hematopoietic deficiency of ddx3x is associated with a reduced proportion and absolute number of B cells, starting at the small pre-B cell stage, and then subsequently immature B cells and mature B cells in bone marrow, relative to wild-type mice. In spleen, the frequency and number of follicular B cells (FOB) were significantly decreased whereas the numbers of marginal zone, germinal center and plasma-blast B cells were sustained. Mixed Vav1ddx3x and wild-type bone marrow chimera mice demonstrated a B-cell intrinsic effect. Surprisingly, despite fewer B cells in the BM and spleen, serum IgM, IgA and IgG titers were elevated in Vav1ddx3x mice, suggesting ddx3x may affect B-cell function as well.

In conclusion, we showed that ddx3x deficiency affects B-cell development resulting in a loss of B cells in bone marrow and the periphery and an increase in serum immunoglobulin levels.

Klinefelter's syndrome (47,XXY) is in excess among men with Sjögren's syndrome

Primary Sjögren's syndrome (pSS) has a strong female bias. We evaluated an X chromosome dose effect by analyzing 47,XXY (Klinefelter's syndrome, 1 in 500 live male births) among subjects with pSS. 47,XXY was determined by examination of fluorescence intensity of single nucleotide polymorphisms from the X and Y chromosomes. Among 136 pSS men there were 4 with 47,XXY. This was significantly different from healthy controls (1 of 1254 had 47,XXY, p=0.0012 by Fisher's exact test) as well men with rheumatoid arthritis (0 of 363 with 47,XXY), but not different compared to men with systemic lupus erythematosus (SLE) (4 of 136 versus 8 of 306, Fisher's exact test p=NS). These results are consistent with the hypothesis that the number of X chromosomes is critical for the female bias of pSS, a property that may be shared with SLE but not RA.

Nanovesicle-targeted Kv1.3 knockdown in memory T cells suppresses CD40L expression and memory phenotype

Ca(2+) signaling controls activation and effector functions of T lymphocytes. Ca(2+) levels also regulate NFAT activation and CD40 ligand (CD40L) expression in T cells. CD40L in activated memory T cells binds to its cognate receptor, CD40, on other cell types resulting in the production of antibodies and pro-inflammatory mediators. The CD40L/CD40 interaction is implicated in the pathogenesis of autoimmune disorders and CD40L is widely recognized as a therapeutic target. Ca(2+) signaling in T cells is regulated by Kv1.3 channels. We have developed lipid nanoparticles that deliver Kv1.3 siRNAs (Kv1.3-NPs) selectively to CD45RO(+) memory T cells and reduce the activation-induced Ca(2+) influx. Herein we report that Kv1.3-NPs reduced NFAT activation and CD40L expression exclusively in CD45RO(+) T cells. Furthermore, Kv1.3-NPs suppressed cytokine release and induced a phenotype switch of T cells from predominantly memory to naïve. These findings indicate that Kv1.3-NPs operate as targeted immune suppressive agents with promising therapeutic potentials.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production

OBJECTIVES

Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants.

METHODS

We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA.

RESULTS

We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10^-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10^-3 and 6.8×10^-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10^-5 and 2.0×10^-4, respectively).

CONCLUSION

We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Focused transcription from the human CR2/CD21 core promoter is regulated by synergistic activity of TATA and Initiator elements in mature B cells

Complement receptor 2 (CR2/CD21) is predominantly expressed on the surface of mature B cells where it forms part of a coreceptor complex that functions, in part, to modulate B-cell receptor signal strength. CR2/CD21 expression is tightly regulated throughout B-cell development such that CR2/CD21 cannot be detected on pre-B or terminally differentiated plasma cells. CR2/CD21 expression is upregulated at B-cell maturation and can be induced by IL-4 and CD40 signaling pathways. We have previously characterized elements in the proximal promoter and first intron of CR2/CD21 that are involved in regulating basal and tissue-specific expression. We now extend these analyses to the CR2/CD21 core promoter. We show that in mature B cells, CR2/CD21 transcription proceeds from a focused TSS regulated by a non-consensus TATA box, an initiator element and a downstream promoter element. Furthermore, occupancy of the general transcriptional machinery in pre-B versus mature B-cell lines correlate with CR2/CD21 expression level and indicate that promoter accessibility must switch from inactive to active during the transitional B-cell window.

ePhenotyping for Abdominal Aortic Aneurysm in the Electronic Medical Records and Genomics (eMERGE) Network: Algorithm Development and Konstanz Information Miner Workflow

BACKGROUND AND OBJECTIVE

We designed an algorithm to identify abdominal aortic aneurysm cases and controls from electronic health records to be shared and executed within the "electronic Medical Records and Genomics" (eMERGE) Network.

MATERIALS AND METHODS

Structured Query Language, was used to script the algorithm utilizing "Current Procedural Terminology" and "International Classification of Diseases" codes, with demographic and encounter data to classify individuals as case, control, or excluded. The algorithm was validated using blinded manual chart review at three eMERGE Network sites and one non-eMERGE Network site. Validation comprised evaluation of an equal number of predicted cases and controls selected at random from the algorithm predictions. After validation at the three eMERGE Network sites, the remaining eMERGE Network sites performed verification only. Finally, the algorithm was implemented as a workflow in the Konstanz Information Miner, which represented the logic graphically while retaining intermediate data for inspection at each node. The algorithm was configured to be independent of specific access to data and was exportable (without data) to other sites.

RESULTS

The algorithm demonstrated positive predictive values (PPV) of 92.8% (CI: 86.8-96.7) and 100% (CI: 97.0-100) for cases and controls, respectively. It performed well also outside the eMERGE Network. Implementation of the transportable executable algorithm as a Konstanz Information Miner workflow required much less effort than implementation from pseudo code, and ensured that the logic was as intended.

DISCUSSION AND CONCLUSION

This ePhenotyping algorithm identifies abdominal aortic aneurysm cases and controls from the electronic health record with high case and control PPV necessary for research purposes, can be disseminated easily, and applied to high-throughput genetic and other studies.

A GWAS Study on Liver Function Test Using eMERGE Network Participants

Introduction

Liver enzyme levels and total serum bilirubin are under genetic control and in recent years genome-wide population-based association studies have identified different susceptibility loci for these traits. We conducted a genome-wide association study in European ancestry participants from the Electronic Medical Records and Genomics (eMERGE) Network dataset of patient medical records with available genotyping data in order to identify genetic contributors to variability in serum bilirubin levels and other liver function tests and to compare the effects between adult and pediatric populations.

Methods

The process of whole genome imputation of eMERGE samples with standard quality control measures have been described previously. After removing missing data and outliers based on principal components (PC) analyses, 3294 samples from European ancestry were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and total serum bilirubin and other liver function tests was tested using linear regression, adjusting for age, gender, site, platform and ancestry principal components (PC).

Results

Consistent with previous results, a strong association signal has been detected for UGT1A gene cluster (best SNP rs887829, beta = 0.15, p = 1.30x10^-118) for total serum bilirubin level. Indeed, in this region more than 176 SNPs (or indels) had p<10⁻⁸ spanning 150Kb on the long arm of chromosome 2q37.1. In addition, we found a similar level of magnitude in a pediatric group (p = 8.26x10^-47, beta = 0.17). Further imputation using sequencing data as a reference panel revealed association of other markers including known TA7 repeat indels (rs8175347) (p = 9.78x10^-117) and rs111741722 (p = 5.41x10^-119) which were in proxy (r2 = 0.99) with rs887829. Among rare variants, two Asian subjects homozygous for coding SNP rs4148323 (G71R) were identified. Additional known effects for total serum bilirubin were also confirmed including organic anion transporters SLCO1B1-SLCO1B3, TDRP and ZMYND8 at FDR<0.05 with no gene-gene interaction effects. Phenome-wide association studies (PheWAS) suggest a protective effect of TA7 repeat against cerebrovascular disease in an adult cohort (OR = 0.75, p = 0.0008). Among other liver function tests, we also confirmed the previous effect of the ABO blood group locus for variation in serum alkaline phosphatase (rs579459, p = 9.44x10^-15).

Conclusions

Taken together, our data present interesting findings with strong confirmation of previous effects by simply using the eMERGE electronic health record phenotyping. In addition, our findings indicate that similar to the adult population, the UGT1A1 is the main locus responsible for normal variation of serum bilirubin in pediatric populations.

Genetic associations of leptin-related polymorphisms with systemic lupus erythematosus

Leptin is abnormally elevated in the plasma of patients with systemic lupus erythematosus (SLE), where it is thought to promote and/or sustain pro-inflammatory responses. Whether this association could reflect an increased genetic susceptibility to develop SLE is not known, and studies of genetic associations with leptin-related polymorphisms in SLE patients have been so far inconclusive. Here we genotyped DNA samples from 15,706 SLE patients and healthy matched controls from four different ancestral groups, to correlate polymorphisms of genes of the leptin pathway to risk for SLE. It was found that although several SNPs showed weak associations, those associations did not remain significant after correction for multiple testing. These data do not support associations between defined leptin-related polymorphisms and increased susceptibility to develop SLE.

Autoantibodies targeting glomerular annexin A2 identify patients with proliferative lupus nephritis

PURPOSE

Patients with systemic lupus erythematosus (SLE) frequently develop lupus nephritis (LN), a complication frequently leading to end stage kidney disease. Immune complex deposition in the glomerulus is central to the development of LN. Using a targeted proteomic approach, we tested the hypothesis that autoantibodies targeting glomerular antigens contribute to the development of LN.

EXPERIMENTAL DESIGN

Human podocyte and glomerular proteins were separated by SDS-PAGE and immunoblotted with sera from SLE patients with and without LN. The regions of those gels corresponding to reactive bands observed with sera from LN patients were analyzed using LC-MS/MS.

RESULTS

LN reactive bands were seen at approximately 50 kDa in podocyte extracts and between 36 and 50 kDa in glomerular extracts. Those bands were analyzed by LC-MS/MS and 102 overlapping proteins were identified. Bioinformatic analysis determined that 36 of those proteins were membrane associated, including a protein previously suggested to contribute to glomerulonephritis and LN, annexin A2. By ELISA, patients with proliferative LN demonstrated significantly increased antibodies against annexin A2.

CONCLUSION AND CLINICAL RELEVANCE

Proteomic approaches identified multiple candidate antigens for autoantibodies in patients with LN. Serum antibodies against annexin A2 were significantly elevated in subjects with proliferative LN, validating those antibodies as potential biomarkers.

Biology-Driven Gene-Gene Interaction Analysis of Age-Related Cataract in the eMERGE Network

Bioinformatics approaches to examine gene-gene models provide a means to discover interactions between multiple genes that underlie complex disease. Extensive computational demands and adjusting for multiple testing make uncovering genetic interactions a challenge. Here, we address these issues using our knowledge-driven filtering method, Biofilter, to identify putative single nucleotide polymorphism (SNP) interaction models for cataract susceptibility, thereby reducing the number of models for analysis. Models were evaluated in 3,377 European Americans (1,185 controls, 2,192 cases) from the Marshfield Clinic, a study site of the Electronic Medical Records and Genomics (eMERGE) Network, using logistic regression. All statistically significant models from the Marshfield Clinic were then evaluated in an independent dataset of 4,311 individuals (742 controls, 3,569 cases), using independent samples from additional study sites in the eMERGE Network: Mayo Clinic, Group Health/University of Washington, Vanderbilt University Medical Center, and Geisinger Health System. Eighty-three SNP-SNP models replicated in the independent dataset at likelihood ratio test P < 0.05. Among the most significant replicating models was rs12597188 (intron of CDH1)-rs11564445 (intron of CTNNB1). These genes are known to be involved in processes that include: cell-to-cell adhesion signaling, cell-cell junction organization, and cell-cell communication. Further Biofilter analysis of all replicating models revealed a number of common functions among the genes harboring the 83 replicating SNP-SNP models, which included signal transduction and PI3K-Akt signaling pathway. These findings demonstrate the utility of Biofilter as a biology-driven method, applicable for any genome-wide association study dataset.

MiR-125a targets effector programs to stabilize Treg-mediated immune homeostasis

Although different autoimmune diseases show discrete clinical features, there are common molecular pathways intimately involved. Here we show that miR-125a is downregulated in peripheral CD4(+) T cells of human autoimmune diseases including systemic lupus erythematosus and Crohn's disease, and relevant autoimmune mouse models. miR-125a stabilizes both the commitment and immunoregulatory capacity of Treg cells. In miR-125a-deficient mice, the balance appears to shift from immune suppression to inflammation, and results in more severe pathogenesis of colitis and experimental autoimmune encephalomyelitis (EAE). The genome-wide target analysis reveals that miR-125a suppresses several effector T-cell factors including Stat3, Ifng and Il13. Using a chemically synthesized miR-125a analogue, we show potential to re-programme the immune homeostasis in EAE models. These findings point to miR-125a as a critical factor that controls autoimmune diseases by stabilizing Treg-mediated immune homeostasis.

Genome-Wide Association Study Identifies Novel Loci Associated With Diisocyanate-Induced Occupational Asthma

Diisocyanates, reactive chemicals used to produce polyurethane products, are the most common causes of occupational asthma. The aim of this study is to identify susceptibility gene variants that could contribute to the pathogenesis of diisocyanate asthma (DA) using a Genome-Wide Association Study (GWAS) approach. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed in 74 diisocyanate-exposed workers with DA and 824 healthy controls using Omni-2.5 and Omni-5 SNP microarrays. We identified 11 SNPs that exceeded genome-wide significance; the strongest association was for the rs12913832 SNP located on chromosome 15, which has been mapped to the HERC2 gene (p = 6.94 × 10(-14)). Strong associations were also found for SNPs near the ODZ3 and CDH17 genes on chromosomes 4 and 8 (rs908084, p = 8.59 × 10(-9) and rs2514805, p = 1.22 × 10(-8), respectively). We also prioritized 38 SNPs with suggestive genome-wide significance (p < 1 × 10(-6)). Among them, 17 SNPs map to the PITPNC1, ACMSD, ZBTB16, ODZ3, and CDH17 gene loci. Functional genomics data indicate that 2 of the suggestive SNPs (rs2446823 and rs2446824) are located within putative binding sites for the CCAAT/Enhancer Binding Protein (CEBP) and Hepatocyte Nuclear Factor 4, Alpha transcription factors (TFs), respectively. This study identified SNPs mapping to the HERC2, CDH17, and ODZ3 genes as potential susceptibility loci for DA. Pathway analysis indicated that these genes are associated with antigen processing and presentation, and other immune pathways. Overlap of 2 suggestive SNPs with likely TF binding sites suggests possible roles in disruption of gene regulation. These results provide new insights into the genetic architecture of DA and serve as a basis for future functional and mechanistic studies.

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 (JIA), resembles familial hemophagocytic lymphohistiocytosis (HLH), a constellation of autosomal-recessive immune disorders resulting from deficiency in cytolytic pathway proteins. We undertook this study to test our hypothesis that MAS predisposition in systemic JIA could be attributed to rare gene sequence variants affecting the cytotolytic pathway.

METHODS

Whole-exome sequencing was used in 14 patients with systemic JIA and MAS and in their parents to identify protein-altering single-nucleotide polymorphisms/indels in known HLH-associated genes. To discover new candidate genes, the entire whole-exome sequencing 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 patients with systemic JIA and MAS (35.7%). This was in contrast to only 4 variants in 4 of 29 patients with systemic JIA without MAS (13.8%). Homozygosity and compound heterozygosity analysis applied to the entire whole-exome sequencing data in systemic JIA/MAS revealed 3 recessive pairs in 3 genes and compound heterozygotes in 73 genes. We also identified 20 heterozygous rare protein-altering variants that occurred in at least 2 patients. Many of the identified genes encoded 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.10 × 10(-5) ).

CONCLUSION

Whole-exome sequencing performed in patients with systemic JIA and MAS identified rare protein-altering variants in known HLH-associated genes as well as in new candidate genes.

Association of systemic lupus erythematosus with uranium exposure in a community living near a uranium-processing plant: a nested case-control study

OBJECTIVE

To explore the hypothesis that cases of systemic lupus erythematosus (SLE) would be found more frequently in community members with high prior uranium exposure in the Fernald Community Cohort (FCC).

METHODS

A nested case-control study was performed using data from the FCC, a volunteer population of individuals who had resided near a uranium ore-processing plant in Fernald, Ohio during the years of plant operation; uranium plant workers were excluded. Members of the FCC were monitored for 18 years. SLE cases were identified using the American College of Rheumatology 1997 revised classification criteria, laboratory testing, and medical record review. Each case was matched to 4 controls by age, race, and sex. Sera from potential cases and controls were screened for autoantibodies. Cumulative exposure to uranium particulates was calculated using a dosimetry model. Logistic regression with covariates was used to calculate the odds ratios (ORs) with 95% confidence intervals (95% CIs) for the probability of an association between uranium exposure and SLE.

RESULTS

The FCC comprised 4,187 individuals with minimal levels of uranium exposure, 1,273 with moderate exposure, and 2,756 with high exposure. The diagnosis of SLE was confirmed in 23 of 31 individuals who had been assigned International Classification of Diseases, Ninth Revision codes for lupus, and was also confirmed in 2 of 43 individuals who had been prescribed hydroxychloroquine. The female to male ratio was 5.25:1. Of the 25 confirmed SLE cases, 12 were in the high exposure group. The presence of SLE was associated with higher levels of uranium exposure (OR 3.92, 95% CI 1.13-13.59; P = 0.031).

CONCLUSION

High uranium exposure is associated with SLE, as compared to matched controls, in this sample of uranium-exposed individuals. Potential explanations for this relationship include possible autoimmune or estrogen effects of uranium, somatic mutation, epigenetic effects, or effects of some other unidentified accompanying exposure.

Lupus risk variants in the PXK locus alter B-cell receptor internalization

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⁻¹⁰, 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.

Genetic association of CD247 (CD3ζ) with SLE in a large-scale multiethnic study

A classic T-cell phenotype in systemic lupus erythematosus (SLE) is the downregulation and replacement of the CD3ζ chain that alters T-cell receptor signaling. However, genetic associations with SLE in the human CD247 locus that encodes CD3ζ are not well established and require replication in independent cohorts. Our aim was therefore to examine, localize and validate CD247-SLE association in a large multiethnic population. We typed 44 contiguous CD247 single-nucleotide polymorphisms (SNPs) in 8922 SLE patients and 8077 controls from four ethnically distinct populations. The strongest associations were found in the Asian population (11 SNPs in intron 1, 4.99 × 10(-4) < P < 4.15 × 10(-2)), where we further identified a five-marker haplotype (rs12141731-rs2949655-rs16859085-rs12144621-rs858554; G-G-A-G-A; P(hap) = 2.12 × 10(-5)) that exceeded the most associated single SNP rs858554 (minor allele frequency in controls = 13%; P = 4.99 × 10(-4), odds ratio = 1.32) in significance. Imputation and subsequent association analysis showed evidence of association (P < 0.05) at 27 additional SNPs within intron 1. Cross-ethnic meta-analysis, assuming an additive genetic model adjusted for population proportions, showed five SNPs with significant P-values (1.40 × 10(-3) < P< 3.97 × 10(-2)), with one (rs704848) remaining significant after Bonferroni correction (P(meta) = 2.66 × 10(-2)). Our study independently confirms and extends the association of SLE with CD247, which is shared by various autoimmune disorders and supports a common T-cell-mediated mechanism.

Phenome-wide association study (PheWAS) in EMR-linked pediatric cohorts, genetically links PLCL1 to speech language development and IL5-IL13 to Eosinophilic Esophagitis

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⁻⁷, 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⁻⁸, 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⁻⁹, 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⁻⁵, 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.

Twin and family studies reveal strong environmental and weaker genetic cues explaining heritability of eosinophilic esophagitis

BACKGROUND

Eosinophilic esophagitis (EoE) is a chronic antigen-driven allergic inflammatory disease, likely involving the interplay of genetic and environmental factors, yet their respective contributions to heritability are unknown.

OBJECTIVE

To quantify the risk associated with genes and environment on familial clustering of EoE.

METHODS

Family history was obtained from a hospital-based cohort of 914 EoE probands (n = 2192 first-degree "Nuclear-Family" relatives) and an international registry of monozygotic and dizygotic twins/triplets (n = 63 EoE "Twins" probands). Frequencies, recurrence risk ratios (RRRs), heritability, and twin concordance were estimated. Environmental exposures were preliminarily examined.

RESULTS

Analysis of the Nuclear-Family-based cohort revealed that the rate of EoE, in first-degree relatives of a proband, was 1.8% (unadjusted) and 2.3% (sex-adjusted). RRRs ranged from 10 to 64, depending on the family relationship, and were higher in brothers (64.0; P = .04), fathers (42.9; P = .004), and males (50.7; P < .001) than in sisters, mothers, and females, respectively. The risk of EoE for other siblings was 2.4%. In the Nuclear-Family cohort, combined gene and common environment heritability was 72.0% ± 2.7% (P < .001). In the Twins cohort, genetic heritability was 14.5% ± 4.0% (P < .001), and common family environment contributed 81.0% ± 4% (P < .001) to phenotypic variance. Probandwise concordance in monozygotic co-twins was 57.9% ± 9.5% compared with 36.4% ± 9.3% in dizygotic co-twins (P = .11). Greater birth weight difference between twins (P = .01), breast-feeding (P = .15), and fall birth season (P = .02) were associated with twin discordance in disease status.

CONCLUSIONS

EoE RRRs are increased 10- to 64-fold compared with the general population. EoE in relatives is 1.8% to 2.4%, depending on relationship and sex. Nuclear-Family heritability appeared to be high (72.0%). However, the Twins cohort analysis revealed a powerful role for common environment (81.0%) compared with additive genetic heritability (14.5%).

Genetic analysis of the pathogenic molecular sub-phenotype interferon-alpha identifies multiple novel loci involved in systemic lupus erythematosus

Systemic Lupus Erythematosus (SLE) is a chronic autoimmune disorder characterized by inflammation of multiple organ systems and dysregulated interferon responses. SLE is both genetically and phenotypically heterogeneous, greatly reducing the power of case-control studies in SLE. Elevated circulating interferon alpha (IFN-α) is a stable, heritable trait in SLE, which has been implicated in primary disease pathogenesis. 40–50% of patients have high IFN-α, and high levels correspond with clinical differences. To study genetic heterogeneity in SLE, we performed a case-case study comparing patients with high vs. low IFN-α in over 1550 SLE cases, including GWAS and replication cohorts. In meta-analysis, the top associations in European ancestry were PRKG1 rs7897633 (P_Meta=2.75 × 10⁻⁸) and PNP rs1049564 (P_Meta=1.24 × 10⁻⁷). We also found evidence for cross-ancestral background associations with the ANKRD44 and PLEKHF2 loci. These loci have not been previously identified in case-control SLE genetic studies. Bioinformatic analyses implicated these loci functionally in dendritic cells and natural killer cells, both of which are involved in IFN-α production in SLE. As case-control studies of heterogeneous diseases reach a limit of feasibility with respect to subject number and detectable effect size, the study of informative pathogenic subphenotypes becomes an attractive strategy for genetic discovery in complex disease.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22 981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10⁻⁸). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.