Function of multiple sclerosis-protective HLA class I alleles revealed by genome-wide protein-quantitative trait loci mapping of interferon signalling.

Interferons (IFNs) are cytokines that are central to the host defence against viruses and other microorganisms. If not properly regulated, IFNs may contribute to the pathogenesis of inflammatory autoimmune, or infectious diseases. To identify genetic polymorphisms regulating the IFN system we performed an unbiased genome-wide protein-quantitative trait loci (pQTL) mapping of cell-type specific type I and type II IFN receptor levels and their responses in immune cells from 303 healthy individuals. Seven genome-wide significant (p < 5.0E-8) pQTLs were identified. Two independent SNPs that tagged the multiple sclerosis (MS)-protective HLA class I alleles A*02/A*68 and B*44, respectively, were associated with increased levels of IFNAR2 in B and T cells, with the most prominent effect in IgD-CD27+ memory B cells. The increased IFNAR2 levels in B cells were replicated in cells from an independent set of healthy individuals and in MS patients. Despite increased IFNAR2 levels, B and T cells carrying the MS-protective alleles displayed a reduced response to type I IFN stimulation. Expression and methylation-QTL analysis demonstrated increased mRNA expression of the pseudogene HLA-J in B cells carrying the MS-protective class I alleles, possibly driven via methylation-dependent transcriptional regulation. Together these data suggest that the MS-protective effects of HLA class I alleles are unrelated to their antigen-presenting function, and propose a previously unappreciated function of type I IFN signalling in B and T cells in MS immune-pathogenesis.

Association of aberrant ASNS imprinting with asparaginase sensitivity and chromosomal abnormality in childhood BCP-ALL.

Karyotype is an important prognostic factor in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL), but the underlying pharmacogenomics remain unknown. Asparaginase is an integral component in current chemotherapy for childhood BCP-ALL. Asparaginase therapy depletes serum asparagine. Normal hematopoietic cells can produce asparagine by asparagine synthetase (ASNS) activity, but ALL cells are unable to synthesize adequate amounts of asparagine. The ASNS gene has a typical CpG island in its promoter. Thus, methylation of the ASNS CpG island could be one of the epigenetic mechanisms for ASNS gene silencing in BCP-ALL. To gain deep insights into the pharmacogenomics of asparaginase therapy, we investigated the association of ASNS methylation status with asparaginase sensitivity. The ASNS CpG island is largely unmethylated in normal hematopoietic cells, but it is allele-specifically methylated in BCP-ALL cells. The ASNS gene is located at 7q21, an evolutionally conserved imprinted gene cluster. ASNS methylation in childhood BCP-ALL is associated with an aberrant methylation of the imprinted gene cluster at 7q21. Aberrant methylation of mouse Asns and a syntenic imprinted gene cluster is also confirmed in leukemic spleen samples from ETV6-RUNX1 knockin mice. In 3 childhood BCP-ALL cohorts, ASNS is highly methylated in BCP-ALL patients with favorable karyotypes but is mostly unmethylated in BCP-ALL patients with poor prognostic karyotypes. Higher ASNS methylation is associated with higher L-asparaginase sensitivity in BCP-ALL through lower ASNS gene and protein expression levels. These observations demonstrate that silencing of the ASNS gene as a result of aberrant imprinting is a pharmacogenetic mechanism for the leukemia-specific activity of asparaginase therapy in BCP-ALL.

Variants in BANK1 are associated with lupus nephritis of European ancestry.

The genetic background of lupus nephritis (LN) has not been completely elucidated. We performed a case-only study of 2886 SLE patients, including 947 (33%) with LN. Renal biopsies were available from 396 patients. The discovery cohort (Sweden, n = 1091) and replication cohort 1 (US, n = 962) were genotyped on the Immunochip and replication cohort 2 (Denmark/Norway, n = 833) on a custom array. Patients with LN, proliferative nephritis, or LN with end-stage renal disease were compared with SLE without nephritis. Six loci were associated with LN (p < 1 × 10-4, NFKBIA, CACNA1S, ITGA1, BANK1, OR2Y, and ACER3) in the discovery cohort. Variants in BANK1 showed the strongest association with LN in replication cohort 1 (p = 9.5 × 10-4) and proliferative nephritis in a meta-analysis of discovery and replication cohort 1. There was a weak association between BANK1 and LN in replication cohort 2 (p = 0.052), and in the meta-analysis of all three cohorts the association was strengthened (p = 2.2 × 10-7). DNA methylation data in 180 LN patients demonstrated methylation quantitative trait loci (meQTL) effects between a CpG site and BANK1 variants. To conclude, we describe genetic variations in BANK1 associated with LN and evidence for genetic regulation of DNA methylation within the BANK1 locus. This indicates a role for BANK1 in LN pathogenesis.

Interaction between the STAT4 rs11889341(T) risk allele and smoking confers increased risk of myocardial infarction and nephritis in patients with systemic lupus erythematosus.

OBJECTIVE: To investigate how genetics influence the risk of smoking-related systemic lupus erythematosus (SLE) manifestations. METHODS: Patients with SLE (ndiscovery cohort=776, nreplication cohort=836) were genotyped using the 200K Immunochip single nucleotide polymorphisms (SNP) Array (Illumina) and a custom array. Sixty SNPs with SLE association (p<5.0×10-8) were analysed. Signal transducer and activator of transcription 4 (STAT4) activation was assessed in in vitro stimulated peripheral blood mononuclear cells from healthy controls (n=45). RESULTS: In the discovery cohort, smoking was associated with myocardial infarction (MI) (OR 1.96 (95% CI 1.09 to 3.55)), with a greater effect in patients carrying any rs11889341 STAT4 risk allele (OR 2.72 (95% CI 1.24 to 6.00)) or two risk alleles (OR 8.27 (95% CI 1.48 to 46.27)).Smokers carrying the risk allele also displayed an increased risk of nephritis (OR 1.47 (95% CI 1.06 to 2.03)). In the replication cohort, the high risk of MI in smokers carrying the risk allele and the association between the STAT4 risk allele and nephritis in smokers were confirmed (OR 6.19 (95% CI 1.29 to 29.79) and 1.84 (95% CI 1.05 to 3.29), respectively).The interaction between smoking and the STAT4 risk allele resulted in further increase in the risk of MI (OR 2.14 (95% CI 1.01 to 4.62)) and nephritis (OR 1.53 (95% CI 1.08 to 2.17)), with 54% (MI) and 34% (nephritis) of the risk attributable to the interaction. Levels of interleukin-12-induced phosphorylation of STAT4 in CD8+ T cells were higher in smokers than in non-smokers (mean geometric fluorescence intensity 1063 vs 565, p=0.0063).Lastly, the IL12A rs564799 risk allele displayed association with MI in both cohorts (OR 1.53 (95% CI 1.01 to 2.31) and 2.15 (95% CI 1.08 to 4.26), respectively). CONCLUSIONS: Smoking in the presence of the STAT4 risk gene variant appears to increase the risk of MI and nephritis in SLE. Our results also highlight the role of the IL12-STAT4 pathway in SLE-cardiovascular morbidity.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing.

OBJECTIVES: Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE. METHODS: We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS). RESULTS: We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes. CONCLUSIONS: Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.

De novo mutations implicate novel genes in systemic lupus erythematosus.

The omnigenic model of complex disease stipulates that the majority of the heritability will be explained by the effects of common variation on genes in the periphery of core disease pathways. Rare variant associations, expected to explain far less of the heritability, may be enriched in core disease genes and thus will be instrumental in the understanding of complex disease pathogenesis and their potential therapeutic targets. Here, using complementary whole-exome sequencing, high-density imputation, and in vitro cellular assays, we identify candidate core genes in the pathogenesis of systemic lupus erythematosus (SLE). Using extreme-phenotype sampling, we sequenced the exomes of 30 SLE parent-affected-offspring trios and identified 14 genes with missense de novo mutations (DNM), none of which are within the >80 SLE susceptibility loci implicated through genome-wide association studies. In a follow-up cohort of 10, 995 individuals of matched European ancestry, we imputed genotype data to the density of the combined UK10K-1000 genomes Phase III reference panel across the 14 candidate genes. Gene-level analyses indicate three functional candidates: DNMT3A, PRKCD, and C1QTNF4. We identify a burden of rare variants across PRKCD associated with SLE risk (P = 0.0028), and across DNMT3A associated with two severe disease prognosis sub-phenotypes (P = 0.0005 and P = 0.0033). We further characterise the TNF-dependent functions of the third candidate gene C1QTNF4 on NF-κB activation and apoptosis, which are inhibited by the p.His198Gln DNM. Our results identify three novel genes in SLE susceptibility and support extreme-phenotype sampling and DNM gene discovery to aid the search for core disease genes implicated through rare variation.

Interferon signature in patients with STAT1 gain-of-function mutation is epigenetically determined.

STAT1 gain-of-function (GOF) variants lead to defective Th17 cell development and chronic mucocutaneous candidiasis (CMC), but frequently also to autoimmunity. Stimulation of cells with STAT1 inducing cytokines like interferons (IFN) result in hyperphosphorylation and delayed dephosphorylation of GOF STAT1. However, the mechanism how the delayed dephosphorylation exactly causes the increased expression of STAT1-dependent genes, and how the intracellular signal transduction from cytokine receptors is affected, remains unknown. In this study we show that the circulating levels of IFN-α were not persistently elevated in STAT1 GOF patients. Nevertheless, the expression of interferon signature genes was evident even in the patient with low or undetectable serum IFN-α levels. Chromatin immunoprecipitation (ChIP) experiments revealed that the active chromatin mark trimethylation of lysine 4 of histone 3 (H3K4me3), was significantly enriched in areas associated with interferon-stimulated genes in STAT1 GOF cells in comparison to cells from healthy donors. This suggests that the chromatin binding of GOF STAT1 variant promotes epigenetic changes compatible with higher gene expression and elevated reactivity to type I interferons, and possibly predisposes for interferon-related autoimmunity. The results also suggest that epigenetic rewiring may be responsible for treatment failure of Janus kinase 1/2 (JAK1/2) inhibitors in certain patients.

Genome-wide association study of angioedema induced by angiotensin-converting enzyme inhibitor and angiotensin receptor blocker treatment.

Angioedema in the mouth or upper airways is a feared adverse reaction to angiotensin-converting enzyme inhibitor (ACEi) and angiotensin receptor blocker (ARB) treatment, which is used for hypertension, heart failure and diabetes complications. This candidate gene and genome-wide association study aimed to identify genetic variants predisposing to angioedema induced by these drugs. The discovery cohort consisted of 173 cases and 4890 controls recruited in Sweden. In the candidate gene analysis, ETV6, BDKRB2, MME, and PRKCQ were nominally associated with angioedema (p < 0.05), but did not pass Bonferroni correction for multiple testing (p < 2.89 × 10-5). In the genome-wide analysis, intronic variants in the calcium-activated potassium channel subunit alpha-1 (KCNMA1) gene on chromosome 10 were significantly associated with angioedema (p < 5 × 10-8). Whilst the top KCNMA1 hit was not significant in the replication cohort (413 cases and 599 ACEi-exposed controls from the US and Northern Europe), a meta-analysis of the replication and discovery cohorts (in total 586 cases and 1944 ACEi-exposed controls) revealed that each variant allele increased the odds of experiencing angioedema 1.62 times (95% confidence interval 1.05-2.50, p = 0.030). Associated KCNMA1 variants are not known to be functional, but are in linkage disequilibrium with variants in transcription factor binding sites active in relevant tissues. In summary, our data suggest that common variation in KCNMA1 is associated with risk of angioedema induced by ACEi or ARB treatment. Future whole exome or genome sequencing studies will show whether rare variants in KCNMA1 or other genes contribute to the risk of ACEi- and ARB-induced angioedema.

High genetic risk score is associated with early disease onset, damage accrual and decreased survival in systemic lupus erythematosus.

OBJECTIVES: To investigate associations between a high genetic disease risk and disease severity in patients with systemic lupus erythematosus (SLE). METHODS: Patients with SLE (n=1001, discovery cohort and n=5524, replication cohort) and healthy controls (n=2802 and n=9859) were genotyped using a 200K Immunochip single nucleotide polymorphism array. A genetic risk score (GRS) was assigned to each individual based on 57 SLE risk loci. RESULTS: SLE was more prevalent in the high, compared with the low, GRS-quartile (OR 12.32 (9.53 to 15.71), p=7.9×10-86 and OR 7.48 (6.73 to 8.32), p=2.2×10-304 for the discovery and the replication cohorts, respectively). In the discovery cohort, patients in the high GRS-quartile had a 6-year earlier mean disease onset (HR 1.47 (1.22 to 1.75), p=4.3×10-5), displayed higher prevalence of damage accrual (OR 1.47 (1.06 to 2.04), p=2.0×10-2), renal disorder (OR 2.22 (1.50 to 3.27), p=5.9×10-5), anti-dsDNA (OR 1.83 (1.19 to 2.81), p=6.1×10-3), end-stage renal disease (ESRD) (OR 5.58 (1.50 to 20.79), p=1.0×10-2), proliferative nephritis (OR 2.42 (1.30 to 4.49), p=5.1×10-3), anti-cardiolipin-IgG (OR 1.89 (1.13 to 3.18), p=1.6×10-2), anti-ß2-glycoprotein-I-IgG (OR 2.29 (1.29 to 4.06), p=4.8×10-3) and positive lupus anticoagulant test (OR 2.12 (1.16 to 3.89), p=1.5×10-2) compared with patients in the low GRS-quartile. Survival analysis showed earlier onset of the first organ damage (HR 1.51 (1.04 to 2.25), p=3.7×10-2), first cardiovascular event (HR 1.65 (1.03 to 2.64), p=2.6×10-2), nephritis (HR 2.53 (1.72 to 3.71), p=9.6×10-7), ESRD (HR 6.78 (1.78 to 26.86), p=6.5×10-3) and decreased overall survival (HR 1.83 (1.02 to 3.30), p=4.3×10-2) in high to low quartile comparison. CONCLUSIONS: A high GRS is associated with increased risk of organ damage, renal dysfunction and all-cause mortality. Our results indicate that genetic profiling may be useful for predicting outcomes in patients with SLE.

Epigenetics in pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.

Whole-genome sequencing identifies complex contributions to genetic risk by variants in genes causing monogenic systemic lupus erythematosus.

Systemic lupus erythematosus (SLE, OMIM 152700) is a systemic autoimmune disease with a complex etiology. The mode of inheritance of the genetic risk beyond familial SLE cases is currently unknown. Additionally, the contribution of heterozygous variants in genes known to cause monogenic SLE is not fully understood. Whole-genome sequencing of DNA samples from 71 Swedish patients with SLE and their healthy biological parents was performed to investigate the general genetic risk of SLE using known SLE GWAS risk loci identified using the ImmunoChip, variants in genes associated to monogenic SLE, and the mode of inheritance of SLE risk alleles in these families. A random forest model for predicting genetic risk for SLE showed that the SLE risk variants were mainly inherited from one of the parents. In the 71 patients, we detected a significant enrichment of ultra-rare ( ≤ 0.1%) missense and nonsense mutations in 22 genes known to cause monogenic forms of SLE. We identified one previously reported homozygous nonsense mutation in the C1QC (Complement C1q C Chain) gene, which explains the immunodeficiency and severe SLE phenotype of that patient. We also identified seven ultra-rare, coding heterozygous variants in five genes (C1S, DNASE1L3, DNASE1, IFIH1, and RNASEH2A) involved in monogenic SLE. Our findings indicate a complex contribution to the overall genetic risk of SLE by rare variants in genes associated with monogenic forms of SLE. The rare variants were inherited from the other parent than the one who passed on the more common risk variants leading to an increased genetic burden for SLE in the child. Higher frequency SLE risk variants are mostly passed from one of the parents to the offspring affected with SLE. In contrast, the other parent, in seven cases, contributed heterozygous rare variants in genes associated with monogenic forms of SLE, suggesting a larger impact of rare variants in SLE than hitherto reported.

Genome-wide repression of eRNA and target gene loci by the ETV6-RUNX1 fusion in acute leukemia.

Approximately 20%-25% of childhood acute lymphoblastic leukemias carry the ETV6-RUNX1 (E/R) fusion gene, a fusion of two central hematopoietic transcription factors, ETV6 (TEL) and RUNX1 (AML1). Despite its prevalence, the exact genomic targets of E/R have remained elusive. We evaluated gene loci and enhancers targeted by E/R genome-wide in precursor B acute leukemia cells using global run-on sequencing (GRO-seq). We show that expression of the E/R fusion leads to widespread repression of RUNX1 motif-containing enhancers at its target gene loci. Moreover, multiple super-enhancers from the CD19+/CD20+-lineage were repressed, implicating a role in impediment of lineage commitment. In effect, the expression of several genes involved in B cell signaling and adhesion was down-regulated, and the repression depended on the wild-type DNA-binding Runt domain of RUNX1. We also identified a number of E/R-regulated annotated and de novo noncoding genes. The results provide a comprehensive genome-wide mapping between E/R-regulated key regulatory elements and genes in precursor B cell leukemia that disrupt normal B lymphopoiesis.

Genetic variations in A20 DUB domain provide a genetic link to citrullination and neutrophil extracellular traps in systemic lupus erythematosus.

OBJECTIVES: Genetic variations in TNFAIP3 (A20) de-ubiquitinase (DUB) domain increase the risk of systemic lupus erythematosus (SLE) and rheumatoid arthritis. A20 is a negative regulator of NF-κB but the role of its DUB domain and related genetic variants remain unclear. We aimed to study the functional effects of A20 DUB-domain alterations in immune cells and understand its link to SLE pathogenesis. METHODS: CRISPR/Cas9 was used to generate human U937 monocytes with A20 DUB-inactivating C103A knock-in (KI) mutation. Whole genome RNA-sequencing was used to identify differentially expressed genes between WT and C103A KI cells. Functional studies were performed in A20 C103A U937 cells and in immune cells from A20 C103A mice and genotyped healthy individuals with A20 DUB polymorphism rs2230926. Neutrophil extracellular trap (NET) formation was addressed ex vivo in neutrophils from A20 C103A mice and SLE-patients with rs2230926. RESULTS: Genetic disruption of A20 DUB domain in human and murine myeloid cells did not give rise to enhanced NF-κB signalling. Instead, cells with C103A mutation or rs2230926 polymorphism presented an upregulated expression of PADI4, an enzyme regulating protein citrullination and NET formation, two key mechanisms in autoimmune pathology. A20 C103A cells exhibited enhanced protein citrullination and extracellular trap formation, which could be suppressed by selective PAD4 inhibition. Moreover, SLE-patients with rs2230926 showed increased NETs and increased frequency of autoantibodies to citrullinated epitopes. CONCLUSIONS: We propose that genetic alterations disrupting the A20 DUB domain mediate increased susceptibility to SLE through the upregulation of PADI4 with resultant protein citrullination and extracellular trap formation.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project--the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Neuronal Expression of Opioid Gene is Controlled by Dual Epigenetic and Transcriptional Mechanism in Human Brain.

Molecular mechanisms that define patterns of neuropeptide expression are essential for the formation and rewiring of neural circuits. The prodynorphin gene (PDYN) gives rise to dynorphin opioid peptides mediating depression and substance dependence. We here demonstrated that PDYN is expressed in neurons in human dorsolateral prefrontal cortex (dlPFC), and identified neuronal differentially methylated region in PDYN locus framed by CCCTC-binding factor binding sites. A short, nucleosome size human-specific promoter CpG island (CGI), a core of this region may serve as a regulatory module, which is hypomethylated in neurons, enriched in 5-hydroxymethylcytosine, and targeted by USF2, a methylation-sensitive E-box transcription factor (TF). USF2 activates PDYN transcription in model systems, and binds to nonmethylated CGI in dlPFC. USF2 and PDYN expression is correlated, and USF2 and PDYN proteins are co-localized in dlPFC. Segregation of activatory TF and repressive CGI methylation may ensure contrasting PDYN expression in neurons and glia in human brain.

SPlinted Ligation Adapter Tagging (SPLAT), a novel library preparation method for whole genome bisulphite sequencing.

Sodium bisulphite treatment of DNA combined with next generation sequencing (NGS) is a powerful combination for the interrogation of genome-wide DNA methylation profiles. Library preparation for whole genome bisulphite sequencing (WGBS) is challenging due to side effects of the bisulphite treatment, which leads to extensive DNA damage. Recently, a new generation of methods for bisulphite sequencing library preparation have been devised. They are based on initial bisulphite treatment of the DNA, followed by adaptor tagging of single stranded DNA fragments, and enable WGBS using low quantities of input DNA. In this study, we present a novel approach for quick and cost effective WGBS library preparation that is based on splinted adaptor tagging (SPLAT) of bisulphite-converted single-stranded DNA. Moreover, we validate SPLAT against three commercially available WGBS library preparation techniques, two of which are based on bisulphite treatment prior to adaptor tagging and one is a conventional WGBS method.

Genome-wide association and Mendelian randomization study of NT-proBNP in patients with acute coronary syndrome.

N-terminal pro-B-type natriuretic peptide (NT-proBNP) is a strong predictor of mortality in coronary artery disease and is widely employed as a prognostic biomarker. However, a causal relationship between NT-proBNP and clinical endpoints has not been established. We have performed a genome-wide association and Mendelian randomization study of NT-proBNP. We used a discovery set of 3740 patients from the PLATelet inhibition and patient Outcomes (PLATO) trial, which enrolled 18 624 patients with acute coronary syndrome (ACS). A further set of 5492 patients, from the same trial, was used for replication. Genetic variants at two novel loci (SLC39A8 and POC1B/GALNT4) were associated with NT-proBNP levels and replicated together with the previously known NPPB locus. The most significant SNP (rs198389, pooled P = 1.07 × 10(-15)) in NPPB interrupts an E-box consensus motif in the gene promoter. The association in SLC39A8 is driven by a deleterious variant (rs13107325, pooled P = 5.99 × 10(-10)), whereas the most significant SNP in POC1B/GALNT4 (rs11105306, pooled P = 1.02 × 10(-16)) is intronic. The SLC39A8 SNP was associated with higher risk of cardiovascular (CV) death (HR = 1.39, 95% CI: 1.08-1.79, P = 0.0095), but the other loci were not associated with clinical endpoints. We have identified two novel loci to be associated with NT-proBNP in patients with ACS. Only the SLC39A8 variant, but not the NPPB variant, was associated with a clinical endpoint. Due to pleotropic effects of SLC39A8, these results do not suggest that NT-proBNP levels have a direct effect on mortality in ACS patients. PLATO Clinical Trial Registration: www.clinicaltrials.gov; NCT00391872.

The STAT4 SLE risk allele rs7574865[T] is associated with increased IL-12-induced IFN-γ production in T cells from patients with SLE.

OBJECTIVES: Genetic variants in the transcription factor STAT4 are associated with increased susceptibility to systemic lupus erythematosus (SLE) and a more severe disease phenotype. This study aimed to clarify how the SLE-associated intronic STAT4 risk allele rs7574865[T] affects the function of immune cells in SLE. METHODS: Peripheral blood mononuclear cells (PBMCs) were isolated from 52 genotyped patients with SLE. Phosphorylation of STAT4 (pSTAT4) and STAT1 (pSTAT1) in response to interferon (IFN)-α, IFN-γ or interleukin (IL)-12, total levels of STAT4, STAT1 and T-bet, and frequency of IFN-γ+ cells on IL-12 stimulation were determined by flow cytometry in subsets of immune cells before and after preactivation of cells with phytohaemagglutinin (PHA) and IL-2. Cellular responses and phenotypes were correlated to STAT4 risk allele carriership. Janus kinase inhibitors (JAKi) selective for TYK2 (TYK2i) or JAK2 (JAK2i) were evaluated for inhibition of IL-12 or IFN-γ-induced activation of SLE PBMCs. RESULTS: In resting PBMCs, the STAT4 risk allele was neither associated with total levels of STAT4 or STAT1, nor cytokine-induced pSTAT4 or pSTAT1. Following PHA/IL-2 activation, CD8+ T cells from STAT4 risk allele carriers displayed increased levels of STAT4 resulting in increased pSTAT4 in response to IL-12 and IFN-α, and an augmented IL-12-induced IFN-γ production in CD8+ and CD4+ T cells. The TYK2i and the JAK2i efficiently blocked IL-12 and IFN-γ-induced activation of PBMCs from STAT4 risk patients, respectively. CONCLUSIONS: T cells from patients with SLE carrying the STAT4 risk allele rs7574865[T] display an augmented response to IL-12 and IFN-α. This subset of patients may benefit from JAKi treatment.

DNA methylation mapping identifies gene regulatory effects in patients with systemic lupus erythematosus.

OBJECTIVES: Systemic lupus erythematosus (SLE) is a chronic autoimmune condition with heterogeneous presentation and complex aetiology where DNA methylation changes are emerging as a contributing factor. In order to discover novel epigenetic associations and investigate their relationship to genetic risk for SLE, we analysed DNA methylation profiles in a large collection of patients with SLE and healthy individuals. METHODS: DNA extracted from blood from 548 patients with SLE and 587 healthy controls were analysed on the Illumina HumanMethylation 450 k BeadChip, which targets 485 000 CpG sites across the genome. Single nucleotide polymorphism (SNP) genotype data for 196 524 SNPs on the Illumina ImmunoChip from the same individuals were utilised for methylation quantitative trait loci (cis-meQTLs) analyses. RESULTS: We identified and replicated differentially methylated CpGs (DMCs) in SLE at 7245 CpG sites in the genome. The largest methylation differences were observed at type I interferon-regulated genes which exhibited decreased methylation in SLE. We mapped cis-meQTLs and identified genetic regulation of methylation levels at 466 of the DMCs in SLE. The meQTLs for DMCs in SLE were enriched for genetic association to SLE, and included seven SLE genome-wide association study (GWAS) loci: PTPRC (CD45), MHC-class III, UHRF1BP1, IRF5, IRF7, IKZF3 and UBE2L3. In addition, we observed association between genotype and variance of methylation at 20 DMCs in SLE, including at the HLA-DQB2 locus. CONCLUSIONS: Our results suggest that several of the genetic risk variants for SLE may exert their influence on the phenotype through alteration of DNA methylation levels at regulatory regions of target genes.

Novel gene variants associated with cardiovascular disease in systemic lupus erythematosus and rheumatoid arthritis.

OBJECTIVES: Patients with systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) have increased risk of cardiovascular disease (CVD). We investigated whether single nucleotide polymorphisms (SNPs) at autoimmunity risk loci were associated with CVD in SLE and RA. METHODS: Patients with SLE (n=1045) were genotyped using the 200K Immunochip SNP array (Illumina). The allele frequency was compared between patients with and without different manifestations of CVD. Results were replicated in a second SLE cohort (n=1043) and in an RA cohort (n=824). We analysed publicly available genetic data from general population, performed electrophoretic mobility shift assays and measured cytokine levels and occurrence of antiphospholipid antibodies (aPLs). RESULTS: We identified two new putative risk loci associated with increased risk for CVD in two SLE populations, which remained after adjustment for traditional CVD risk factors. An IL19 risk allele, rs17581834(T) was associated with stroke/myocardial infarction (MI) in SLE (OR 2.3 (1.5 to 3.4), P=8.5×10-5) and RA (OR 2.8 (1.4 to 5.6), P=3.8×10-3), meta-analysis (OR 2.5 (2.0 to 2.9), P=3.5×10-7), but not in population controls. The IL19 risk allele affected protein binding, and SLE patients with the risk allele had increased levels of plasma-IL10 (P=0.004) and aPL (P=0.01). An SRP54-AS1 risk allele, rs799454(G) was associated with stroke/transient ischaemic attack in SLE (OR 1.7 (1.3 to 2.2), P=2.5×10-5) but not in RA. The SRP54-AS1 risk allele is an expression quantitative trait locus for four genes. CONCLUSIONS: The IL19 risk allele was associated with stroke/MI in SLE and RA, but not in the general population, indicating that shared immune pathways may be involved in the CVD pathogenesis in inflammatory rheumatic diseases.