Anthropometric traits and risk of multiple myeloma: differences by race, sex and diagnostic clinical features.

BACKGROUND: Obesity is an established modifiable risk factor for multiple myeloma (MM). However, associations of obesity and MM risk in Black populations, for whom obesity and MM are more common, is less clear. METHODS: Using participants enrolled in the Integrative Molecular And Genetic Epidemiology study, we evaluated the association of anthropometric traits with MM risk overall, stratified by race and sex. Among cases, we assessed the association of BMI with the presence of myeloma-defining events. RESULTS: We observed an 18% increase in MM risk for every 5 kg/m2 increase in usual adult BMI. Participants with severe obesity (BMI ≥ 40 kg/m2) had the highest risk compared to those with a normal usual adult BMI (18.5-24.9 kg/m2; OR = 1.87, 95% CI 1.25-2.80), particularly among Black men (OR = 3.94, 95% CI 0.90-17.36). Furthermore, MM cases with overweight/obesity (BMI ≥ 25 kg/m2) were more likely to present at diagnosis with low renal function (OR = 1.62, 95% CI 1.09-2.40), deletion 13q (OR = 1.73, 95% CI 1.08-2.76) and lytic lesions or compression fractures (OR = 2.39, 95% CI 0.82-7.01) and less likely to present with severe diffuse osteopenia (OR = 0.51, 95% CI 0.31-0.81). CONCLUSIONS: Findings underscore the importance of obesity as a modifiable risk factor for MM, particularly in high-risk populations, and for the clinical presentation of disease.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk.

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

Genome-wide DNA methylation analysis implicates enrichment of interferon pathway in African American patients with Systemic Lupus Erythematosus and European Americans with lupus nephritis.

Systemic Lupus Erythematosus (SLE) is a chronic, multisystem, inflammatory autoimmune disease that disproportionately affects women. Trends in SLE prevalence and clinical course differ by ancestry, with those of African American ancestry presenting with more active, severe and rapidly progressive disease than European Americans. Previous research established altered epigenetic signatures in SLE patients compared to controls. However, the contribution of aberrant DNA methylation (DNAm) to the risk of SLE by ancestry and differences among patients with SLE-associated Lupus Nephritis (LN) has not been well described. We evaluated the DNA methylomes of 87 individuals including 41 SLE patients, with and without LN, and 46 controls enrolled in an ancestry diverse, well-characterized cohort study of established SLE (41 SLE patients [20 SLE-LN+, 21 SLE-LN-] and 46 sex-, race- and age-matched controls; 55% African American, 45% European American). Participants were genotyped using the Infinium Global Diversity Array (GDA), and genetic ancestry was estimated using principal components. Genome-wide DNA methylation was initially measured using the Illumina MethylationEPIC 850K Beadchip array followed by methylation-specific qPCR to validate the methylation status at putative loci. Differentially Methylated Positions (DMP) were identified using a case-control approach adjusted for ancestry. We identified a total of 51 DMPs in CpGs among SLE patients compared to controls. Genes proximal to these CpGs were highly enriched for involvement in type I interferon signaling. DMPs among European American SLE patients with LN were similar to African American SLE patients with and without LN. Our findings were validated using an orthogonal, methyl-specific PCR for three SLE-associated DMPs near or proximal to MX1, USP18, and IFITM1. Our study confirms previous reports that DMPs in CpGs associated with SLE are enriched in type I interferon genes. However, we show that European American SLE patients with LN have similar DNAm patterns to African American SLE patients irrespective of LN, suggesting that aberrant DNAm alters activity of type I interferon pathway leading to more severe disease independent of ancestry.

Expression quantitative trait loci of genes predicting outcome are associated with survival of multiple myeloma patients.

Gene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype-tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10-7 either in EBV-transformed B-lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression-free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4-rs1992292, TBRG4-rs2287535 and ENTPD1-rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta-analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4-rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04-1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

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.

Infectious mononucleosis, immune genotypes, and non-Hodgkin lymphoma (NHL): an InterLymph Consortium study.

PURPOSE: We explored the interaction between non-Hodgkin lymphoma (NHL), infectious mononucleosis (IM) history, and immune-related genotypes in a pooled case-control analysis. METHODS: A total of 7,926 NHL patients and 10,018 controls from 12 case-control studies were included. Studies were conducted during various time periods between 1988 and 2008, and participants were 17-96 years of age at the time of ascertainment/recruitment. Self-reported IM history and immune response genotypes were provided by the InterLymph Data Coordinating Center at Mayo Clinic. Odds ratios (OR) were estimated using multivariate logistic regression, and interactions were estimated using the empirical Bayes method. PACT was used to account for multiple comparisons. RESULTS: There was evidence of an interaction effect between IM history and two variants on T-cell lymphoma (TCL) risk: rs1143627 in interleukin-1B (IL1B) (pinteraction = 0.04, ORinteraction = 0.09, 95% confidence interval [CI] 0.01, 0.87) and rs1800797 in interleukin-6 (IL6) (pinteraction = 0.03, ORinteraction = 0.08, 95% CI 0.01, 0.80). Neither interaction effect withstood adjustment for multiple comparisons. There were no statistically significant interactions between immune response genotypes and IM on other NHL subtypes. CONCLUSIONS: Genetic risk variants in IL1B and IL6 may affect the association between IM and TCL, possibly by influencing T-cell activation, growth, and differentiation in the presence of IM, thereby decreasing risk of immune cell proliferation.

Correction to: Infectious mononucleosis, immune genotypes, and non-Hodgkin lymphoma (NHL): an InterLymph Consortium study.

Unfortunately, the word "Group" is missed in the article title of the original publication. It has been corrected by this erratum.

Pristine Graphene Microspheres by the Spreading and Trapping of Graphene at an Interface.

The interfacial spreading and exfoliation of graphene was used to create low-density, hollow microspheres defined by a thin shell of graphene. The spheres were templated by a thermodynamically driven self-assembly process in which graphite spontaneously exfoliated and spread at the high-energy interfaces of a water-in-oil emulsion. Graphene thus acted as a 2D surfactant to stabilize the dispersed water droplets utilized as polymerization templates. Using a mixture of organic solvent and monomer as the emulsion oil phase, polystyrene-coated hollow graphene microspheres were created. These spheres were characterized by optical and electron microscopy, thermo-gravimetric analysis, nanoindentation, and particle sizing. The mechanism leading to the microsphere surface morphology and shape is discussed, with the oil phase composition shown to play a critical role.

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.

Fibronectin on the Surface of Myeloma Cell-derived Exosomes Mediates Exosome-Cell Interactions.

Exosomes regulate cell behavior by binding to and delivering their cargo to target cells; however, the mechanisms mediating exosome-cell interactions are poorly understood. Heparan sulfates on target cell surfaces can act as receptors for exosome uptake, but the ligand for heparan sulfate on exosomes has not been identified. Using exosomes isolated from myeloma cell lines and from myeloma patients, we identify exosomal fibronectin as a key heparan sulfate-binding ligand and mediator of exosome-cell interactions. We discovered that heparan sulfate plays a dual role in exosome-cell interaction; heparan sulfate on exosomes captures fibronectin, and on target cells it acts as a receptor for fibronectin. Removal of heparan sulfate from the exosome surface releases fibronectin and dramatically inhibits exosome-target cell interaction. Antibody specific for the Hep-II heparin-binding domain of fibronectin blocks exosome interaction with tumor cells or with marrow stromal cells. Regarding exosome function, fibronectin-mediated binding of exosomes to myeloma cells activated p38 and pERK signaling and expression of downstream target genes DKK1 and MMP-9, two molecules that promote myeloma progression. Antibody against fibronectin inhibited the ability of myeloma-derived exosomes to stimulate endothelial cell invasion. Heparin or heparin mimetics including Roneparstat, a modified heparin in phase I trials in myeloma patients, significantly inhibited exosome-cell interactions. These studies provide the first evidence that fibronectin binding to heparan sulfate mediates exosome-cell interactions, revealing a fundamental mechanism important for exosome-mediated cross-talk within tumor microenvironments. Moreover, these results imply that therapeutic disruption of fibronectin-heparan sulfate interactions will negatively impact myeloma tumor growth and progression.

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.

Impact of marital status, insurance status, income, and race/ethnicity on the survival of younger patients diagnosed with multiple myeloma in the United States.

BACKGROUND: Recent advances in the treatment of multiple myeloma (MM) have been associated with improved survival, predominantly among young and white patients. The authors hypothesized that sociodemographic factors, adjusted for race/ethnicity, influence the survival of younger patients with MM. METHODS: Overall survival (OS) data were obtained for individuals included in the Surveillance, Epidemiology, and End Results (SEER-18) program who were diagnosed with MM before the age of 65 years between 2007 and 2012. The sociodemographic variables addressed were marital status, insurance status, median household income, and educational achievement in the county of residence. Race/ethnicity was defined as a self-reported construct including Hispanic (regardless of race), non-Hispanic black, non-Hispanic white, and other. RESULTS: There were 10,161 cases of MM included with a median follow-up of 27 months (range, 0-71 months; 22,179 person-years). Using multivariable Cox proportional hazards analysis, SEER registry; age; male sex; and 3 sociodemographic factors including marital status (other than married), insurance status (uninsured or Medicaid), and county-level income (lowest 2 quartiles), but not race/ethnicity, were found to be associated with an increased risk of death. The 4-year estimated OS rate was 71.1%, 63.2%, 53.4%, and 46.5% (P<.001), respectively, for patients with 0, 1, 2, or 3 adverse sociodemographic factors. Hispanic and non-Hispanic black individuals were found to have more adverse sociodemographic factors and worse OS. However, when the population was stratified by the cumulative number of sociodemographic factors, no consistent association between race/ethnicity and OS was observed after adjustment for confounders. CONCLUSIONS: Sociodemographic factors that potentially affect care, but not race/ethnicity, were found to influence the survival of younger patients with MM. Cancer 2016;122:3183-90. © 2016 American Cancer Society.

Multiple myeloma and family history of lymphohaematopoietic cancers: Results from the International Multiple Myeloma Consortium.

Family clusters of multiple myeloma (MM) suggest disease heritability. Nevertheless, patterns of inheritance and the importance of genetic versus environmental risk factors in MM aetiology remain unclear. We pooled data from eleven case-control studies from the International Multiple Myeloma Consortium to characterize the association of MM risk with having a first-degree relative with a history of a lympho-haematapoietic cancer. Unconditional logistic regression models, adjusted for study, sex, age and education level, were used to estimate associations between MM risk and having a first-degree relative with a history of non-Hodgkin lymphoma, Hodgkin lymphoma, leukaemia or MM. Sex, African American race/ethnicity and age were explored as effect modifiers. A total of 2843 cases and 11 470 controls were included. MM risk was elevated in association with having a first-degree relative with any lympho-haematapoietic cancer (Odds Ratio (OR) = 1·29, 95% Confidence Interval (CI): 1·08-1·55). The association was particularly strong for having a first-degree relative with MM (OR = 1·90, 95% CI: 1·26-2·87), especially among men (OR = 4·13, 95% CI: 2·17-7·85) and African Americans (OR = 5·52, 95% CI: 1·87-16·27).These results support the hypothesis that genetic inheritance plays a role in MM aetiology. Future studies are warranted to characterize interactions of genetic markers with environmental exposures.

Family history of hematologic malignancies and risk of multiple myeloma: differences by race and clinical features.

PURPOSE: Multiple myeloma (MM) is the most common hematologic malignancy affecting Blacks in the USA, with standardized incidence rates that are twofold to threefold higher than Whites. The rationale for the disparity is unclear. METHODS: Using participants enrolled in the Molecular And Genetic Epidemiology study of myeloma (259 MM cases; 461 controls), we examined the risk of MM associated with family history of cancer, differences by race and among cases, defining clinical features. Risk estimates were calculated using odds ratios and corresponding 95% confidence intervals from logistic regression adjusted for confounders. RESULTS: Overall, MM risk in cases with relatives affected with any hematologic malignancy was significantly elevated compared to controls (OR 1.89, 95% CI 1.25-2.86). Myeloma risk associated with a family history of MM was higher than the risk associated with any hematologic malignancy (OR 3.75, 95% CI 1.75-8.05), and the effect was greater for Blacks (OR 20.9, 95% CI 2.59-168) than Whites (OR 2.04, 95% 0.83-5.04), among cases with early onset (≤60 years; OR 4.58, 95% CI 1.21-17.3) and with increasing numbers of affected relatives (p trend = 0.001). Overall, frequencies of end organ damage differed in cases with relatives affected with any hematologic malignancy and significantly more cases exhibited κ light chain restriction (OR 3.23, 95% CI 1.13-9.26). CONCLUSIONS: The excess risk of MM observed in Blacks and the variation in clinical features observed in MM patients according to family history of hematologic malignancy may be attributed to a shared germline and environmental susceptibility.

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.

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.

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.

MicroRNA-3148 modulates allelic expression of toll-like receptor 7 variant associated with systemic lupus erythematosus.

We previously reported that the G allele of rs3853839 at 3'untranslated region (UTR) of Toll-like receptor 7 (TLR7) was associated with elevated transcript expression and increased risk for systemic lupus erythematosus (SLE) in 9,274 Eastern Asians [P = 6.5×10(-10), odds ratio (OR) (95%CI) = 1.27 (1.17-1.36)]. Here, we conducted trans-ancestral fine-mapping in 13,339 subjects including European Americans, African Americans, and Amerindian/Hispanics and confirmed rs3853839 as the only variant within the TLR7-TLR8 region exhibiting consistent and independent association with SLE (Pmeta = 7.5×10(-11), OR = 1.24 [1.18-1.34]). The risk G allele was associated with significantly increased levels of TLR7 mRNA and protein in peripheral blood mononuclear cells (PBMCs) and elevated luciferase activity of reporter gene in transfected cells. TLR7 3'UTR sequence bearing the non-risk C allele of rs3853839 matches a predicted binding site of microRNA-3148 (miR-3148), suggesting that this microRNA may regulate TLR7 expression. Indeed, miR-3148 levels were inversely correlated with TLR7 transcript levels in PBMCs from SLE patients and controls (R(2) = 0.255, P = 0.001). Overexpression of miR-3148 in HEK-293 cells led to significant dose-dependent decrease in luciferase activity for construct driven by TLR7 3'UTR segment bearing the C allele (P = 0.0003). Compared with the G-allele construct, the C-allele construct showed greater than two-fold reduction of luciferase activity in the presence of miR-3148. Reduced modulation by miR-3148 conferred slower degradation of the risk G-allele containing TLR7 transcripts, resulting in elevated levels of gene products. These data establish rs3853839 of TLR7 as a shared risk variant of SLE in 22,613 subjects of Asian, EA, AA, and Amerindian/Hispanic ancestries (Pmeta  = 2.0×10(-19), OR = 1.25 [1.20-1.32]), which confers allelic effect on transcript turnover via differential binding to the epigenetic factor miR-3148.

Preferential binding to Elk-1 by SLE-associated IL10 risk allele upregulates IL10 expression.

Immunoregulatory cytokine interleukin-10 (IL-10) is elevated in sera from patients with systemic lupus erythematosus (SLE) correlating with disease activity. The established association of IL10 with SLE and other autoimmune diseases led us to fine map causal variant(s) and to explore underlying mechanisms. We assessed 19 tag SNPs, covering the IL10 gene cluster including IL19, IL20 and IL24, for association with SLE in 15,533 case and control subjects from four ancestries. The previously reported IL10 variant, rs3024505 located at 1 kb downstream of IL10, exhibited the strongest association signal and was confirmed for association with SLE in European American (EA) (P = 2.7×10⁻8, OR = 1.30), but not in non-EA ancestries. SNP imputation conducted in EA dataset identified three additional SLE-associated SNPs tagged by rs3024505 (rs3122605, rs3024493 and rs3024495 located at 9.2 kb upstream, intron 3 and 4 of IL10, respectively), and SLE-risk alleles of these SNPs were dose-dependently associated with elevated levels of IL10 mRNA in PBMCs and circulating IL-10 protein in SLE patients and controls. Using nuclear extracts of peripheral blood cells from SLE patients for electrophoretic mobility shift assays, we identified specific binding of transcription factor Elk-1 to oligodeoxynucleotides containing the risk (G) allele of rs3122605, suggesting rs3122605 as the most likely causal variant regulating IL10 expression. Elk-1 is known to be activated by phosphorylation and nuclear localization to induce transcription. Of interest, phosphorylated Elk-1 (p-Elk-1) detected only in nuclear extracts of SLE PBMCs appeared to increase with disease activity. Co-expression levels of p-Elk-1 and IL-10 were elevated in SLE T, B cells and monocytes, associated with increased disease activity in SLE B cells, and were best downregulated by ERK inhibitor. Taken together, our data suggest that preferential binding of activated Elk-1 to the IL10 rs3122605-G allele upregulates IL10 expression and confers increased risk for SLE in European Americans.