Parsing the Interferon Transcriptional Network and Its Disease Associations.

Type 1 interferon (IFN) is a key mediator of organismal responses to pathogens, eliciting prototypical "interferon signature genes" that encode antiviral and inflammatory mediators. For a global view of IFN signatures and regulatory pathways, we performed gene expression and chromatin analyses of the IFN-induced response across a range of immunocyte lineages. These distinguished ISGs by cell-type specificity, kinetics, and sensitivity to tonic IFN and revealed underlying changes in chromatin configuration. We combined 1,398 human and mouse datasets to computationally infer ISG modules and their regulators, validated by genetic analysis in both species. Some ISGs are controlled by Stat1/2 and Irf9 and the ISRE DNA motif, but others appeared dependent on non-canonical factors. This regulatory framework helped to interpret JAK1 blockade pharmacology, different clusters being affected under tonic or IFN-stimulated conditions, and the IFN signatures previously associated with human diseases, revealing unrecognized subtleties in disease footprints, as affected by human ancestry.

TDP-43 loss and ALS-risk SNPs drive mis-splicing and depletion of UNC13A.

Variants of UNC13A, a critical gene for synapse function, increase the risk of amyotrophic lateral sclerosis and frontotemporal dementia1-3, two related neurodegenerative diseases defined by mislocalization of the RNA-binding protein TDP-434,5. Here we show that TDP-43 depletion induces robust inclusion of a cryptic exon in UNC13A, resulting in nonsense-mediated decay and loss of UNC13A protein. Two common intronic UNC13A polymorphisms strongly associated with amyotrophic lateral sclerosis and frontotemporal dementia risk overlap with TDP-43 binding sites. These polymorphisms potentiate cryptic exon inclusion, both in cultured cells and in brains and spinal cords from patients with these conditions. Our findings, which demonstrate a genetic link between loss of nuclear TDP-43 function and disease, reveal the mechanism by which UNC13A variants exacerbate the effects of decreased TDP-43 function. They further provide a promising therapeutic target for TDP-43 proteinopathies.

Fine-mapping of Parkinson's disease susceptibility loci identifies putative causal variants.

Recent genome-wide association studies have identified 78 loci associated with Parkinson's disease susceptibility but the underlying mechanisms remain largely unclear. To identify likely causal variants for disease risk, we fine-mapped these Parkinson's-associated loci using four different fine-mapping methods. We then integrated multi-assay cell type-specific epigenomic profiles to pinpoint the likely mechanism of action of each variant, allowing us to identify Consensus single nucleotide polymorphism (SNPs) that disrupt LRRK2 and FCGR2A regulatory elements in microglia, an MBNL2 enhancer in oligodendrocytes, and a DYRK1A enhancer in neurons. This genome-wide functional fine-mapping investigation of Parkinson's disease substantially advances our understanding of the causal mechanisms underlying this complex disease while avoiding focus on spurious, non-causal mechanisms. Together, these results provide a robust, comprehensive list of the likely causal variants, genes and cell-types underlying Parkinson's disease risk as demonstrated by consistently greater enrichment of our fine-mapped SNPs relative to lead GWAS SNPs across independent functional impact annotations. In addition, our approach prioritized an average of 3/85 variants per locus as putatively causal, making downstream experimental studies both more tractable and more likely to yield disease-relevant, actionable results. Large-scale studies comparing individuals with Parkinson's disease to age-matched controls have identified many regions of the genome associated with the disease. However, there is widespread correlation between different parts of the genome, making it difficult to tell which genetic variants cause Parkinson's and which are simply co-inherited with causal variants. We therefore applied a suite of statistical models to identify the most likely causal genetic variants (i.e. fine-mapping). We then linked these genetic variants with epigenomic and gene expression signatures across a wide variety of tissues and cell types to identify how these variants cause disease. Therefore, this study provides a comprehensive and robust list of cellular and molecular mechanisms that may serve as targets in the development of more effective Parkinson's therapeutics.

echolocatoR: an automated end-to-end statistical and functional genomic fine-mapping pipeline.

SUMMARY: echolocatoR integrates a diverse suite of statistical and functional fine-mapping tools to identify, test enrichment in, and visualize high-confidence causal consensus variants in any phenotype. It requires minimal input from users (a summary statistics file), can be run in a single R function, and provides extensive access to relevant datasets (e.g. reference linkage disequilibrium panels, quantitative trait loci, genome-wide annotations, cell-type-specific epigenomics), thereby enabling rapid, robust and scalable end-to-end fine-mapping investigations. AVAILABILITY AND IMPLEMENTATION: echolocatoR is an open-source R package available through GitHub under the GNU General Public License (Version 3) license: https://github.com/RajLabMSSM/echolocatoR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Brain DNA methylomic analysis of frontotemporal lobar degeneration reveals OTUD4 in shared dysregulated signatures across pathological subtypes.

Frontotemporal lobar degeneration (FTLD) is an umbrella term describing the neuropathology of a clinically, genetically and pathologically heterogeneous group of diseases, including frontotemporal dementia (FTD) and progressive supranuclear palsy (PSP). Among the major FTLD pathological subgroups, FTLD with TDP-43 positive inclusions (FTLD-TDP) and FTLD with tau-positive inclusions (FTLD-tau) are the most common, representing about 90% of the cases. Although alterations in DNA methylation have been consistently associated with neurodegenerative diseases, including Alzheimer's disease and Parkinson's disease, little is known for FTLD and its heterogeneous subgroups and subtypes. The main goal of this study was to investigate DNA methylation variation in FTLD-TDP and FTLD-tau. We used frontal cortex genome-wide DNA methylation profiles from three FTLD cohorts (142 FTLD cases and 92 controls), generated using the Illumina 450K or EPIC microarrays. We performed epigenome-wide association studies (EWAS) for each cohort followed by meta-analysis to identify shared differentially methylated loci across FTLD subgroups/subtypes. In addition, we used weighted gene correlation network analysis to identify co-methylation signatures associated with FTLD and other disease-related traits. Wherever possible, we also incorporated relevant gene/protein expression data. After accounting for a conservative Bonferroni multiple testing correction, the EWAS meta-analysis revealed two differentially methylated loci in FTLD, one annotated to OTUD4 (5'UTR-shore) and the other to NFATC1 (gene body-island). Of these loci, OTUD4 showed consistent upregulation of mRNA and protein expression in FTLD. In addition, in the three independent co-methylation networks, OTUD4-containing modules were enriched for EWAS meta-analysis top loci and were strongly associated with the FTLD status. These co-methylation modules were enriched for genes implicated in the ubiquitin system, RNA/stress granule formation and glutamatergic synaptic signalling. Altogether, our findings identified novel FTLD-associated loci, and support a role for DNA methylation as a mechanism involved in the dysregulation of biological processes relevant to FTLD, highlighting novel potential avenues for therapeutic development.

Tensor decomposition of stimulated monocyte and macrophage gene expression profiles identifies neurodegenerative disease-specific trans-eQTLs.

Recent human genetic studies suggest that cells of the innate immune system have a primary role in the pathogenesis of neurodegenerative diseases. However, the results from these studies often do not elucidate how the genetic variants affect the biology of these cells to modulate disease risk. Here, we applied a tensor decomposition method to uncover disease associated gene networks linked to distal genetic variation in stimulated human monocyte and macrophage gene expression profiles. We report robust evidence that some disease associated genetic variants affect the expression of multiple genes in trans. These include a Parkinson's disease locus influencing the expression of genes mediated by a protease that controls lysosomal function, and Alzheimer's disease loci influencing the expression of genes involved in type 1 interferon signaling, myeloid phagocytosis, and complement cascade pathways. Overall, we uncover gene networks in induced innate immune cells linked to disease associated genetic variants, which may help elucidate the underlying biology of disease.

Multi-omic insights into Parkinson's Disease: From genetic associations to functional mechanisms.

Genome-Wide Association Studies (GWAS) have elucidated the genetic components of Parkinson's Disease (PD). However, because the vast majority of GWAS association signals fall within non-coding regions, translating these results into an interpretable, mechanistic understanding of the disease etiology remains a major challenge in the field. In this review, we provide an overview of the approaches to prioritize putative causal variants and genes as well as summarise the primary findings of previous studies. We then discuss recent efforts to integrate multi-omics data to identify likely pathogenic cell types and biological pathways implicated in PD pathogenesis. We have compiled full summary statistics of cell-type, tissue, and phentoype enrichment analyses from multiple studies of PD GWAS and provided them in a standardized format as a resource for the research community (https://github.com/RajLabMSSM/PD_omics_review). Finally, we discuss the experimental, computational, and conceptual advances that will be necessary to fully elucidate the effects of functional variants and genes on cellular dysregulation and disease risk.

Transcriptomic analysis of frontotemporal lobar degeneration with TDP-43 pathology reveals cellular alterations across multiple brain regions.

Frontotemporal lobar degeneration (FTLD) is a group of heterogeneous neurodegenerative disorders affecting the frontal and temporal lobes of the brain. Nuclear loss and cytoplasmic aggregation of the RNA-binding protein TDP-43 represents the major FTLD pathology, known as FTLD-TDP. To date, there is no effective treatment for FTLD-TDP due to an incomplete understanding of the molecular mechanisms underlying disease development. Here we compared postmortem tissue RNA-seq transcriptomes from the frontal cortex, temporal cortex, and cerebellum between 28 controls and 30 FTLD-TDP patients to profile changes in cell-type composition, gene expression and transcript usage. We observed downregulation of neuronal markers in all three regions of the brain, accompanied by upregulation of microglia, astrocytes, and oligodendrocytes, as well as endothelial cells and pericytes, suggesting shifts in both immune activation and within the vasculature. We validate our estimates of neuronal loss using neuropathological atrophy scores and show that neuronal loss in the cortex can be mainly attributed to excitatory neurons, and that increases in microglial and endothelial cell expression are highly correlated with neuronal loss. All our analyses identified a strong involvement of the cerebellum in the neurodegenerative process of FTLD-TDP. Altogether, our data provides a detailed landscape of gene expression alterations to help unravel relevant disease mechanisms in FTLD.

Dysregulated coordination of MAPT exon 2 and exon 10 splicing underlies different tau pathologies in PSP and AD.

Understanding regulation of MAPT splicing is important to the etiology of many nerurodegenerative diseases, including Alzheimer disease (AD) and progressive supranuclear palsy (PSP), in which different tau isoforms accumulate in pathologic inclusions. MAPT, the gene encoding the tau protein, undergoes complex alternative pre-mRNA splicing to generate six isoforms. Tauopathies can be categorized by the presence of tau aggregates containing either 3 (3R) or 4 (4R) microtubule-binding domain repeats (determined by inclusion/exclusion of exon 10), but the role of the N-terminal domain of the protein, determined by inclusion/exclusion of exons 2 and 3 has been less well studied. Using a correlational screen in human brain tissue, we observed coordination of MAPT exons 2 and 10 splicing. Expressions of exon 2 splicing regulators and subsequently exon 2 inclusion are differentially disrupted in PSP and AD brain, resulting in the accumulation of 1N4R isoforms in PSP and 0N isoforms in AD temporal cortex. Furthermore, we identified different N-terminal isoforms of tau present in neurofibrillary tangles, dystrophic neurites and tufted astrocytes, indicating a role for differential N-terminal splicing in the development of disparate tau neuropathologies. We conclude that N-terminal splicing and combinatorial regulation with exon 10 inclusion/exclusion is likely to be important to our understanding of tauopathies.

Distinct non-inflammatory signature of microglia in post-mortem brain tissue of patients with major depressive disorder.

Findings from epidemiological studies, biomarker measurements and animal experiments suggest a role for aberrant immune processes in the pathogenesis of major depressive disorder (MDD). Microglia, the resident immune cells of the brain, are likely to play a key role in these processes. Previous post-mortem studies reported conflicting findings regarding microglial activation and an in-depth profiling of those cells in MDD is lacking. The aim of this study was therefore to characterize the phenotype and function of microglia in MDD. We isolated microglia from post-mortem brain tissue of patients with MDD (n = 13-19) and control donors (n = 12-25). Using flow cytometry and quantitative Polymerase Chain Reaction (qPCR), we measured protein and mRNA levels of a panel of microglial markers across four different brain regions (medial frontal gyrus, superior temporal gyrus, thalamus, and subventricular zone). In MDD cases, we found a significant upregulation of CX3CR1 and TMEM119 mRNA expression and a downregulation of CD163 mRNA expression and CD14 protein expression across the four brain regions. Expression levels of microglial activation markers, such as HLA-DRA, IL6, and IL1ß, as well as the inflammatory responses to lipopolysaccharide and dexamethasone were unchanged. Our findings suggest that microglia enhance homeostatic functions in MDD but are not immune activated.

Landscape of Conditional eQTL in Dorsolateral Prefrontal Cortex and Co-localization with Schizophrenia GWAS.

Causal genes and variants within genome-wide association study (GWAS) loci can be identified by integrating GWAS statistics with expression quantitative trait loci (eQTL) and determining which variants underlie both GWAS and eQTL signals. Most analyses, however, consider only the marginal eQTL signal, rather than dissect this signal into multiple conditionally independent signals for each gene. Here we show that analyzing conditional eQTL signatures, which could be important under specific cellular or temporal contexts, leads to improved fine mapping of GWAS associations. Using genotypes and gene expression levels from post-mortem human brain samples (n = 467) reported by the CommonMind Consortium (CMC), we find that conditional eQTL are widespread; 63% of genes with primary eQTL also have conditional eQTL. In addition, genomic features associated with conditional eQTL are consistent with context-specific (e.g., tissue-, cell type-, or developmental time point-specific) regulation of gene expression. Integrating the 2014 Psychiatric Genomics Consortium schizophrenia (SCZ) GWAS and CMC primary and conditional eQTL data reveals 40 loci with strong evidence for co-localization (posterior probability > 0.8), including six loci with co-localization of conditional eQTL. Our co-localization analyses support previously reported genes, identify novel genes associated with schizophrenia risk, and provide specific hypotheses for their functional follow-up.

Collagenous Colitis Is Associated With HLA Signature and Shares Genetic Risks With Other Immune-Mediated Diseases.

BACKGROUND & AIMS: Collagenous colitis (CC) is an inflammatory bowel disorder with unknown etiopathogenesis involving HLA-related immune-mediated responses and environmental and genetic risk factors. We carried out an array-based genetic association study in a cohort of patients with CC and investigated the common genetic basis between CC and Crohn's disease (CD), ulcerative colitis (UC), and celiac disease. METHODS: DNA from 804 CC formalin-fixed, paraffin-embedded tissue samples was genotyped with Illumina Immunochip. Matching genotype data on control samples and CD, UC, and celiac disease cases were provided by the respective consortia. A discovery association study followed by meta-analysis with an independent cohort, polygenic risk score calculation, and cross-phenotype analyses were performed. Enrichment of regulatory expression quantitative trait loci among the CC variants was assessed in hemopoietic and intestinal cells. RESULTS: Three HLA alleles (HLA-B∗08:01, HLA-DRB1∗03:01, and HLA-DQB1∗02:01), related to the ancestral haplotype 8.1, were significantly associated with increased CC risk. We also identified an independent protective effect of HLA-DRB1∗04:01 on CC risk. Polygenic risk score quantifying the risk across multiple susceptibility loci was strongly associated with CC risk. An enrichment of expression quantitative trait loci was detected among the CC-susceptibility variants in various cell types. The cross-phenotype analysis identified a complex pattern of polygenic pleiotropy between CC and other immune-mediated diseases. CONCLUSIONS: In this largest genetic study of CC to date with histologically confirmed diagnosis, we strongly implicated the HLA locus and proposed potential non-HLA mechanisms in disease pathogenesis. We also detected a shared genetic risk between CC, celiac disease, CD, and UC, which supports clinical observations of comorbidity.

Genetic analysis of isoform usage in the human anti-viral response reveals influenza-specific regulation of ERAP2 transcripts under balancing selection.

While genetic variants are known to be associated with overall gene abundance in stimulated immune cells, less is known about their effects on alternative isoform usage. By analyzing RNA-seq profiles of monocyte-derived dendritic cells from 243 individuals, we uncovered thousands of unannotated isoforms synthesized in response to influenza infection and type 1 interferon stimulation. We identified more than a thousand quantitative trait loci (QTLs) associated with alternate isoform usage (isoQTLs), many of which are independent of expression QTLs (eQTLs) for the same gene. Compared with eQTLs, isoQTLs are enriched for splice sites and untranslated regions, but depleted of sequences upstream of annotated transcription start sites. Both eQTLs and isoQTLs explain a significant proportion of the disease heritability attributed to common genetic variants. At the ERAP2 locus, we shed light on the function of the gene and how two frequent, highly differentiated haplotypes with intermediate frequencies could be maintained by balancing selection. At baseline and following type 1 interferon stimulation, the major haplotype is associated with low ERAP2 expression caused by nonsense-mediated decay, while the minor haplotype, known to increase Crohn's disease risk, is associated with high ERAP2 expression. In response to influenza infection, we found two uncharacterized isoforms expressed from the major haplotype, likely the result of multiple perfectly linked variants affecting the transcription and splicing at the locus. Thus, genetic variants at a single locus could modulate independent gene regulatory processes in innate immune responses and, in the case of ERAP2, may confer a historical fitness advantage in response to virus.

HnRNP K mislocalisation is a novel protein pathology of frontotemporal lobar degeneration and ageing and leads to cryptic splicing.

Heterogeneous nuclear ribonucleoproteins (HnRNPs) are a group of ubiquitously expressed RNA-binding proteins implicated in the regulation of all aspects of nucleic acid metabolism. HnRNP K is a member of this highly versatile hnRNP family. Pathological redistribution of hnRNP K to the cytoplasm has been linked to the pathogenesis of several malignancies but, until now, has been underexplored in the context of neurodegenerative disease. Here we show hnRNP K mislocalisation in pyramidal neurons of the frontal cortex to be a novel neuropathological feature that is associated with both frontotemporal lobar degeneration and ageing. HnRNP K mislocalisation is mutually exclusive to TDP-43 and tau pathological inclusions in neurons and was not observed to colocalise with mitochondrial, autophagosomal or stress granule markers. De-repression of cryptic exons in RNA targets following TDP-43 nuclear depletion is an emerging mechanism of potential neurotoxicity in frontotemporal lobar degeneration and the mechanistically overlapping disorder amyotrophic lateral sclerosis. We silenced hnRNP K in neuronal cells to identify the transcriptomic consequences of hnRNP K nuclear depletion. Intriguingly, by performing RNA-seq analysis we find that depletion of hnRNP K induces 101 novel cryptic exon events. We validated cryptic exon inclusion in an SH-SY5Y hnRNP K knockdown and in FTLD brain exhibiting hnRNP K nuclear depletion. We, therefore, present evidence for hnRNP K mislocalisation to be associated with FTLD and for this to induce widespread changes in splicing.

Inferring Relevant Cell Types for Complex Traits by Using Single-Cell Gene Expression.

Previous studies have prioritized trait-relevant cell types by looking for an enrichment of genome-wide association study (GWAS) signal within functional regions. However, these studies are limited in cell resolution by the lack of functional annotations from difficult-to-characterize or rare cell populations. Measurement of single-cell gene expression has become a popular method for characterizing novel cell types, and yet limited work has linked single-cell RNA sequencing (RNA-seq) to phenotypes of interest. To address this deficiency, we present RolyPoly, a regression-based polygenic model that can prioritize trait-relevant cell types and genes from GWAS summary statistics and gene expression data. RolyPoly is designed to use expression data from either bulk tissue or single-cell RNA-seq. In this study, we demonstrated RolyPoly's accuracy through simulation and validated previously known tissue-trait associations. We discovered a significant association between microglia and late-onset Alzheimer disease and an association between schizophrenia and oligodendrocytes and replicating fetal cortical cells. Additionally, RolyPoly computes a trait-relevance score for each gene to reflect the importance of expression specific to a cell type. We found that differentially expressed genes in the prefrontal cortex of individuals with Alzheimer disease were significantly enriched with genes ranked highly by RolyPoly gene scores. Overall, our method represents a powerful framework for understanding the effect of common variants on cell types contributing to complex traits.

Large-Scale trans-eQTLs Affect Hundreds of Transcripts and Mediate Patterns of Transcriptional Co-regulation.

Efforts to decipher the causal relationships between differences in gene regulation and corresponding differences in phenotype have been stymied by several basic technical challenges. Although detecting local, cis-eQTLs is now routine, trans-eQTLs, which are distant from the genes of origin, are far more difficult to find because millions of SNPs must currently be compared to thousands of transcripts. Here, we demonstrate an alternative approach: we looked for SNPs associated with the expression of many genes simultaneously and found that hundreds of trans-eQTLs each affect hundreds of transcripts in lymphoblastoid cell lines across three African populations. These trans-eQTLs target the same genes across the three populations and show the same direction of effect. We discovered that target transcripts of a high-confidence set of trans-eQTLs encode proteins that interact more frequently than expected by chance, are bound by the same transcription factors, and are enriched for pathway annotations indicative of roles in basic cell homeostasis. We thus demonstrate that our approach can uncover trans-acting transcriptional control circuits that affect co-regulated groups of genes: a key to understanding how cellular pathways and processes are orchestrated.

A characterization of the molecular phenotype and inflammatory response of schizophrenia patient-derived microglia-like cells.

Different lines of evidence support a causal role for microglia in the pathogenesis of schizophrenia. However, how schizophrenia patient-derived microglia are affected at the phenotypic and functional level is still largely unknown. We used a recently described model to induce patient-derived microglia-like cells and used this to analyze changes in the molecular phenotype and function of myeloid cells in schizophrenia. We isolated monocytes from twenty recent-onset schizophrenia patients and twenty non-psychiatric controls. We cultured the cells towards an induced microglia-like phenotype (iMG), analyzed the phenotype of the cells by RNA sequencing and mass cytometry, and their response to LPS. Mass cytometry showed a high heterogeneity of iMG in cells derived from patients as well as controls. The prevalence of two iMG clusters was significantly higher in schizophrenia patients (adjusted p-value < 0.001). These subsets are characterized by expression of ApoE, Ccr2, CD18, CD44, and CD95, as well as IRF8, P2Y12, Cx3cr1 and HLA-DR. In addition, we found that patient-derived iMG show an enhanced response to LPS, with increased secretion of TNF-α. Further studies are needed to replicate these findings, to determine whether similar subclusters are present in schizophrenia patients in vivo, and to address how these subclusters are related to the increased response to LPS, as well as other microglial functions.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Rheumatoid arthritis-associated RBPJ polymorphism alters memory CD4+ T cells.

Notch signaling has recently emerged as an important regulator of immune responses in autoimmune diseases. The recombination signal-binding protein for immunoglobulin kappa J region (RBPJ) is a transcriptional repressor, but converts into a transcriptional activator upon activation of the canonical Notch pathway. Genome-wide association studies of rheumatoid arthritis (RA) identified a susceptibility locus, rs874040(CC), which implicated the RBPJ gene. Here, chromatin state mapping generated using the chromHMM algorithm reveals strong enhancer regions containing DNase I hypersensitive sites overlapping the rs874040 linkage disequilibrium block in human memory, but not in naïve CD4(+) T cells. The rs874040 overlapping this chromatin state was associated with increased RBPJ expression in stimulated memory CD4(+) T cells from healthy subjects homozygous for the risk allele (CC) compared with memory CD4(+) T cells bearing the protective allele (GG). Transcriptomic analysis of rs874040(CC) memory T cells showed a repression of canonical Notch target genes IL (interleukin)-9, IL-17 and interferon (IFN)γ in the basal state. Interestingly, activation of the Notch pathway using soluble Notch ligand, Jagged2-Fc, induced IL-9 and IL-17A while delta-like 4Fc, another Notch ligand, induced higher IFNγ expression in the rs874040(CC) memory CD4(+) T cells compared with their rs874040(GG) counterparts. In RA, RBPJ expression is elevated in memory T cells from RA patients compared with control subjects, and this was associated with induced inflammatory cytokines IL-9, IL-17A and IFNγ in response to Notch ligation in vitro. These findings demonstrate that the rs874040(CC) allele skews memory T cells toward a pro-inflammatory phenotype involving Notch signaling, thus increasing the susceptibility to develop RA.