Complement genes contribute sex-biased vulnerability in diverse disorders.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Harmonized Multisite MRI-Based Quantification of Human Liver Fat and Stiffness: A Pilot Study.

BACKGROUND: Nonalcoholic fatty liver disease (NAFLD) is a leading cause of end-stage liver disease. NAFLD diagnosis and follow-up relies on a combination of clinical data, liver imaging, and/or liver biopsy. However, intersite imaging differences impede diagnostic consistency and reduce the repeatability of the multisite clinical trials necessary to develop effective treatments. PURPOSE/HYPOTHESIS: The goal of this pilot study was to harmonize commercially available 3 T magnetic resonance imaging (MRI) measurements of liver fat and stiffness in human participants across academic sites and MRI vendors. STUDY TYPE: Cohort. SUBJECTS: Four community-dwelling adults with obesity. FIELD STRENGTH/SEQUENCE: 1.5 and 3 T, multiecho 3D imaging, PRESS, and GRE. ASSESSMENT: Harmonized proton density fat fraction (PDFF) and magnetic resonance spectroscopy (MRS) protocols were used to quantify the FF of synthetic phantoms and human participants with obesity using standard acquisition parameters at four sites that had four different 3 T MRI instruments. In addition, a harmonized magnetic resonance elastography (MRE) protocol was used to quantify liver stiffness among participants at two different sites at 1.5 and 3 T field strengths. Data were sent to a single data coordinating site for postprocessing. STATISTICAL TESTS: Linear regression in MATLAB, ICC analyses using SAS 9.4, one-sided 95% confidence intervals for the ICC. RESULTS: PDFF and MRS FF measurements were highly repeatable among sites in both humans and phantoms. MRE measurements of liver stiffness in three individuals at two sites using one 1.5 T and one 3 T instrument showed repeatability that was high although lower than that of MRS and PDFF. CONCLUSIONS: We demonstrated harmonization of PDFF, MRS, and MRE-based quantification of liver fat and stiffness through synthetic phantoms, traveling participants, and standardization of postprocessing analysis. Multisite MRI harmonization could contribute to multisite clinical trials assessing the efficacy of interventions and therapy for NAFLD. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 2.

Skewed allelic expression on X chromosome associated with aberrant expression of XIST on systemic lupus erythematosus lymphocytes.

A common feature of autoimmune diseases, including systemic lupus erythematosus (SLE), is an increased prevalence in women. However, the molecular basis for sex disparity in SLE remains poorly understood. To examine the role of X-linked transcription in SLE adaptive immune cells, we performed RNA-seq in T cell and B cell subsets from either healthy donors or patients with SLE. Analyses of allelic expression (AE) profiles identified a pattern of increased allelic imbalance across the entire X chromosome in SLE lymphocytes. X-linked genes exhibiting AE in SLE had an extensive overlap with genes known to escape X chromosome inactivation (XCI). XIST RNA was overexpressed in SLE patients. Differential XIST expression correlated with AE profiles more positively at X-linked genes than the genome-wide background. Analysis of three independent RNA-seq data verified the XIST-associated skewed AE on X chromosome in SLE. Integrative analyses of DNA methylation profiles showed an increased variability of DNA methylation levels at these AE-related X-linked genes. In cultured lymphoblastic cells, knockdown of XIST specifically altered allelic imbalance patterns between X chromosomes. Our study provides genetic evidence that upregulation of XIST accompanied with more skewed allelic expression on X chromosome is associated with the pathogenesis of SLE and may provide mechanistic insights into the increased incidence of SLE in females.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans.

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

MRLR: unraveling high-resolution meiotic recombination by linked reads.

MOTIVATION: Meiotic recombination facilitates the transmission of exchanged genetic material between homologous chromosomes and plays a crucial role in increasing the genetic variations in eukaryotic organisms. In humans, thousands of crossover events have been identified by genotyping related family members. However, most of these crossover regions span tens to hundreds of kb, which is not sufficient resolution to accurately identify the crossover breakpoints in a typical trio family. RESULTS: We have developed MRLR, a software using 10X linked reads to identify crossover events at a high resolution. By reconstructing the gamete genome, MRLR only requires a trio family dataset and can efficiently discover the crossover events. Using MRLR, we revealed a fine-scale pattern of crossover regions in six human families. From the two closest heterozygous alleles around the crossovers, we determined that MRLR achieved a median resolution 4.5 kb. This method can delineate a genome-wide landscape of crossover events at a precise scale, which is important for both functional and genomic features analysis of meiotic recombination. AVAILABILITY AND IMPLEMENTATION: MRLR is freely available at https://github.com/ChongLab/MRLR, implemented in Perl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Genetic variants at the 16p13 locus confer risk for eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a chronic inflammatory disease of the esophagus triggered by immune hypersensitivity to food. Herein, we tested whether genetic risk factors for known, non-allergic, immune-mediated diseases, particularly those involving autoimmunity, were associated with EoE risk. We used the high-density Immunochip platform, encoding 200,000 genetic variants for major auto-immune disease. Accordingly, 1214 subjects with EoE of European ancestry and 3734 population controls were genotyped and assessed using data directly generated or imputed from the previously published GWAS. We found lack of association of EoE with the genetic variants in the major histocompatibility complex (MHC) class I, II, and III genes and nearly all other loci using a highly powered study design with dense genotyping throughout the locus. Importantly, we identified an EoE risk locus at 16p13 with genome-wide significance (Pcombined=2.05 × 10-9, odds ratio = 0.76-0.81). This region is known to encode for the genes CLEC16A, DEXI, and CIITI, which are expressed in immune cells and esophageal epithelial cells. Suggestive EoE risk were also seen 5q23 (intergenic) and 7p15 (JAZF1). Overall, we have identified an additional EoE risk locus at 16p13 and highlight a shared and unique genetic etiology of EoE with a spectrum of immune-associated diseases.

Individualized decision aid for diverse women with lupus nephritis (IDEA-WON): A randomized controlled trial.

BACKGROUND: Treatment decision-making regarding immunosuppressive therapy is challenging for individuals with lupus. We assessed the effectiveness of a decision aid for immunosuppressive therapy in lupus nephritis. METHODS AND FINDINGS: In a United States multicenter, open-label, randomized controlled trial (RCT), adult women with lupus nephritis, mostly from racial/ethnic minority backgrounds with low socioeconomic status (SES), seen in in- or outpatient settings, were randomized to an individualized, culturally tailored, computerized decision aid versus American College of Rheumatology (ACR) lupus pamphlet (1:1 ratio), using computer-generated randomization. We hypothesized that the co-primary outcomes of decisional conflict and informed choice regarding immunosuppressive medications would improve more in the decision aid group. Of 301 randomized women, 298 were analyzed; 47% were African-American, 26% Hispanic, and 15% white. Mean age (standard deviation [SD]) was 37 (12) years, 57% had annual income of <$40,000, and 36% had a high school education or less. Compared with the provision of the ACR lupus pamphlet (n = 147), participants randomized to the decision aid (n = 151) had (1) a clinically meaningful and statistically significant reduction in decisional conflict, 21.8 (standard error [SE], 2.5) versus 12.7 (SE, 2.0; p = 0.005) and (2) no difference in informed choice in the main analysis, 41% versus 31% (p = 0.08), but clinically meaningful and statistically significant difference in sensitivity analysis (net values for immunosuppressives positive [in favor] versus negative [against]), 50% versus 35% (p = 0.006). Unresolved decisional conflict was lower in the decision aid versus pamphlet groups, 22% versus 44% (p < 0.001). Significantly more patients in the decision aid versus pamphlet group rated information to be excellent for understanding lupus nephritis (49% versus 33%), risk factors (43% versus 27%), medication options (50% versus 33%; p ≤ 0.003 for all); and the ease of use of materials was higher in the decision aid versus pamphlet groups (51% versus 38%; p = 0.006). Key study limitations were the exclusion of men, short follow-up, and the lack of clinical outcomes, including medication adherence. CONCLUSIONS: An individualized decision aid was more effective than usual care in reducing decisional conflict for choice of immunosuppressive medications in women with lupus nephritis. TRIAL REGISTRATION: Clinicaltrials.gov, NCT02319525.

Association of response to TNF inhibitors in rheumatoid arthritis with quantitative trait loci for CD40 and CD39.

OBJECTIVES: We sought to investigate whether genetic effects on response to TNF inhibitors (TNFi) in rheumatoid arthritis (RA) could be localised by considering known genetic susceptibility loci for relevant traits and to evaluate the usefulness of these genetic loci for stratifying drug response. METHODS: We studied the relation of TNFi response, quantified by change in swollen joint counts ( Δ SJC) and erythrocyte sedimentation rate ( Δ ESR) with locus-specific scores constructed from genome-wide assocation study summary statistics in 2938 genotyped individuals: 37 scores for RA; scores for 19 immune cell traits; scores for expression or methylation of 93 genes with previously reported associations between transcript level and drug response. Multivariate associations were evaluated in penalised regression models by cross-validation. RESULTS: We detected a statistically significant association between Δ SJC and the RA score at the CD40 locus (p=0.0004) and an inverse association between Δ SJC and the score for expression of CD39 on CD4 T cells (p=0.00005). A previously reported association between CD39 expression on regulatory T cells and response to methotrexate was in the opposite direction. In stratified analysis by concomitant methotrexate treatment, the inverse association was stronger in the combination therapy group and dissipated in the TNFi monotherapy group. Overall, ability to predict TNFi response from genotypic scores was limited, with models explaining less than 1% of phenotypic variance. CONCLUSIONS: The association with the CD39 trait is difficult to interpret because patients with RA are often prescribed TNFi after failing to respond to methotrexate. The CD39 and CD40 pathways could be relevant for targeting drug therapy.

Genome-wide association studies suggest that APOL1-environment interactions more likely trigger kidney disease in African Americans with nondiabetic nephropathy than strong APOL1-second gene interactions.

African Americans carrying two apolipoprotein L1 gene (APOL1) renal risk variants have a high risk for nephropathy. However, only a minority develops end-stage renal disease (ESRD). Hence, modifying factors likely contribute to initiation of kidney disease such as endogenous (HIV infection) or exogenous (interferon treatment) environmental modifiers. In this report, genome-wide association studies and a meta-analysis were performed to identify novel loci for nondiabetic ESRD in African Americans and to detect genetic modifiers in APOL1-associated nephropathy. Two African American cohorts were analyzed, 1749 nondiabetic ESRD cases and 1136 controls from Wake Forest and 901 lupus nephritis (LN)-ESRD cases and 520 controls with systemic lupus erythematosus but lacking nephropathy from the LN-ESRD Consortium. Association analyses adjusting for APOL1 G1/G2 renal-risk variants were completed and stratified by APOL1 risk genotype status. Individual genome-wide association studies and meta-analysis results of all 2650 ESRD cases and 1656 controls did not detect significant genome-wide associations with ESRD beyond APOL1. Similarly, no single nucleotide polymorphism showed significant genome-wide evidence of an interaction with APOL1 risk variants. Thus, although variants with small individual effects cannot be ruled out and are likely to exist, our results suggest that APOL1-environment interactions may be of greater clinical importance in triggering nephropathy in African Americans than APOL1 interactions with other single nucleotide polymorphisms.

Lupus Risk Variant Increases pSTAT1 Binding and Decreases ETS1 Expression.

Genetic variants at chromosomal region 11q23.3, near the gene ETS1, have been associated with systemic lupus erythematosus (SLE), or lupus, in independent cohorts of Asian ancestry. Several recent studies have implicated ETS1 as a critical driver of immune cell function and differentiation, and mice deficient in ETS1 develop an SLE-like autoimmunity. We performed a fine-mapping study of 14,551 subjects from multi-ancestral cohorts by starting with genotyped variants and imputing to all common variants spanning ETS1. By constructing genetic models via frequentist and Bayesian association methods, we identified 16 variants that are statistically likely to be causal. We functionally assessed each of these variants on the basis of their likelihood of affecting transcription factor binding, miRNA binding, or chromatin state. Of the four variants that we experimentally examined, only rs6590330 differentially binds lysate from B cells. Using mass spectrometry, we found more binding of the transcription factor signal transducer and activator of transcription 1 (STAT1) to DNA near the risk allele of rs6590330 than near the non-risk allele. Immunoblot analysis and chromatin immunoprecipitation of pSTAT1 in B cells heterozygous for rs6590330 confirmed that the risk allele increased binding to the active form of STAT1. Analysis with expression quantitative trait loci indicated that the risk allele of rs6590330 is associated with decreased ETS1 expression in Han Chinese, but not other ancestral cohorts. We propose a model in which the risk allele of rs6590330 is associated with decreased ETS1 expression and increases SLE risk by enhancing the binding of pSTAT1.

Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription.

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Blocking of targeted microRNAs from next-generation sequencing libraries.

Highly abundant microRNAs (miRNAs) in small RNA sequencing libraries make it difficult to obtain efficient measurements of more lowly expressed species. We present a new method that allows for the selective blocking of specific, abundant miRNAs during preparation of sequencing libraries. This technique is specific with little off-target effects and has no impact on the reproducibility of the measurement of non-targeted species. In human plasma samples, we demonstrate that blocking of highly abundant hsa-miR-16-5p leads to improved detection of lowly expressed miRNAs and more precise measurement of differential expression overall. Furthermore, we establish the ability to target a second abundant miRNA and to multiplex the blocking of two miRNAs simultaneously. For small RNA sequencing, this technique could fill a similar role as do ribosomal or globin removal technologies in messenger RNA sequencing.

Characterization of genetic loci that affect susceptibility to inflammatory bowel diseases in African Americans.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) has familial aggregation in African Americans (AAs), but little is known about the molecular genetic susceptibility. Mapping studies using the Immunochip genotyping array expand the number of susceptibility loci for IBD in Caucasians to 163, but the contribution of the 163 loci and European admixture to IBD risk in AAs is unclear. We performed a genetic mapping study using the Immunochip to determine whether IBD susceptibility loci in Caucasians also affect risk in AAs and identify new associated loci. METHODS: We recruited AAs with IBD and without IBD (controls) from 34 IBD centers in the United States; additional controls were collected from 4 other Immunochip studies. Association and admixture loci were mapped for 1088 patients with Crohn's disease, 361 with ulcerative colitis, 62 with IBD type unknown, and 1797 controls; 130,241 autosomal single-nucleotide polymorphisms (SNPs) were analyzed. RESULTS: The strongest associations were observed between ulcerative colitis and HLA rs9271366 (P = 7.5 × 10(-6)), Crohn's disease and 5p13.1 rs4286721 (P = 3.5 × 10(-6)), and IBD and KAT2A rs730086 (P = 2.3 × 10(-6)). Additional suggestive associations (P < 4.2 × 10(-5)) were observed between Crohn's disease and IBD and African-specific SNPs in STAT5A and STAT3; between IBD and SNPs in IL23R, IL12B, and C2orf43; and between ulcerative colitis and SNPs near HDAC11 and near LINC00994. The latter 3 loci have not been previously associated with IBD, but require replication. Established Caucasian associations were replicated in AAs (P < 3.1 × 10(-4)) at NOD2, IL23R, 5p15.3, and IKZF3. Significant admixture (P < 3.9 × 10(-4)) was observed for 17q12-17q21.31 (IZKF3 through STAT3), 10q11.23-10q21.2, 15q22.2-15q23, and 16p12.2-16p12.1. Network analyses showed significant enrichment (false discovery rate <1 × 10(-5)) in genes that encode members of the JAK-STAT, cytokine, and chemokine signaling pathways, as well those involved in pathogenesis of measles. CONCLUSIONS: In a genetic analysis of 3308 AA IBD cases and controls, we found that many variants associated with IBD in Caucasians also showed association evidence with these diseases in AAs; we also found evidence for variants and loci not previously associated with IBD. The complex genetic factors that determine risk for or protection against IBD in different populations require further study.

Preferential association of a functional variant in complement receptor 2 with antibodies to double-stranded DNA.

OBJECTIVES: Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association. METHODS: Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR. RESULTS: The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10(-4), OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10(-7), OR 0.71; case-only pmeta=1.9×10(-4), OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR. CONCLUSIONS: These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.

Decreased SMG7 expression associates with lupus-risk variants and elevated antinuclear antibody production.

OBJECTIVES: Following up the systemic lupus erythematosus (SLE) genome-wide association studies (GWAS) identification of NMNAT2 at rs2022013, we fine-mapped its 150 kb flanking regions containing NMNAT2 and SMG7 in a 15 292 case-control multi-ancestry population and tested functions of identified variants. METHODS: We performed genotyping using custom array, imputation by IMPUTE 2.1.2 and allele specific functions using quantitative real-time PCR and luciferase reporter transfections. SLE peripheral blood mononuclear cells (PBMCs) were cultured with small interfering RNAs to measure antinuclear antibody (ANA) and cyto/chemokine levels in supernatants using ELISA. RESULTS: We confirmed association at NMNAT2 in European American (EA) and Amerindian/Hispanic ancestries, and identified independent signal at SMG7 tagged by rs2702178 in EA only (p=2.4×10-8, OR=1.23 (95% CI 1.14 to 1.32)). In complete linkage disequilibrium with rs2702178, rs2275675 in the promoter region robustly associated with SMG7 mRNA levels in multiple expression quantitative trait locus (eQTL) datasets. Its risk allele was dose-dependently associated with decreased SMG7 mRNA levels in PBMCs of 86 patients with SLE and 119 controls (p=1.1×10-3 and 6.8×10-8, respectively) and conferred reduced transcription activity in transfected HEK-293 (human embryonic kidney cell line) and Raji cells (p=0.0035 and 0.0037, respectively). As a critical component in the nonsense-mediated mRNA decay pathway, SMG7 could regulate autoantigens including ribonucleoprotein (RNP) and Smith (Sm). We showed SMG7 mRNA levels in PBMCs correlated inversely with ANA titres of patients with SLE (r=-0.31, p=0.01), and SMG7 knockdown increased levels of ANA IgG and chemokine (C-C motif) ligand 19 in SLE PBMCs (p=2.0×10-5 and 2.0×10-4, respectively). CONCLUSION: We confirmed NMNAT2 and identified independent SMG7 association with SLE. The inverse relationship between levels of the risk allele-associated SMG7 mRNAs and ANA suggested the novel contribution of mRNA surveillance pathway to SLE pathogenesis.

Genome-wide DNA methylation analysis of systemic lupus erythematosus reveals persistent hypomethylation of interferon genes and compositional changes to CD4+ T-cell populations.

Systemic lupus erythematosus (SLE) is an autoimmune disease with known genetic, epigenetic, and environmental risk factors. To assess the role of DNA methylation in SLE, we collected CD4+ T-cells, CD19+ B-cells, and CD14+ monocytes from 49 SLE patients and 58 controls, and performed genome-wide DNA methylation analysis with Illumina Methylation 450 microarrays. We identified 166 CpGs in B-cells, 97 CpGs in monocytes, and 1,033 CpGs in T-cells with highly significant changes in DNA methylation levels (p < 1 × 10(-8)) among SLE patients. Common to all three cell-types were widespread and severe hypomethylation events near genes involved in interferon signaling (type I). These interferon-related changes were apparent in patients collected during active and quiescent stages of the disease, suggesting that epigenetically-mediated hypersensitivity to interferon persists beyond acute stages of the disease and is independent of circulating interferon levels. This interferon hypersensitivity was apparent in memory, naïve and regulatory T-cells, suggesting that this epigenetic state in lupus patients is established in progenitor cell populations. We also identified a widespread, but lower amplitude shift in methylation in CD4+ T-cells (> 16,000 CpGs at FDR < 1%) near genes involved in cell division and MAPK signaling. These cell type-specific effects are consistent with disease-specific changes in the composition of the CD4+ population and suggest that shifts in the proportion of CD4+ subtypes can be monitored at CpGs with subtype-specific DNA methylation patterns.

Admixture mapping in lupus identifies multiple functional variants within IFIH1 associated with apoptosis, inflammation, and autoantibody production.

Systemic lupus erythematosus (SLE) is an inflammatory autoimmune disease with a strong genetic component. African-Americans (AA) are at increased risk of SLE, but the genetic basis of this risk is largely unknown. To identify causal variants in SLE loci in AA, we performed admixture mapping followed by fine mapping in AA and European-Americans (EA). Through genome-wide admixture mapping in AA, we identified a strong SLE susceptibility locus at 2q22-24 (LOD=6.28), and the admixture signal is associated with the European ancestry (ancestry risk ratio ~1.5). Large-scale genotypic analysis on 19,726 individuals of African and European ancestry revealed three independently associated variants in the IFIH1 gene: an intronic variant, rs13023380 [P(meta) = 5.20×10(-14); odds ratio, 95% confidence interval = 0.82 (0.78-0.87)], and two missense variants, rs1990760 (Ala946Thr) [P(meta) = 3.08×10(-7); 0.88 (0.84-0.93)] and rs10930046 (Arg460His) [P(dom) = 1.16×10(-8); 0.70 (0.62-0.79)]. Both missense variants produced dramatic phenotypic changes in apoptosis and inflammation-related gene expression. We experimentally validated function of the intronic SNP by DNA electrophoresis, protein identification, and in vitro protein binding assays. DNA carrying the intronic risk allele rs13023380 showed reduced binding efficiency to a cellular protein complex including nucleolin and lupus autoantigen Ku70/80, and showed reduced transcriptional activity in vivo. Thus, in SLE patients, genetic susceptibility could create a biochemical imbalance that dysregulates nucleolin, Ku70/80, or other nucleic acid regulatory proteins. This could promote antibody hypermutation and auto-antibody generation, further destabilizing the cellular network. Together with molecular modeling, our results establish a distinct role for IFIH1 in apoptosis, inflammation, and autoantibody production, and explain the molecular basis of these three risk alleles for SLE pathogenesis.

Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8 × 10(-8)), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3' UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1 × 10(-11) in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.