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.

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.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes.

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Does a Multiple Myeloma Polygenic Risk Score Predict Overall Survival of Patients with Myeloma?

BACKGROUND: Genome-wide association studies (GWAS) of multiple myeloma in populations of European ancestry (EA) identified and confirmed 24 susceptibility loci. For other cancers (e.g., colorectum and melanoma), risk loci have also been associated with patient survival. METHODS: We explored the possible association of all the known risk variants and their polygenic risk score (PRS) with multiple myeloma overall survival (OS) in multiple populations of EA [the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the International Lymphoma Epidemiology consortium, CoMMpass, and the German GWAS] for a total of 3,748 multiple myeloma cases. Cox proportional hazards regression was used to assess the association between each risk SNP with OS under the allelic and codominant models of inheritance. All analyses were adjusted for age, sex, country of origin (for IMMEnSE) or principal components (for the others) and disease stage (ISS). SNP associations were meta-analyzed. RESULTS: SNP associations were meta-analyzed. From the meta-analysis, two multiple myeloma risk SNPs were associated with OS (P < 0.05), specifically POT1-AS1-rs2170352 [HR = 1.37; 95% confidence interval (CI) = 1.09-1.73; P = 0.007] and TNFRSF13B-rs4273077 (HR = 1.19; 95% CI = 1.01-1.41; P = 0.04). The association between the combined 24 SNP MM-PRS and OS, however, was not significant. CONCLUSIONS: Overall, our results did not support an association between the majority of multiple myeloma risk SNPs and OS. IMPACT: This is the first study to investigate the association between multiple myeloma PRS and OS in multiple myeloma.

B-Cell NHL Subtype Risk Associated with Autoimmune Conditions and PRS.

BACKGROUND: A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non-Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes. METHODS: In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell-mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression. RESULTS: We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08-1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59-2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell-mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction. CONCLUSIONS: Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions. IMPACT: Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.

Distinct phenotype of neutrophil, monocyte, and eosinophil populations indicates altered myelopoiesis in a subset of patients with multiple myeloma.

Hematologic malignancies, including multiple myeloma (MM), promote systemic immune dysregulation resulting in an alteration and increased plasticity of myeloid cell subsets. To determine the heterogeneity of the myeloid cell compartment in the peripheral blood of patients with MM, we performed a detailed investigation of the phenotype and function of myeloid subpopulations. We report that a subset of MM patients exhibits a specific myeloid cell phenotype indicative of altered myelopoiesis characterized by significant changes in the properties of circulating granulocytic, monocytic, and eosinophilic populations. The subset, referred to as MM2, is defined by a markedly elevated level of CD64 (FcγRI) on the surface of circulating neutrophils. Compared to healthy controls or MM1 patients displaying intermediate levels of CD64, neutrophils from MM2 patients exhibit a less differentiated phenotype, low levels of CD10 and CXC chemokine receptor 2 (CXCR2), increased capacity for the production of mitochondrial reactive oxygen species, and an expansion of CD16neg immature neutrophil subset. Classical and patrolling monocytes from MM2 patients express elevated levels of CD64 and activation markers. MM2 eosinophils display lower levels of C-C Chemokine receptor 3 (CCR3), Toll-like receptor 4 (TLR4, CD284), and tissue factor (TF, CD142). The MM2 (CD64high) phenotype is independent of age, race, sex, and treatment type. Characteristic features of the MM2 (CD64high) phenotype are associated with myeloma-defining events including elevated involved/uninvolved immunoglobulin free light chain (FLC) ratio at diagnosis. Detailed characterization of the altered myeloid phenotype in multiple myeloma will likely facilitate the identification of patients with an increased risk of disease progression and open new avenues for the rational design of novel therapeutic approaches.

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.

PolyHIPE foams from pristine graphene: Strong, porous, and electrically conductive materials templated by a 2D surfactant.

Graphene is attractive as a functional 2D surfactant for polymerized high internal phase emulsions (polyHIPEs) due to its remarkable mechanical and electrical properties. We have developed polyHIPEs stabilized by pristine, unoxidized graphene via the spontaneous exfoliation of graphite at high-energy aqueous/organic interfaces. The exfoliated graphene self-assembles into a percolating network and incorporates into the polyHIPE cell walls as verified by TEM. The resulting composites showed compressive strengths of 7.0 MPa at densities of 0.22 g/cm3 and conductivities up to 0.36 S/m. Systematically reducing the concentration of monomer in the oil phase by dilution with a porogenic-acting solvent increased the porosity and lowered the density of the polyHIPEs. Characterization of these composites indicated that graphene's high compressive strength and modulus was transferred to the polyHIPEs and provided mechanical reinforcement even at low polymer content. SEM showed that the morphology of the polymer changed with decreasing monomer content while the graphene lined cells retained their shape. Moreover, we show that the polyHIPEs contain a continuous graphene percolating network resulting in electrically conductive materials at low graphene loading.

Coinherited genetics of multiple myeloma and its precursor, monoclonal gammopathy of undetermined significance.

So far, 23 germline susceptibility loci have been associated with multiple myeloma (MM) risk. It is unclear whether the genetic variation associated with MM susceptibility also predisposes to its precursor, monoclonal gammopathy of undetermined significance (MGUS). Leveraging 2434 MM cases, 754 MGUS cases, and 2 independent sets of controls (2567/879), we investigated potential shared genetic susceptibility of MM and MGUS by (1) performing MM and MGUS genome-wide association studies (GWAS); (2) validating the association of a polygenic risk score (PRS) based on 23 established MM loci (MM-PRS) with risk of MM, and for the first time with MGUS; and (3) examining genetic correlation of MM and MGUS. Heritability and genetic estimates yielded 17% (standard error [SE] ±0.04) and 15% (SE ±0.11) for MM and MGUS risk, respectively, and a 55% (SE ±0.30) genetic correlation. The MM-PRS was associated with risk of MM when assessed continuously (odds ratio [OR], 1.17 per SD; 95% confidence interval [CI], 1.13-1.21) or categorically (OR, 1.70; 95% CI, 1.38-2.09 for highest; OR, 0.71; 95% CI, 0.55-0.90 for lowest compared with middle quintile). The MM-PRS was similarly associated with MGUS (OR, 1.19 per SD; 95% CI, 1.14-1.26 as a continuous measure, OR, 1.77, 95%CI: 1.29-2.43 for highest and OR, 0.70, 95%CI: 0.50-0.98 for lowest compared with middle quintile). MM and MGUS associations did not differ by age, sex, or MM immunoglobulin isotype. We validated a 23-SNP MM-PRS in an independent series of MM cases and provide evidence for its association with MGUS. Our results suggest shared common genetic susceptibility to MM and MGUS.

Long-term survival of 1338 MM patients treated with tandem autologous vs. autologous-allogeneic transplantation.

Contrary to tandem autologous transplant (auto-auto), autologous followed by reduced intensity conditioning allogenic transplantation (auto-allo) offers graft-versus-myeloma (GVM) effect but with higher toxicity. Trials comparing these two strategies relied on availability of HLA-matched sibling donors for arm allocation (biological randomization) and yielded conflicting results. A pooled analysis of multiple trials with extended follow up provides an opportunity to compare these strategies. We obtained individual patient data from participants of four trials comparing auto-auto vs. auto-allo after induction therapy. There were 899 patients in auto-auto and 439 in auto-allo. Median follow up of survivors was 118.5 months. Median overall survival (OS) was 78.0 months in auto-auto and 98.3 months in auto-allo (HR = 0.84, P = 0.02). OS was 36.4% vs. 44.1% at 10 years (P = 0.01) for auto-auto and auto-allo, respectively. Progression-free survival was also improved in auto-allo (HR = 0.84, P = 0.004). Risk of non-relapse mortality was higher in auto-allo (10 year 8.3% vs. 19.7%, P < 0.001), while risk of disease progression was higher in auto-auto (10 year 77.2% vs. 61.6%, P < 0.001). Median post relapse survival was 41.5 months in auto-auto and 62.3 months in auto-allo (HR = 0.71, P < 0.001). This supports the existence of durable GVM effect enhancing myeloma control with subsequent therapies.

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.

Serum free light chain level at diagnosis in myeloma cast nephropathy-a multicentre study.

Myeloma cast nephropathy (MCN) is a common cause of severe renal impairment in multiple myeloma (MM). The level of free light chain (FLC) that causes MCN varies substantially and there is uncertainty about the threshold level that should be used to inform clinical practice. In a multicentre cohort study of 103 patients with a diagnosis of MM and biopsy-confirmed MCN made between 2002-2014, we report prospectively measured levels of serum FLC at diagnosis obtained using a single nephelometric assay (Freelite®) and we explore the relationship between serum FLC level at diagnosis with renal outcome and patient survival. Using a landmark approach, overall survival (OS) was compared between patients who achieved independence from dialysis compared to those who remained dialysis dependent at 3-month, 6-month, 9-month, and 12-month time points. The median serum FLC level at diagnosis was 7531 mg/L (range 107-114600). Serum creatinine was 535 µmol/L (range 168-2993) and eGFR 7 ml/min/1.73 m2 (range 1-34). Six patients (5.8%) had an FLC level <1500 mg/L, which is the International Myeloma Working Group threshold for MCN and two patients were below the International Kidney and Monoclonal Gammopathy working group threshold of 500 mg/L; one was hypercalcaemic, and one had high-normal serum calcium level and had received a non-steroidal anti-inflammatory agent. Sixty-nine (67%) patients required haemodialysis treatment of whom 36 (52.1%) recovered independent renal function. Sixty-six (64%) patients died with a median OS of 2.5 years (95% CI 1.8-3.3). A landmark analysis revealed that independence from dialysis was associated with improved survival at 3-months (P = 0.003), 6-months (P = 0.035) and 9-months (P = 0.014); there was no survival benefit observed beyond 12 months (P = 0.146). Serum FLC level at diagnosis was neither associated with renal function recovery nor with OS. This is the largest reported cohort of patients with biopsy-confirmed MCN and prospectively measured serum FLC levels. These results indicate that a serum monoclonal FLC > 500 mg/L should be considered the threshold level associated with the development of MCN.

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.

Heparanase promotes myeloma stemness and in vivo tumorigenesis.

Heparanase is known to enhance the progression of many cancer types and is associated with poor patient prognosis. We recently reported that after patients with multiple myeloma were treated with high dose chemotherapy, the tumor cells that emerged upon relapse expressed a much higher level of heparanase than was present prior to therapy. Because tumor cells having stemness properties are thought to seed tumor relapse, we investigated whether heparanase had a role in promoting myeloma stemness. When plated at low density and grown in serum-free conditions that support survival and expansion of stem-like cells, myeloma cells expressing a low level of heparanase formed tumor spheroids poorly. In contrast, cells expressing a high level of heparanase formed significantly more and larger spheroids than did the heparanase low cells. Importantly, heparanase-low expressing cells exhibited plasticity and were induced to exhibit stemness properties when exposed to recombinant heparanase or to exosomes that contained a high level of heparanase cargo. The spheroid-forming heparanase-high cells had elevated expression of GLI1, SOX2 and ALDH1A1, three genes known to be associated with myeloma stemness. Inhibitors that block the heparan sulfate degrading activity of heparanase significantly diminished spheroid formation and expression of stemness genes implying a direct role of the enzyme in regulating stemness. Blocking the NF-κB pathway inhibited spheroid formation and expression of stemness genes demonstrating a role for NF-κB in heparanase-mediated stemness. Myeloma cells made deficient in heparanase exhibited decreased stemness properties in vitro and when injected into mice they formed tumors poorly compared to the robust tumorigenic capacity of cells expressing higher levels of heparanase. These studies reveal for the first time a role for heparanase in promoting cancer stemness and provide new insight into its function in driving tumor progression and its association with poor prognosis in cancer patients.

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.

Two high-risk susceptibility loci at 6p25.3 and 14q32.13 for Waldenström macroglobulinemia.

Waldenström macroglobulinemia (WM)/lymphoplasmacytic lymphoma (LPL) is a rare, chronic B-cell lymphoma with high heritability. We conduct a two-stage genome-wide association study of WM/LPL in 530 unrelated cases and 4362 controls of European ancestry and identify two high-risk loci associated with WM/LPL at 6p25.3 (rs116446171, near EXOC2 and IRF4; OR = 21.14, 95% CI: 14.40-31.03, P = 1.36 × 10-54) and 14q32.13 (rs117410836, near TCL1; OR = 4.90, 95% CI: 3.45-6.96, P = 8.75 × 10-19). Both risk alleles are observed at a low frequency among controls (~2-3%) and occur in excess in affected cases within families. In silico data suggest that rs116446171 may have functional importance, and in functional studies, we demonstrate increased reporter transcription and proliferation in cells transduced with the 6p25.3 risk allele. Although further studies are needed to fully elucidate underlying biological mechanisms, together these loci explain 4% of the familial risk and provide insights into genetic susceptibility to this malignancy.

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.

HLA Class I and II Diversity Contributes to the Etiologic Heterogeneity of Non-Hodgkin Lymphoma Subtypes.

A growing number of loci within the human leukocyte antigen (HLA) region have been implicated in non-Hodgkin lymphoma (NHL) etiology. Here, we test a complementary hypothesis of "heterozygote advantage" regarding the role of HLA and NHL, whereby HLA diversity is beneficial and homozygous HLA loci are associated with increased disease risk. HLA alleles at class I and II loci were imputed from genome-wide association studies (GWAS) using SNP2HLA for 3,617 diffuse large B-cell lymphomas (DLBCL), 2,686 follicular lymphomas (FL), 2,878 chronic lymphocytic leukemia/small lymphocytic lymphomas (CLL/SLL), 741 marginal zone lymphomas (MZL), and 8,753 controls of European descent. Both DLBCL and MZL risk were elevated with homozygosity at class I HLA-B and -C loci (OR DLBCL = 1.31, 95% CI = 1.06-1.60; OR MZL = 1.45, 95% CI = 1.12-1.89) and class II HLA-DRB1 locus (OR DLBCL = 2.10, 95% CI = 1.24-3.55; OR MZL = 2.10, 95% CI = 0.99-4.45). Increased FL risk was observed with the overall increase in number of homozygous HLA class II loci (P trend < 0.0001, FDR = 0.0005). These results support a role for HLA zygosity in NHL etiology and suggests that distinct immune pathways may underly the etiology of the different NHL subtypes.Significance: HLA gene diversity reduces risk for non-Hodgkin lymphoma. Cancer Res; 78(14); 4086-96. ©2018 AACR.