High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis

A novel method to identify hub pathways of rheumatoid arthritis based on differential pathway networks

The aim of the current study was to identify hub pathways of rheumatoid arthritis (RA) using a novel method based on differential pathway network (DPN) analysis. The present study proposed a DPN where protein-protein interaction (PPI) network was integrated with pathway-pathway interactions. Pathway data was obtained from background PPI network and the Reactome pathway database. Subsequently, pathway interactions were extracted from the pathway data by building randomized gene-gene interactions and a weight value was assigned to each pathway interaction using Spearman correlation coefficient (SCC) to identify differential pathway interactions. Differential pathway interactions were visualized using Cytoscape to construct a DPN. Topological analysis was conducted to identify hub pathways that possessed the top 5% degree distribution of DPN. Modules of DPN were mined according to ClusterONE. A total of 855 pathways were selected to build pathway interactions. By filtrating pathway interactions of weight values >0.7, a DPN with 312 nodes and 791 edges was obtained. Topological degree analysis revealed 15 hub pathways, such as heparan sulfate/heparin-glycosaminoglycan (HS-GAG) degradation, HS-GAG metabolism and keratan sulfate degradation for RA based on DPN. Furthermore, hub pathways were also important in modules, which validated the significance of hub pathways. In conclusion, the proposed method is a computationally efficient way to identify hub pathways of RA, which identified 15 hub pathways that may be potential biomarkers and provide insight to future investigation and treatment of RA.

Genotypic variability based association identifies novel non-additive loci DHCR7 and IRF4 in sero-negative rheumatoid arthritis

Sero-negative rheumatoid arthritis (RA) is a highly heterogeneous disorder with only a few additive loci identified to date. We report a genotypic variability-based genome-wide association study (vGWAS) of six cohorts of sero-negative RA recruited in Europe and the US that were genotyped with the Immunochip. A two-stage approach was used: (1) a mixed model to partition dichotomous phenotypes into an additive component and non-additive residuals on the liability scale and (2) the Levene’s test to assess equality of the residual variances across genotype groups. The vGWAS identified rs2852853 (P = 1.3e-08, DHCR7) and rs62389423 (P = 1.8e-05, near IRF4) in addition to two previously identified loci (HLA-DQB1 and ANKRD55), which were all statistically validated using cross validation. DHCR7 encodes an enzyme important in cutaneous synthesis of vitamin D and DHCR7 mutations are believed to be important for early humans to adapt to Northern Europe where residents have reduced ultraviolet-B exposure and tend to have light skin color. IRF4 is a key locus responsible for skin color, with a vitamin D receptor-binding interval. These vGWAS results together suggest that vitamin D deficiency is potentially causal of sero-negative RA and provide new insights into the pathogenesis of the disorder.

Rheumatoid Arthritis Naive T Cells Share Hypermethylation Sites With Synoviocytes

Objective

To determine whether differentially methylated CpGs in synovium‐derived fibroblast‐like synoviocytes (FLS) of patients with rheumatoid arthritis (RA) were also differentially methylated in RA peripheral blood (PB) samples.

Methods

For this study, 371 genome‐wide DNA methylation profiles were measured using Illumina HumanMethylation450 BeadChips in PB samples from 63 patients with RA and 31 unaffected control subjects, specifically in the cell subsets of CD14+ monocytes, CD19+ B cells, CD4+ memory T cells, and CD4+ naive T cells.

Results

Of 5,532 hypermethylated FLS candidate CpGs, 1,056 were hypermethylated in CD4+ naive T cells from RA PB compared to control PB. In analyses of a second set of CpG candidates based on single‐nucleotide polymorphisms from a genome‐wide association study of RA, 1 significantly hypermethylated CpG in CD4+ memory T cells and 18 significant CpGs (6 hypomethylated, 12 hypermethylated) in CD4+ naive T cells were found. A prediction score based on the hypermethylated FLS candidates had an area under the curve of 0.73 for association with RA case status, which compared favorably to the association of RA with the HLA–DRB1 shared epitope risk allele and with a validated RA genetic risk score.

Conclusion

FLS‐representative DNA methylation signatures derived from the PB may prove to be valuable biomarkers for the risk of RA or for disease status.

Association Between Genetic Variation in FOXO3 and Reductions in Inflammation and Disease Activity in Inflammatory Polyarthritis

Objective

Genetic variation in FOXO3 (tagged by rs12212067) has been associated with a milder course of rheumatoid arthritis (RA) and shown to limit monocyte‐driven inflammation through a transforming growth factor β1–dependent pathway. This genetic association, however, has not been consistently observed in other RA cohorts. We sought to clarify the contribution of FOXO3 to prognosis in RA by combining detailed analysis of nonradiographic disease severity measures with an in vivo model of arthritis.

Methods

Collagen‐induced arthritis, the most commonly used mouse model of RA, was used to assess how Foxo3 contributes to arthritis severity. Using clinical, serologic, and biochemical methods, the arthritis that developed in mice carrying a loss‐of‐function mutation in Foxo3 was compared with that which occurred in littermate controls. The association of rs12212067 with nonradiographic measures of RA severity, including the C‐reactive protein level, the swollen joint count, the tender joint count, the Disease Activity Score in 28 joints, and the Health Assessment Questionnaire score, were modeled longitudinally in a large prospective cohort of patients with early RA.

Results

Loss of Foxo3 function resulted in more severe arthritis in vivo (both clinically and histologically) and was associated with higher titers of anticollagen antibodies and interleukin‐6 in the blood. Similarly, rs12212067 (a single‐nucleotide polymorphism that increases FOXO3 transcription) was associated with reduced inflammation, both biochemically and clinically, and with lower RA activity scores.

Conclusion

Consistent with its known role in restraining inflammatory responses, FOXO3 limits the severity of in vivo arthritis and, through genetic variation that increases its transcription, is associated with reduced inflammation and disease activity in RA patients, effects that result in less radiographic damage.

HLA class II and rheumatoid arthritis: the bumpy road of revelation

Rheumatoid arthritis (RA) is a chronic auto-immune disease primarily targeting the joints. Approximately 1% of the population is affected by RA, and despite the improvements in therapeutic interventions, elucidation of the disease pathogenesis is still in its infancy. RA patients can be subdivided on basis of the presence of autoantibodies, especially anti-citrullinated protein antibodies (ACPA). ACPA⁺ and ACPA^- disease most likely differ in aetiology, as different genetic and environmental risk factors are associated with these two disease entities. For ACPA⁺ RA disease, the genetic factors associating with disease mainly comprised of human leukocyte antigen (HLA) class II molecules. The predisposing HLA-DR alleles have been depicted as the 'HLA Shared Epitope (SE) alleles', as these alleles encode a similar sequence, the shared epitope sequence, within the beta chain of the HLA-DR molecule. In addition to the involvement of the HLA-SE alleles in the development of ACPA⁺ RA disease, other HLA-DR molecules have been shown to confer protection against this disease entity. The protective HLA molecules have, instead of the SE-motif, a different but shared sequence at the same location in the beta chain of HLA-DR molecules, consisting of the amino acid residues DERAA. The possible contributions of the predisposing and protective HLA molecules in association with ACPA-positive RA are discussed in this review.

Replication of Associations of Genetic Loci Outside the HLA Region With Susceptibility to Anti-Cyclic Citrullinated Peptide-Negative Rheumatoid Arthritis

Objective

Genetic polymorphisms within the HLA region explain only a modest proportion of anti–cyclic citrullinated peptide (anti‐CCP)–negative rheumatoid arthritis (RA) heritability. However, few non‐HLA markers have been identified so far. This study was undertaken to replicate the associations of anti‐CCP–negative RA with non‐HLA genetic polymorphisms demonstrated in a previous study.

Methods

The Rheumatoid Arthritis Consortium International densely genotyped 186 autoimmune‐related regions in3,339 anti‐CCP–negative RA patients and 15,870 controls across 6 different populations using the Illumina ImmunoChip array. We performed a case–control replication study of the anti‐CCP–negative markers with the strongest associations in that discovery study, in an independent cohort of anti‐CCP–negative UK RA patients. Individuals from the arcOGENConsortium and Wellcome Trust Case Control Consortium were used as controls. Genotyping in cases was performed using Sequenom MassArray technology. Genome‐wide data from controls were imputed using the 1000 Genomes Phase I integrated variant call set release version 3 as a reference panel.

Results

After genotyping and imputation quality control procedures, data were available for 15 non‐HLA single‐nucleotide polymorphisms in 1,024 cases and 6,348 controls. We confirmed the known markers ANKRD55 (meta‐analysis odds ratio [OR] 0.80; P = 2.8 × 10⁻¹³)and BLK (OR 1.13; P = 7.0 × 10⁻⁶) and identified new and specific markers of anti‐CCP–negative RA (prolactin [PRL] [OR 1.13; P = 2.1 × 10⁻⁶] and NFIA [OR 0.85; P = 2.5 × 10⁻⁶]). Neither of these loci is associated with other common, complex autoimmune diseases.

Conclusion

Anti‐CCP–negative RA and anti‐CCP–positive RA are genetically different disease subsets that only partially share susceptibility factors. Genetic polymorphisms located near the PRL and NFIA genes represent examples of genetic susceptibility factors specific for anti‐CCP–negative RA.

Reduced DICER1 Expression Bestows Rheumatoid Arthritis Synoviocytes Proinflammatory Properties and Resistance to Apoptotic Stimuli

OBJECTIVE

While the regulatory role of individual microRNAs (miRNAs) in rheumatoid arthritis (RA) is well established, the role of DICER1 in the pathogenesis of the disease has not yet been investigated. The purpose of this study was to analyze the expression of factors involved in miRNA biogenesis in fibroblast-like synoviocytes (FLS) from RA patients and to monitor the arthritis triggered by K/BxN serum transfer in mice deficient in the Dicer gene (Dicer(d/d) ).

METHODS

The expression of genes and precursor miRNAs was quantified by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). MicroRNA macroarray profiling was monitored by qRT-PCR. Cytokines were quantified by enzyme-linked immunosorbent assay. Experimental arthritis in mice was achieved by the transfer of serum from K/BxN donors. Apoptosis was quantified using an enzyme-linked immunosorbent assay.

RESULTS

We found decreased DICER1 and mature miRNA expression in synovial fibroblasts from RA patients. These cells were hyperresponsive to lipopolysaccharide, as evidenced by their increased interleukin-6 secretion upon stimulation. Experimental serum-transfer arthritis in Dicer(d/d) mice confirmed that an unbalanced biogenesis of miRNAs correlated with an enhanced inflammatory response. Synoviocytes from both RA patients and Dicer(d/d) mice exhibited increased resistance to apoptotic stimuli.

CONCLUSION

The findings of this study further substantiate the important role of DICER1 in the maintenance of homeostasis and the regulation of inflammatory responses.

Dense Genotyping of Immune-Related Regions Identifies Loci for Rheumatoid Arthritis Risk and Damage in African Americans

Over 100 risk loci for rheumatoid arthritis (RA) have been identified in individuals of European and Asian descent, but the genetic basis for RA in African Americans is less well understood. We genotyped 610 African Americans with autoantibody positive RA and 933 African American controls on the ImmunoChip (iChip) array. Using multivariable regression we evaluated the association between iChip markers and the risk of RA and radiographic severity. The single nucleotide polymorphism (SNP) rs1964995 (OR = 1.97, p = 1.28 × 10^-15) near HLA-DRB1 was the most strongly associated risk SNP for RA susceptibility; SNPs in AFF3, TNFSF11, and TNFSF18 loci were suggestively associated (10^-4 < p < 3.1 × 10^-6). Trans-ethnic fine mapping of AFF3 identified a 90% credible set containing previously studied variants including rs9653442, rs7608424, and rs6712515 as well as the novel candidate variant rs11681966; several of these likely influence AFF3 gene expression level. Variants in TNFRSF9, CTLA4, IL2RA, C5/TRAF1, and ETS1 - but no variants within the major histocompatibility complex - were associated with RA radiographic severity. Conditional regression and pairwise linkage disequilibrium (LD) analyses suggest that additional pathogenic variants may be found in ETS1 and IL2RA beyond those found in other ethnicities. In summary, we use the dense genotyping of the iChip array and unique LD structure of African Americans to validate known risk loci for RA susceptibility and radiographic severity, and to better characterize the associations of AFF3, ETS1, and IL2RA.

High-Density Genetic Mapping Identifies New Susceptibility Variants in Sarcoidosis Phenotypes and Shows Genomic-driven Phenotypic Differences

RATIONALE

Sarcoidosis is a multisystem disease of unknown cause. Löfgren's syndrome (LS) is a characteristic subgroup of sarcoidosis that is associated with a good prognosis in sarcoidosis. However, little is known about its genetic architecture or its broader phenotype, non-LS sarcoidosis.

OBJECTIVES

To address the genetic architecture of sarcoidosis phenotypes, LS and non-LS.

METHODS

An association study in a white Swedish cohort of 384 LS, 664 non-LS, and 2,086 control subjects, totaling 3,134 subjects using a fine-mapping genotyping platform was conducted. Replication was performed in four independent cohorts, three of white European descent (Germany, n = 4,975; the Netherlands, n = 613; and Czech Republic, n = 521), and one of black African descent (United States, n = 1,657), totaling 7,766 subjects.

MEASUREMENTS AND MAIN RESULTS

A total of 727 LS-associated variants expanding throughout the extended major histocompatibility complex (MHC) region and 68 non-LS-associated variants located in the MHC class II region were identified and confirmed. A shared overlap between LS and non-LS defined by 17 variants located in the MHC class II region was found. Outside the MHC region, two LS-associated loci, in ADCY3 and between CSMD1 and MCPH1, were observed and replicated.

CONCLUSIONS

Comprehensive and integrative analyses of genetics, transcription, and pathway modeling on LS and non-LS indicates that these sarcoidosis phenotypes have different genetic susceptibility, genomic distributions, and cellular activities, suggesting distinct molecular mechanisms in pathways related to immune response with a common region.

The genetics revolution in rheumatology: large scale genomic arrays and genetic mapping

Susceptibility to rheumatic diseases, such as osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, systemic lupus erythematosus, juvenile idiopathic arthritis and psoriatic arthritis, includes a large genetic component. Understanding how an individual's genetic background influences disease onset and outcome can lead to a better understanding of disease biology, improved diagnosis and treatment, and, ultimately, to disease prevention or cure. The past decade has seen great progress in the identification of genetic variants that influence the risk of rheumatic diseases. The challenging task of unravelling the function of these variants is ongoing. In this Review, the major insights from genetic studies, gained from advances in technology, bioinformatics and study design, are discussed in the context of rheumatic disease. In addition, pivotal genetic studies in the main rheumatic diseases are highlighted, with insights into how these studies have changed the way we view these conditions in terms of disease overlap, pathways of disease and potential new therapeutic targets. Finally, the limitations of genetic studies, gaps in our knowledge and ways in which current genetic knowledge can be fully translated into clinical benefit are examined.

Genetics of rheumatoid arthritis susceptibility, severity, and treatment response

A decade after the first genome-wide association study in rheumatoid arthritis (RA), a plethora of genetic association studies have been published on RA and its clinical or serological subtypes. We review the major milestones in the study of the genetic architecture of RA susceptibility, severity, and response to treatment. We set the scientific context necessary for non-geneticists to understand the potential clinical applications of human genetics and its significance for a stratified approach to the management of RA in the future.

Phenome-wide scanning identifies multiple diseases and disease severity phenotypes associated with HLA variants

Although many phenotypes have been associated with variants in human leukocyte antigen (HLA) genes, the full phenotypic impact of HLA variants across all diseases is unknown. We imputed HLA genomic variation from two populations of 28,839 and 8431 European ancestry individuals and tested association of HLA variation with 1368 phenotypes. A total of 104 four-digit and 92 two-digit HLA allele phenotype associations were significant in both discovery and replication cohorts, the strongest being HLA-DQB1*03:02 and type 1 diabetes. Four previously unidentified associations were identified across the spectrum of disease with two- and four-digit HLA alleles and 10 with nonsynonymous variants. Some conditions associated with multiple HLA variants and stronger associations with more severe disease manifestations were identified. A comprehensive, publicly available catalog of clinical phenotypes associated with HLA variation is provided. Examining HLA variant disease associations in this large data set allows comprehensive definition of disease associations to drive further mechanistic insights.

Immune-Array Analysis in Sporadic Inclusion Body Myositis Reveals HLA-DRB1 Amino Acid Heterogeneity Across the Myositis Spectrum

OBJECTIVE

Inclusion body myositis (IBM) is characterized by a combination of inflammatory and degenerative changes affecting muscle. While the primary cause of IBM is unknown, genetic factors may influence disease susceptibility. To determine genetic factors contributing to the etiology of IBM, we conducted the largest genetic association study of the disease to date, investigating immune-related genes using the Immunochip.

METHODS

A total of 252 Caucasian patients with IBM were recruited from 11 countries through the Myositis Genetics Consortium and compared with 1,008 ethnically matched controls. Classic HLA alleles and amino acids were imputed using SNP2HLA.

RESULTS

The HLA region was confirmed as the most strongly associated region in IBM (P = 3.58 × 10-33 ). HLA imputation identified 3 independent associations (with HLA-DRB1*03:01, DRB1*01:01, and DRB1*13:01), although the strongest association was with amino acid positions 26 and 11 of the HLA-DRB1 molecule. No association with anti-cytosolic 5'-nucleotidase 1A-positive status was found independent of HLA-DRB1*03:01. There was no association of HLA genotypes with age at onset of IBM. Three non-HLA regions reached suggestive significance, including the chromosome 3 p21.31 region, an established risk locus for autoimmune disease, where a frameshift mutation in CCR5 is thought to be the causal variant.

CONCLUSION

This is the largest, most comprehensive genetic association study to date in IBM. The data confirm that HLA is the most strongly associated region and identifies novel amino acid associations that may explain the risk in this locus. These amino acid associations differentiate IBM from polymyositis and dermatomyositis and may determine properties of the peptide-binding groove, allowing it to preferentially bind autoantigenic peptides. A novel suggestive association within the chromosome 3 p21.31 region suggests a role for CCR5.

Systematic analysis of coronary artery disease datasets revealed the potential biomarker and treatment target

Coronary artery disease caused about 1 of every 7 deaths in the United States and early prevention was potential to decrease the incidence and mortality. We aimed to figure the genes involving in the coronary artery disease using meta-anlaysis. Five datasets of coronary heart disease from GEO series were retrieved and data preprocessing and quality control were carried out. Moderated t-test was used to decide the differentially expressed genes for a single dataset. And the combined p-value using systematic-analysis methods were conducted using MetaDE. The pathway enrichment was carried out using Reactome database. Protein-protein interactions of the identified differentially expressed genes were also analyzed using STRING v10.0 online tool. After removing unidentified or intermediate samples and a total of 238 cases and 189 matched or partially matched control from five microarray datasets were retrieved from GEO. Six different quality control measures were calculated and PCA biplots were plotted in order to visualize the quantitative measure. The first two PCs captured 91% of the variance and we decided to include all of the datasets for systematic analysis. Using the FDR cut-off as 0.1, nine genes, including LFNG, ID3, PLA2G7, FOLR3, PADI4, ARG1, IL1R2, NFIL3 and MGAM, were differentially expressed according to maxP. Their protein-protein interactions showed that they were closely connected and 24 Reactome pathways were related to coronary artery disease. We concluded that pathways related to immune responses, especially neutrophil degranulation, were associated with coronary heart disease.

Neutrophil-mediated IL-6 receptor trans-signaling and the risk of chronic obstructive pulmonary disease and asthma

The Asp358Ala variant in the interleukin-6 receptor (IL-6R) gene has been implicated in asthma, autoimmune and cardiovascular disorders, but its role in other respiratory conditions such as chronic obstructive pulmonary disease (COPD) has not been investigated. The aims of this study were to evaluate whether there is an association between Asp358Ala and COPD or asthma risk, and to explore the role of the Asp358Ala variant in sIL-6R shedding from neutrophils and its pro-inflammatory effects in the lung. We undertook logistic regression using data from the UK Biobank and the ECLIPSE COPD cohort. Results were meta-analyzed with summary data from a further three COPD cohorts (7,519 total cases and 35,653 total controls), showing no association between Asp358Ala and COPD (OR = 1.02 [95% CI: 0.96, 1.07]). Data from the UK Biobank showed a positive association between the Asp358Ala variant and atopic asthma (OR = 1.07 [1.01, 1.13]). In a series of in vitro studies using blood samples from 37 participants, we found that shedding of sIL-6R from neutrophils was greater in carriers of the Asp358Ala minor allele than in non-carriers. Human pulmonary artery endothelial cells cultured with serum from homozygous carriers showed an increase in MCP-1 release in carriers of the minor allele, with the difference eliminated upon addition of tocilizumab. In conclusion, there is evidence that neutrophils may be an important source of sIL-6R in the lungs, and the Asp358Ala variant may have pro-inflammatory effects in lung cells. However, we were unable to identify evidence for an association between Asp358Ala and COPD.

Micro RNA-146a But Not IRAK1 is Associated with Rheumatoid Arthritis in the Tunisian Population

BACKGROUND

Rheumatoid arthritis (RA) is characterized by the production of an array of proinflammatory cytokines through the nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling pathway. Interleukin-1 receptor (IL-1R) and Toll-like receptors contain a common cytoplasmic motif the Toll/IL-1R (TIR) homology domain. This motif is required for NF-κB activation. IL-1R-associated kinase 1 (IRAK1) is a key adapter molecule recruited during the signaling cascade of the TIR. Its gene expression is regulated by the micro-RNA (miR)-146a.

OBJECTIVE

We investigated the role of IRAK1 single-nucleotide polymorphism (SNP) rs3027898 (IRAK1 rs3027898) and miR-146a SNP rs2910164 (miR-146a rs2910164) in Tunisian patients with RA and their association with C reactive protein (CRP), rheumatoid factor (RF), anticyclic citrullinated peptide (anti-CCP) antibodies, and erosion.

PATIENTS AND METHODS

In a cohort of 172 adult RA patients and 224 matched controls, IRAK1 rs3027898 genotyping was determined by mutagenically separated polymerase chain reaction (MS-PCR) with newly designed primers, and miR-146a rs2910164 genotyping was determined by fragment length polymorphism PCR-restriction (RFLP-PCR).

RESULTS

The IRAK1 rs3027898 A allele was detected in 67% of RA patients and 70% of controls indicating that it is not associated with RA in codominant, dominant, or recessive models even after stratification by age and gender. The miR-146a rs2910164 G allele was detected in 76% of RA patients and 68% of controls, thus the C allele confers some protection based on a dominant model [CC+GC (odds ratio (95% confidence interval) = 0.6 (0.3-0.9), p = 0.03)]. No association with CRP, RF, anti-CCP, or erosion was found for either SNPs.

CONCLUSION

The IRAK1 rs3027898 was not associated with RA, whereas C allele of miR-146a rs2910164 was found to be protective. Functional studies are required to investigate the exact role of miR-146a rs2910164 during RA.

Polymorphisms of RAD51B are associated with rheumatoid arthritis and erosion in rheumatoid arthritis patients

Rheumatoid arthritis (RA) is a common, chronic autoimmune disease affecting 0.5–1.0% of adults worldwide, including approximately 4.5–5.0 million patients in China. The genetic etiology and pathogenesis of RA have not yet been fully elucidated. Recently, one new RA susceptibility gene (RAD51B) has been identified in Korean and European populations. In this study, we designed a two-stage case-control study to further assess the relationship of common variants in the RAD51B gene with increased risk of RA in a total of 965 RA patients and 2,511 unrelated healthy controls of Han Chinese ancestry. We successfully identified a common variant, rs911263, as being significantly associated with the disease status of RA (P = 4.8 × 10⁻⁵, OR = 0.64). In addition, this SNP was shown to be related to erosion, a clinical assessment of disease severity in RA (P = 2.89 × 10⁻⁵, OR = 0.52). These findings shed light on the role of RAD51B in the onset and severity of RA. More research in the future is needed to clarify the underlying functional link between rs911263 and the disease.

Genetic pleiotropy between age-related macular degeneration and 16 complex diseases and traits

BACKGROUND

Age-related macular degeneration (AMD) is a common condition of vision loss with disease development strongly influenced by environmental and genetic factors. Recently, 34 loci were associated with AMD at genome-wide significance. So far, little is known about a genetic overlap between AMD and other complex diseases or disease-relevant traits.

METHODS

For each of 60 complex diseases/traits with publicly available genome-wide significant association data, the lead genetic variant per independent locus was extracted and a genetic score was calculated for each disease/trait as the weighted sum of risk alleles. The association with AMD was estimated based on 16,144 AMD cases and 17,832 controls using logistic regression.

RESULTS

Of the respective disease/trait variance, the 60 genetic scores explained on average 4.8% (0.27-20.69%) and 16 of them were found to be significantly associated with AMD (Q-values < 0.01, p values from < 1.0 × 10^-16 to 1.9 × 10^-3). Notably, an increased risk for AMD was associated with reduced risk for cardiovascular diseases, increased risk for autoimmune diseases, higher HDL and lower LDL levels in serum, lower bone-mineral density as well as an increased risk for skin cancer. By restricting the analysis to 1824 variants initially used to compute the 60 genetic scores, we identified 28 novel AMD risk variants (Q-values < 0.01, p values from 1.1 × 10^-7 to 3.0 × 10^-4), known to be involved in cardiovascular disorders, lipid metabolism, autoimmune diseases, anthropomorphic traits, ocular disorders, and neurological diseases. The latter variants represent 20 novel AMD-associated, pleiotropic loci. Genes in the novel loci reinforce previous findings strongly implicating the complement system in AMD pathogenesis.

CONCLUSIONS

We demonstrate a substantial overlap of the genetics of several complex diseases/traits with AMD and provide statistically significant evidence for an additional 20 loci associated with AMD. This highlights the possibility that so far unrelated pathologies may have disease pathways in common.

Interplay between genetic and epigenetic mechanisms in rheumatoid arthritis

Genetic and environmental factors contribute to the risk for rheumatoid arthritis (RA), with epigenetics serving as a possible interface through which risk factors contribute to RA. High-throughput technologies for interrogating genome and epigenome, and the availability of genetic and epigenetic datasets across a diversity of cell types, enable the identification of candidate causal genetic variants for RA to study their function in core RA processes. To date, RA risk variants were studied in the immune cells but not joint resident cells, for example, synovial fibroblasts. Synovial fibroblasts from different joints are distinct, anatomically specialized cells, defined by joint-specific transcriptomes, epigenomes and phenotypes. Cell type-specific analysis of epigenetic changes, together with genetic fine mapping and interrogation of chromatin 3D interactions may identify new disease relevant pathways, potential therapeutic targets and biomarkers for RA progression or therapy response.

Human IL-6RhiTIGIT- CD4+CD127lowCD25+ T cells display potent in vitro suppressive capacity and a distinct Th17 profile

To date many clinical studies aim to increase the number and/or fitness of CD4⁺ CD127^lowCD25⁺ regulatory T cells (Tregs) in vivo to harness their regulatory potential in the context of treating autoimmune disease. Here, we sought to define the phenotype and function of Tregs expressing the highest levels of IL-6 receptor (IL-6R). We have identified a population of CD4⁺ CD127^lowCD25⁺ TIGIT⁻ T cells distinguished by their elevated IL-6R expression that lacked expression of HELIOS, showed higher CTLA-4 expression, and displayed increased suppressive capacity compared to IL-6R^hiTIGIT⁺ Tregs. IL-6R^hiTIGIT⁻ CD127^lowCD25⁺ T cells contained a majority of cells demethylated at FOXP3 and displayed a Th17 transcriptional signature, including RORC (RORγt) and the capacity of producing both pro- and anti-inflammatory cytokines, such as IL-17, IL-22 and IL-10. We propose that in vivo, in the presence of IL-6-associated inflammation, the suppressive function of CD4⁺ CD127^lowCD25⁺ FOXP3⁺ IL-6R^hiTIGIT⁻ T cells is temporarily disarmed allowing further activation of the effector functions and potential pathogenic tissue damage.

•

IL-6R is highly expressed in certain Treg subsets.

•

IL-6R^hiTIGIT⁻ CD127^lowCD25⁺ T cells contain a subset of antigen-experienced Tregs with potent suppression capacity.

•

IL-6R^hiTIGIT⁻ Tregs display a Th17 transcriptional profile ex vivo, and the capacity to migrate to the gut.

•

IL-2 treatment in humans elicits the trafficking and expansion of Tregs in circulation.

Dense genotyping of immune-related loci identifies HLA variants associated with increased risk of collagenous colitis

OBJECTIVE

Collagenous colitis (CC) is a major cause of chronic non-bloody diarrhoea, particularly in the elderly female population. The aetiology of CC is unknown, and still poor is the understanding of its pathogenesis. This possibly involves dysregulated inflammation and immune-mediated reactions in genetically predisposed individuals, but the contribution of genetic factors to CC is underinvestigated. We systematically tested immune-related genes known to impact the risk of several autoimmune diseases for their potential CC-predisposing role.

DESIGN

Three independent cohorts of histologically confirmed CC cases (N=314) and controls (N=4299) from Sweden and Germany were included in a 2-step association analysis. Immunochip and targeted single nucleotide polymorphism (SNP) genotype data were produced, respectively, for discovery and replication purposes. Classical human leucocyte antigen (HLA) variants at 2-digit and 4-digit resolution were obtained via imputation from single marker genotypes. SNPs and HLA variants passing quality control filters were tested for association with CC with logistic regression adjusting for age, sex and country of origin.

RESULTS

Forty-two markers gave rise to genome-wide significant association signals, all contained within the HLA region on chromosome 6 (best p=4.2×10^-10 for SNP rs4143332). Among the HLA variants, most pronounced risk effects were observed for 8.1 haplotype alleles including DQ2.5, which was targeted and confirmed in the replication data set (p=2.3×10^-11; OR=2.06; 95% CI (1.67 to 2.55) in the combined analysis).

CONCLUSIONS

HLA genotype associates with CC, thus implicating HLA-related immune mechanisms in its pathogenesis.

MHC class II alleles associated with Th1 rather than Th17 type immunity drive the onset of early arthritis in a rat model of rheumatoid arthritis

Polymorphisms in the MHC class II (MHCII) genes are strongly associated with rheumatoid arthritis, supporting the importance of autoreactive T helper (Th) cells for the development of this disease. Here, we used pristane-induced arthritis (PIA), induced by the non-antigenic hydrocarbon pristane, to study the impact of different MHCII alleles on T-cell activation and differentiation. In MHCII-congenic rats with disease-promoting MHCII alleles, pristane primarily induced activation of Th1 cells, whereas activated T cells were Th17 biased in rats with protective MHCII alleles. Neutralization of IFN-γ during T-cell activation abrogated the development of disease, suggesting that Th1 immunity is important for disease induction. Neutralization of IL-17, by contrast, suppressed arthritis only when performed in rats with established disease. Adoptive T-cell transfers showed that T cells acquired arthritogenic capacity earlier in strains with a prevailing Th1 response. Moreover, upon pristane injection, these strains exhibited more Ag-primed OX40+ and proliferating T cells of polyclonal origin. These data show that T cells are polarized upon the first encounter with peptide-MHCII complexes in an allele-dependent fashion. In PIA, the polyclonal expansion of autoreactive Th1 cells was necessary for the onset of arthritis, while IL-17 mediated immunity contributed to the progression to chronic disease.

Receptor activator of nuclear factor kappa-B ligand (RANKL) but not sclerostin or gene polymorphisms is related to joint destruction in early rheumatoid arthritis

The aim of this study was to analyze relationships between receptor activator of nuclear factor kappa-B (RANKL), sclerostin and their gene polymorphisms with radiological progression in patients with early rheumatoid arthritis (RA). Patients with early RA (n = 407, symptomatic <1 year) (ARA criteria) examined radiologically at inclusion and after 24 months were consecutively included. Disease activity score and C-reactive protein were regularly recorded. Sclerostin, RANKL, and anti-CCP2 antibodies were analyzed in plasma at baseline using ELISAs. Data on gene polymorphism for sclerostin and RANKL were extracted from Immunochip analysis. Sex- and age-matched controls (n = 71) were identified from the Medical Biobank of Northern Sweden. The concentration of RANKL was significantly higher in patients compared with controls, median (IQR) 0.56 (0.9) nmol/L and 0.20 (0.25) nmol/L (p < 0.001), and in anti-CCP2-positive patients compared with sero-negative individuals. Sclerostin was significantly increased in female patients 0.59 (0.47–0.65) ng/mL compared with female controls 0.49 (0.4–0.65) ng/mL (p < 0.02). RANKL concentration was related to the Larsen score at baseline (p < 0.01), after 24 months (p < 0.001), and to radiological progression at 24 months (p < 0.001). Positivity of RANKL and anti-CCP2 yielded significant risk for progression with negativity for both as reference. No single nucleotide polymorphism encoding TNFSF11 or SOST was associated with increased concentrations of the factors. The concentration of RANKL was related to the Larsen score at baseline, at 24 months, and radiological progression at 24 months particularly in anti-CCP2-positive patients, while the concentration of sclerostin was unrelated to radiological findings.

TMEM187-IRAK1 Polymorphisms Associated with Rheumatoid Arthritis Susceptibility in Tunisian and French Female Populations: Influence of Geographic Origin

Polymorphisms have been identified in the Xq28 locus as risk loci for rheumatoid arthritis (RA). Here, we investigated the association between three polymorphisms in the Xq28 region containing TMEM187 and IRAK1 (rs13397, rs1059703, and rs1059702) in two unstudied populations: Tunisian and French. The rs13397 G and rs1059703 T major alleles were significantly increased in RA patients (n = 408) compared with age-matched controls (n = 471) in both Tunisian and French women. These results were confirmed by a meta-analysis replication study including two independent Greek and Korean cohorts. The rs1059702 C major allele was significantly associated with RA, only with French women. In the French population, the GTC haplotype displayed a protective effect against RA, while the ATC, GCC, and GTT haplotypes conferred significant risk for RA. No association for these haplotypes was found in the Tunisian population. Our results replicated for the first time the association of the three Xq28 polymorphisms with RA risk in Tunisian and French populations and suggested that RA susceptibility is associated with TMEM187-IRAK1 polymorphisms in women. Our data further support the involvement of X chromosome in RA susceptibility and evidence ethnicities differences that might be explained by differences in the frequencies of SE HLA-DRB1 alleles between both populations.

Discovery of new candidate genes for rheumatoid arthritis through integration of genetic association data with expression pathway analysis

Background

Here we integrate verified signals from previous genetic association studies with gene expression and pathway analysis for discovery of new candidate genes and signaling networks, relevant for rheumatoid arthritis (RA).

Method

RNA-sequencing-(RNA-seq)-based expression analysis of 377 genes from previously verified RA-associated loci was performed in blood cells from 5 newly diagnosed, non-treated patients with RA, 7 patients with treated RA and 12 healthy controls. Differentially expressed genes sharing a similar expression pattern in treated and untreated RA sub-groups were selected for pathway analysis. A set of “connector” genes derived from pathway analysis was tested for differential expression in the initial discovery cohort and validated in blood cells from 73 patients with RA and in 35 healthy controls.

Results

There were 11 qualifying genes selected for pathway analysis and these were grouped into two evidence-based functional networks, containing 29 and 27 additional connector molecules. The expression of genes, corresponding to connector molecules was then tested in the initial RNA-seq data. Differences in the expression of ERBB2, TP53 and THOP1 were similar in both treated and non-treated patients with RA and an additional nine genes were differentially expressed in at least one group of patients compared to healthy controls. The ERBB2, TP53. THOP1 expression profile was successfully replicated in RNA-seq data from peripheral blood mononuclear cells from healthy controls and non-treated patients with RA, in an independent collection of samples.

Conclusion

Integration of RNA-seq data with findings from association studies, and consequent pathway analysis implicate new candidate genes, ERBB2, TP53 and THOP1 in the pathogenesis of RA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1220-5) contains supplementary material, which is available to authorized users.

Genetic and environmental risk factors for rheumatoid arthritis

Multiple genetic and environmental factors have been associated with an increased risk for rheumatoid arthritis (RA). Of these, the strongest associations have been seen with female sex, a family history of RA, the genetic factor the "shared epitope," and exposure to tobacco smoke. There is also renewed interest in mucosal inflammation and microbial factors as contributors to the development of RA. However, the identification of a "preclinical" period of RA that can be defined as local or systemic autoimmunity as measured by autoantibodies and other biomarkers prior to the development of clinically apparent synovitis suggests that the risk factors for RA are acting long prior to first clinical evidence of IA. As such, a major challenge to the field will be to investigate the full spectrum of the development of RA, from initiation and propagation of autoimmunity during preclinical RA and transition to clinically apparent synovitis and classifiable RA, to determine which genetic and environmental factors are important at each stage of disease development. Understanding the exact role and timing of action of risk factors for RA is especially important given the advent of prevention trials in RA, and the hope that a full understanding of genetic and environmental factors in RA could lead to effective preventive interventions.

VAV1 regulates experimental autoimmune arthritis and is associated with anti-CCP negative rheumatoid arthritis

Rheumatoid arthritis (RA) patients can be stratified into two subgroups defined by the presence or absence of antibodies against citrullinated circular peptides (anti-CCP) with most of the genetic association found in anti-CCP positive RA. Here we addressed the role of VAV1, previously associated to multiple sclerosis (MS), in the pathogenesis of RA in experimental models and in a genetic association study. Experimental arthritis triggered by pristane or collagen type II was induced in DA rats and in the DA.BN-R25 congenic line that carries a polymorphism in Vav1. Difference in arthritis severity was observed only after immunization with pristane. In a case-control study, 34 SNPs from VAV1 locus were analyzed by Immunochip genotyping in 11475 RA patients (7573 anti-CCP positive and 3902 negative) and 15,870 controls in six cohorts of European Caucasians. A combination of the previous MS-associated haplotype and two additional SNPs was associated with anti-CCP negative RA (alleles G-G-A-A of rs682626-rs2546133-rs2617822-rs12979659, OR=1.13, P=1.27 × 10^-5). The same markers also contributed to activity of RA at baseline with the strongest association in the anti-CCP negative group for the rs682626-rs12979659 G-A haplotype (β=-0.283, P=0.0048). Our study suggests a role for VAV1 and T-cell signaling in the pathology of anti-CCP-negative RA.

Immunochip analysis identifies novel susceptibility loci in the human leukocyte antigen region for acquired thrombotic thrombocytopenic purpura

Essentials Genetic predisposition to acquired thrombotic thrombocytopenic purpura (aTTP) is mainly unknown. Genetic risk factors for aTTP were studied by Immunochip analysis and replication study. Human leukocyte antigen (HLA) variant rs6903608 conferred a 2.5-fold higher risk of developing aTTP. rs6903608 and HLA-DQB1*05:03 may explain most of the HLA association signal in aTTP. Click to hear Dr Cataland's presentation on acquired thrombotic thrombocytopenic purpura SUMMARY: Background Acquired thrombotic thrombocytopenic purpura (TTP) is a rare, life-threatening thrombotic microangiopathy associated with the development of autoantibodies against the von Willebrand factor-cleaving protease ADAMTS-13. Similarly to what has been found for other autoimmune disorders, there is evidence of a genetic contribution, including the association of the human leukocyte antigen (HLA) class II complex with disease risk. Objective To identify novel genetic risk factors in acquired TTP. Patients/Methods We undertook a case-control genetic association study in 190 European-origin TTP patients and 1255 Italian healthy controls by using the Illumina Immunochip. Replication analysis in 88 Italian cases and 456 controls was performed with single-nucleotide polymorphism (SNP) TaqMan assays. Results and conclusion We identified one common variant (rs6903608) located within the HLA class II locus that was independently associated with acquired TTP at genome-wide significance and conferred a 2.6-fold increased risk of developing a TTP episode (95% confidence interval [CI] 2.02-3.27, P = 1.64 × 10-14 ). We also found five non-HLA variants mapping to chromosomes 2, 6, 8 and X that were suggestively associated with the disease: rs9490550, rs115265285, rs5927472, rs7823314, and rs1334768 (nominal P-values ranging from 1.59 × 10-5 to 7.60 × 10-5 ). Replication analysis confirmed the association of HLA variant rs6903608 with acquired TTP (pooled P = 3.95 × 10-19 ). Imputation of classic HLA genes followed by stepwise conditional analysis revealed that the combination of rs6903608 and HLA-DQB1*05:03 may explain most of the HLA association signal in acquired TTP. Our results refined the association of the HLA class II locus with acquired TTP, confirming its importance in the etiology of this autoimmune disease.

Gene-Based Genome-Wide Association Analysis in European and Asian Populations Identified Novel Genes for Rheumatoid Arthritis

Objective

Rheumatoid arthritis (RA) is a complex autoimmune disease. Using a gene-based association research strategy, the present study aims to detect unknown susceptibility to RA and to address the ethnic differences in genetic susceptibility to RA between European and Asian populations.

Methods

Gene-based association analyses were performed with KGG 2.5 by using publicly available large RA datasets (14,361 RA cases and 43,923 controls of European subjects, 4,873 RA cases and 17,642 controls of Asian Subjects). For the newly identified RA-associated genes, gene set enrichment analyses and protein-protein interactions analyses were carried out with DAVID and STRING version 10.0, respectively. Differential expression verification was conducted using 4 GEO datasets. The expression levels of three selected ‘highly verified’ genes were measured by ELISA among our in-house RA cases and controls.

Results

A total of 221 RA-associated genes were newly identified by gene-based association study, including 71‘overlapped’, 76 ‘European-specific’ and 74 ‘Asian-specific’ genes. Among them, 105 genes had significant differential expressions between RA patients and health controls at least in one dataset, especially for 20 genes including 11 ‘overlapped’ (ABCF1, FLOT1, HLA-F, IER3, TUBB, ZKSCAN4, BTN3A3, HSP90AB1, CUTA, BRD2, HLA-DMA), 5 ‘European-specific’ (PHTF1, RPS18, BAK1, TNFRSF14, SUOX) and 4 ‘Asian-specific’ (RNASET2, HFE, BTN2A2, MAPK13) genes whose differential expressions were significant at least in three datasets. The protein expressions of two selected genes FLOT1 (P value = 1.70E-02) and HLA-DMA (P value = 4.70E-02) in plasma were significantly different in our in-house samples.

Conclusion

Our study identified 221 novel RA-associated genes and especially highlighted the importance of 20 candidate genes on RA. The results addressed ethnic genetic background differences for RA susceptibility between European and Asian populations and detected a long list of overlapped or ethnic specific RA genes. The study not only greatly increases our understanding of genetic susceptibility to RA, but also provides important insights into the ethno-genetic homogeneity and heterogeneity of RA in both ethnicities.

Clinical genomic profiling identifies TYK2 mutation and overexpression in patients with neurofibromatosis type 1-associated malignant peripheral nerve sheath tumors

BACKGROUND

Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive sarcomas that arise at an estimated frequency of 8% to 13% in individuals with neurofibromatosis type 1 (NF1). Compared with their sporadic counterparts, NF1-associated MPNSTs (NF1-MPNSTs) develop in young adults, frequently recur (approximately 50% of cases), and carry a dismal prognosis. As such, most individuals affected with NF1-MPNSTs die within 5 years of diagnosis, despite surgical resection combined with radiotherapy and chemotherapy.

METHODS

Clinical genomic profiling was performed using 1000 ng of DNA from 7 cases of NF1-MPNST, and bioinformatic analyses were conducted to identify genes with actionable mutations.

RESULTS

A total of 3 women and 4 men with NF1-MPNST were identified (median age, 38 years). Nonsynonymous mutations were discovered in 4 genes (neurofibromatosis type 1 [NF1], ROS proto-oncogene 1 [ROS1], tumor protein p53 [TP53], and tyrosine kinase 2 [TYK2]), which in addition were mutated in other MPNST cases in this sample set. Consistent with their occurrence in individuals with NF1, all tumors had at least 1 mutation in the NF1 gene. Whereas TP53 gene mutations are frequently observed in other cancers, ROS1 mutations are common in melanoma (15%-35%), another neural crest-derived malignancy. In contrast, TYK2 mutations are uncommon in other malignancies (<7%). In the current series, recurrent TYK2 mutations were identified in 2 cases of NF1-MPNST (30% of cases), whereas TYK2 protein overexpression was observed in 60% of MPNST cases using an independently generated tissue microarray, regardless of NF1 status.

CONCLUSIONS

Clinical genomic analysis of the current series of NF1-MPNST cases found that TYK2 is a new gene mutated in MPNST. Future work will focus on examining the utility of TYK2 expression as a biomarker and therapeutic target for these cancers. Cancer 2017;123:1194-1201. © 2016 American Cancer Society.

Autoimmune risk variants in ERAP2 are associated with gene-expression levels in thymus

Genetic polymorphisms in the endoplasmic reticulum aminopeptidase (ERAP)1 and ERAP2 genes have been associated with several autoimmune diseases (AIDs) at a genome-wide significance level. In this study, we performed a cis expression quantitative trait locus (eQTL) screen to investigate whether seven fine-mapped AID single-nucleotide polymorphisms (SNPs) in the ERAP-region influence the gene-expression levels of ERAP1 and ERAP2 in thymus. After quality control, we identified six significant eQTLs. We further assessed the peak eQTL signals, and both genes showed highly significant and independent thymic eQTL signals (P=2.16 × 10^-15 and P=8.22 × 10^-23, respectively). Interestingly, the peak eQTL signal overlapped with the AID risk loci in ERAP2 (r²>0.94), but were distinct in ERAP1 (r²<0.4). Finally, among the SNPs showing the most significant eQTL associations with ERAP2 (P<3.4 × 10^-20), six were located within transcription factor motifs in an enhancer region in thymus. Our study therefore reveals the fine-mapped AID risk variants that act as eQTLs with ERAP2 in thymus, and highlights the potential causal regulatory variants.

Capture Hi-C identifies a novel causal gene, IL20RA, in the pan-autoimmune genetic susceptibility region 6q23

BACKGROUND

The identification of causal genes from genome-wide association studies (GWAS) is the next important step for the translation of genetic findings into biologically meaningful mechanisms of disease and potential therapeutic targets. Using novel chromatin interaction detection techniques and allele specific assays in T and B cell lines, we provide compelling evidence that redefines causal genes at the 6q23 locus, one of the most important loci that confers autoimmunity risk.

RESULTS

Although the function of disease-associated non-coding single nucleotide polymorphisms (SNPs) at 6q23 is unknown, the association is generally assigned to TNFAIP3, the closest gene. However, the DNA fragment containing the associated SNPs interacts through chromatin looping not only with TNFAIP3, but also with IL20RA, located 680 kb upstream. The risk allele of the most likely causal SNP, rs6927172, is correlated with both a higher frequency of interactions and increased expression of IL20RA, along with a stronger binding of both the NFκB transcription factor and chromatin marks characteristic of active enhancers in T-cells.

CONCLUSIONS

Our results highlight the importance of gene assignment for translating GWAS findings into biologically meaningful mechanisms of disease and potential therapeutic targets; indeed, monoclonal antibody therapy targeting IL-20 is effective in the treatment of rheumatoid arthritis and psoriasis, both with strong GWAS associations to this region.

TcoF-DB v2: update of the database of human and mouse transcription co-factors and transcription factor interactions

Transcription factors (TFs) play a pivotal role in transcriptional regulation, making them crucial for cell survival and important biological functions. For the regulation of transcription, interactions of different regulatory proteins known as transcription co-factors (TcoFs) and TFs are essential in forming necessary protein complexes. Although TcoFs themselves do not bind DNA directly, their influence on transcriptional regulation and initiation, although indirect, has been shown to be significant, with the functionality of TFs strongly influenced by the presence of TcoFs. In the TcoF-DB v2 database, we collect information on TcoFs. In this article, we describe updates and improvements implemented in TcoF-DB v2. TcoF-DB v2 provides several new features that enables exploration of the roles of TcoFs. The content of the database has significantly expanded, and is enriched with information from Gene Ontology, biological pathways, diseases and molecular signatures. TcoF-DB v2 now includes many more TFs; has substantially increased the number of human TcoFs to 958, and now includes information on mouse (418 new TcoFs). TcoF-DB v2 enables the exploration of information on TcoFs and allows investigations into their influence on transcriptional regulation in humans and mice. TcoF-DB v2 can be accessed at http://tcofdb.org/.

Rheumatoid Arthritis and Coronary Artery Disease: Genetic Analyses Do Not Support a Causal Relation

OBJECTIVE

Inflammatory diseases, specifically rheumatoid arthritis (RA), are assumed to increase the risk of coronary artery disease (CAD). More recently, multiple single-nucleotide polymorphisms (SNP) associated with RA risk were identified. If causal mechanisms affecting risks of RA and CAD are overlapping, risk alleles for RA might also increase the risk of CAD.

METHODS

Sixty-one SNP associating with RA in genome-wide significant analyses were tested for association with CAD in CARDIoGRAM (Coronary ARtery DIsease Genome wide Replication and Meta-analysis), a metaanalysis including genome-wide association data (22,233 CAD cases, 64,762 controls). In parallel, a set of SNP being associated with low-density lipoprotein cholesterol (LDL-C) was tested as a positive control.

RESULTS

Twenty-nine RA-associated SNP displayed a directionality-consistent association with CAD (OR range 1.002-1.073), whereas 32 RA-associated SNP were not associated with CAD (OR range 0.96-0.99 per RA risk-increasing allele). The proportion (48%) of directionality-consistent associated SNP equaled the proportion expected by chance (50%, p = 0.09). Of only 5 RA-associated SNP showing p values for CAD < 0.05, 4 loci (C5orf30, IL-6R, PTPN22, and RAD51B) showed directionality-consistent effects on CAD, and 1 (rs10774624, locus SH2B3) reached study-wide significance (p = 7.29E-06). By contrast, and as a proof of concept, 46 (74%) out of 62 LDL-C-associated SNP displayed a directionality-consistent association with CAD, a proportion that was significantly different from 50% (p = 5.9E-05).

CONCLUSION

We found no evidence that RA-associated SNP as a group are associated with CAD. Even though we were not able to study potential effects of all genetic variants individually, shared nongenetic factors may more plausibly explain the observed coincidence of the 2 conditions.

Metabolic profile of Guge Fengtong tablet in rat urine and bile after oral administration, using high-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry

In the present study, we developed and validated a rapid analytical method using high performance liquid chromatography coupled with quadrupole time of flight mass spectrometry (HPLC-Q-TOF/MS) to investigate the metabolic profile of Guge Fengtong tablet (GGFTT), a traditional Chinese medicine. The urine and bile samples were collected with 24 h after oral administration of GGFTT. A total of 34 compounds, including 11 parent compounds and 23 metabolites were unambiguously or tentatively identified. Our results indicated that glucuronidation, oxidation and methylation were the major metabolic pathways of the constituents in GGFTT. In addition, the results of this work also demonstrated the feasibility of HPLC-ESI-Q-TOF/MS for reliable characterization of the in vivo metabolites from herbal preparations.

A Comparative Analysis of the Peptide Repertoires of HLA-DR Molecules Differentially Associated With Rheumatoid Arthritis

OBJECTIVE

To evaluate similarity of the peptide repertoires bound to HLA-DR molecules that are differentially associated with rheumatoid arthritis (RA), and to define structural features of the shared peptides.

METHODS

Peptide pools bound to HLA-DRB1*01:01, HLA-DRB1*04:01, and HLA-DRB1*10:01 (RA associated) and those bound to HLA-DRB1*15:01 (non-RA-associated) were purified and analyzed by liquid chromatography (LC) matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MS) and LC-ion-trap MS. Peptide pools from each allotype were compared in terms of size, protein origin, composition, and affinity (both theoretical and experimental with some peptides). Finally, 1 peptide sequenced from DR1, DR4, and DR10, but not from DR15, was modeled in complex with all 4 HLA-DRB1 molecules and HLA-DRB5*01:01.

RESULTS

A total of 6,309 masses and 962 unique peptide sequences were compared. DR10 shared 29 peptides with DR1, 9 with DR4, and 1 with DR15; DR1 shared 6 peptides with DR4 and 9 with DR15; and DR4 and DR15 shared 4 peptides. The direct identification of peptide ligands indicated that DR1 and DR10 were the most similar molecules regarding the peptides that they could share. The peptides common to these molecules contained a high proportion of Leu at P4 and basic residues at P8 binding core positions.

CONCLUSION

The degree of overlap between peptide repertoires associated with different HLA-DR molecules is low. The repertoires associated with DR1 and DR10 have the highest similarity among the molecules analyzed (∼10% overlap). Among the peptides shared between DR1 and DR10, a high proportion contained Leu(4) and basic residues at the P8 position of the binding core.

Genetic analysis of emerging risk factors in coronary artery disease

BACKGROUND AND AIMS

Type 2 diabetes (T2D), low-density lipoprotein-cholesterol (LDL-c), body mass index (BMI), blood pressure and smoking are established risk factors that play a causal role in coronary artery disease (CAD). Numerous common genetic variants associating with these and other risk factors have been identified, but their association with CAD has not been comprehensively examined in a single study. Our goal was to comprehensively evaluate the associations of established and emerging risk factors with CAD using genetic variants identified from Genome-wide Association Studies (GWAS).

METHODS

We tested the effect of 60 traditional and putative risk factors with CAD, using summary statistics obtained in GWAS. We approximated the regression of a response variable onto an additive multi-SNP genetic risk score in the Coronary Artery DIsease Genomewide Replication And Meta-analysis (CARDIoGRAM) consortium dataset weighted by the effect of the SNP on the risk factors.

RESULTS

The strongest association with risk of CAD was for LDL-c SNPs (p = 3.96E-34). For non-established CAD risk factors, we found significant CAD associations for coronary artery calcification (CAC), Lp(a), LP-PLA2 activity, plaque, vWF and FVIII. In an attempt to identify independent associations between risk factors and CAD, only SNPs with an effect on the target trait were included. This identified CAD associations for Lp(a)(p = 1.77E-21), LDL-c (p = 4.16E-06), triglycerides (TG) (p = 1.94E-05), height (p = 2.06E-05), CAC (p = 3.13E-23) and carotid plaque (p = 2.08E-05).

CONCLUSIONS

We identified SNPs associated with the emerging risk factors Lp(a), TG, plaque, height and CAC to be independently associated with risk of CAD. This provides further support for-ongoing clinical trials of Lp(a) and TG, and suggests that CAC and plaque could be used as surrogate markers for CAD in clinical trials.

Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells

Background

The observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells – “Immunoseq”.

Results

We performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations.

Conclusions

We show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0220-7) contains supplementary material, which is available to authorized users.

Genetic Analysis with the Immunochip Platform in Behçet Disease. Identification of Residues Associated in the HLA Class I Region and New Susceptibility Loci

Behcet's disease (BD) is an immuno-mediated vasculitis in which knowledge of its etiology and genetic basis is limited. To improve the current knowledge, a genetic analysis performed with the Immunochip platform was carried out in a population from Spain. A discovery cohort comprising 278 BD cases and 1,517 unaffected controls were genotyped using the Immunochip platform. The validation step was performed on an independent replication cohort composed of 130 BD cases and 600 additional controls. The strongest association signals were observed in the HLA class I region, being HLA-B*51 the highest peak (overall P = 6.82E-32, OR = 3.82). A step-wise conditional logistic regression with classical alleles identified HLA-B*57 and HLA-A*03 as additional independent markers. The amino acid model that best explained the association, includes the position 97 of the HLA-B molecule and the position 66 of the HLA-A. Among the non-HLA loci, the most significant in the discovery analysis were: IL23R (rs10889664: P = 3.81E-12, OR = 2.00), the JRKL/CNTN5 region (rs2848479: P = 5.00E-08, OR = 1.68) and IL12A (rs1874886: P = 6.67E-08, OR = 1.72), which were confirmed in the validation phase (JRKL/CNTN5 rs2848479: P = 3.29E-10, OR = 1.66; IL12A rs1874886: P = 1.62E-08, OR = 1.61). Our results confirm HLA-B*51 as a primary-association marker in predisposition to BD and suggest additional independent signals within the class I region, specifically in the genes HLA-A and HLA-B. Regarding the non-HLA genes, in addition to IL-23R, previously reported in our population; IL12A, described in other populations, was found to be a BD susceptibility factor also in Spaniards; finally, a new associated locus was found in the JRKL/CNTN5 region.

Immune recognition of citrullinated epitopes

Conversion of arginine into citrulline is a post-translational modification that is observed in normal physiological processes. However, abnormal citrullination can provoke autoimmunity by generating altered self-epitopes that are specifically targeted by autoantibodies and T cells. In this review we discuss the recognition of citrullinated antigens in human autoimmune diseases and the role that this modification plays in increasing antigenic diversity and circumventing tolerance mechanisms. Early published work demonstrated that citrullinated proteins are specifically targeted by autoantibodies in rheumatoid arthritis and that citrullinated peptides are more readily presented to T cells by arthritis-susceptible HLA class II 'shared epitope' proteins. Emerging data support the relevance of citrullinated epitopes in other autoimmune diseases, including type 1 diabetes and multiple sclerosis, whose susceptible HLA haplotypes also preferentially present citrullinated peptides. In these settings, autoimmune patients have been shown to have elevated responses to citrullinated epitopes derived from tissue-specific antigens. Contrasting evidence implicates autophagy or perforin and complement-mediated membrane attack as inducers of ectopic citrullination. In either case, the peptidyl deiminases responsible for citrullination are activated in response to inflammation or insult, providing a mechanistic link between this post-translational modification and interactions with the environment and infection. As such, it is likely that immune recognition of citrullinated epitopes also plays a role in pathogen clearance. Indeed, our recent data suggest that responses to citrullinated peptides facilitate recognition of novel influenza strains. Therefore, increased understanding of responses to citrullinated epitopes may provide important insights about the initiation of autoimmunity and recognition of heterologous viruses.

Rheumatoid arthritis onset in postmenopausal women: Does the ACPA seropositive subset result from genetic effects, estrogen deficiency, skewed profile of CD4(+) T-cells, and their interactions?

Rheumatoid arthritis (RA) incidence displays a differentiated age-dependent female-to-male ratio in which women outnumber men. Evidence that the peak incidence of RA in women coincides with menopause age, suggests a potential estrogenic role to disease etiology. Estrogens exert physiologically both stimulatory and inhibitory effects on the immune system. Epidemiologic and animal model studies with estrogen deprivation or supplementation suggested estrogens as to play, mainly, a protective role in RA immunopathology. In this review, we propose that some yet unidentified disturbances associated with estrogen circulating levels, differentiated by the menopausal status, play a major role in women's RA susceptibility. We focus on the interaction between estrogen deprivation and genetic risk alleles for anti-citrullinated protein antibodies (ACPA) seropositive RA, as a major driving force for increased immune reactivity and RA susceptibility, in postmenopausal women. This opens up new fields for research concerning the association among different irregular estrogenic conditions, the cytokine milieu, and age/menopausal status bias in RA.

Influence of TYK2 in systemic sclerosis susceptibility: a new locus in the IL-12 pathway

OBJECTIVES

TYK2 is a common genetic risk factor for several autoimmune diseases. This gene encodes a protein kinase involved in interleukin 12 (IL-12) pathway, which is a well-known player in the pathogenesis of systemic sclerosis (SSc). Therefore, we aimed to assess the possible role of this locus in SSc.

METHODS

This study comprised a total of 7103 patients with SSc and 12 220 healthy controls of European ancestry from Spain, USA, Germany, the Netherlands, Italy and the UK. Four TYK2 single-nucleotide polymorphisms (V362F (rs2304256), P1104A (rs34536443), I684S (rs12720356) and A928V (rs35018800)) were selected for follow-up based on the results of an Immunochip screening phase of the locus. Association and dependence analyses were performed by the means of logistic regression and conditional logistic regression. Meta-analyses were performed using the inverse variance method.

RESULTS

Genome-wide significance level was reached for TYK2 V362F common variant in our pooled analysis (p=3.08×10(-13), OR=0.83), while the association of P1104A, A928V and I684S rare and low-frequency missense variants remained significant with nominal signals (p=2.28×10(-3), OR=0.80; p=1.27×10(-3), OR=0.59; p=2.63×10(-5), OR=0.83, respectively). Interestingly, dependence and allelic combination analyses showed that the strong association observed for V362F with SSc, corresponded to a synthetic association dependent on the effect of the three previously mentioned TYK2 missense variants.

CONCLUSIONS

We report for the first time the association of TYK2 with SSc and reinforce the relevance of the IL-12 pathway in SSc pathophysiology.

Dense genotyping of immune-related loci in idiopathic inflammatory myopathies confirms HLA alleles as the strongest genetic risk factor and suggests different genetic background for major clinical subgroups

OBJECTIVES

The idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of rare autoimmune diseases characterised by muscle weakness and extramuscular manifestations such as skin rashes and interstitial lung disease. We genotyped 2566 IIM cases of Caucasian descent using the Immunochip; a custom array covering 186 established autoimmune susceptibility loci. The cohort was predominantly comprised of patients with dermatomyositis (DM, n=879), juvenile DM (JDM, n=481), polymyositis (PM, n=931) and inclusion body myositis (n=252) collected from 14 countries through the Myositis Genetics Consortium.

RESULTS

The human leucocyte antigen (HLA) and PTPN22 regions reached genome-wide significance (p<5×10(-8)). Nine regions were associated at a significance level of p<2.25×10(-5), including UBE2L3, CD28 and TRAF6, with evidence of independent effects within STAT4. Analysis of clinical subgroups revealed distinct differences between PM, and DM and JDM. PTPN22 was associated at genome-wide significance with PM, but not DM and JDM, suggesting this effect is driven by PM. Additional suggestive associations including IL18R1 and RGS1 in PM and GSDMB in DM were identified. HLA imputation confirmed that alleles HLA-DRB1*03:01 and HLA-B*08:01 of the 8.1 ancestral haplotype (8.1AH) are most strongly associated with IIM, and provides evidence that amino acids within the HLA, such as HLA-DQB1 position 57 in DM, may explain part of the risk in this locus. Associations with alleles outside the 8.1AH reveal differences between PM, DM and JDM.

CONCLUSIONS

This work represents the largest IIM genetic study to date, reveals new insights into the genetic architecture of these rare diseases and suggests different predominating pathophysiology in different clinical subgroups.

Investigation of the genetic overlap between rheumatoid arthritis and psoriatic arthritis in a Greek population

OBJECTIVES

Several rheumatoid arthritis (RA) susceptibility loci have also been found to be associated with psoriatic arthritis (PsA), demonstrating that there is a degree of genetic overlap between various autoimmune diseases. We sought to investigate whether single nucleotide polymorphisms (SNPs) mapping to previously reported RA and/or PsA susceptibility loci, including PLCL2, CCL21, REL, STAT4, CD226, PTPN22, and TYK2, are associated with risk for the two diseases in a genetically homogeneous Greek population.

METHOD

This study included 392 RA patients, 126 PsA patients, and 521 healthy age- and sex-matched controls from Greece. Genotyping of the SNPs was performed with Taqman primer/probe sets. Bioinformatic analysis was performed using BlastP, PyMOL, and Maestro and Desmond.

RESULTS

A significant association was detected between the GC genotype of rs34536443 (TYK2) in both the PsA and RA cohorts. The C allele of this SNP was associated with PsA only. Evidence for association with PsA was also found for the GG genotype and G allele of the rs10181656 SNP of STAT4. The TC genotype of the rs763361 SNP of CD226 was associated with PsA only.

CONCLUSIONS

Genetic overlap between PsA and RA was detected for the rs34536443 SNP of the TYK2 gene within a Greek population. An association of STAT4 (rs10181656) with PsA was confirmed whereas CD226 (rs763361) was associated with PsA but not with RA, in contrast to previous reports. The different findings of this study compared to previous ones highlights the importance of comparative studies that include various ethnic or racial populations.

RBP-J-Regulated miR-182 Promotes TNF-α-Induced Osteoclastogenesis

Increased osteoclastogenesis is responsible for osteolysis, which is a severe consequence of inflammatory diseases associated with bone destruction, such as rheumatoid arthritis and periodontitis. The mechanisms that limit osteoclastogenesis under inflammatory conditions are largely unknown. We previously identified transcription factor RBP-J as a key negative regulator that restrains TNF-α-induced osteoclastogenesis and inflammatory bone resorption. In this study, we tested whether RBP-J suppresses inflammatory osteoclastogenesis by regulating the expression of microRNAs (miRNAs) important for this process. Using high-throughput sequencing of miRNAs, we obtained the first, to our knowledge, genome-wide profile of miRNA expression induced by TNF-α in mouse bone marrow-derived macrophages/osteoclast precursors during inflammatory osteoclastogenesis. Furthermore, we identified miR-182 as a novel miRNA that promotes inflammatory osteoclastogenesis driven by TNF-α and whose expression is suppressed by RBP-J. Downregulation of miR-182 dramatically suppressed the enhanced osteoclastogenesis program induced by TNF-α in RBP-J-deficient cells. Complementary loss- and gain-of-function approaches showed that miR-182 is a positive regulator of osteoclastogenic transcription factors NFATc1 and B lymphocyte-induced maturation protein-1. Moreover, we identified that direct miR-182 targets, Foxo3 and Maml1, play important inhibitory roles in TNF-α-mediated osteoclastogenesis. Thus, RBP-J-regulated miR-182 promotes TNF-α-induced osteoclastogenesis via inhibition of Foxo3 and Maml1. Suppression of miR-182 by RBP-J serves as an important mechanism that restrains TNF-α-induced osteoclastogenesis. Our results provide a novel miRNA-mediated mechanism by which RBP-J inhibits osteoclastogenesis and suggest that targeting of the newly described RBP-J-miR-182-Foxo3/Maml1 axis may represent an effective therapeutic approach to suppress inflammatory osteoclastogenesis and bone resorption.