Sex-specific effects of injury and beta-adrenergic activation on metabolic and inflammatory mediators in a murine model of post-traumatic osteoarthritis

OBJECTIVE

Metabolic processes are intricately linked to the resolution of innate inflammation and tissue repair, two critical steps for treating post-traumatic osteoarthritis (PTOA). Based on lipolytic and immunoregulatory actions of norepinephrine, we hypothesized that intra-articular β-adrenergic receptor (βAR) stimulation would suppress PTOA-associated inflammation in the infrapatellar fat pad (IFP) and synovium.

DESIGN

We used the βAR agonist isoproterenol to perturb intra-articular metabolism 3.5 weeks after applying a non-invasive single-load compression injury to knees of 12-week-old male and female mice. We examined the acute effects of intra-articular isoproterenol treatment relative to saline on IFP histology, multiplex gene expression of synovium-IFP tissue, synovial fluid metabolomics, and mechanical allodynia.

RESULTS

Injured knees developed PTOA pathology characterized by heterotopic ossification, articular cartilage loss, and IFP atrophy and fibrosis. Isoproterenol suppressed the upregulation of pro-fibrotic genes and downregulated the expression of adipose genes and pro-inflammatory genes (Adam17, Cd14, Icam1, Csf1r, and Casp1) in injured joints of female (but not male) mice. Analysis of published single-cell RNA-seq data identified elevated catecholamine-associated gene expression in resident-like synovial-IFP macrophages after injury. Injury substantially altered synovial fluid metabolites by increasing amino acids, peptides, sphingolipids, phospholipids, bile acids, and dicarboxylic acids, but these changes were not appreciably altered by isoproterenol. Intra-articular injection of either isoproterenol or saline increased mechanical allodynia in female mice, whereas neither substance affected male mice.

CONCLUSIONS

Acute βAR activation altered synovial-IFP transcription in a sex and injury-dependent manner, suggesting that women with PTOA may be more sensitive than men to treatments targeting sympathetic neural signaling pathways.

Corrigendum: Sex differences in the genetics of sarcoidosis across European and African ancestry populations

[This corrects the article DOI: 10.3389/fmed.2023.1132799.].

Examination of eQTL Polymorphisms Associated with Increased Risk of Progressive Complicated Sarcoidosis in European and African Descent Subjects

Background

A limited pool of SNPs are linked to the development and severity of sarcoidosis, a systemic granulomatous inflammatory disease. By integrating genome-wide association studies (GWAS) data and expression quantitative trait loci (eQTL) single nuclear polymorphisms (SNPs), we aimed to identify novel sarcoidosis SNPs potentially influencing the development of complicated sarcoidosis.

Methods

A GWAS (Affymetrix 6.0) involving 209 African-American (AA) and 193 European-American (EA, 75 and 51 complicated cases respectively) and publicly-available GWAS controls (GAIN) was utilized. Annotation of multi-tissue eQTL SNPs present on the GWAS created a pool of ~46,000 eQTL SNPs examined for association with sarcoidosis risk and severity (Logistic Model, Plink). The most significant EA/AA eQTL SNPs were genotyped in a sarcoidosis validation cohort (n=1034) and cross-validated in two independent GWAS cohorts.

Results

No single GWAS SNP achieved significance (p<1x10-8), however, analysis of the eQTL/GWAS SNP pool yielded 621 eQTL SNPs (p<10-4) associated with 730 genes that highlighted innate immunity, MHC Class II, and allograft rejection pathways with multiple SNPs validated in an independent sarcoidosis cohort (105 SNPs analyzed) (NOTCH4, IL27RA, BTNL2, ANXA11, HLA-DRB1). These studies confirm significant association of eQTL/GWAS SNPs in EAs and AAs with sarcoidosis risk and severity (complicated sarcoidosis) involving HLA region and innate immunity.

Conclusion

Despite the challenge of deciphering the genetic basis for sarcoidosis risk/severity, these results suggest that integrated eQTL/GWAS approaches may identify novel variants/genes and support the contribution of dysregulated innate immune responses to sarcoidosis severity.

Identification of Environmental Exposures Associated with Risk of Sarcoidosis in African Americans

Rationale: Sarcoidosis is a racially disparate granulomatous disease likely caused by environmental exposures, genes, and their interactions. Despite increased risk in African Americans, few environmental risk factor studies in this susceptible population exist. Objectives: To identify environmental exposures associated with the risk of sarcoidosis in African Americans and those that differ in effect by self-identified race and genetic ancestry. Methods: The study sample comprised 2,096 African Americans (1,205 with and 891 without sarcoidosis) compiled from three component studies. Unsupervised clustering and multiple correspondence analyses were used to identify underlying clusters of environmental exposures. Mixed-effects logistic regression was used to evaluate the association of these exposure clusters and the 51 single-component exposures with risk of sarcoidosis. A comparison case-control sample of 762 European Americans (388 with and 374 without sarcoidosis) was used to assess heterogeneity in exposure risk by race. Results: Seven exposure clusters were identified, five of which were associated with risk. The exposure cluster with the strongest risk association was composed of metals (P < 0.001), and within this cluster, exposure to aluminum had the highest risk (odds ratio, 3.30; 95% confidence interval [95% CI], 2.23-4.09; P < 0.001). This effect also differed by race (P < 0.001), with European Americans having no significant association with exposure (odds ratio, 0.86; 95% CI, 0.56-1.33). Within African Americans, the increased risk was dependent on genetic African ancestry (P = 0.047). Conclusions: Our findings support African Americans having sarcoidosis environmental exposure risk profiles that differ from those of European Americans. These differences may underlie racially disparate incidence rates that are partially explained by genetic variation differing by African ancestry.

Multiple Correspondence Analysis and HLA-Associations of Organ Involvement in a Large Cohort of African-American and European-American Patients with Sarcoidosis

Sarcoidosis is a systemic granulomatous disease with predominant pulmonary involvement and vast heterogeneity of clinical manifestations and disease outcomes. African American (AA) patients suffer greater morbidity and mortality. Using Multiple Correspondence Analysis, we identified seven clusters of organ involvement in European American (EA; n = 385) patients which were similar to those previously described in a Pan-European (GenPhenReSa) and a Spanish cohort (SARCOGEAS). In contrast, AA (n = 987) had six, less well-defined and overlapping clusters with little similarity to the cluster identified in the EA cohort evaluated at the same U.S. institutions. Association of cluster membership with two-digit HLA-DRB1 alleles demonstrated ancestry-specific patterns of association and replicated known HLA effects.These results further support the notion that genetically influenced immune risk profiles, which differ based on ancestry, play a role in phenotypic heterogeneity. Dissecting such risk profiles will move us closer to personalized medicine for this complex disease.

Sex differences in the genetics of sarcoidosis across European and African ancestry populations

Background

Sex differences in the susceptibility of sarcoidosis are unknown. The study aims to identify sex-dependent genetic variations in two clinical sarcoidosis phenotypes: Löfgren's syndrome (LS) and non-Löfgren's syndrome (non-LS).

Methods

A meta-analysis of genome-wide association studies was conducted on Europeans and African Americans, totaling 10,103 individuals from three population-based cohorts, Sweden (n = 3,843), Germany (n = 3,342), and the United States (n = 2,918), followed by an SNP lookup in the UK Biobank (UKB, n = 387,945). A genome-wide association study based on Immunochip data consisting of 141,000 single nucleotide polymorphisms (SNPs) was conducted in the sex groups. The association test was based on logistic regression using the additive model in LS and non-LS sex groups independently. Additionally, gene-based analysis, gene expression, expression quantitative trait loci (eQTL) mapping, and pathway analysis were performed to discover functionally relevant mechanisms related to sarcoidosis and biological sex.

Results

We identified sex-dependent genetic variations in LS and non-LS sex groups. Genetic findings in LS sex groups were explicitly located in the extended Major Histocompatibility Complex (xMHC). In non-LS, genetic differences in the sex groups were primarily located in the MHC class II subregion and ANXA11. Gene-based analysis and eQTL enrichment revealed distinct sex-specific gene expression patterns in various tissues and immune cell types. In LS sex groups, a pathway map related to antigen presentation machinery by IFN-gamma. In non-LS, pathway maps related to immune response lectin-induced complement pathway in males and related to maturation and migration of dendritic cells in skin sensitization in females were identified.

Conclusion

Our findings provide new evidence for a sex bias underlying sarcoidosis genetic architecture, particularly in clinical phenotypes LS and non-LS. Biological sex likely plays a role in disease mechanisms in sarcoidosis.

Low Gut Microbial Diversity Augments Estrogen-Driven Pulmonary Fibrosis in Female-Predominant Interstitial Lung Disease

Although profibrotic cytokines, such as IL-17A and TGF-β1, have been implicated in the pathogenesis of interstitial lung disease (ILD), the interactions between gut dysbiosis, gonadotrophic hormones and molecular mediators of profibrotic cytokine expression, such as the phosphorylation of STAT3, have not been defined. Here, through chromatin immunoprecipitation sequencing (ChIP-seq) analysis of primary human CD4+ T cells, we show that regions within the STAT3 locus are significantly enriched for binding by the transcription factor estrogen receptor alpha (ERa). Using the murine model of bleomycin-induced pulmonary fibrosis, we found significantly increased regulatory T cells compared to Th17 cells in the female lung. The genetic absence of ESR1 or ovariectomy in mice significantly increased pSTAT3 and IL-17A expression in pulmonary CD4+ T cells, which was reduced after the repletion of female hormones. Remarkably, there was no significant reduction in lung fibrosis under either condition, suggesting that factors outside of ovarian hormones also contribute. An assessment of lung fibrosis among menstruating females in different rearing environments revealed that environments favoring gut dysbiosis augment fibrosis. Furthermore, hormone repletion following ovariectomy further augmented lung fibrosis, suggesting pathologic interactions between gonadal hormones and gut microbiota in relation to lung fibrosis severity. An analysis of female sarcoidosis patients revealed a significant reduction in pSTAT3 and IL-17A levels and a concomitant increase in TGF-β1 levels in CD4+ T cells compared to male sarcoidosis patients. These studies reveal that estrogen is profibrotic in females and that gut dysbiosis in menstruating females augments lung fibrosis severity, supporting a critical interaction between gonadal hormones and gut flora in lung fibrosis pathogenesis.

Low Gut Microbial Diversity Augments Estrogen-driven Pulmonary Fibrosis in Female-Predominant Interstitial Lung Disease

Although profibrotic cytokines such as IL-17A and TGF-β1 have been implicated in interstitial lung disease (ILD) pathogenesis, interactions between gut dysbiosis, gonadotrophic hormones and molecular mediators of profibrotic cytokine expression, such as phosphorylation of STAT3, have not been defined. Here we show by chromatin immunoprecipitation sequencing (ChIP-seq) analysis of primary human CD4+ T cells that regions within the STAT3 locus are significantly enriched for binding by the transcription factor estrogen receptor alpha (ERa). Using the murine model of bleomycin-induced pulmonary fibrosis, we found significantly increased regulatory T cells compared to Th17 cells in the female lung. Genetic absence of ESR1 or ovariectomy in mice significantly increased pSTAT3 and IL-17A expression in pulmonary CD4+ T cells, which was reduced after repletion of female hormones. Remarkably, there was no significant reduction in lung fibrosis under either condition, suggesting that factors outside of ovarian hormones also contribute. Assessment of lung fibrosis among menstruating females in different rearing environments revealed that environments favoring gut dysbiosis augment fibrosis. Furthermore, hormone repletion following ovariectomy further augmented lung fibrosis, suggesting pathologic interactions between gonadal hormones and gut microbiota on lung fibrosis severity. Analysis in female sarcoidosis patients revealed a significant reduction in pSTAT3 and IL-17A levels and a concomitant increase in TGF-β1 levels in CD4+ T cells, compared to male sarcoidosis patients. These studies reveal that estrogen is profibrotic in females and that gut dysbiosis in menstruating females augments lung fibrosis severity, supporting a critical interaction between gonadal hormones and gut flora in lung fibrosis pathogenesis.

Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses

OBJECTIVE

Sjögren's disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised.

METHODS

Whole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo-); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays.

RESULTS

Whole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene, RSAD2, in SjDRo+ (r≥0.65 or ≤-0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells.

CONCLUSION

LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Rare genetic variants explain missing heritability in smoking

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

Novel HLA associations with outcomes of Mycobacterium tuberculosis exposure and sarcoidosis in individuals of African ancestry using nearest-neighbor feature selection

Tuberculosis and sarcoidosis are inflammatory diseases characterized by granulomas that may occur in any organ but are often found in the lung. The panoply of classical human leukocyte antigen (HLA) alleles associated with occurrence and/or severity of both diseases varies considerably across studies. This heterogeneity of results, due to variation in factors like ancestry and disease subphenotype, as well as the use of simple modeling strategies to elucidate likely complex relationships, has made conclusions about underlying commonalities difficult. Here we perform HLA association analyses in individuals of African ancestry, using a greater resolution to include subphenotypes of disease and employing more comprehensive analytical techniques. Using a novel application of nearest-neighbor feature selection to score allelic importance, we investigated HLA allele association with Mycobacterium tuberculosis exposure outcomes in the first analysis of both latent Mycobacterium tuberculosis infection and active disease compared with those who, despite long-term exposure to active index cases, have neither positive diagnostic tests nor display clinical symptoms. We also compared persistent to resolved sarcoidosis. This led to the identification of novel HLA associations and evidence of main effects and interaction effects. We found strikingly similar main effects and interaction effects at HLA-DRB1, -DQB1, and -DPB1 in those resistant to tuberculosis (either latent or active) and persistent sarcoidosis.

Systemic immune response to vimentin and granuloma formation in a model of pulmonary sarcoidosis

A characteristic feature of sarcoidosis is a dysregulated immune response to persistent stimuli, often leading to the formation of non-necrotizing granulomas in various organs. Although genetic susceptibility is an essential factor in disease development, the etiology of sarcoidosis is not fully understood. Specifically, whether autoimmunity contributes to the initiation or progression of the disease is uncertain. In this study, we investigated systemic autoimmunity to vimentin in sarcoidosis. IgG antibodies to human vimentin were measured in sera from sarcoidosis patients and healthy controls. Mice immunized with recombinant murine vimentin were challenged intravenously with vimentin-coated beads to mimic pulmonary sarcoidosis. Lungs from treated mice were studied for cellular infiltration, granuloma formation, and gene expression. Immune cells in the bronchoalveolar lavage fluid were evaluated by flow cytometry. Compared to healthy controls, sarcoidosis patients had a higher frequency and levels of circulating anti-vimentin IgG. Vimentin-immunized mice developed lung granulomas following intravenous challenge with vimentin-coated beads. These sarcoidosis-like granulomas showed the presence of Langhans and foreign body multinucleated giant cells, CD4 T cells, and a heterogeneous collection of MHC II positive and arginase 1-expressing macrophages. The lungs showed upregulated pro-inflammatory gene expression, including Ifng, Il17, and Tnfa, reflecting TH1/TH17 responses typical of sarcoidosis. In addition, genes in the TH2 canonical pathway were also upregulated, congruent with increased numbers of ILC2 in the bronchoalveolar lavage. Overall, these results further validate vimentin as an autoantigen in sarcoidosis and provide evidence for an anti-vimentin immune response in disease pathogenesis. Our study also highlights the possible role of ILC2-driven TH2-like responses in the formation of lung granulomas in sarcoidosis.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries

Traditional Hardy-Weinberg equilibrium (HWE) tests (the χ2 test and the exact test) have long been used as a metric for evaluating genotype quality, as technical artifacts leading to incorrect genotype calls often can be identified as deviations from HWE. However, in data sets composed of individuals from diverse ancestries, HWE can be violated even without genotyping error, complicating the use of HWE testing to assess genotype data quality. In this manuscript, we present the Robust Unified Test for HWE (RUTH) to test for HWE while accounting for population structure and genotype uncertainty, and to evaluate the impact of population heterogeneity and genotype uncertainty on the standard HWE tests and alternative methods using simulated and real sequence data sets. Our results demonstrate that ignoring population structure or genotype uncertainty in HWE tests can inflate false-positive rates by many orders of magnitude. Our evaluations demonstrate different tradeoffs between false positives and statistical power across the methods, with RUTH consistently among the best across all evaluations. RUTH is implemented as a practical and scalable software tool to rapidly perform HWE tests across millions of markers and hundreds of thousands of individuals while supporting standard VCF/BCF formats. RUTH is publicly available at https://www.github.com/statgen/ruth.

Genome-Wide Association Study of Ocular Sarcoidosis Confirms HLA Associations and Implicates Barrier Function and Autoimmunity in African Americans

Purpose: Identify genes associated with ocular sarcoidosis (OS).Methods: We genotyped 1.1 million genetic variants to identify significant OS associations, defined as those that achieved p < 5 × 10-8 in a genome-wide comparison of OS cases to healthy controls in our European- or African-American cohorts (EA, AA). Potential functional roles of all associated variants were assessed.Results: Eight significant non-HLA variants were found in AA OS cases compared to healthy controls and confirmed as at least suggestive when comparing OS to non-OS cases. Seven of these were within MAGI1 and include transcription factor binding sites and expression quantitative trait loci. Our EA cohort, while showing similar effect sizes at variants within MAGI1, had no significant variants. Association analysis of HLA-DRB1 alleles confirmed association to OS in EA to *04:01.Conclusion: Our results support organ-specific genetic risk in OS in a compelling candidate, MAGI1, known to be associated with barrier function and autoimmunity.

Single Cell Transcriptomics Implicate Novel Monocyte and T Cell Immune Dysregulation in Sarcoidosis

Sarcoidosis is a systemic inflammatory disease characterized by infiltration of immune cells into granulomas. Previous gene expression studies using heterogeneous cell mixtures lack insight into cell-type-specific immune dysregulation. We performed the first single-cell RNA-sequencing study of sarcoidosis in peripheral immune cells in 48 patients and controls. Following unbiased clustering, differentially expressed genes were identified for 18 cell types and bioinformatically assessed for function and pathway enrichment. Our results reveal persistent activation of circulating classical monocytes with subsequent upregulation of trafficking molecules. Specifically, classical monocytes upregulated distinct markers of activation including adhesion molecules, pattern recognition receptors, and chemokine receptors, as well as enrichment of immunoregulatory pathways HMGB1, mTOR, and ephrin receptor signaling. Predictive modeling implicated TGFβ and mTOR signaling as drivers of persistent monocyte activation. Additionally, sarcoidosis T cell subsets displayed patterns of dysregulation. CD4 naïve T cells were enriched for markers of apoptosis and Th17/T_reg differentiation, while effector T cells showed enrichment of anergy-related pathways. Differentially expressed genes in regulatory T cells suggested dysfunctional p53, cell death, and TNFR2 signaling. Using more sensitive technology and more precise units of measure, we identify cell-type specific, novel inflammatory and regulatory pathways. Based on our findings, we suggest a novel model involving four convergent arms of dysregulation: persistent hyperactivation of innate and adaptive immunity via classical monocytes and CD4 naïve T cells, regulatory T cell dysfunction, and effector T cell anergy. We further our understanding of the immunopathology of sarcoidosis and point to novel therapeutic targets.

Recent advances in sarcoidosis genomics: epigenetics, gene expression, and gene by environment (G × E) interaction studies

PURPOSE OF REVIEW

We aim to review the most recent findings in genomics of sarcoidosis and highlight the gaps in the field.

RECENT FINDINGS

Original explorations of sarcoidosis subphenotypes, including cases associated with the World Trade Center and ocular sarcoidosis, have identified novel risk loci. Innovative gene--environment interaction studies utilizing modern analytical techniques have discovered risk loci associated with smoking and insecticide exposure. The application of whole-exome sequencing has identified genetic variants associated with persistent sarcoidosis and rare functional variations. A single epigenomics study has provided background knowledge of DNA methylation mechanisms in comparison with gene expression data. The application of machine-learning techniques has suggested new drug repositioning for the treatment of sarcoidosis. Several gene expression studies have identified prominent inflammatory pathways enriched in the affected tissue.

SUMMARY

Certainly, sarcoidosis research has recently advanced in the exploration of disease subphenotypes, utilizing novel analytical techniques, and including measures of clinical variation. Nevertheless, large-scale and diverse cohorts investigated with advanced sequencing methods, such as whole-genome and single-cell RNA sequencing, epigenomics, and meta-analysis coupled with cutting-edge analytic approaches, when employed, will broaden and translate genomics findings into clinical applications, and ultimately open venues for personalized medicine.

Nearest-Neighbor Projected Distance Regression for Epistasis Detection in GWAS With Population Structure Correction

Nearest-neighbor Projected-Distance Regression (NPDR) is a feature selection technique that uses nearest-neighbors in high dimensional data to detect complex multivariate effects including epistasis. NPDR uses a regression formalism that allows statistical significance testing and efficient control for multiple testing. In addition, the regression formalism provides a mechanism for NPDR to adjust for population structure, which we apply to a GWAS of systemic lupus erythematosus (SLE). We also test NPDR on benchmark simulated genetic variant data with epistatic effects, main effects, imbalanced data for case-control design and continuous outcomes. NPDR identifies potential interactions in an epistasis network that influences the SLE disorder.

American Indians Have a Higher Risk of Sjögren's Syndrome and More Disease Activity Than European Americans and African Americans

OBJECTIVE

To describe the clinical and serologic manifestations of Sjögren's syndrome (SS) in ethnic groups of the US.

METHODS

This was a cross-sectional study of 648 patients with primary SS: 20 African American (AA), 164 American Indian (AI), 426 European American (EA), and 38 patients of other races evaluated in a multidisciplinary Sjögren's International Collaborative Clinical Alliance research clinic.

RESULTS

AA subjects comprised 3.1% of the SS cohort, much lower than the percentage of AA in the state of Oklahoma (P = 3.01 × E - 05), the US (P = 2.24E - 13), or a systemic lupus erythematosus (SLE) cohort at the same institution (P = 4.26 × 10E - 27). In contrast, the percentage of AI in the SS cohort (25.3%) was much higher than expected (P = 3.17E - 09 versus SLE cohort, P = 6.36 - 26 versus Oklahoma, and P = 8.14E - 96 versus US population). The SS classification criteria were similar between AA and EA, but subjects of AI ancestry had lower rates of abnormal tear and salivary flow, as well as anti-Ro/SSA and anti-La/SSB antibodies. Paradoxically, AI had higher levels of disease activity (mean ± SD European League Against Rheumatism Sjögren's Syndrome Disease Activity Index score 3.77 ± 4.78) in comparison to EA (2.90 ± 4.12; P = 0.011) and more extraglandular manifestations affecting mainly the articular and glandular domains. Meanwhile, AA patients were characterized by higher rates of hypergammaglobulinemia (odds ratio [OR] 1.39 [95% confidence interval (95% CI) 1.39-8.65]; P = 0.01), elevated erythrocyte sedimentation rate (ESR) (OR 3.95 [95% CI 1.46-9.95]; P = 0.009), and parotid enlargement (OR 4.40 [95% CI 1.49-13.07]; P = 0.02).

CONCLUSION

AI are affected at high rates with SS but present with few classical features, potentially preventing timely diagnosis. In contrast to SLE, SS is infrequent and not more severe among AA, but the triad of hypergammaglobulinemia, increased ESR, and parotid enlargement warrants extra vigilance for lymphomagenesis.

PD-1 up-regulation on CD4+ T cells promotes pulmonary fibrosis through STAT3-mediated IL-17A and TGF-β1 production

Pulmonary fibrosis is a progressive inflammatory disease with high mortality and limited therapeutic options. Previous genetic and immunologic investigations suggest common intersections between idiopathic pulmonary fibrosis (IPF), sarcoidosis, and murine models of pulmonary fibrosis. To identify immune responses that precede collagen deposition, we conducted molecular, immunohistochemical, and flow cytometric analysis of human and murine specimens. Immunohistochemistry revealed programmed cell death-1 (PD-1) up-regulation on IPF lymphocytes. PD-1⁺CD4⁺ T cells with reduced proliferative capacity and increased transforming growth factor-β (TGF-β)/interleukin-17A (IL-17A) expression were detected in IPF, sarcoidosis, and bleomycin CD4⁺ T cells. PD-1⁺ T helper 17 cells are the predominant CD4⁺ T cell subset expressing TGF-β. Coculture of PD-1⁺CD4⁺ T cells with human lung fibroblasts induced collagen-1 production. Strikingly, ex vivo PD-1 pathway blockade resulted in reductions in TGF-β and IL-17A expression from CD4⁺ T cells, with concomitant declines in collagen-1 production from fibroblasts. Molecular analysis demonstrated PD-1 regulation of the transcription factor STAT3 (signal transducer and activator of transcription 3). Chemical blockade of STAT3, using the inhibitor STATTIC, inhibited collagen-1 production. Both bleomycin administration to PD-1 null mice or use of antibody against programmed cell death ligand 1 (PD-L1) demonstrated significantly reduced fibrosis compared to controls. This work identifies a critical, previously unrecognized role for PD-1⁺CD4⁺ T cells in pulmonary fibrosis, supporting the use of readily available therapeutics that directly address interstitial lung disease pathophysiology.

P6288Autoantibodies in heart failure associate with disease severity and differentially expressed genes in apoptotic, fibrotic, and hypoxia pathways in cardiomyocytes

Background

Previous studies suggest that autoantibodies against cardiac myosin lead to dilated cardiomyopathy (DCM). Anti-cardiac myosin antibodies cross-react with the beta adrenergic receptor (βAR) and signal cAMP-dependent protein kinase A (PKA) in cardiomyocytes leading to apoptosis, fibrosis, dilated cardiomyopathy and arrhythmias.

Purpose

To determine if cross-reactive anti-cardiac myosin/anti-βAR autoantibodies which signal cardiomyocytes through PKA might play a role to establish DCM by promoting remodeling, apoptosis, and fibrosis.

Methods

Forty-one adults with DCM were enrolled <6 months from symptom onset and followed for 12 months. Serum and myocarditis/DCM-derived human mAb were analyzed by ELISA for autoantibodies, and a PKA assay measured anti-HCM/βAR antibody-mediated signaling of cardiomyocytes (ATCC primary heart cell line H9c2). The top 50 genes differentially expressed in the cardiomyocytes treated with sera or human mAb were identified and compared to genes differentially expressed in blood of DCM patients to identify shared disease-specific genes.

Results

Anti-HCM autoantibodies including autoantibody responses against 32 overlapping synthetic peptides of the S2 fragment of HCM were significantly elevated in patients whose ejection fraction did not improve over 1-year compared to those with improved ejection fraction. The human mAb confirmed our results with HCM, βAR, specific HCM peptides, and PKA signaling. RNA sequencing revealed differentially expressed genes in serum/mAb-treated cardiomyocytes compared to genes identified after RNA sequencing of peripheral blood of patients (n=10) with DCM for >1 year from onset. A primary heart cell line (H9c2-ATCC) treated with myocarditis/DCM patient sera or human mAb revealed differentially expressed genes associated with cardiac hypertrophy and heart failure, and included inflammasome component NLRP3 and complement factor H. Ingenuity Pathway Analyses revealed 27, 7, and 1 differentially expressed genes related to apoptosis, fibrosis, and hypoxia, respectively. Gene expression of CASZ1, a transcription factor important in protection against DCM, was strongly correlated with PKA signaling (r=0.89). The KDM6B gene for lysine demethylase associated with hypoxia and apoptosis pathways and was shared between cardiomyocyte and peripheral blood analysis of DCM patients. Overall, 5 genes were shared in heart failure vs in vitro Ab-treated cardiomyocyte RNA sequencing analysis: CYP4F3, KDM6B, MBOAT7, SMAP2, and DDIT4, which affects phosphorylation of mTOR to promote autophagy and cell death, cardiac hypertrophy and dysfunction.

Conclusions

Significantly higher responses to cardiac myosin in patients with DCM were related to lack of left ventricular function improvement and to differential expression of genes promoting apoptosis, fibrosis and disease severity. These studies identify autoantibody-directed gene signaling as a potential novel therapeutic target in DCM.

Acknowledgement/Funding

National Heart, Lung, and Blood Institute, Bethesda, MD, USA

Reproducibility of Ocular Surface Staining in the Assessment of Sjögren Syndrome-Related Keratoconjunctivitis Sicca: Implications on Disease Classification

Objective

The objective of this study was to assess the performance and reproducibility of the two currently used ocular surface staining scores in the assessment of keratoconjunctivitis sicca in Sjögren syndrome (SS) research classification.

Methods

In a multidisciplinary clinic for the evaluation of sicca, we performed all tests for the American European Consensus Group (AECG) and the American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) classification criteria, including the van Bijsterveld score (vBS) and the Ocular Staining Score (OSS), in 994 participants with SS or with non-SS sicca. We analyzed the concordance between the scores, the diagnostic accuracy and correlation with clinical variables, and interrater and intrasubject reproducibility.

Results

A total of 308 (31.1%) participants had a discordant vBS and OSS that was due to extra corneal staining points in the OSS. The presence of one or more of the additional points was highly predictive of SS classification (odds ratio = 3.66; P = 1.65 × 10e-20) and was associated with abnormal results of all measures of autoimmunity and glandular dysfunction. Receiver operating characteristic curves showed optimal cutoff values of four for the vBS (sensitivity = 0.62; specificity = 0.71; Youden's J = 0.33) and five for the OSS (sensitivity = 0.56; specificity = 0.75; Youden's J = 0.31). Notably, there was very poor consistency in interobserver mean scores and distributions (P < 0.0001) and in intrasubject scores after a median of 5.5 years (35% changed status of the ocular criterion).

Conclusion

Ocular surface staining scores are useful for SS research classification; however, they are subject to significant interrater and intrasubject variability, which could result in changes in classification in 5%-10% of all subjects. These results highlight the need for objective and reproducible markers of disease that have thus far remained elusive for SS.

Microbial exposure and autoimmunity in the etiopathogenesis of Sarcoidosis

Sarcoidosis is an enigmatic inflammatory disease affecting multiple organs such as the lungs, lymph nodes, brain, heart, and skin. It is characterized by the occurrence of multiple non-caseating granulomas comprised of macrophages, giant cells, and CD4 T cells. Although genetic and environmental factors increase the risk for sarcoidosis, the etiology remains unknown. Recent studies have shown that HLA-DR3+ sarcoidosis patients harbor CD4+ T cells reactive with a vimentin peptide 424–443. However, whether autoimmunity is a cause or consequence of the disease is not known.

This study was undertaken to investigate the mechanisms responsible for activation of vimentin reactive T cells and determine their pathogenic role in sarcoidosis. For this purpose, we immunized the NOD-DR3 mice with the vimentin peptide 424–443 and induced a robust T cell response against the peptide. Fine epitope mapping using T cell hybridomas reactive with the peptide showed that both N-terminal and C-terminal amino acids were critical for T cell activation. Homology searches of microbial protein database identified a 100% homologous sequence in a protein from Klebsiella pneumoniae. Immunization of mice with this peptide, followed by intravenous transfer of peptide-coupled agarose beads induced granulomatous lesions in the lungs. Immunohistochemical analysis of these lesions showed the presence of CD4+ T cells, macrophages, and an upregulated MHC-II expression.

In summary, our data demonstrate for the first time that autoimmunity is involved in the pathogenesis of sarcoidosis. Furthermore, the demonstration of molecular mimicry between a lung infecting microbe and vimentin might explain the role of microbial exposures as a risk factor for sarcoidosis.

Extended methods for gene-environment-wide interaction scans in studies of admixed individuals with varying degrees of relationships

The etiology of many complex diseases involves both environmental exposures and inherited genetic predisposition as well as interactions between them. Gene-environment-wide interaction studies (GEWIS) provide a means to identify the interactions between genetic variation and environmental exposures that underlie disease risk. However, current GEWIS methods lack the capability to adjust for the potentially complex correlations in studies with varying degrees of relationships (both known and unknown) among individuals in admixed populations. We developed novel generalized estimating equation (GEE) based methods-GEE-adaptive and GEE-joint-to account for phenotypic correlations due to kinship while accounting for covariates, including, measures of genome-wide ancestry. In simulation studies of admixed individuals, both methods controlled family-wise error rates, an advantage over the case-only approach. They demonstrated higher power than traditional case-control methods across a wide range of underlying alternative hypotheses, especially where both marginal and interaction effects were present. We applied the proposed method to conduct a GEWIS of a known sarcoidosis risk factor (insecticide exposure) and risk of sarcoidosis in African Americans and identified two novel loci with suggestive evidence of G × E interaction.

ARID3a gene profiles are strongly associated with human interferon alpha production

Type I interferons (IFN) causes inflammatory responses to pathogens, and can be elevated in autoimmune diseases such as systemic lupus erythematosus (SLE). We previously reported unexpected associations of increased numbers of B lymphocytes expressing the DNA-binding protein ARID3a with both IFN alpha (IFNα) expression and increased disease activity in SLE. Here, we determined that IFNα producing low density neutrophils (LDNs) and plasmacytoid dendritic cells (pDCs) from SLE patients exhibit strong associations between ARID3a protein expression and IFNα production. Moreover, SLE disease activity indices correlate most strongly with percentages of ARID3a⁺ LDNs, but were also associated, less significantly, with IFNα expression in LDNs and pDCs. Hierarchical clustering and transcriptome analyses of LDNs and pDCs revealed SLE patients with low ARID3a expression cluster with healthy controls and identified gene profiles associated with increased proportions of ARID3a- and IFNα-expressing cells of each type. These data identify ARID3a as a potential transcription regulator of IFNα-related inflammatory responses and other pathways important for SLE disease activity.

Molecular mimicry and signaling by human monoclonal autoantibody derived from human myocarditis and heart failure may contribute to fibrosis and remodeling in cardiomyopathy

Previous studies in animal models of myocarditis suggest that immune responses against cardiac myosin may lead to remodeling in the heart and dilated cardiomyopathy. In human myocarditis and dilated cardiomyopathy, anti-cardiac myosin autoantibodies significantly associated with myocarditis in cohorts of adults and children and reacted significantly with both human cardiac myosin (HCM) and the β-adrenergic receptors (βAR). We produced a human monoclonal antibody (mAb) derived from human myocarditis which reacted with HCM and specific peptides from pathogenic regions of HCM. Autoantibodies against HCM and βAR from myocarditis and dilated cardiomyopathy including the human myocarditis-derived mAb targeted heart cells by signaling protein kinase A (PKA) in a primary heart cell line. To further define potential pathogenic mechanisms of the human mAb and serum autoantibodies from myocarditis and cardiomyopathy and to test our hypothesis that autoantibodies may alter gene expression related to extracellular matrix, remodeling and fibrosis, we investigated genes differentially expressed by RNA sequence analysis of a primary heart cell line after treatment with the human mAb and sera from both children and adults with myocarditis. The overall heatmap of differentially expressed genes in heart cells demonstrated a significant difference between the control sera group and both human mAb and myocarditis sera treatment groups. We identified upregulated genes related to fibrosis, apoptosis, inflammation and cardiomyopathy, suggesting that the myocarditis- and cardiomyopathy-derived autoantibodies may induce signaling in primary heart cells and potentially contribute to fibrosis and remodeling in the heart.

Genetic, Immunologic, and Environmental Basis of Sarcoidosis

Sarcoidosis is a multisystem disease with tremendous heterogeneity in disease manifestations, severity, and clinical course that varies among different ethnic and racial groups. To better understand this disease and to improve the outcomes of patients, a National Heart, Lung, and Blood Institute workshop was convened to assess the current state of knowledge, gaps, and research needs across the clinical, genetic, environmental, and immunologic arenas. We also explored to what extent the interplay of the genetic, environmental, and immunologic factors could explain the different phenotypes and outcomes of patients with sarcoidosis, including the chronic phenotypes that have the greatest healthcare burden. The potential use of current genetic, epigenetic, and immunologic tools along with study approaches that integrate environmental exposures and precise clinical phenotyping were also explored. Finally, we made expert panel-based consensus recommendations for research approaches and priorities to improve our understanding of the effect of these factors on the health outcomes in sarcoidosis.

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.

Downregulation of E Protein Activity Augments an ILC2 Differentiation Program in the Thymus

Innate lymphoid cells (ILCs) are important regulators in various immune responses. The current paradigm states that all newly made ILCs originate from common lymphoid progenitors in the bone marrow. Id2, an inhibitor of E protein transcription factors, is indispensable for ILC differentiation. Unexpectedly, we found that ectopically expressing Id1 or deleting two E protein genes in the thymus drastically increased ILC2 counts in the thymus and other organs where ILC2 normally reside. Further evidence suggests a thymic origin of these mutant ILC2s. The mutant mice exhibit augmented spontaneous infiltration of eosinophils and heightened responses to papain in the lung and increased ability to expulse the helminth parasite, Nippostrongylus brasiliensis These results prompt the questions of whether the thymus naturally has the capacity to produce ILC2s and whether E proteins restrain such a potential. The abundance of ILC2s in Id1 transgenic mice also offers a unique opportunity for testing the biological functions of ILC2s.

Association of IFIH1 and pro-inflammatory mediators: Potential new clues in SLE-associated pathogenesis

Antiviral defenses are inappropriately activated in systemic lupus erythematosus (SLE) and association between SLE and the antiviral helicase gene, IFIH1, is well established. We sought to extend the previously reported association of pathogenic soluble mediators and autoantibodies with mouse Mda5 to its human ortholog, IFIH1. To better understand the role this gene plays in human lupus, we assessed association of IFIH1 variants with soluble mediators and autoantibodies in 357 European-American SLE patients, first-degree relatives, and unrelated, unaffected healthy controls. Association between each of 135 genotyped SNPs in IFIH1 and four lupus-associated plasma mediators, IL-6, TNF-α, IFN-β, and IP-10, were investigated via linear regression. No significant associations were found to SNPs orthologous to those identified in exon 13 of the mouse. However, outside of this region there were significant associations between IL-6 and rs76162067 (p = 0.008), as well as IP-10 and rs79711023 (p = 0.003), located in a region of IFIH1 previously shown to directly influence MDA-5 mediated IP-10 and IL-6 secretion. SLE patients and FDRs carrying the minor allele for rs79711023 demonstrated lower levels of IP-10, while only FDRs carrying the minor allele for rs76162067 demonstrated an increased level of IL-6. This would suggest that the change in IP-10 is genotypically driven, while the change in IL-6 may be reflective of SLE transition status. These data suggest that IFIH1 may contribute to SLE pathogenesis via altered inflammatory mechanisms.

Single-cell analysis of glandular T cell receptors in Sjögren's syndrome

CD4⁺ T cells predominate in salivary gland (SG) inflammatory lesions in Sjögren's syndrome (SS). However, their antigen specificity, degree of clonal expansion, and relationship to clinical disease features remain unknown. We used multiplex reverse-transcriptase PCR to amplify paired T cell receptor α (TCRα) and β transcripts of single CD4⁺CD45RA^- T cells from SG and peripheral blood (PB) of 10 individuals with primary SS, 9 of whom shared the HLA DR3/DQ2 risk haplotype. TCRα and β sequences were obtained from a median of 91 SG and 107 PB cells per subject. The degree of clonal expansion and frequency of cells expressing two productively rearranged α genes were increased in SG versus PB. Expanded clones from SG exhibited complementary-determining region 3 (CDR3) sequence similarity both within and among subjects, suggesting antigenic selection and shared antigen recognition. CDR3 similarities were shared among expanded clones from individuals discordant for canonical Ro and La autoantibodies, suggesting recognition of alternative SG antigen(s). The extent of SG clonal expansion correlated with reduced saliva production and increased SG fibrosis, linking expanded SG T cells with glandular dysfunction. Knowledge of paired TCRα and β sequences enables further work toward identification of target antigens and development of novel therapies.

Development and validation of a simple lupus severity index using ACR criteria for classification of SLE

Objective

To develop a simple systemic lupus erythematosus (SLE) severity index that requires knowledge of only American College of Rheumatology (ACR) criteria and subcriteria.

Methods

This study used demographic, mortality and medical records data of 1915 patients with lupus from the Lupus Family Registry and Repository. The data were randomly split (2:1 ratio) into independent training and validation sets. A logistic regression with ridge penalty was used to model the probability of being prescribed major immunosuppressive drugs—a surrogate indicator of lupus severity. ACR criteria and subcriteria were used as predictor variables in this model, and the resulting regression coefficient estimates obtained from the training data were used as item weightings to construct the severity index.

Results

The resulting index was tested on the independent validation dataset and was found to have high predictive accuracy for immunosuppressive use and early mortality. The index was also found to be strongly correlated with a previously existing severity score for lupus. In addition, demographic factors known to influence lupus severity (eg, age of onset, gender and ethnicity) all showed robust associations with our severity index that were consistent with observed clinical trends.

Conclusions

This new index can be easily computed using ACR criteria, which may be among the most readily available data elements from patient medical records. This tool may be useful in lupus research, especially large dataset analyses to stratify patients by disease severity, an important prognostic indicator in SLE.

Association of HLA-DRB1 with Sarcoidosis Susceptibility and Progression in African Americans

HLA-DRB1 is a sarcoidosis risk gene, and the *03:01 allele is strongly associated with disease resolution in European sarcoidosis cases. Whereas the HLA-DRB1 variation is associated with sarcoidosis susceptibility in African Americans, DRB1 risk alleles are not as well defined, and associations with disease resolution have not been studied. Associations between genotyped and imputed HLA-DRB1 alleles and disease susceptibility/resolution were evaluated in a sample of 1,277 African-American patients with sarcoidosis and 1,467 control subjects. In silico binding assays were performed to assess the functional significance of the associated alleles. Increased disease susceptibility was associated with the HLA-DRB1 alleles *12:01 (odds ratio [OR], 2.11; 95% confidence interval [CI], 1.65-2.69; P = 3.2 × 10(-9)) and *11:01 (OR, 1.69; 95% CI, 1.42-2.01; P = 3.0 × 10(-9)). The strongest protective association was found with *03:01 (OR, 0.56; 95% CI, 0.44-0.73; P = 1.0 × 10(-5)). The African-derived allele *03:02 was associated with decreased risk of persistent radiographic disease (OR, 0.52; 95% CI, 0.37-0.72; P = 1.3 × 10(-4)), a finding consistent across the three component studies comprising the analytic sample. The DRB1*03:01 association with disease persistence was dependent upon local ancestry, with carriers of at least one European allele at DRB1 at a decreased risk of persistent disease (OR, 0.36; 95% CI, 0.14-0.94; P = 0.037). Results of in silico binding analyses showed that DRB1*03:01 consistently demonstrated the highest binding affinities for six bacterial peptides previously found in sarcoidosis granulomas, whereas *12:01 displayed the lowest binding affinities. This study has identified DRB1*03:01 and *03:02 as novel alleles associated with disease susceptibility and course in African Americans. Further investigation of DRB1*03 alleles may uncover immunologic factors that favor sarcoidosis protection and resolution among African Americans.

dcVar: a method for identifying common variants that modulate differential correlation structures in gene expression data

Recent studies have implicated the role of differential co-expression or correlation structure in gene expression data to help explain phenotypic differences. However, few attempts have been made to characterize the function of variants based on their role in regulating differential co-expression. Here, we describe a statistical methodology that identifies pairs of transcripts that display differential correlation structure conditioned on genotypes of variants that regulate co-expression. Additionally, we present a user-friendly, computationally efficient tool, dcVar, that can be applied to expression quantitative trait loci (eQTL) or RNA-Seq datasets to infer differential co-expression variants (dcVars). We apply dcVar to the HapMap3 eQTL dataset and demonstrate the utility of this methodology at uncovering novel function of variants of interest with examples from a height genome-wide association and cancer drug resistance. We provide evidence that differential correlation structure is a valuable intermediate molecular phenotype for further characterizing the function of variants identified in GWAS and related studies.

Identification of Immune-Relevant Factors Conferring Sarcoidosis Genetic Risk

RATIONALE

Genetic variation plays a significant role in the etiology of sarcoidosis. However, only a small fraction of its heritability has been explained so far.

OBJECTIVES

To define further genetic risk loci for sarcoidosis, we used the Immunochip for a candidate gene association study of immune-associated loci.

METHODS

Altogether the study population comprised over 19,000 individuals. In a two-stage design, 1,726 German sarcoidosis cases and 5,482 control subjects were genotyped for 128,705 single-nucleotide polymorphisms using the Illumina Immunochip for the screening step. The remaining 3,955 cases, 7,514 control subjects, and 684 parents of affected offspring were used for validation and replication of 44 candidate and two established risk single-nucleotide polymorphisms.

MEASUREMENTS AND MAIN RESULTS

Four novel susceptibility loci were identified with genome-wide significance in the European case-control populations, located on chromosomes 12q24.12 (rs653178; ATXN2/SH2B3), 5q33.3 (rs4921492; IL12B), 4q24 (rs223498; MANBA/NFKB1), and 2q33.2 (rs6748088; FAM117B). We further defined three independent association signals in the HLA region with genome-wide significance, peaking in the BTNL2 promoter region (rs5007259), at HLA-B (rs4143332/HLA-B*0801) and at HLA-DPB1 (rs9277542), and found another novel independent signal near IL23R (rs12069782) on chromosome 1p31.3.

CONCLUSIONS

Functional predictions and protein network analyses suggest a prominent role of the drug-targetable IL23/Th17 signaling pathway in the genetic etiology of sarcoidosis. Our findings reveal a substantial genetic overlap of sarcoidosis with diverse immune-mediated inflammatory disorders, which could be of relevance for the clinical application of modern therapeutics.

Fine mapping of chromosome 15q25 implicates ZNF592 in neurosarcoidosis patients

Neurosarcoidosis is a clinical subtype of sarcoidosis characterized by the presence of granulomas in the nervous system. Here, we report a highly significant association with a variant (rs75652600, P = 3.12 × 10⁻⁸, odds ratios = 4.34) within a zinc finger gene, ZNF592, from an imputation-based fine-mapping study of the chromosomal region 15q25 in African-Americans with neurosarcoidosis. We validate the association with ZNF592, a gene previously shown to cause cerebellar ataxia, in a cohort of European-Americans with neurosarcoidosis by uncovering low-frequency variants with a similar risk effect size (chr15:85309284, P = 0.0021, odds ratios = 5.36).

The IRF5-TNPO3 association with systemic lupus erythematosus has two components that other autoimmune disorders variably share

Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5-TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5-TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10(-49); OR = 1.38-1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10(-27)-10(-32), OR = 1.7-1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5-TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5-TNPO3.

Performance of HLA allele prediction methods in African Americans for class II genes HLA-DRB1, -DQB1, and -DPB1

Background

The expense of human leukocyte antigen (HLA) allele genotyping has motivated the development of imputation methods that use dense single nucleotide polymorphism (SNP) genotype data and the region’s haplotype structure, but the performance of these methods in admixed populations (such as African Americans) has not been adequately evaluated. We compared genotype-based—derived from both genome-wide genotyping and targeted sequencing—imputation results to existing allele data for HLA–DRB1, −DQB1, and –DPB1.

Results

In European Americans, the newly-developed HLA Genotype Imputation with Attribute Bagging (HIBAG) method outperformed HLA*IMP:02. In African Americans, HLA*IMP:02 performed marginally better than HIBAG pre-built models, but HIBAG models constructed using a portion of our African American sample with both SNP genotyping and four-digit HLA class II allele typing had consistently higher accuracy than HLA*IMP:02. However, HIBAG was significantly less accurate in individuals heterozygous for local ancestry (p ≤0.04). Accuracy improved in models with equal numbers of African and European chromosomes. Variants added by targeted sequencing and SNP imputation further improved both imputation accuracy and the proportion of high quality calls.

Conclusion

Combining the HIBAG approach with local ancestry and dense variant data can produce highly-accurate HLA class II allele imputation in African Americans.

Efficient generalized least squares method for mixed population and family-based samples in genome-wide association studies

Genome-wide association studies (GWAS) that draw samples from multiple studies with a mixture of relationship structures are becoming more common. Analytical methods exist for using mixed-sample data, but few methods have been proposed for the analysis of genotype-by-environment (G×E) interactions. Using GWAS data from a study of sarcoidosis susceptibility genes in related and unrelated African Americans, we explored the current analytic options for genotype association testing in studies using both unrelated and family-based designs. We propose a novel method-generalized least squares (GLX)-to estimate both SNP and G×E interaction effects for categorical environmental covariates and compared this method to generalized estimating equations (GEE), logistic regression, the Cochran-Armitage trend test, and the WQLS and MQLS methods. We used simulation to demonstrate that the GLX method reduces type I error under a variety of pedigree structures. We also demonstrate its superior power to detect SNP effects while offering computational advantages and comparable power to detect G×E interactions versus GEE. Using this method, we found two novel SNPs that demonstrate a significant genome-wide interaction with insecticide exposure-rs10499003 and rs7745248, located in the intronic and 3' UTR regions of the FUT9 gene on chromosome 6q16.1.

Inactivating mutation in the prostaglandin transporter gene, SLCO2A1, associated with familial digital clubbing, colon neoplasia, and NSAID resistance

HPGDand SLCO2A1 genes encode components of the prostaglandin catabolic pathway, with HPGD encoding the degradative enzyme 15-hydroxyprostaglandin dehydrogenase (15-PGDH), and SLCO2A1 encoding the prostaglandin transporter PGT that brings substrate to 15-PGDH. HPGD-null mice show increased prostaglandin E2 (PGE2), marked susceptibility to developing colon tumors, and resistance to colon tumor prevention by nonsteroidal anti-inflammatory drugs (NSAID). But in humans, HPGD and SLCO2A1 mutations have only been associated with familial digital clubbing. We, here, characterize a family with digital clubbing and early-onset colon neoplasia. Whole-exome sequencing identified a heterozygous nonsense mutation (G104X) in the SLCO2A1 gene segregating in 3 males with digital clubbing. Two of these males further demonstrated notably early-onset colon neoplasia, 1 with an early-onset colon cancer and another with an early-onset sessile serrated colon adenoma. Two females also carried the mutation, and both these women developed sessile serrated colon adenomas without any digital clubbing. Males with clubbing also showed marked elevations in the levels of urinary prostaglandin E2 metabolite, PGE-M, whereas, female mutation carriers were in the normal range. Furthermore, in the male proband, urinary PGE-M remained markedly elevated during NSAID treatment with either celecoxib or sulindac. Thus, in this human kindred, a null SLCO2A1 allele mimics the phenotype of the related HPGD-null mouse, with increased prostaglandin levels that cannot be normalized by NSAID therapy, plus with increased colon neoplasia. The development of early-onset colon neoplasia in male and female human SLCO2A1 mutation carriers suggests that disordered prostaglandin catabolism can mediate inherited susceptibility to colon neoplasia in man.

Comparison of the American-European Consensus Group Sjogren's syndrome classification criteria to newly proposed American College of Rheumatology criteria in a large, carefully characterised sicca cohort

OBJECTIVE

To compare the performance of the American-European Consensus Group (AECG) and the newly proposed American College of Rheumatology (ACR) classification criteria for Sjögren's Syndrome (SS) in a well-characterised sicca cohort, given ongoing efforts to resolve discrepancies and weaknesses in the systems.

METHODS

In a multidisciplinary clinic for the evaluation of sicca, we assessed features of salivary and lacrimal gland dysfunction and autoimmunity as defined by tests of both AECG and ACR criteria in 646 participants. Global gene expression profiles were compared in a subset of 180 participants.

RESULTS

Application of the AECG and ACR criteria resulted in classification of 279 and 268 participants with SS, respectively. Both criteria were met by 244 participants (81%). In 26 of the 35 AECG+/ACR participants, the minor salivary gland biopsy focal score was ≥1 (74%), while nine had positive anti-Ro/La (26%). There were 24 AECG-/ACR+ who met ACR criteria mainly due to differences in the scoring of corneal staining. All patients with SS, regardless of classification, had similar gene expression profiles, which were distinct from the healthy controls.

CONCLUSIONS

The two sets of classification criteria yield concordant results in the majority of cases and gene expression profiling suggests that patients meeting either set of criteria are more similar to other SS participants than to healthy controls. Thus, there is no clear evidence for increased value of the new ACR criteria over the old AECG criteria from the clinical or biological perspective. It is our contention, supported by this report, that improvements in diagnostic acumen will require a more fundamental understanding of the pathogenic mechanisms than is at present available.

Variants at multiple loci implicated in both innate and adaptive immune responses are associated with Sjögren's syndrome

Sjögren’s syndrome is a common autoimmune disease (~0.7% of European Americans) typically presenting as keratoconjunctivitis sicca and xerostomia. In addition to strong association within the HLA region at 6p21 (P_meta=7.65×10⁻¹¹⁴), we establish associations with IRF5-TNPO3 (P_meta=2.73×10⁻¹⁹), STAT4 (P_meta=6.80×10⁻¹⁵), IL12A (P_meta =1.17×10⁻¹⁰), FAM167A-BLK (P_meta=4.97×10⁻¹⁰), DDX6-CXCR5 (P_meta=1.10×10⁻⁸), and TNIP1 (P_meta=3.30×10⁻⁸). Suggestive associations with P_meta<5×10⁻⁵ were observed with 29 regions including TNFAIP3, PTTG1, PRDM1, DGKQ, FCGR2A, IRAK1BP1, ITSN2, and PHIP amongst others. These results highlight the importance of genes involved in both innate and adaptive immunity in Sjögren’s syndrome.

Encore: Genetic Association Interaction Network centrality pipeline and application to SLE exome data

Open source tools are needed to facilitate the construction, analysis, and visualization of gene-gene interaction networks for sequencing data. To address this need, we present Encore, an open source network analysis pipeline for genome-wide association studies and rare variant data. Encore constructs Genetic Association Interaction Networks or epistasis networks using two optional approaches: our previous information-theory method or a generalized linear model approach. Additionally, Encore includes multiple data filtering options, including Random Forest/Random Jungle for main effect enrichment and Evaporative Cooling and Relief-F filters for enrichment of interaction effects. Encore implements SNPrank network centrality for identifying susceptibility hubs (nodes containing a large amount of disease susceptibility information through the combination of multivariate main effects and multiple gene-gene interactions in the network), and it provides appropriate files for interactive visualization of a network using tools from our online Galaxy instance. We implemented these algorithms in C++ using OpenMP for shared-memory parallel analysis on a server or desktop. To demonstrate Encore's utility in analysis of genetic sequencing data, we present an analysis of exome resequencing data from healthy individuals and those with Systemic Lupus Erythematous (SLE). Our results verify the importance of the previously associated SLE genes HLA-DRB and NCF2, and these two genes had the highest gene-gene interaction degrees among the susceptibility hubs. An additional 14 genes previously associated with SLE emerged in our epistasis network model of the exome data, and three novel candidate genes, ST8SIA4, CMTM4, and C2CD4B, were implicated in the model. In summary, we present a comprehensive tool for epistasis network analysis and the first such analysis of exome data from a genetic study of SLE.

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway.

METHODS

We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively.

RESULTS

We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression.

CONCLUSION

Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

Association of ANXA11 genetic variation with sarcoidosis in African Americans and European Americans

A recent genome-wide association study in a German population and two subsequent studies in European populations found that a non-synonymous single-nucleotide polymorphism (SNP), rs1049550, within the annexin A11 (ANXA11) gene was associated with susceptibility to sarcoidosis. We sought to identify additional ANXA11 variants independently associated with sarcoidosis, determine whether any sarcoidosis-associated ANXA11 variants were associated with chest radiographic phenotypes, and explore human leukocyte antigen (HLA) SNP-SNP interactions with ANXA11. A total of 209 SNPs spanning 100 kb including the 5' promoter, coding, and 3' untranslated regions of ANXA11 were genotyped for 1689 sarcoidosis cases and 1252 controls. After adjustment for rs1049550, two additional novel ANXA11 sarcoidosis associations were identified only in African Americans--rs61860052 (odds ratio (OR)=0.62; 95% confidence interval (CI)=0.40-0.97) and rs4377299 (OR=1.31; 95% CI=1.06-1.63). These associations were more pronounced in radiologically-classified Scadding stage IV sarcoidosis cases. We also identified a significant SNP-SNP interaction between rs1049550 and a sarcoidosis risk SNP (rs9268839) near the HLA-DRA locus. This further genetic dissection of ANXA11 may provide additional insight into the immune dysregulation characteristic of sarcoidosis pathophysiology.

Genome-wide association study of African and European Americans implicates multiple shared and ethnic specific loci in sarcoidosis susceptibility

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, P(AA-meta) = 6.51 × 10(-10)) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.

Genetics of Sjögren's syndrome in the genome-wide association era

While Sjögren's syndrome (SS) is more common than related autoimmune disorders, such as systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), scientific and medical research in SS has lagged behind significantly. This is especially true in the field of SS genetics, where efforts to date have relied heavily on candidate gene approaches. Within the last decade, the advent of the genome-wide association (GWA) scan has altered our understanding of disease pathogenesis in hundreds of disorders through the successful identification of novel risk loci. With strong evidence for a genetic component in SS as evidenced by familial aggregation of SS as well as similarities between SS and SLE and RA, the application of GWA approaches would likely yield numerous novel risk loci in SS. Here we review the fundamental scientific principles employed in GWA scans as well as the limitations of this tool, and we discuss the application of GWA scans in determining genetic variants at play in complex disease. We also examine the successful application of GWA scans in SLE, which now has more than 40 confirmed risk loci, and consider the possibility for a similar trajectory of SS genetic discovery in the era of GWA scans. Ultimately, the GWA studies that will be performed in SS have the potential to identify a myriad of novel genetic loci that will allow scientists to begin filling in the gaps in our understanding of the SS pathogenesis.

The genomics of autoimmune disease in the era of genome-wide association studies and beyond

Recent advances in the field of genetics have dramatically changed our understanding of autoimmune disease. Candidate gene and, more recently, genome-wide association (GWA) studies have led to an explosion in the number of loci and pathways known to contribute to autoimmune phenotypes. Since the 1970s, researchers have known that several alleles in the MHC region play a role in the pathogenesis of many autoimmune diseases. More recent work has identified numerous risk loci involving both the innate and adaptive immune responses. However, much remains to be learned about the heritability of autoimmune conditions. Most regions found through GWA scans have yet to isolate the association to the causal allele(s) responsible for conferring disease risk. A role for rare variants (allele frequencies of <1%) has begun to emerge. Future research will use next-generation sequencing (NGS) technology to comprehensively evaluate the human genome for risk variants. Whole-transcriptome sequencing is now possible, which will provide much more detailed gene expression data. The dramatic drop in the cost and time required to sequence the entire human genome will ultimately make it possible for this technology to be used as a clinical diagnostic tool.