Association with HLA-DRß1 position 37 distinguishes juvenile dermatomyositis from adult-onset myositis.

Juvenile dermatomyositis (JDM) is a rare, severe autoimmune disease and the most common idiopathic inflammatory myopathy of children. JDM and adult-onset dermatomyositis (DM) have similar clinical, biological and serological features, although these features differ in prevalence between childhood-onset and adult-onset disease, suggesting that age of disease onset may influence pathogenesis. Therefore, a JDM-focused genetic analysis was performed using the largest collection of JDM samples to date. Caucasian JDM samples (n = 952) obtained via international collaboration were genotyped using the Illumina HumanCoreExome chip. Additional non-assayed human leukocyte antigen (HLA) loci and genome-wide single-nucleotide polymorphisms (SNPs) were imputed. HLA-DRB1*03:01 was confirmed as the classical HLA allele most strongly associated with JDM [odds ratio (OR) 1.66; 95% confidence interval (CI) 1.46, 1.89; P = 1.4 × 10-14], with an independent association at HLA-C*02:02 (OR = 1.74; 95% CI 1.42, 2.13, P = 7.13 × 10-8). Analyses of amino acid positions within HLA-DRB1 indicated that the strongest association was at position 37 (omnibus P = 3.3 × 10-19), with suggestive evidence this association was independent of position 74 (omnibus P = 5.1 × 10-5), the position most strongly associated with adult-onset DM. Conditional analyses also suggested that the association at position 37 of HLA-DRB1 was independent of some alleles of the Caucasian HLA 8.1 ancestral haplotype (AH8.1) such as HLA-DQB1*02:01 (OR = 1.62; 95% CI 1.36, 1.93; P = 8.70 × 10-8), but not HLA-DRB1*03:01 (OR = 1.49; 95% CR 1.24, 1.80; P = 2.24 × 10-5). No associations outside the HLA region were identified. Our findings confirm previous associations with AH8.1 and HLA-DRB1*03:01, HLA-C*02:02 and identify a novel association with amino acid position 37 within HLA-DRB1, which may distinguish JDM from adult DM.

Pre-treatment calprotectin (MRP8/14) provides no added value to testing CRP alone in terms of predicting response to TNF inhibitors in rheumatoid arthritis in a post hoc analysis.

OBJECTIVES: The inflammatory protein calprotectin (MRP8/14) has been identified as a promising biomarker of treatment response in rheumatoid arthritis (RA). Our aim was to test MRP8/14 as a biomarker of response to tumour necrosis factor (TNF)-inhibitors in the largest RA cohort to date and to compare with C-reactive protein (CRP). METHODS: Serum MRP8/14 was measured in 470 patients with RA about to commence treatment with adalimumab (n=196) or etanercept (n=274). Additionally, MRP8/14 was measured in the 3-month sera of 179 adalimumab-treated patients. Response was determined using European League against Rheumatism (EULAR) response criteria calculated using the traditional 4-component (4C) DAS28-CRP and alternate validated versions using 3-component (3C) and 2-component (2C), clinical disease activity index (CDAI) improvement criteria and change in individual outcome measures. Logistic/linear regression models were fitted for response outcome. RESULTS: In the 3C and 2C models, patients with RA were 1.92 (CI: 1.04 to 3.54) and 2.03 (CI: 1.09 to 3.78) times more likely to be classified as EULAR responders if they had high (75th quartile) pre-treatment levels of MRP8/14 compared with low (25th quartile). No significant associations were observed for the 4C model. When only using CRP as a predictor, in the 3C and 2C analyses, patients above the 75th quartile were 3.79 (CI: 1.81 to 7.93) and 3.58 (CI: 1.74 to 7.35) times more likely to be EULAR responders and addition of MRP8/14 did not significantly improve model fit (p values=0.62 and 0.80, respectively). No significant associations were observed in the 4C analysis. Exclusion of CRP from the outcome measure (CDAI) did not result in any significant associations with MRP8/14 (OR 1.00 (CI: 0.99 to 1.01), suggesting that the associations were due to the correlation with CRP and that there is no additional utility of MRP8/14 beyond use of CRP in patients with RA starting TNFi therapy. CONCLUSION: Beyond correlation with CRP, we found no evidence to suggest that MRP8/14 explains additional variability in response to TNFi in patients with RA over and above CRP alone.

Take your PICS: moving from GWAS to immune function.

Many of the hits identified through genome-wide association studies are located outside protein-coding regions, making it difficult to define mechanism. In Nature, Farh et al., (2014) describe an approach to identify causal variants in autoimmune disease as first step to assigning function.

No evidence that genetic predictors of susceptibility predict changes in core outcomes in JIA.

OBJECTIVES: The clinical progression of JIA is unpredictable. Knowing who will develop severe disease could facilitate rapid intensification of therapies. We use genetic variants conferring susceptibility to JIA to predict disease outcome measures. METHODS: A total of 713 JIA patients with genotype data and core outcome variables (COVs) at diagnosis (baseline) and 1 year follow-up were identified from the Childhood Arthritis Prospective Study (CAPS). A weighted genetic risk score (GRS) was generated, including all single nucleotide polymorphisms (SNPs) previously associated with JIA susceptibility (P-value < 5×10-08). We used multivariable linear regression to test the GRS for association with COVS (limited joint count, active joint count, physician global assessment, parent/patient general evaluation, childhood HAQ and ESR) at baseline and change in COVS from baseline to 1 year, adjusting for baseline COV and International League of Associations of Rheumatology (ILAR) category. The GRS was split into quintiles to identify high (quintile 5) and low (quintile 1) risk groups. RESULTS: Patients in the high-risk group for the GRS had a younger age at presentation (median low risk 7.79, median high risk 3.51). No association was observed between the GRS and any outcome measures at 1 year follow-up or baseline. CONCLUSION: For the first time we have used all known JIA genetic susceptibility loci (P=<5×10-08) in a GRS to predict changes in disease outcome measured over time. Genetic susceptibility variants are poor predictors of changes in core outcome measures, it is likely that genetic factors predicting disease outcome are independent to those predicting susceptibility. The next step will be to conduct a genome-wide association analysis of JIA outcome.

Germline intergenic duplications at Xq26.1 underlie Bazex-Dupré-Christol basal cell carcinoma susceptibility syndrome.

BACKGROUND: Bazex-Dupré-Christol syndrome (BDCS; MIM301845) is a rare X-linked dominant genodermatosis characterized by follicular atrophoderma, congenital hypotrichosis and multiple basal cell carcinomas (BCCs). Previous studies have linked BDCS to an 11·4-Mb interval on chromosome Xq25-q27.1. However, the genetic mechanism of BDCS remains an open question. OBJECTIVES: To investigate the genetic aetiology and molecular mechanisms underlying BDCS. METHODS: We ascertained multiple individuals from eight unrelated families affected with BDCS (F1-F8). Whole-exome (F1 and F2) and genome sequencing (F3) were performed to identify putative disease-causing variants within the linkage region. Array comparative genomic hybridization and quantitative polymerase chain reaction (PCR) were used to explore copy number variations, followed by long-range gap PCR and Sanger sequencing to amplify the duplication junctions and to define the head-tail junctions. Hi-C was performed on dermal fibroblasts from two affected individuals with BDCS and one control. Public datasets and tools were used to identify regulatory elements and transcription factor binding sites within the minimal duplicated region. Immunofluorescence was performed in hair follicles, BCCs and trichoepitheliomas from patients with BDCS and sporadic BCCs. The ACTRT1 variant c.547dup (p.Met183Asnfs*17), previously proposed to cause BDCS, was evaluated with t allele frequency calculator. RESULTS: In eight families with BDCS, we identified overlapping 18-135-kb duplications (six inherited and two de novo) at Xq26.1, flanked by ARHGAP36 and IGSF1. Hi-C showed that the duplications did not affect the topologically associated domain, but may alter the interactions between flanking genes and putative enhancers located in the minimal duplicated region. We detected ARHGAP36 expression near the control hair follicular stem cell compartment, and found increased ARHGAP36 levels in hair follicles in telogen, in BCCs and in trichoepitheliomas from patients with BDCS. ARHGAP36 was also detected in sporadic BCCs from individuals without BDCS. Our modelling showed the predicted maximum tolerated minor allele frequency of ACTRT1 variants in control populations to be orders of magnitude higher than expected for a high-penetrant ultra-rare disorder, suggesting loss of function of ACTRT1 variants to be an unlikely cause for BDCS. CONCLUSIONS: Noncoding Xq26.1 duplications cause BDCS. The BDCS duplications most likely lead to dysregulation of ARHGAP36. ARHGAP36 is a potential therapeutic target for both inherited and sporadic BCCs. What is already known about this topic? Bazex-Dupré-Christol syndrome (BDCS) is a rare X-linked basal cell carcinoma susceptibility syndrome linked to an 11·4-Mb interval on chromosome Xq25-q27.1. Loss-of-function variants in ACTRT1 and its regulatory elements were suggested to cause BDCS. What does this study add? BDCS is caused by small tandem noncoding intergenic duplications at chromosome Xq26.1. The Xq26.1 BDCS duplications likely dysregulate ARHGAP36, the flanking centromeric gene. ACTRT1 loss-of-function variants are unlikely to cause BDCS. What is the translational message? This study provides the basis for accurate genetic testing for BDCS, which will aid precise diagnosis and appropriate surveillance and clinical management. ARHGAP36 may be a novel therapeutic target for all forms of sporadic basal cell carcinomas.

Combined genetic analysis of juvenile idiopathic arthritis clinical subtypes identifies novel risk loci, target genes and key regulatory mechanisms.

OBJECTIVES: Juvenile idiopathic arthritis (JIA) is the most prevalent form of juvenile rheumatic disease. Our understanding of the genetic risk factors for this disease is limited due to low disease prevalence and extensive clinical heterogeneity. The objective of this research is to identify novel JIA susceptibility variants and link these variants to target genes, which is essential to facilitate the translation of genetic discoveries to clinical benefit. METHODS: We performed a genome-wide association study (GWAS) in 3305 patients and 9196 healthy controls, and used a Bayesian model selection approach to systematically investigate specificity and sharing of associated loci across JIA clinical subtypes. Suggestive signals were followed-up for meta-analysis with a previous GWAS (2751 cases/15 886 controls). We tested for enrichment of association signals in a broad range of functional annotations, and integrated statistical fine-mapping and experimental data to identify target genes. RESULTS: Our analysis provides evidence to support joint analysis of all JIA subtypes with the identification of five novel significant loci. Fine-mapping nominated causal single nucleotide polymorphisms with posterior inclusion probabilities ≥50% in five JIA loci. Enrichment analysis identified RELA and EBF1 as key transcription factors contributing to disease risk. Our integrative approach provided compelling evidence to prioritise target genes at six loci, highlighting mechanistic insights for the disease biology and IL6ST as a potential drug target. CONCLUSIONS: In a large JIA GWAS, we identify five novel risk loci and describe potential function of JIA association signals that will be informative for future experimental works and therapeutic strategies.

Lymphocyte DNA methylation mediates genetic risk at shared immune-mediated disease loci.

BACKGROUND: Defining regulatory mechanisms through which noncoding risk variants influence the cell-mediated pathogenesis of immune-mediated disease (IMD) has emerged as a priority in the post-genome-wide association study era. OBJECTIVES: With a focus on rheumatoid arthritis, we sought new insight into genetic mechanisms of adaptive immune dysregulation to help prioritize molecular pathways for targeting in this and related immune pathologies. METHODS: Whole-genome methylation and transcriptional data from isolated CD4+ T cells and B cells of more than 100 genotyped and phenotyped patients with inflammatory arthritis, all of whom were naive to immunomodulatory treatments, were obtained. Analysis integrated these comprehensive data with genome-wide association study findings across IMDs and other publicly available resources. RESULTS: We provide strong evidence that disease-associated DNA variants regulate cis-CpG methylation in CD4+ T and/or B cells at 37% RA loci. Using paired, cell-specific transcriptomic data and causal inference testing, we identify examples where site-specific DNA methylation in turn mediates gene expression, including FCRL3 in both cell types and ORMDL3/GSDMB, IL6ST/ANKRD55, and JAZF1 in CD4+ T cells. A number of genes regulated in this way highlight mechanisms common to RA and other IMDs including multiple sclerosis and asthma, in turn distinguishing them from osteoarthritis, a primarily degenerative disease. Finally, we corroborate the observed effects experimentally. CONCLUSIONS: Our observations highlight important mechanisms of genetic risk in RA and the wider context of immune dysregulation. They confirm the utility of DNA methylation profiling as a tool for causal gene prioritization and, potentially, therapeutic targeting in complex IMD.

Genetics of rheumatoid arthritis: 2018 status.

Study of the genetics of rheumatoid arthritis (RA) began about four decades ago with the discovery of HLA-DRB1 Since the beginning of this century, a number of non-HLA risk loci have been identified through genome-wide association studies (GWAS). We now know that over 100 loci are associated with RA risk. Because genetic information implies a clear causal relationship to the disease, research into the pathogenesis of RA should be promoted. However, only 20% of GWAS loci contain coding variants, with the remaining variants occurring in non-coding regions, and therefore, the majority of causal genes and causal variants remain to be identified. The use of epigenetic studies, high-resolution mapping of open chromatin, chromosomal conformation technologies and other approaches could identify many of the missing links between genetic risk variants and causal genetic components, thus expanding our understanding of RA genetics.

A large-scale genetic analysis reveals a strong contribution of the HLA class II region to giant cell arteritis susceptibility.

We conducted a large-scale genetic analysis on giant cell arteritis (GCA), a polygenic immune-mediated vasculitis. A case-control cohort, comprising 1,651 case subjects with GCA and 15,306 unrelated control subjects from six different countries of European ancestry, was genotyped by the Immunochip array. We also imputed HLA data with a previously validated imputation method to perform a more comprehensive analysis of this genomic region. The strongest association signals were observed in the HLA region, with rs477515 representing the highest peak (p = 4.05 × 10(-40), OR = 1.73). A multivariate model including class II amino acids of HLA-DRß1 and HLA-DQα1 and one class I amino acid of HLA-B explained most of the HLA association with GCA, consistent with previously reported associations of classical HLA alleles like HLA-DRB1(∗)04. An omnibus test on polymorphic amino acid positions highlighted DRß1 13 (p = 4.08 × 10(-43)) and HLA-DQα1 47 (p = 4.02 × 10(-46)), 56, and 76 (both p = 1.84 × 10(-45)) as relevant positions for disease susceptibility. Outside the HLA region, the most significant loci included PTPN22 (rs2476601, p = 1.73 × 10(-6), OR = 1.38), LRRC32 (rs10160518, p = 4.39 × 10(-6), OR = 1.20), and REL (rs115674477, p = 1.10 × 10(-5), OR = 1.63). Our study provides evidence of a strong contribution of HLA class I and II molecules to susceptibility to GCA. In the non-HLA region, we confirmed a key role for the functional PTPN22 rs2476601 variant and proposed other putative risk loci for GCA involved in Th1, Th17, and Treg cell function.

Optimisation of methods for bacterial skin microbiome investigation: primer selection and comparison of the 454 versus MiSeq platform.

BACKGROUND: The composition of the skin microbiome is predicted to play a role in the development of conditions such as atopic eczema and psoriasis. 16S rRNA gene sequencing allows the investigation of bacterial microbiota. A significant challenge in this field is development of cost effective high throughput methodologies for the robust interrogation of the skin microbiota, where biomass is low. Here we describe validation of methodologies for 16S rRNA (ribosomal ribonucleic acid) gene sequencing from the skin microbiome, using the Illumina MiSeq platform, the selection of primer to amplify regions for sequencing and we compare results with the current standard protocols.. METHODS: DNA was obtained from two low density mock communities of 11 diverse bacterial strains (with and without human DNA supplementation) and from swabs taken from the skin of healthy volunteers. This was amplified using primer pairs covering hypervariable regions of the 16S rRNA gene: primers 63F and 519R (V1-V3); and 347F and 803R (V3-V4). The resultant libraries were indexed for the MiSeq and Roche454 and sequenced. Both data sets were denoised, cleaned of chimeras and analysed using QIIME. RESULTS: There was no significant difference in the diversity indices at the phylum and the genus level observed between the platforms. The capture of diversity using the low density mock community samples demonstrated that the primer pair spanning the V3-V4 hypervariable region had better capture when compared to the primer pair for the V1-V3 region and was robust to spiking with human DNA. The pilot data generated using the V3-V4 region from the skin of healthy volunteers was consistent with these results, even at the genus level (Staphylococcus, Propionibacterium, Corynebacterium, Paracoccus, Micrococcus, Enhydrobacter and Deinococcus identified at similar abundances on both platforms). CONCLUSIONS: The results suggest that the bacterial community diversity captured using the V3-V4 16S rRNA hypervariable region from sequencing using the MiSeq platform is comparable to the Roche454 GS Junior platform. These findings provide evidence that the optimised method can be used in human clinical samples of low bacterial biomass such as the investigation of the skin microbiota.

The predictive value of serum S100A9 and response to etanercept is not confirmed in a large UK rheumatoid arthritis cohort.

Objective: The aim was to correlate protein concentrations of S100A9 in pretreatment serum samples with response to the tumour-necrosis factor (TNF) inhibitor drugs etanercept in a large UK replication cohort. Methods: Pretreatment serum samples from patients with RA (n = 236) about to commence treatment with etanercept had S100A9 serum concentration measured using an ELISA. Following the experimental procedure, S100A9 concentrations were analysed with respect to EULAR response. Results: No evidence of association between S100A9 concentration and EULAR response to the TNF-inhibitor biologic drug etanercept was observed following multinomial logistic regression analysis (non-responder vs moderate responder, P = 0.957; and non-responder vs good responder, P = 0.316). Furthermore, no significant associations were observed when correlating pretreatment S100A9 concentrations with clinical parameters of disease activity (P > 0.05). Conclusion: In the largest replication cohort conducted to date, no evidence for association was observed to support the use of S100A9 as a clinical biomarker predictive of response to the TNF-inhibitor biologic drug etanercept.

The skin microbiome in psoriatic arthritis: methodology development and pilot data.

BACKGROUND: Skin microbiota are likely to be important in the development of conditions such as psoriatic arthritis. Profiling the bacterial community in the psosriatic plaques will contribute to our understanding of the role of the skin microbiome in these conditions. The aim of this work was to determine the optimum study design for work on the skin microbiome with use of the MiSeq platform. The objectives were to compare data generated from two platforms for two primer pairs in a low density mock bacterial community. METHODS: DNA was obtained from two low density mock communities of 11 diverse bacterial strains (with and without human DNA supplementation) and from swabs taken from the skin of four healthy volunteers. The DNA was amplified with primer pairs covering hypervariable regions of the 16S rRNA gene: primers 63F and 519R (V1-V3), and 347F and 803R (V3-V4). The resultant libraries were indexed for the MiSeq and Roche454 platforms and sequenced. Both datasets were de-noised, cleaned of chimeras, and analysed by use of QIIME software (version 1.8.0). FINDINGS: No significant difference in the diversity indices at the phylum and the genus level between the platforms was seen. Comparison of the diversity indices for the mock community data for the two primer pairs demonstrated that the V3-V4 hypervariable region had significantly better capture of bacterial diversity than did the V1-V3 region. Amplification with the same primer pairs showed strong concordance within each platform (98·9-99·8%), with negligible effect of spiked human DNA contamination. Comparison at the family level classification between samples processed on the MiSeq and Roche454 platforms using the V3-V4 hypervariable region also showed a high level of concordance (87%), although less so for the V1-V3 primers (10%). The pilot data from healthy volunteers were similar. INTERPRETATION: Results obtained from the V3-V4 16S rRNA hypervariable region, sequencing on the MiSeq and Roche454 platforms, were concordant between replicates, and between each other. These findings suggest that the MiSeq platform, and these primers, is a comparable method for determining skin microbiota to the widely used Roche454 methodology. FUNDING: NIHR Manchester Musculoskeletal Biomedical Research Unit.

Predicting the risk of rheumatoid arthritis and its age of onset through modelling genetic risk variants with smoking.

The improved characterisation of risk factors for rheumatoid arthritis (RA) suggests they could be combined to identify individuals at increased disease risks in whom preventive strategies may be evaluated. We aimed to develop an RA prediction model capable of generating clinically relevant predictive data and to determine if it better predicted younger onset RA (YORA). Our novel modelling approach combined odds ratios for 15 four-digit/10 two-digit HLA-DRB1 alleles, 31 single nucleotide polymorphisms (SNPs) and ever-smoking status in males to determine risk using computer simulation and confidence interval based risk categorisation. Only males were evaluated in our models incorporating smoking as ever-smoking is a significant risk factor for RA in men but not women. We developed multiple models to evaluate each risk factor's impact on prediction. Each model's ability to discriminate anti-citrullinated protein antibody (ACPA)-positive RA from controls was evaluated in two cohorts: Wellcome Trust Case Control Consortium (WTCCC: 1,516 cases; 1,647 controls); UK RA Genetics Group Consortium (UKRAGG: 2,623 cases; 1,500 controls). HLA and smoking provided strongest prediction with good discrimination evidenced by an HLA-smoking model area under the curve (AUC) value of 0.813 in both WTCCC and UKRAGG. SNPs provided minimal prediction (AUC 0.660 WTCCC/0.617 UKRAGG). Whilst high individual risks were identified, with some cases having estimated lifetime risks of 86%, only a minority overall had substantially increased odds for RA. High risks from the HLA model were associated with YORA (P<0.0001); ever-smoking associated with older onset disease. This latter finding suggests smoking's impact on RA risk manifests later in life. Our modelling demonstrates that combining risk factors provides clinically informative RA prediction; additionally HLA and smoking status can be used to predict the risk of younger and older onset RA, respectively.

The bacterial skin microbiome in psoriatic arthritis, an unexplored link in pathogenesis: challenges and opportunities offered by recent technological advances.

The resident microbial community, harboured by humans in sites such as the skin and gastrointestinal tract, is enormous, representing a candidate environmental factor affecting susceptibility to complex diseases, where both genetic and environmental risk factors are important. The potential of microorganisms to influence the human immune system is considerable, given their ubiquity. The impact of the host-gene-microbe interaction on the maintenance of health and the development of disease has not yet been assessed robustly in chronic inflammatory conditions. PsA represents a model inflammatory disease to explore the role of the microbiome because skin involvement and overlap with IBD implicates both the skin and gastrointestinal tract as sources of microbial triggers for PsA. In parallel with genetic studies, characterization of the host microbiota may benefit our understanding of the microbial contribution to disease pathogenesis-knowledge that may eventually inform the development of novel therapeutics.

Identification of BACH2 and RAD51B as rheumatoid arthritis susceptibility loci in a meta-analysis of genome-wide data.

OBJECTIVE: A recent high-density fine-mapping (ImmunoChip) study of genetic associations in rheumatoid arthritis (RA) identified 14 risk loci with validated genome-wide significance, as well as a number of loci showing associations suggestive of significance (P = 5 × 10(-5) < 5 × 10(-8)), but these have yet to be replicated. The aim of this study was to determine whether these potentially significant loci are involved in the pathogenesis of RA, and to explore whether any of the loci are associated with a specific RA serotype. METHODS: A total of 16 single-nucleotide polymorphisms (SNPs) were selected for genotyping and association analyses in 2 independent validation cohorts, comprising 6,106 RA cases and 4,290 controls. A meta-analysis of the data from the original ImmunoChip discovery cohort and from both validation cohorts was carried out, for a combined total of 17,581 RA cases and 20,160 controls. In addition, stratified analysis of patient subsets, defined according to their anti-cyclic citrullinated peptide (anti-CCP) antibody status, was performed. RESULTS: A significant association with RA risk (P < 0.05) was replicated for 6 of the SNPs assessed in the validation cohorts. All SNPs in the validation study had odds ratios (ORs) for RA susceptibility in the same direction as those in the ImmunoChip discovery study. One SNP, rs72928038, mapping to an intron of BACH2, achieved genome-wide significance in the meta-analysis (P = 1.2 × 10(-8), OR 1.12), and a second SNP, rs911263, mapping to an intron of RAD51B, was significantly associated in the anti-CCP-positive RA subgroup (P = 4 × 10(-8), OR 0.89), confirming that both are RA susceptibility loci. CONCLUSION: This study provides robust evidence for an association of RA susceptibility with genes involved in B cell differentiation (BACH2) and DNA repair (RAD51B). The finding that the RAD51B gene exhibited different associations based on serologic subtype adds to the expanding knowledge base in defining subgroups of RA.

Functional Genomics and Insights into the Pathogenesis and Treatment of Psoriasis.

Psoriasis is a lifelong, systemic, immune mediated inflammatory skin condition, affecting 1-3% of the world's population, with an impact on quality of life similar to diseases like cancer or diabetes. Genetics are the single largest risk factor in psoriasis, with Genome-Wide Association (GWAS) studies showing that many psoriasis risk genes lie along the IL-23/Th17 axis. Potential psoriasis risk genes determined through GWAS can be annotated and characterised using functional genomics, allowing the identification of novel drug targets and the repurposing of existing drugs. This review is focused on the IL-23/Th17 axis, providing an insight into key cell types, cytokines, and intracellular signaling pathways involved. This includes examination of currently available biological treatments, time to relapse post drug withdrawal, and rates of primary/secondary drug failure, showing the need for greater understanding of the underlying genetic mechanisms of psoriasis and how they can impact treatment. This could allow for patient stratification towards the treatment most likely to reduce the burden of disease for the longest period possible.

The potential use of expression profiling: implications for predicting treatment response in rheumatoid arthritis.

Whole genome expression profiling, or transcriptomics, is a high throughput technology with the potential for major impacts in both clinical settings and drug discovery and diagnostics. In particular, there is much interest in this technique as a mechanism for predicting treatment response. Gene expression profiling entails the quantitative measurement of messenger RNA levels for thousands of genes simultaneously with the inherent possibility of identifying biomarkers of response to a particular therapy or by singling out those at risk of serious adverse events. This technology should contribute to the era of stratified medicine, in which patient specific populations are matched to potentially beneficial drugs via clinical tests. Indeed, in the oncology field, gene expression testing is already recommended to allow rational use of therapies to treat breast cancer. However, there are still many issues surrounding the use of the various testing platforms available and the statistical analysis associated with the interpretation of results generated. This review will discuss the implications this promising technology has in predicting treatment response and outline the various advantages and pitfalls associated with its use.

Identification of differences in CD4+ T-cell gene expression between people with asthma and healthy controls.

Functional enrichment analysis of genome-wide association study (GWAS)-summary statistics has suggested that CD4+ T-cells play an important role in asthma pathogenesis. Despite this, CD4+ T-cells are under-represented in asthma transcriptome studies. To fill the gap, 3'-RNA-Seq was used to generate gene expression data on CD4+ T-cells (isolated within 2 h from collection) from peripheral blood from participants with well-controlled asthma (n = 32) and healthy controls (n = 11). Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify sets of co-expressed genes (modules) associated with the asthma phenotype. We identified three modules associated with asthma, which are strongly enriched for GWAS-identified asthma genes, antigen processing/presentation and immune response to viral infections. Through integration of publicly available eQTL and GWAS summary statistics (colocalisation), and protein-protein interaction (PPI) data, we identified PTPRC, a potential druggable target, as a putative master regulator of the asthma gene-expression profiles. Using a co-expression network approach, with integration of external genetic and PPI data, we showed that CD4+ T-cells from peripheral blood from asthmatics have different expression profiles, albeit small in magnitude, compared to healthy controls, for sets of genes involved in immune response to viral infections (upregulated) and antigen processing/presentation (downregulated).

Identification of a regulatory pathway governing TRAF1 via an arthritis-associated non-coding variant.

TRAF1/C5 was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of TRAF1. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by FOSL2, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.

ASSIMILATOR: a new tool to inform selection of associated genetic variants for functional studies.

MOTIVATION: Fine-mapping experiments from genome-wide association studies (GWAS) are underway for many complex diseases. These are likely to identify a number of putative causal variants, which cannot be separated further in terms of strength of genetic association due to linkage disequilibrium. The challenge will be selecting which variant to prioritize for subsequent expensive functional studies. A wealth of functional information generated from wet lab experiments now exists but cannot be easily interrogated by the user. Here, we describe a program designed to quickly assimilate this data called ASSIMILATOR and validate the method by interrogating two regions to show its effectiveness. AVAILABILITY: http://www.medicine.manchester.ac.uk/musculoskeletal/research/arc/genetics/bioinformatics/assimilator/.