Variation at 3p24.1 and 6q23.3 influences the risk of Hodgkin's lymphoma

Integration of transcriptome-wide association study and gene-based association analysis identifies candidate genes for Hodgkin lymphoma

BACKGROUND

Genome-wide association studies (GWASs) have pinpointed many susceptibility loci for Hodgkin Lymphoma (HL), but their underlying biological mechanisms remain unclear.

METHODS

Utilizing GWAS data from the UK Biobank and FinnGen, along with expression quantitative trait loci (eQTL) statistics from the Genotype-Tissue Expression (GTEx) and the eQTL Catalogue, we carried out a large-scale gene-level association study using Omnibus Transcriptome Test with Expression Reference Summary data (OTTERS), and gene-based analysis with eQTL Multi-marker Analysis of Genomic Annotation (E-MAGMA).

RESULTS

We identified sixteen susceptibility genes for HL (FDR < 0.01), primarily immune-related, including HLA-DQA1, HLA-DQA2, HLA-DQB1, HLA-DRB1, HLA-DRB5, HLA-DMA, and HLA-DPB1, alongside genes involved in apoptosis, RNA processing, transcriptional regulation, and signal transduction. We identified five novel plausible genes, including HLA-DMA, HLA-DPB1, LSM2, AAR2, and NOTCH4.

CONCLUSION

These findings highlight the role of the exogenous antigen presentation pathway in HL, shedding light on potential mechanisms.

Establishing a robust triangulation framework to explore the relationship between hearing loss and Parkinson's disease

The relationship between hearing loss (HL) and Parkinson's disease (PD) remains unclear. Using individual-level and summary-level data from the UK Biobank and the largest genome-wide association studies, we examined this link through observational, Mendelian randomization and genetic pleiotropy analyses. Among 158,229 participants, PD risk rose with HL severity especially in elder and males, and hearing aids significantly reduced PD risk in males. Although our results did not support a causal association, genetic correlation analysis suggested a localized genetic overlap (17q21.31). We identified 1545 SNPs and 63 genes with pleiotropic effects on HL and PD, including 79 novel SNPs across 6 loci, with 3 showing strong co-localization. These loci were enriched in key tissues like brain, heart, liver and pancreas, linked to the dihydrolipoyl dehydrogenase complex pathway, and targeted by drugs such as Warfarin and Phenprocoumon. Overall, this study reveals the risk association, genetic basis, and pleiotropic loci connecting HL and PD.

Genetic drivers and cellular selection of female mosaic X chromosome loss

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Deciphering the Genetic Complexity of Classical Hodgkin Lymphoma: Insights and Effective Strategies

Understanding the genetics of susceptibility to classical Hodgkin lymphoma (cHL) is considerably limited compared to other cancers due to the rare Hodgkin and Reed-Sternberg (HRS) tumor cells, which coexist with the predominant non-malignant microenvironment. This article offers insights into genetic abnormalities in cHL, as well as nucleotide variants and their associated target genes, elucidated through recent technological advancements. Oncogenomes in HRS cells highlight the survival and proliferation of these cells through hyperactive signaling in specific pathways (e.g., NF-kB) and their interplay with microenvironmental cells (e.g., CD4+ T cells). In contrast, the susceptibility genes identified from genome-wide association studies and expression quantitative trait locus analyses only vaguely implicate their potential roles in susceptibility to more general cancers. To pave the way for the era of precision oncology, more intensive efforts are imperative, employing the following strategies: exploring genetic heterogeneity by gender and cHL subtype, investigating colocalization with various types of expression quantitative trait loci, and leveraging single-cell analysis. These approaches provide valuable perspectives for unraveling the genetic complexities of cHL.

Classic Hodgkin lymphoma in young people

Classic Hodgkin lymphoma (CHL) is a unique form of lymphoid cancer featuring a heterogeneous tumor microenvironment and a relative paucity of malignant Hodgkin and Reed-Sternberg (HRS) cells with characteristic phenotype. Younger individuals (children, adolescents and young adults) are affected as often as the elderly, producing a peculiar bimodal age-incidence profile that has generated immense interest in this disease and its origins. Decades of epidemiological investigations have documented the populations most susceptible and identified multiple risk factors that can be broadly categorized as either biological or environmental in nature. Most risk factors result in overt immunodeficiency or confer more subtle alterations to baseline health, physiology or immune function. Epstein Barr virus, however, is both a risk factor and well-established driver of lymphomagenesis in a significant subset of cases. Epigenetic changes, along with the accumulation of somatic driver mutations and cytogenetic abnormalities are required for the malignant transformation of germinal center-experienced HRS cell precursors. Chromosomal instability and the influence of endogenous mutational processes are critical in this regard, by impacting genes involved in key signaling pathways that promote the survival and proliferation of HRS cells and their escape from immune destruction. Here we review the principal features, known risk factors and lymphomagenic mechanisms relevant to newly diagnosed CHL, with an emphasis on those most applicable to young people.

Improve the model of disease subtype heterogeneity by leveraging external summary data

Researchers are often interested in understanding the disease subtype heterogeneity by testing whether a risk exposure has the same level of effect on different disease subtypes. The polytomous logistic regression (PLR) model provides a flexible tool for such an evaluation. Disease subtype heterogeneity can also be investigated with a case-only study that uses a case-case comparison procedure to directly assess the difference between risk effects on two disease subtypes. Motivated by a large consortium project on the genetic basis of non-Hodgkin lymphoma (NHL) subtypes, we develop PolyGIM, a procedure to fit the PLR model by integrating individual-level data with summary data extracted from multiple studies under different designs. The summary data consist of coefficient estimates from working logistic regression models established by external studies. Examples of the working model include the case-case comparison model and the case-control comparison model, which compares the control group with a subtype group or a broad disease group formed by merging several subtypes. PolyGIM efficiently evaluates risk effects and provides a powerful test for disease subtype heterogeneity in situations when only summary data, instead of individual-level data, is available from external studies due to various informatics and privacy constraints. We investigate the theoretic properties of PolyGIM and use simulation studies to demonstrate its advantages. Using data from eight genome-wide association studies within the NHL consortium, we apply it to study the effect of the polygenic risk score defined by a lymphoid malignancy on the risks of four NHL subtypes. These results show that PolyGIM can be a valuable tool for pooling data from multiple sources for a more coherent evaluation of disease subtype heterogeneity.

Identification and Interpretation of eQTL and eGenes for Hodgkin Lymphoma Susceptibility

Genome-wide association studies (GWAS) have revealed approximately 100 genomic signals associated with Hodgkin lymphoma (HL); however, their target genes and underlying mechanisms causing HL susceptibility remain unclear. In this study, transcriptome-wide analysis of expression quantitative trait loci (eQTL) was conducted to identify target genes associated with HL GWAS signals. A mixed model, which explains polygenic regulatory effects by the genomic covariance among individuals, was implemented to discover expression genes (eGenes) using genotype data from 462 European/African individuals. Overall, 80 eGenes were identified to be associated with 20 HL GWAS signals. Enrichment analysis identified apoptosis, immune responses, and cytoskeletal processes as functions of these eGenes. The eGene of rs27524 encodes ERAP1 that can cleave peptides attached to human leukocyte antigen in immune responses; its minor allele may help Reed-Sternberg cells to escape the immune response. The eGene of rs7745098 encodes ALDH8A1 that can oxidize the precursor of acetyl-CoA for the production of ATP; its minor allele may increase oxidization activity to evade apoptosis of pre-apoptotic germinal center B cells. Thus, these minor alleles may be genetic risk factors for HL susceptibility. Experimental studies on genetic risk factors are needed to elucidate the underlying mechanisms of HL susceptibility and improve the accuracy of precision oncology.

An update on molecular features and therapeutic perspectives of pediatric classical Hodgkin Lymphoma. What the clinician needs to know?

Our understanding of Hodgkin lymphoma (HL) molecular biology has been radically transformed over recent years due to the advent and the spreading of the new generation sequencing approaches. These advances offer new insights about genetic predisposition to HL in children and are currently being translated into promising and more selective drugs (brentuximab and checkpoint inhibitors) offering the perspective to reduce treatment-related toxicity. Thus, as more than 90% of pediatric patients are cured after the first line treatment, a major emphasis is placed on survivorship by reducing treatment intensity, in particular, the use of radiotherapy and chemotherapy associated with long-term toxicities. The purposes of this review are to summarize the recent advances performed in the field of molecular biology of HL, in particular the promising development of liquid biopsies. We also provide an update review of immunodeficiencies associated to HL in children recently identified. Finally, we report the recent studies supporting the efficacy of new targeted therapeutics in adult and pediatric cHL (anti-CD30 and anti-PD1).

Genetic associations with lymphomas in Polish patients: A pooled-DNA genome-wide association analysis

Several single nucleotide polymorphisms (SNPs) associated with susceptibility to Hodgkin lymphoma (HL) and diffuse large B-cell lymphoma (DLBCL) have been identified. The aim of this study was to identify susceptibility loci for HL and DLBCL in Polish patients. Altogether, DLBCL (n = 218 and HL patients (n = 224) and healthy individuals (n = 1181) were recruited. Lymphoma diagnosis was based on standard criteria. Genome-wide association study (GWAS) was performed using pooled-DNA samples on llumina Infinium Omni2.5 Exome-8 v1.3, and selected loci were replicated by TaqMan SNP genotyping of individuals. GWAS detected thirteen and seven SNPs associated with DLBCL and HL, respectively. In the replication study, six and seven SNPs reached significance after correction for multiple testing in the DLBCL and HL cohorts, respectively. One and four SNPs associated with DLBCL and HL, respectively, were localized within, and two SNPs-near the major histocompatibility complex (MHC) region. In conclusion, the majority of loci associated with HL and DLBCL aetiology in previous studies have potential roles in immune function. Our pooled-DNA GWAS enabled the identification of several susceptibility loci for DLBCL and HL in the Polish population; some of them were mapped within or adjacent to the MHC, and other associated SNPs were located outside the MHC.

Association of Single-Nucleotide Variants in the Human Leukocyte Antigen and Other Loci With Childhood Hodgkin Lymphoma

IMPORTANCE

Studies focusing on genetic susceptibility of childhood Hodgkin lymphoma (HL) are limited.

OBJECTIVES

To identify genetic variants associated with childhood-onset HL vs adult-onset HL.

DESIGN, SETTING, AND PARTICIPANTS

This genetic association study was performed with 3 cohorts: the St Jude Lifetime Cohort Study (SJLIFE), initiated in 2007 with ongoing follow-up, and the original and expansion cohorts of the Childhood Cancer Survivor Study (CCSS), initiated in the 1990s with ongoing follow-up. Results of these genome-wide association studies (GWASs) were combined via meta-analysis. Data were analyzed from June 2021 to June 2022.

MAIN OUTCOMES AND MEASURES

Childhood HL was the focused outcome. Single-nucleotide variant (SNV, formerly single-nucleotide polymorphism) array genotyping and imputation were conducted for the CCSS original cohort, and whole-genome sequencing was performed for the SJLIFE and CCSS expansion cohort.

RESULTS

A total of 1286 HL cases (mean diagnosis [SD] age, 14.6 [3.9] years), 6193 non-HL childhood cancer cases, and 369 noncancer controls, all of European ancestry, were included in the analysis. Using step-wise conditional logistic regression, the odds ratios (ORs) for each of the 3 independent SNVs identified in the human leukocyte antigen (HLA) locus were 1.80 (95% CI, 1.59-2.03; P = 2.14 × 10-21) for rs28383311, 1.53 (95% CI, 1.37-1.70; P = 2.05 × 10-14) for rs3129198, and 1.51 (95% CI, 1.35-1.69; P = 6.21 × 10-13) for rs3129890. Further HLA imputation revealed 9 alleles and 55 amino acid changes that potentially conferred HL susceptibility. In addition, 5 non-HLA loci were identified: (1) rs1432297 (OR, 1.29; 95% CI, 1.18-1.41; P = 2.5 × 10-8; r2 = 0.55; D' = 0.75 with previously reported rs1432295, REL); (2) rs2757647 (OR, 1.30; 95% CI, 1.18-1.42; P = 3.5 × 10-8; r2 = 0.59; D' = 0.83 with previously reported rs6928977, AHI1); (3) rs13279159 (OR, 1.33; 95% CI, 1.20-1.47; P = 1.7 × 10-8; r2 = 0.75; D' = 1.00 with previously reported rs2019960, PVT1); (4) rs3824662 (OR, 1.52; 95% CI, 1.33-1.73; P = 3.9 × 10-10; r2 = 0.91; D' = 1.00 with previously reported rs3781093, GATA3); and (5) rs117953624 (OR, 1.98; 95% CI, 1.56-2.51; P = 1.5 × 10-8; minor allele frequency, 0.02), a novel uncommon SNV mapped to PDGFD. Twelve of 18 previously reported genome-wide significant non-HLA SNVs (67%) were replicated with statistically significant results.

CONCLUSIONS AND RELEVANCE

In this genetic association study, a predominantly common and potentially unique genetic etiology was found between childhood-onset and adulthood-onset HL.

CD4+ Cytotoxic T Cells Involved in the Development of EBV-Associated Diseases

Activated cytotoxic CD4 T cells (HLA-DR+) play an important role in the control of EBV infection, especially in cells with latency I (EBNA-1). One of the evasion mechanisms of these latency cells is generated by gp42, which, via peripherally binding to the β1 domain of the β chain of MHC class II (HLA-DQ, -DR, and -DP) of the infected B lymphocyte, can block/alter the HLA class II/T-cell receptor (TCR) interaction, and confer an increased level of susceptibility towards the development of EBV-associated autoimmune diseases or cancer in genetically predisposed individuals (HLA-DRB1* and DQB1* alleles). The main developments predisposing the factors of these diseases are: EBV infection; HLA class II risk alleles; sex; and tissue that is infiltrated with EBV-latent cells, forming ectopic lymphoid structures. Therefore, there is a need to identify treatments for eliminating cells with EBV latency, because the current treatments (e.g., antivirals and rituximab) are ineffective.

HLA Expression in Relation to HLA Type in Classic Hodgkin Lymphoma Patients

Simple Summary

Classic Hodgkin lymphoma (cHL) is a B-cell malignancy with involvement of Epstein–Barr virus (EBV) in about 30% of the European population. The risk to develop cHL is strongly linked to genetic variants in the human leukocyte antigen (HLA) genomic region and to certain HLA alleles. This may be caused by the function of HLA alleles, or by genetic linkage to non-HLA genes. HLA can present EBV-derived and tumour-cell specific antigens and this may lead to anti-tumour immune responses. However, the tumour cells downregulate HLA expression in a proportion of the cases, which may result in immune escape. In this study, we tested whether the loss of HLA expression is related to the presence of certain protective HLA alleles. We found that loss and retention of HLA expression is indeed associated with presence of known susceptibility HLA alleles. These findings suggest that HLA itself is involved in development of cHL.

Abstract

Several human leukocyte antigen (HLA) alleles are strongly associated with susceptibility to classic Hodgkin lymphoma (cHL), also in subgroups stratified for presence of the Epstein–Barr virus (EBV). We tested the hypothesis that the pressure on cHL tumour cells to lose HLA expression is associated with HLA susceptibility alleles. A meta-analysis was carried out to identify consistent protective and risk HLA alleles in a combined cohort of 839 cHL patients from the Netherlands and the United Kingdom. Tumour cell HLA expression was studied in 338 cHL cases from these two cohorts and correlated to the presence of specific susceptibility HLA alleles. Carriers of the HLA-DRB1*07 protective allele frequently lost HLA class II expression in cHL overall. Patients carrying the HLA-DRB1*15/16 (DR2) risk allele retained HLA class II expression in EBV− cHL and patients with the HLA-B*37 risk allele retained HLA class I expression more frequently than non-carriers in EBV+ cHL. The other susceptibility alleles showed no significant differences in expression. Thus, HLA expression by tumour cells is associated with a subset of the protective and risk alleles. This strongly suggests that HLA associations in cHL are related to peptide binding capacities of specific HLA alleles.

Genetic variation near CXCL12 is associated with susceptibility to HIV-related non-Hodgkin lymphoma

Human immunodeficiency virus (HIV) infection is associated with an increased risk of non-Hodgkin lymphoma (NHL). Even in the era of suppressive antiretroviral treatment, HIV-infected individuals remain at higher risk of developing NHL compared to the general population. In order to identify potential genetic risk loci, we performed case-control genome-wide association studies and a meta-analysis across three cohorts of HIV-infected patients of European ancestry, including a total of 278 cases and 1,924 matched controls. We observed a significant association with NHL susceptibility in the C-X-C motif chemokine ligand 12 (CXCL12) region on chromosome 10. A fine mapping analysis identified rs7919208 as the most likely causal variant (P=4.77e-11), with the G>A polymorphism creating a new transcription factor binding site for BATF and JUND. These results suggest a modulatory role of CXCL12 regulation in the increased susceptibility to NHL observed in the HIV-infected population.

EOMES and IL-10 regulate antitumor activity of T regulatory type 1 CD4+ T cells in chronic lymphocytic leukemia

The transcription factor eomesodermin (EOMES) promotes interleukin (IL)-10 expression in CD4+ T cells, which has been linked to immunosuppressive and cytotoxic activities. We detected cytotoxic, programmed cell death protein-1 (PD-1) and EOMES co-expressing CD4+ T cells in lymph nodes (LNs) of patients with chronic lymphocytic leukemia (CLL) or diffuse large B-cell lymphoma. Transcriptome and flow cytometry analyses revealed that EOMES does not only drive IL-10 expression, but rather controls a unique transcriptional signature in CD4+ T cells, that is enriched in genes typical for T regulatory type 1 (TR1) cells. The TR1 cell identity of these CD4+ T cells was supported by their expression of interferon gamma and IL-10, as well as inhibitory receptors including PD-1. TR1 cells with cytotoxic capacity accumulate also in Eµ-TCL1 mice that develop CLL-like disease. Whereas wild-type CD4+ T cells control TCL1 leukemia development after adoptive transfer in leukopenic Rag2-/- mice, EOMES-deficient CD4+ T cells failed to do so. We further show that TR1 cell-mediated control of TCL1 leukemia requires IL-10 receptor (IL-10R) signaling, as Il10rb-deficient CD4+ T cells showed impaired antileukemia activity. Altogether, our data demonstrate that EOMES is indispensable for the development of IL-10-expressing, cytotoxic TR1 cells, which accumulate in LNs of CLL patients and control TCL1 leukemia in mice in an IL-10R-dependent manner.

Functional multigenic variations associated with hodgkin lymphoma

INTRODUCTION

The current study aimed to describe genotypes associated with Hodgkin lymphoma (HL) in a cohort of Saudi and non-Saudi patients and discuss their possible susceptibility to HL.

METHODS

We studied clinical, histopathological, and laboratory findings of HL patients admitted over 12 years duration, at King Fahd University Hospital, KSA. The genomic DNAs of HL patients (n = 61) and normal control subjects (n = 36) were extracted, and genotyping was performed using the Illumina human exome bead chip. Set of HL patients and set of normal controls were included in this study.

RESULTS

A total of 35 DNA variants were found to be highly significant with the P-value <9.90 × 10^-11 among 243 345 exonic biomarkers and obeying the Hardy-Weinberg equilibrium. Nine, MEGF11-rs150945752 (P-value 1.20 × 10^-12 ), CACNA1I- s58055559 (P-value 1.93 × 10^-12 ), DECR2-rs146760080 (P-value 2.19 × 10^-12 ), STAB1-rs143894786 (P-value 2.45 × 10^-12 ), ZNF526-rs144433879 (P-value 2.76 × 10^-12 ), CPLANE1-rs200612080 (P-value 3.77 × 10^-12 ), DLK1-rs1058009 (P-value 5.95 × 10^-12 ), RTN4RL2-rs61745214 (P-value 7.71 × 10^-12 ), and PGRMC1-rs145582672 (P-value 8.56 × 10^-12 ), exonic variants on chromosomes 15, 22, and 16 were highly associated with HL cases. THE HIGHLY SIGNIFICANT HAPLOTYPES AT CHROMOSOME 3: rs143894786G; rs149982219G with P-value = 3.43 × 10^-14 was found to be the risk haplotype for the HL patients. The opposite alleles at chromosome 3: rs143894786A; rs149982219G is protective with P-value = 2.46 × 10^-12 . Maximum number of SNPs at the chromosome 19: rs144433879C; rs181265966G; rs201144421C; rs145591797G; rs200560875G; rs77270337G (risk P-value = 2.24 × 10^-12 ) and its opposite allele rs144433879A; rs181265966A; rs201144421T; rs145591797A; rs200560875A; rs77270337A (protective P-value = 2.60 × 10^-9 ) were found to be associated haplotype with the HL and controls, respectively, in Saudi population.

CONCLUSION

Our study concludes that the HL is genetically heterogeneous with multigene causation.

Sequencing at lymphoid neoplasm susceptibility loci maps six myeloma risk genes

Inherited genetic risk factors play a role in multiple myeloma (MM), yet considerable missing heritability exists. Rare risk variants at genome-wide association study (GWAS) loci are a new avenue to explore. Pleiotropy between lymphoid neoplasms (LNs) has been suggested in family history and genetic studies, but no studies have interrogated sequencing for pleiotropic genes or rare risk variants. Sequencing genetically enriched cases can help discover rarer variants. We analyzed exome sequencing in familial or early-onset MM cases to identify rare, functionally relevant variants near GWAS loci for a range of LNs. A total of 149 distinct and significant LN GWAS loci have been published. We identified six recurrent, rare, potentially deleterious variants within 5 kb of significant GWAS single nucleotide polymorphisms in 75 MM cases. Mutations were observed in BTNL2, EOMES, TNFRSF13B, IRF8, ACOXL and TSPAN32. All six genes replicated in an independent set of 255 early-onset MM or familial MM or precursor cases. Expansion of our analyses to the full length of these six genes resulted in a list of 39 rare and deleterious variants, seven of which segregated in MM families. Three genes also had significant rare variant burden in 733 sporadic MM cases compared with 935 control individuals: IRF8 (P = 1.0 × 10-6), EOMES (P = 6.0 × 10-6) and BTNL2 (P = 2.1 × 10-3). Together, our results implicate six genes in MM risk, provide support for genetic pleiotropy between LN subtypes and demonstrate the utility of sequencing genetically enriched cases to identify functionally relevant variants near GWAS loci.

EOMES is essential for antitumor activity of CD8+ T cells in chronic lymphocytic leukemia

Genome-wide association studies identified a single-nucleotide polymorphism (SNP) affecting the transcription factor Eomesodermin (EOMES) associated with a significantly increased risk to develop chronic lymphocytic leukemia (CLL). Epigenetic analyses, RNA sequencing, and flow cytometry revealed that EOMES is not expressed in CLL cells, but in CD8⁺ T cells for which EOMES is a known master regulator. We thus hypothesized that the increased CLL risk associated with the EOMES SNP might be explained by its negative impact on CD8⁺ T-cell-mediated immune control of CLL. Flow cytometry analyses revealed a higher EOMES expression in CD8⁺ T cells of CLL patients compared to healthy individuals, and an accumulation of PD-1⁺ EOMES⁺ CD8⁺ T cells in lymph nodes rather than blood or bone marrow in CLL. This was in line with an observed expansion of EOMES⁺ CD8⁺ T cells in the spleen of leukemic Eµ-TCL1 mice. As EOMES expression was highest in CD8⁺ T cells that express inhibitory receptors, an involvement of EOMES in T-cell exhaustion and dysfunction seems likely. Interestingly, Eomes-deficiency in CD8⁺ T cells resulted in their impaired expansion associated with decreased CLL control in mice. Overall, these observations suggest that EOMES is essential for CD8⁺ T-cell expansion and/or maintenance, and therefore involved in adaptive immune control of CLL.

Interaction between ERAP Alleles and HLA Class I Types Support a Role of Antigen Presentation in Hodgkin Lymphoma Development

Simple Summary

Hodgkin lymphoma (HL) is a common lymphoma in young adults derived from B cells. Emerging evidence suggests that antigen presentation by the malignant B cells is critically involved in HL pathogenesis. In fact, genetic variants of the antigen presenting Human Leukocyte Antigens (HLA) are strongly associated with HL susceptibility. Interestingly, the endoplasmic reticulum aminopeptidase (ERAP)1 and ERAP2 genes, that code for enzymes that process antigens, also appear to be associated. In this study, we show that genetic variants of ERAP genes strongly affect expression levels of ERAP1 and ERAP2. In addition, we find that certain ERAP variants interact with specific HLA class I types in HL patients. This suggests that mechanisms that determine the repertoire of antigens that are presented to the immune system, affect the chance of developing HL. Our findings therefore support a prominent role of antigen presentation in HL susceptibility.

Abstract

Genetic variants in the HLA region are the strongest risk factors for developing Hodgkin lymphoma (HL), suggesting an important role for antigen presentation. This is supported by another HL-associated genomic region which contains the loci of two enzymes that process endogenous proteins to peptides to be presented by HLA class I, i.e., endoplasmic reticulum aminopeptidase 1 (ERAP1) and ERAP2. We hypothesized that ERAP and HLA class I type interact in HL susceptibility, as shown previously for several autoimmune diseases. We detected ERAP1 and ERAP2 expression in tumor cells and cells in the microenvironment in primary HL tissue samples. Seven ERAP SNPs and ERAP1 haplotypes showed strong associations with RNA and protein levels of ERAP1 and ERAP2 in LCLs and HL cell lines. Analysis of HLA class I types, ERAP SNPs and ERAP haplotypes by direct genotyping or imputation from genome-wide association data in 390 HL patients revealed significant interactions between HLA-A11, rs27038 and the rs27038 associated ERAP haplotype, as well as between HLA-Cw2 and rs26618. In conclusion, our results show that ERAP and HLA class I interact in genetic susceptibility to HL, providing further evidence that antigen presentation is an important process in HL susceptibility and pathogenesis.

Identification of Familial Hodgkin Lymphoma Predisposing Genes Using Whole Genome Sequencing

Hodgkin lymphoma (HL) is a lymphoproliferative malignancy of B-cell origin that accounts for 10% of all lymphomas. Despite evidence suggesting strong familial clustering of HL, there is no clear understanding of the contribution of genes predisposing to HL. In this study, whole genome sequencing (WGS) was performed on 7 affected and 9 unaffected family members from three HL-prone families and variants were prioritized using our Familial Cancer Variant Prioritization Pipeline (FCVPPv2). WGS identified a total of 98,564, 170,550, and 113,654 variants which were reduced by pedigree-based filtering to 18,158, 465, and 26,465 in families I, II, and III, respectively. In addition to variants affecting amino acid sequences, variants in promoters, enhancers, transcription factors binding sites, and microRNA seed sequences were identified from upstream, downstream, 5′ and 3′ untranslated regions. A panel of 565 cancer predisposing and other cancer-related genes and of 2,383 potential candidate HL genes were also screened in these families to aid further prioritization. Pathway analysis of segregating genes with Combined Annotation Dependent Depletion Tool (CADD) scores >20 was performed using Ingenuity Pathway Analysis software which implicated several candidate genes in pathways involved in B-cell activation and proliferation and in the network of “Cancer, Hematological disease and Immunological Disease.” We used the FCVPPv2 for further in silico analyses and prioritized 45 coding and 79 non-coding variants from the three families. Further literature-based analysis allowed us to constrict this list to one rare germline variant each in families I and II and two in family III. Functional studies were conducted on the candidate from family I in a previous study, resulting in the identification and functional validation of a novel heterozygous missense variant in the tumor suppressor gene DICER1 as potential HL predisposition factor. We aim to identify the individual genes responsible for predisposition in the remaining two families and will functionally validate these in further studies.

Risk of cancer among multiple sclerosis patients, siblings, and population controls: A prospective cohort study

Background:

Risk of cancer in multiple sclerosis (MS) patients compared to their siblings is unknown.

Objective:

The objective was to prospectively investigate the risk of cancer among MS patients compared to siblings without MS and to population controls.

Methods:

We retrieved data on MS patients born between 1930 and 1979 from the Norwegian Multiple Sclerosis Registry and population studies and on cancer diagnosis from the Cancer Registry of Norway. We used adjusted Cox proportional hazard regression to estimate cancer risk among 6883 MS patients, 8918 siblings without MS, and 37,919 population controls.

Results:

During 65 years of follow-up, cancer risk among MS patients was higher than that among population controls (hazard ratio (HR) = 1.14, 95% confidence interval (CI): 1.05–1.23) in respiratory organs (HR = 1.66, 95% CI: 1.26–2.19), urinary organs (HR = 1.51, 95% CI: 1.12–2.04), and the central nervous system (HR = 1.52, 95% CI: 1.11–2. 09). Siblings had higher risk of hematological cancers compared with MS patients (HR = 1.82, 95% CI: 1.21–2.73) and population controls (HR = 1.72, 95% CI: 1.36–2.18).

Conclusion:

MS patients were associated with increased risk of cancer compared to population controls. Siblings had increased risk of hematological cancer. This indicates that MS and hematological cancer could share a common etiology.

The Unsolved Puzzle of c-Rel in B Cell Lymphoma

Aberrant constitutive activation of Rel/NF-κB transcription factors is a hallmark of numerous cancers. Of the five Rel family members, c-Rel has the strongest direct links to tumorigenesis. c-Rel is the only member that can malignantly transform lymphoid cells in vitro. Furthermore, c-Rel is implicated in human B cell lymphoma through the frequent occurrence of REL gene locus gains and amplifications. In normal physiology, high c-Rel expression predominates in the hematopoietic lineage and a diverse range of stimuli can trigger enhanced expression and activation of c-Rel. Both expression and activation of c-Rel are tightly regulated on multiple levels, indicating the necessity to keep its functions under control. In this review we meta-analyze and integrate studies reporting gene locus aberrations to provide an overview on the frequency of REL gains in human B cell lymphoma subtypes, namely follicular lymphoma, diffuse large B cell lymphoma, primary mediastinal B cell lymphoma, and classical Hodgkin lymphoma. We also summarize current knowledge on c-Rel expression and protein localization in these human B cell lymphomas and discuss the co-amplification of BCL11A with REL. In addition, we highlight and illustrate key pathways of c-Rel activation and regulation with a specific focus on B cell biology.

Genome-wide association study implicates immune dysfunction in the development of Hodgkin lymphoma

To further our understanding of inherited susceptibility to Hodgkin lymphoma (HL), we performed a meta-analysis of 7 genome-wide association studies totaling 5325 HL cases and 22 423 control patients. We identify 5 new HL risk loci at 6p21.31 (rs649775; P = 2.11 × 10-10), 6q23.3 (rs1002658; P = 2.97 × 10-8), 11q23.1 (rs7111520; P = 1.44 × 10-11), 16p11.2 (rs6565176; P = 4.00 × 10-8), and 20q13.12 (rs2425752; P = 2.01 × 10-8). Integration of gene expression, histone modification, and in situ promoter capture Hi-C data at the 5 new and 13 known risk loci implicates dysfunction of the germinal center reaction, disrupted T-cell differentiation and function, and constitutive NF-κB activation as mechanisms of predisposition. These data provide further insights into the genetic susceptibility and biology of HL.

Genetic Architecture of Adaptive Immune System Identifies Key Immune Regulators

The immune system is highly diverse, but characterization of its genetic architecture has lagged behind the vast progress made by genome-wide association studies (GWASs) of emergent diseases. Our GWAS for 54 functionally relevant phenotypes of the adaptive immune system in 489 healthy individuals identifies eight genome-wide significant associations explaining 6%-20% of variance. Coding and splicing variants in PTPRC and COMMD10 are involved in memory T cell differentiation. Genetic variation controlling disease-relevant T helper cell subsets includes RICTOR and STON2 associated with Th2 and Th17, respectively, and the interferon-lambda locus controlling regulatory T cell proliferation. Early and memory B cell differentiation stages are associated with variation in LARP1B and SP4. Finally, the latrophilin family member ADGRL2 correlates with baseline pro-inflammatory interleukin-6 levels. Suggestive associations reveal mechanisms of autoimmune disease associations, in particular related to pro-inflammatory cytokine production. Pinpointing these key human immune regulators offers attractive therapeutic perspectives.

Are EBV-related and EBV-unrelated Hodgkin lymphomas different with regard to susceptibility to checkpoint blockade?

Epstein-Barr virus (EBV)–related and EBV-unrelated classical Hodgkin lymphomas (cHLs) are morphologically and phenotypically indistinguishable. However, the tumor microenvironment of EBV-related cHLs contains higher numbers of macrophages and higher expression levels of PD-L1 than that of EBV-unrelated cHLs. Moreover, viral oncoprotein LMP1 may sustain an immunosuppressive microenvironment by inducing/enhancing production of immunosuppressive cytokines and the expression of PD-1. The presence of enhanced immunosuppressive features in EBV-related cHL should make EBV-related cHL patients more susceptible to checkpoint blockade.

Combined linkage and association analysis of classical Hodgkin lymphoma

The heritability of classical Hodgkin lymphoma (cHL) has yet to be fully deciphered. We report a family with five members diagnosed with nodular sclerosis cHL. Genetic analysis of the family provided evidence of linkage at chromosomes 2q35-37, 3p14-22 and 21q22, with logarithm of odds score >2. We excluded the possibility of common genetic variation influencing cHL risk at regions of linkage, by analysing GWAS data from 2,201 cHL cases and 12,460 controls. Whole exome sequencing of affected family members identified the shared missense mutations p.(Arg76Gln) in FAM107A and p.(Thr220Ala) in SLC26A6 at 3p21 as being predicted to impact on protein function. FAM107A expression was shown to be low or absent in lymphoblastoid cell lines and SLC26A6 expression lower in lymphoblastoid cell lines derived from p.(Thr220Ala) mutation carriers. Expression of FAM107A and SLC26A6 was low or absent in Hodgkin Reed-Sternberg (HRS) cell lines and in HRS cells in Hodgkin lymphoma tissue. No sequence variants were detected in KLHDC8B, a gene previously suggested as a cause of familial cHL linked to 3p21. Our findings provide evidence for candidate gene susceptibility to familial cHL.

Genome-wide association study of classical Hodgkin lymphoma identifies key regulators of disease susceptibility

Several susceptibility loci for classical Hodgkin lymphoma have been reported. However, much of the heritable risk is unknown. Here, we perform a meta-analysis of two existing genome-wide association studies, a new genome-wide association study, and replication totalling 5,314 cases and 16,749 controls. We identify risk loci for all classical Hodgkin lymphoma at 6q22.33 (rs9482849, P = 1.52 × 10-8) and for nodular sclerosis Hodgkin lymphoma at 3q28 (rs4459895, P = 9.43 × 10-17), 6q23.3 (rs6928977, P = 4.62 × 10-11), 10p14 (rs3781093, P = 9.49 × 10-13), 13q34 (rs112998813, P = 4.58 × 10-8) and 16p13.13 (rs34972832, P = 2.12 × 10-8). Additionally, independent loci within the HLA region are observed for nodular sclerosis Hodgkin lymphoma (rs9269081, HLA-DPB1*03:01, Val86 in HLA-DRB1) and mixed cellularity Hodgkin lymphoma (rs1633096, rs13196329, Val86 in HLA-DRB1). The new and established risk loci localise to areas of active chromatin and show an over-representation of transcription factor binding for determinants of B-cell development and immune response.

Genome-wide association studies of cancer: current insights and future perspectives

Genome-wide association studies (GWAS) provide an agnostic approach for investigating the genetic basis of complex diseases. In oncology, GWAS of nearly all common malignancies have been performed, and over 450 genetic variants associated with increased risks have been identified. As well as revealing novel pathways important in carcinogenesis, these studies have shown that common genetic variation contributes substantially to the heritable risk of many common cancers. The clinical application of GWAS is starting to provide opportunities for drug discovery and repositioning as well as for cancer prevention. However, deciphering the functional and biological basis of associations is challenging and is in part a barrier to fully unlocking the potential of GWAS.

Genome-Wide Association Studies of Cancer in Diverse Populations

Genome-wide association studies (GWAS) of cancer have identified more than 700 risk loci, of which approximately 80% were first discovered in European ancestry populations, approximately 15% in East Asians, 3% in multiethnic scans, and less than 1% in African and Latin American populations. These percentages closely mirror the distribution of samples included in the discovery phase of cancer GWAS to date (84% European, 11% East Asian, 4% African, and 1% Latin American ancestry). GWAS in non-European ancestry populations have provided insight into ancestry-specific variation in cancer and have pointed to regions of susceptibility that are of particular importance in certain populations. Uncovering and characterizing cancer risk loci in diverse populations is critical for understanding underlying biological mechanisms and developing future genetic risk prediction models in non-European ancestry populations. New GWAS and continued collaborations will be required to eliminate population inequalities in the number of studies, sample sizes, and variant content on GWAS arrays, and to better align genetic research in cancer to the global distribution of race/ethnicity Cancer Epidemiol Biomarkers Prev; 27(4); 405-17. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Germline variation in the 3'-untranslated region of the POU2AF1 gene is associated with susceptibility to lymphoma

Genetic variations in certain genes may alter the susceptibility to lymphoma. We searched electronic databases and selected candidate single-nucleotide polymorphisms (SNPs) located within 3'-untranslated regions (3'-UTRs) that might affect miRNA-binding ability in the 50 most dysregulated genes in lymphoma for further study. We found that rs1042752-located in the 3'-UTR of POU2AF1, which plays a vital role in lymphomagenesis-was significantly associated with lymphoma risk in a case-control study with 793 patients and 939 controls. Compared with individuals with the rs1042752TT genotype, those with the rs1042752CC genotype had a higher risk of lymphoma (OR = 2.14, 95% CI: 1.55-2.95, P < 0.001), even in stratified analysis for non-Hodgkin lymphoma (OR = 4.58, 95% CI: 2.38-8.81, P < 0.001), B-cell lymphoma (OR = 4.89, 95% CI: 2.46-9.73, P < 0.001), T-cell lymphoma (OR = 4.20, 95% CI: 1.76-10.05, P = 0.001), and Hodgkin lymphoma (OR = 3.62, 95% CI: 1.25-10.46, P = 0.018). Similar results were also observed in a recessive genetic model. Database findings suggested that rs1042752 might affect the interaction of POU2AF1 mRNA with hsa-miR-633. Functional assays confirmed that rs1042752C altered the binding site of hsa-miR-633 and increased POU2AF1 expression in Ramos, HuT 102, and Jurkat E6-1 cell lines. These findings demonstrate for the first time that functional polymorphism in the 3'-UTR of POU2AF1 is associated with susceptibility, and that SNP interaction with hsa-miR-633 affects gene expression and increases the risk of lymphoma.

Epigenetic and genetic alterations and their influence on gene regulation in chronic lymphocytic leukemia

BACKGROUND

To understand the changes of gene regulation in carcinogenesis, we explored signals of DNA methylation - a stable epigenetic mark of gene regulatory elements - and designed a computational model to profile loss and gain of regulatory elements (REs) during carcinogenesis. We also utilized sequencing data to analyze the allele frequency of single nucleotide polymorphisms (SNPs) and detected the cancer-associated SNPs, i.e., the SNPs displaying the significant allele frequency difference between cancer and normal samples.

RESULTS

After applying this model to chronic lymphocytic leukemia (CLL) data, we identified REs differentially activated (dREs) between normal and CLL cells, consisting of 6,802 dREs gained and 4,606 dREs lost in CLL. The identified regulatory perturbations coincide with changes in the expression of target genes. In particular, the genes encoding DNA methyltransferases harbor multiple lost-in-cancer dREs and zero gained-in-cancer dREs, indicating that the damaged regulation of these genes might be one of the key causes of tumor formation. dREs display a significantly elevated density of the genome-wide association study (GWAS) SNPs associated with CLL and CLL-related traits. We observed that most of dRE GWAS SNPs associated with CLL and CLL-related traits (83%) display a significant haplotype association among the identified cancer-associated alleles and the risk alleles that have been reported in GWAS. Also dREs are enriched for the binding sites of the well-established B-cell and CLL transcription factors (TFs) NF-kB, AP2, P53, E2F1, PAX5, and SP1. We also identified CLL-associated SNPs and demonstrated that the mutations at these SNPs change the binding sites of key TFs much more frequently than expected.

CONCLUSIONS

Through exploring sequencing data measuring DNA methylation, we identified the epigenetic alterations (more specifically, DNA methylation) and genetic mutations along non-coding genomic regions CLL, and demonstrated that these changes play a critical role in carcinogenesis through damaging the regulation of key genes and alternating the binding of key TFs in B and CLL cells.

Candidate gene association studies and risk of Hodgkin lymphoma: a systematic review and meta-analysis

To evaluate the contribution of association studies of candidate polymorphisms to inherited predisposition to Hodgkin lymphoma (HL), we conducted a systematic review and meta-analysis of published case-control studies. Of the variants examined more than once in candidate gene association studies, we identified 21 studies that reported on 12 polymorphic variants in 10 genes. Data were also extracted from a published genome wide association study to allow analysis of an additional 47 variants in a further 30 genes. Promising associations were seen in nine of the variants (p < 0.05). Given that the estimated false positive report probabilities (FPRPs) for all associations are high (i.e. FPRP > 0.2), these findings should be interpreted with caution. While studies of candidate polymorphisms may be an attractive means of identifying risk factors for HL, future studies should employ sample sizes adequately powered to identify variants having only modest effects on HL risk. Furthermore, because of aetiological heterogeneity within HL, stratification of genotyping according to age, tumour Epstein-Barr virus status and histology is essential. Copyright © 2015 John Wiley & Sons, Ltd.

Genome-wide association analysis of chronic lymphocytic leukaemia, Hodgkin lymphoma and multiple myeloma identifies pleiotropic risk loci

B-cell malignancies (BCM) originate from the same cell of origin, but at different maturation stages and have distinct clinical phenotypes. Although genetic risk variants for individual BCMs have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. We explored genome-wide association studies of chronic lymphocytic leukaemia (CLL, N = 1,842), Hodgkin lymphoma (HL, N = 1,465) and multiple myeloma (MM, N = 3,790). We identified a novel pleiotropic risk locus at 3q22.2 (NCK1, rs11715604, P = 1.60 × 10-9) with opposing effects between CLL (P = 1.97 × 10-8) and HL (P = 3.31 × 10-3). Eight established non-HLA risk loci showed pleiotropic associations. Within the HLA region, Ser37 + Phe37 in HLA-DRB1 (P = 1.84 × 10-12) was associated with increased CLL and HL risk (P = 4.68 × 10-12), and reduced MM risk (P = 1.12 × 10-2), and Gly70 in HLA-DQB1 (P = 3.15 × 10-10) showed opposing effects between CLL (P = 3.52 × 10-3) and HL (P = 3.41 × 10-9). By integrating eQTL, Hi-C and ChIP-seq data, we show that the pleiotropic risk loci are enriched for B-cell regulatory elements, as well as an over-representation of binding of key B-cell transcription factors. These data identify shared biological pathways influencing the development of CLL, HL and MM. The identification of these risk loci furthers our understanding of the aetiological basis of BCMs.

Immunochip analysis identifies novel susceptibility loci in the human leukocyte antigen region for acquired thrombotic thrombocytopenic purpura

Essentials Genetic predisposition to acquired thrombotic thrombocytopenic purpura (aTTP) is mainly unknown. Genetic risk factors for aTTP were studied by Immunochip analysis and replication study. Human leukocyte antigen (HLA) variant rs6903608 conferred a 2.5-fold higher risk of developing aTTP. rs6903608 and HLA-DQB1*05:03 may explain most of the HLA association signal in aTTP. Click to hear Dr Cataland's presentation on acquired thrombotic thrombocytopenic purpura SUMMARY: Background Acquired thrombotic thrombocytopenic purpura (TTP) is a rare, life-threatening thrombotic microangiopathy associated with the development of autoantibodies against the von Willebrand factor-cleaving protease ADAMTS-13. Similarly to what has been found for other autoimmune disorders, there is evidence of a genetic contribution, including the association of the human leukocyte antigen (HLA) class II complex with disease risk. Objective To identify novel genetic risk factors in acquired TTP. Patients/Methods We undertook a case-control genetic association study in 190 European-origin TTP patients and 1255 Italian healthy controls by using the Illumina Immunochip. Replication analysis in 88 Italian cases and 456 controls was performed with single-nucleotide polymorphism (SNP) TaqMan assays. Results and conclusion We identified one common variant (rs6903608) located within the HLA class II locus that was independently associated with acquired TTP at genome-wide significance and conferred a 2.6-fold increased risk of developing a TTP episode (95% confidence interval [CI] 2.02-3.27, P = 1.64 × 10-14 ). We also found five non-HLA variants mapping to chromosomes 2, 6, 8 and X that were suggestively associated with the disease: rs9490550, rs115265285, rs5927472, rs7823314, and rs1334768 (nominal P-values ranging from 1.59 × 10-5 to 7.60 × 10-5 ). Replication analysis confirmed the association of HLA variant rs6903608 with acquired TTP (pooled P = 3.95 × 10-19 ). Imputation of classic HLA genes followed by stepwise conditional analysis revealed that the combination of rs6903608 and HLA-DQB1*05:03 may explain most of the HLA association signal in acquired TTP. Our results refined the association of the HLA class II locus with acquired TTP, confirming its importance in the etiology of this autoimmune disease.

Whole exome sequencing in families at high risk for Hodgkin lymphoma: identification of a predisposing mutation in the KDR gene

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

Evidence of Inbreeding in Hodgkin Lymphoma

Genome-wide association studies (GWASs) have identified several, mainly co-dominantly acting, single-nucleotide polymorphisms (SNPs) associated with Hodgkin lymphoma (HL). We searched for recessively acting disease loci by performing an analysis of runs of homozygosity (ROH) based on windows of homozygous SNP-blocks and by calculating genomic inbreeding coefficients on a SNP-wise basis. We used data from a previous GWAS with 906 cases and 1217 controls from a population with a long history of no matings between relatives. Ten recurrent ROHs were identified among 25 055 ROHs across all individuals but their association with HL was not genome-wide significant. All recurrent ROHs showed significant evidence for natural selection. As a novel finding genomic inbreeding among cases was significantly higher than among controls (P = 2.11*10-14) even after correcting for covariates. Higher inbreeding among the cases was mainly based on a group of individuals with a higher average length of ROHs per person. This result suggests a correlation of higher levels of inbreeding with higher cancer incidence and might reflect the existence of recessive alleles causing HL. Genomic inbreeding may result in a higher expression of deleterious recessive genes within a population.

Human leukocyte antigens and genetic susceptibility to lymphoma

Familial aggregation, coupled with ethnic variation in incidence, suggests that inherited susceptibility plays a role in the development of lymphoma, and the search for genetic risk factors has highlighted the contribution of the human leukocyte antigen (HLA) complex. In a landmark study published almost 50 years ago, Hodgkin lymphoma (HL) was the first disease to be associated with HLA variation. It is now clear that Epstein-Barr virus (EBV)-positive and -negative HL are strongly associated with specific HLA polymorphisms but these differ by EBV status of the tumours. HLA class I alleles are consistently associated with EBV-positive HL while a polymorphism in HLA class II is the strongest predictor of risk of EBV-negative HL. Recent investigations, particularly genome-wide association studies (GWAS), have also revealed associations between HLA and common types of non-Hodgkin lymphoma (NHL). Follicular lymphoma is strongly associated with two distinct haplotypes in HLA class II whereas diffuse large B-cell lymphoma is most strongly associated with HLA-B*08. Although chronic lymphocytic leukaemia is associated with variation in HLA class II, the strongest signals in GWAS are from non-HLA polymorphisms, suggesting that inherited susceptibility is explained by co-inheritance of multiple low risk variants. Associations between B-cell derived lymphoma and HLA variation suggest that antigen presentation, or lack of, plays an important role in disease pathogenesis but the precise mechanisms have yet to be elucidated.

Meta-analysis of genome-wide association studies discovers multiple loci for chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is a common lymphoid malignancy with strong heritability. To further understand the genetic susceptibility for CLL and identify common loci associated with risk, we conducted a meta-analysis of four genome-wide association studies (GWAS) composed of 3,100 cases and 7,667 controls with follow-up replication in 1,958 cases and 5,530 controls. Here we report three new loci at 3p24.1 (rs9880772, EOMES, P=2.55 × 10(-11)), 6p25.2 (rs73718779, SERPINB6, P=1.97 × 10(-8)) and 3q28 (rs9815073, LPP, P=3.62 × 10(-8)), as well as a new independent SNP at the known 2q13 locus (rs9308731, BCL2L11, P=1.00 × 10(-11)) in the combined analysis. We find suggestive evidence (P<5 × 10(-7)) for two additional new loci at 4q24 (rs10028805, BANK1, P=7.19 × 10(-8)) and 3p22.2 (rs1274963, CSRNP1, P=2.12 × 10(-7)). Pathway analyses of new and known CLL loci consistently show a strong role for apoptosis, providing further evidence for the importance of this biological pathway in CLL susceptibility.

A Novel Risk Locus at 6p21.3 for Epstein-Barr Virus-Positive Hodgkin Lymphoma

BACKGROUND

A proportion of the genetic variants involved in susceptibility to Hodgkin lymphoma differ by the tumor's Epstein-Barr virus (EBV) status, particularly within the MHC region.

METHODS

We have conducted an SNP imputation study of the MHC region, considering tumor EBV status in 1,200 classical Hodgkin lymphoma (cHL) cases and 5,726 control subjects of European origin. Notable findings were genotyped in an independent study population of 468 cHL cases and 551 controls.

RESULTS

We identified and subsequently replicated a novel association between a common genetic variant rs6457715 and cHL. Although strongly associated with EBV-positive cHL [OR, 2.33; 95% confidence interval (CI), 1.83-2.97; P = 7 × 10(-12)], there was little evidence for association between rs6457715 and the EBV-negative subgroup of cHL (OR, 1.06; 95% CI, 0.92-1.21), indicating that this association was specific to the EBV-positive subgroup (Phet < P = 10(-8)). Furthermore, the association was limited to EBV-positive cHL subgroups within mixed cell (MCHL) and nodular sclerosis subtypes (NSHL), suggesting that the association is independent of histologic subtype of cHL.

CONCLUSIONS

rs6457715, located near the HLA-DPB1 gene, is associated with EBV-positive cHL and suggests this region as a novel susceptibility locus for cHL.

IMPACT

This expands the number of genetic variants that are associated with cHL and provides additional evidence for a critical and specific role of EBV in the etiology of this disease.

Response

Genome-wide association studies (GWASs) have identified several single-nucleotide polymorphisms (SNPs) influencing the risk of Hodgkin's lymphoma (HL) and demonstrated the association of common genetic variation for this type of cancer. Such evidence for inherited genetic risk is also provided by the family history and the very high concordance between monozygotic twins. However, little is known about the genetic and environmental contributions. A common measure for describing the phenotypic variation due to genetics is the heritability. Using GWAS data on 906 HL cases by considering all typed SNPs simultaneously, we have calculated that the common variance explained by SNPs accounts for >35% of the total variation on the liability scale in HL (95% confidence interval 6–62%). These findings are consistent with similar heritability estimates of ∼0.40 (95% confidence interval 0.17–0.58) based on Swedish population data. Our estimates support the underlying polygenic basis for susceptibility to HL, and show that heritability based on the population data is somehow larger than heritability based on the genomic data because of the possibility of some missing heritability in the GWAS data. Besides that there is still major evidence for multiple loci causing HL on chromosomes other than chromosome 6 that need to be detected. Because of limited findings in prior GWASs, it seems worth checking for more loci causing susceptibility to HL.

Familial predisposition and genetic risk factors for lymphoma

Our understanding of familial predisposition to lymphoma (collectively defined as non-Hodgkin lymphoma [NHL], Hodgkin lymphoma [HL], and chronic lymphocytic leukemia [CLL]) outside of rare hereditary syndromes has progressed rapidly during the last decade. First-degree relatives of NHL, HL, and CLL patients have an ∼1.7-fold, 3.1-fold, and 8.5-fold elevated risk of developing NHL, HL, and CLL, respectively. These familial risks are elevated for multiple lymphoma subtypes and do not appear to be confounded by nongenetic risk factors, suggesting at least some shared genetic etiology across the lymphoma subtypes. However, a family history of a specific subtype is most strongly associated with risk for that subtype, supporting subtype-specific genetic factors. Although candidate gene studies have had limited success in identifying susceptibility loci, genome-wide association studies (GWAS) have successfully identified 67 single nucleotide polymorphisms from 41 loci, predominately associated with specific subtypes. In general, these GWAS-discovered loci are common (minor allele frequency >5%), have small effect sizes (odds ratios, 0.60-2.0), and are of largely unknown function. The relatively low incidence of lymphoma, modest familial risk, and the lack of a screening test and associated intervention, all argue against active clinical surveillance for lymphoma in affected families at this time.

Genome-wide homozygosity signature and risk of Hodgkin lymphoma

Recent studies have reported that regions of homozygosity (ROH) in the genome are detectable in outbred populations and can be associated with an increased risk of malignancy. To examine whether homozygosity is associated with an increased risk of developing Hodgkin lymphoma (HL) we analysed 589 HL cases and 5,199 controls genotyped for 484,072 tag single nucleotide polymorphisms (SNPs). Across the genome the cumulative distribution of ROH was not significantly different between cases and controls. Seven ROH at 4q22.3, 4q32.2, 7p12.3-14.1, 7p22.2, 10p11.22-23, 19q13.12-2 and 19p13.2 were associated with HL risk at P < 0.01. Intriguingly 4q22.3 harbours an ROH to which the nuclear factor NF-kappa-B p105 subunit (NFKB1) maps (P = 0.002). The ROH at 19q13.12-2 has previously been implicated in B-cell precursor acute lymphoblastic leukaemia. Aside from these observations which require validation, it is unlikely that levels of measured homozygosity caused by autozygosity, uniparental isodisomy or hemizygosity play a major role in defining HL risk in predominantly outbred populations.

Risk of familial classical Hodgkin lymphoma by relationship, histology, age, and sex: a joint study from five Nordic countries

We aimed to provide the familial risk of classical Hodgkin lymphoma (HL) by relationship, histology, age at diagnosis, and sex. A cohort of 57,475 first-degree relatives of 13,922 HL patients diagnosed between 1955 and 2009 in 5 European countries was observed for HL incidence. The overall lifetime cumulative risk (CR) of HL in first-degree relatives of a patient with HL was 0.6%, which represents a threefold (standardized incidence ratio [SIR], 3.3; 95% confidence interval [CI], 2.8-3.9) increased risk over the general population risk. The risk in siblings (6.0-fold; 95% CI, 4.8- to 7.4-fold) was significantly higher than in parents and/or children (2.1-fold; 95% CI, 1.6- to 2.6-fold). Very high lifetime risk of HL was found for those with multiple affected first-degree relatives (13-fold; 95% CI, 2.8- to 39-fold) and for same-sex twins (57-fold; 95% CI, 21- to 125-fold). We found high familial risks between some concordant histologic subtypes of HL such as lymphocyte-rich (81-fold; 95% CI, 30- to 177-fold) and nodular sclerosis (4.6-fold; 95% CI, 2.9- to 7.0-fold) and also between some discordant subtypes. The familial risk in sisters (9.4-fold; 95% CI, 5.9- to 14-fold) was higher than in brothers (4.5-fold; 95% CI, 2.9- to 6.7-fold) or unlike-sex siblings (5.9-fold; 95% CI, 4.3- to 8.1-fold). The lifetime risk of HL was higher when first-degree relatives were diagnosed at early ages (before age 30 years). This study provides tangible absolute risk estimates for relatives of HL patients, which can be used as a sex-, age-, and family history-based risk calculator for classical HL by oncologists and genetic counselors.

Genetic variation at MECOM, TERT, JAK2 and HBS1L-MYB predisposes to myeloproliferative neoplasms

Clonal proliferation in myeloproliferative neoplasms (MPN) is driven by somatic mutations in JAK2, CALR or MPL, but the contribution of inherited factors is poorly characterized. Using a three-stage genome-wide association study of 3,437 MPN cases and 10,083 controls, we identify two SNPs with genome-wide significance in JAK2(V617F)-negative MPN: rs12339666 (JAK2; meta-analysis P=1.27 × 10(-10)) and rs2201862 (MECOM; meta-analysis P=1.96 × 10(-9)). Two additional SNPs, rs2736100 (TERT) and rs9376092 (HBS1L/MYB), achieve genome-wide significance when including JAK2(V617F)-positive cases. rs9376092 has a stronger effect in JAK2(V617F)-negative cases with CALR and/or MPL mutations (Breslow-Day P=4.5 × 10(-7)), whereas in JAK2(V617F)-positive cases rs9376092 associates with essential thrombocythemia (ET) rather than polycythemia vera (allelic χ(2) P=7.3 × 10(-7)). Reduced MYB expression, previously linked to development of an ET-like disease in model systems, associates with rs9376092 in normal myeloid cells. These findings demonstrate that multiple germline variants predispose to MPN and link constitutional differences in MYB expression to disease phenotype.

Modeling HLA associations with EBV-positive and -negative Hodgkin lymphoma suggests distinct mechanisms in disease pathogenesis

HLA genotyping and genome wide association studies provide strong evidence for associations between Human Leukocyte Antigen (HLA) alleles and classical Hodgkin lymphoma (cHL). Analysis of these associations is complicated by the extensive linkage disequilibrium within the major histocompatibility region and recent data suggesting that associations with EBV‐positive and EBV‐negative cHL are largely distinct. To distinguish independent and therefore potentially causal associations from associations confounded by linkage disequilibrium, we applied a variable selection regression modeling procedure to directly typed HLA class I and II genes and selected SNPs from EBV‐stratified patient subgroups. In final models, HLA‐A*01:01 and B*37:01 were associated with an increased risk of EBV‐positive cHL whereas DRB1*15:01 and DPB1*01:01 were associated with decreased risk. Effects were independent of a prior history of infectious mononucleosis. For EBV‐negative cHL the class II SNP rs6903608 remained the strongest predictor of disease risk after adjusting for the effects of common HLA alleles. Associations with “all cHL” and differences by case EBV status reflected the subgroup analysis. In conclusion, this study extends previous findings by identifying novel HLA associations with EBV‐stratified subgroups of cHL, highlighting those alleles likely to be biologically relevant and strengthening evidence implicating genetic variation associated with the SNP rs6903608.

What's new?

Strong evidence exists for associations between HLA alleles and classical Hodgkin lymphoma (cHL). Analysis is however complicated by the linkage disequilibrium within the MHC region and data suggesting that associations with Epstein‐Barr virus (EBV)‐positive and negative cHL are distinct. In the largest study to date to investigate associations between EBV‐stratified cHL subgroups and directly typed HLA alleles, the authors extend associations with EBV‐positive cHL to novel HLA class II alleles, which are associated with decreased disease risk. For EBV‐negative disease, the class II SNP rs6903608 remains the strongest predictor of risk after adjusting for the effects of common HLA alleles.

Genetic associations in classical hodgkin lymphoma: a systematic review and insights into susceptibility mechanisms

Both targeted and genome-wide studies have revealed genetic associations for susceptibility, prognosis, and treatment-induced secondary malignancies and toxicities in classical Hodgkin lymphoma (cHL). This review gives a systematic and comprehensive overview of significant associations and places them into a biologic context. The strongest susceptibility polymorphisms have been found for the human leukocyte antigen (HLA) genes. These associations are specific for cHL overall or for subgroups based on tumor cell Epstein-Barr virus (EBV) status. These findings strongly suggest that EBV-specific immune responses influence cHL susceptibility in EBV(+) cHL and that immune responses targeting other tumor-associated antigens are important in EBV(-) cHL. Accordingly, most of the numerous other susceptibility loci map to genes that affect functionality of the immune system, underscoring the crucial role of the immune system in cHL development. The number of association studies on cHL prognosis is limited with one consistent association for the drug-metabolizing UGT1A1 gene. PRDM1 is associated with radiation-induced secondary malignancies and a small number of genes are associated with treatment-related toxicities. In conclusion, most loci showing genetic associations in cHL harbor genes with a potential functional relevance for cHL susceptibility.

A meta-analysis of Hodgkin lymphoma reveals 19p13.3 TCF3 as a novel susceptibility locus

Recent genome-wide association studies (GWAS) of Hodgkin lymphoma (HL) have identified associations with genetic variation at both HLA and non-HLA loci; however, much of heritable HL susceptibility remains unexplained. Here we perform a meta-analysis of three HL GWAS totaling 1,816 cases and 7,877 controls followed by replication in an independent set of 1,281 cases and 3,218 controls to find novel risk loci. We identify a novel variant at 19p13.3 associated with HL (rs1860661; odds ratio (OR)=0.81, 95% confidence interval (95% CI) = 0.76-0.86, P(combined) = 3.5 × 10(-10)), located in intron 2 of TCF3 (also known as E2A), a regulator of B- and T-cell lineage commitment known to be involved in HL pathogenesis. This meta-analysis also notes associations between previously published loci at 2p16, 5q31, 6p31, 8q24 and 10p14 and HL subtypes. We conclude that our data suggest a link between the 19p13.3 locus, including TCF3, and HL risk.