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

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.

Design and quality control of large-scale two-sample Mendelian randomization studies

Background

Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors.

Methods

We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set.

Results

We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls.

Conclusions

In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

Complementary Alu sequences mediate enhancer-promoter selectivity

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters^1-4. However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs^5,6. Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions.

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.

Cross-cancer pleiotropic analysis identifies three novel genetic risk loci for colorectal cancer

BACKGROUND

To understand the shared genetic basis between colorectal cancer (CRC) and other cancers and identify potential pleiotropic loci for compensating the missing genetic heritability of CRC.

METHODS

We conducted a systematic genome-wide pleiotropy scan to appraise associations between cancer-related genetic variants and CRC risk among European populations. Single nucleotide polymorphism (SNP)-set analysis was performed using data from the UK Biobank and the Study of Colorectal Cancer in Scotland (10 039 CRC cases and 30 277 controls) to evaluate the overlapped genetic regions for susceptibility of CRC and other cancers. The variant-level pleiotropic associations between CRC and other cancers were examined by CRC genome-wide association study meta-analysis and the pleiotropic analysis under composite null hypothesis (PLACO) pleiotropy test. Gene-based, co-expression and pathway enrichment analyses were performed to explore potential shared biological pathways. The interaction between novel genetic variants and common environmental factors was further examined for their effects on CRC.

RESULTS

Genome-wide pleiotropic analysis identified three novel SNPs (rs2230469, rs9277378 and rs143190905) and three mapped genes (PIP4K2A, HLA-DPB1 and RTEL1) to be associated with CRC. These genetic variants were significant expressions quantitative trait loci in colon tissue, influencing the expression of their mapped genes. Significant interactions of PIP4K2A and HLA-DPB1 with environmental factors, including smoking and alcohol drinking, were observed. All mapped genes and their co-expressed genes were significantly enriched in pathways involved in carcinogenesis.

CONCLUSION

Our findings provide an important insight into the shared genetic basis between CRC and other cancers. We revealed several novel CRC susceptibility loci to help understand the genetic architecture of CRC.

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.

New molecular targets in Hodgkin and Reed-Sternberg cells

Recent discoveries shed light on molecular mechanisms responsible for classical Hodgkin lymphoma (HL) development and progression, along with features of Hodgkin - Reed and Sternberg cells (HRS). Here, we summarize current knowledge on characteristic molecular alterations in HL, as well as existing targeted therapies and potential novel treatments for this disease. We discuss the importance of cluster of differentiation molecule 30 (CD30) and the programmed cell death-1 protein (PD-1) and ligands (PD-L1/2), and other molecules involved in immune modulation in HL. We highlight emerging evidence indicating that the altered function of SWI/SNF-type chromatin remodeling complexes, PRC2, and other epigenetic modifiers, contribute to variations in chromatin status, which are typical for HL. We postulate that despite of the existence of plentiful molecular data, the understanding of HL development remains incomplete. We therefore propose research directions involving analysis of reverse signaling in the PD-1/PD-L1 mechanism, chromatin remodeling, and epigenetics-related alterations, in order to identify HL features at the molecular level. Such attempts may lead to the identification of new molecular targets, and thus will likely substantially contribute to the future development of more effective targeted therapies.

Discovery of novel predisposing coding and noncoding variants in familial Hodgkin lymphoma

Familial aggregation of Hodgkin lymphoma (HL) has been demonstrated in large population studies, pointing to genetic predisposition to this hematological malignancy. To understand the genetic variants associated with the development of HL, we performed whole genome sequencing on 234 individuals with and without HL from 36 pedigrees that had 2 or more first-degree relatives with HL. Our pedigree selection criteria also required at least 1 affected individual aged <21 years, with the median age at diagnosis of 21.98 years (3-55 years). Family-based segregation analysis was performed for the identification of coding and noncoding variants using linkage and filtering approaches. Using our tiered variant prioritization algorithm, we identified 44 HL-risk variants in 28 pedigrees, of which 33 are coding and 11 are noncoding. The top 4 recurrent risk variants are a coding variant in KDR (rs56302315), a 5' untranslated region variant in KLHDC8B (rs387906223), a noncoding variant in an intron of PAX5 (rs147081110), and another noncoding variant in an intron of GATA3 (rs3824666). A newly identified splice variant in KDR (c.3849-2A>C) was observed for 1 pedigree, and high-confidence stop-gain variants affecting IRF7 (p.W238∗) and EEF2KMT (p.K116∗) were also observed. Multiple truncating variants in POLR1E were found in 3 independent pedigrees as well. Whereas KDR and KLHDC8B have previously been reported, PAX5, GATA3, IRF7, EEF2KMT, and POLR1E represent novel observations. Although there may be environmental factors influencing lymphomagenesis, we observed segregation of candidate germline variants likely to predispose HL in most of the pedigrees studied.

Learning from the nexus of autoimmunity and cancer

The immune system plays critical roles in both autoimmunity and cancer, diseases at opposite ends of the immune spectrum. Autoimmunity arises from loss of T cell tolerance against self, while in cancer, poor immunity against transformed self fails to control tumor growth. Blockade of pathways that preserve self-tolerance is being leveraged to unleash immunity against many tumors; however, widespread success is hindered by the autoimmune-like toxicities that arise in treated patients. Knowledge gained from the treatment of autoimmunity can be leveraged to treat these toxicities in patients. Further, the understanding of how T cell dysfunction arises in cancer can be leveraged to induce a similar state in autoreactive T cells. Here, we review what is known about the T cell response in autoimmunity and cancer and highlight ways in which we can learn from the nexus of these two diseases to improve the application, efficacy, and management of immunotherapies.

A Genome-Wide Association Study Identified Novel Genetic Susceptibility Loci for Oral Cancer in Taiwan

Taiwan has the highest incidence rate of oral cancer in the world. Although oral cancer is mostly an environmentally induced cancer, genetic factors also play an important role in its etiology. Genome-wide association studies (GWAS) have identified nine susceptibility regions for oral cancers in populations of European descent. In this study, we performed the first GWAS of oral cancer in Taiwan with 1529 cases and 44,572 controls. We confirmed two previously reported loci on the 6p21.33 (HLA-B) and 6p21.32 (HLA-DQ gene cluster) loci, highlighting the importance of the human leukocyte antigen and, hence, the immunologic mechanisms in oral carcinogenesis. The TERT-CLMPT1L locus on 5p15.33, the 4q23 ADH1B locus, and the LAMC3 locus on 9q34.12 were also consistent in the Taiwanese. We found two new independent loci on 6p21.32, rs401775 in SKIV2L gene and rs9267798 in TNXB gene. We also found two suggestive novel Taiwanese-specific loci near the TPRS1 gene on 8q23.3 and in the TMED3 gene on 15q25.1. This study identified both common and unique oral cancer susceptibility loci in the Taiwanese as compared to populations of European descent and shed significant light on the etiology of oral cancer in Taiwan.

NFkB Pathway and Hodgkin Lymphoma

The tumor cells that drive classical Hodgkin lymphoma (cHL), namely, Hodgkin and Reed-Sternberg (HRS) cells, display hallmark features that include their rareness in contrast with an extensive and rich reactive microenvironment, their loss of B-cell phenotype markers, their immune escape capacity, and the activation of several key biological pathways, including the constitutive activation of the NFkB pathway. Both canonical and alternative pathways are deregulated by genetic alterations of their components or regulators, EBV infection and interaction with the microenvironment through multiple receptors, including CD30, CD40, BAFF, RANK and BCMA. Therefore, NFkB target genes are involved in apoptosis, cell proliferation, JAK/STAT pathway activation, B-cell marker expression loss, cellular interaction and a positive NFkB feedback loop. Targeting this complex pathway directly (NIK inhibitors) or indirectly (PIM, BTK or NOTCH) remains a challenge with potential therapeutic relevance. Nodular predominant HL (NLPHL), a distinct and rare HL subtype, shows a strong NFkB activity signature because of mechanisms that differ from those observed in cHL, which is discussed in this review.

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.

Genetic association and single-cell transcriptome analyses reveal distinct features connecting autoimmunity with cancers

Autoimmune diseases (ADs) are at a significantly higher risk of cancers with unclear mechanism. By searching GWAS catalog database and Medline, susceptible genes for five common ADs, including systemic lupus erythematosus (SLE), rheumatoid arthritis, Sjögren syndrome, systemic sclerosis, and idiopathic inflammatory myopathies, were collected and then were overlapped with cancer driver genes. Single-cell transcriptome analysis was performed in the comparation between SLE and related cancer. We identified 45 carcinogenic autoimmune disease risk (CAD) genes, which were mainly enriched in T cell signaling pathway and B cell signaling pathway. Integrated single-cell analysis revealed immune cell signaling was significantly downregulated in renal cancer compared with SLE, while stemness signature was significantly enriched in both renal cancer or lymphoma and SLE in specific subpopulations. Drugs targeting CAD genes were shared between ADs and cancer. Our study highlights the common and specific features between ADs and related cancers, and sheds light on a new discovery of treatments.

Genome-wide association study identifies susceptibility loci for acute myeloid leukemia

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

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.

Ultrasound-Guided Anterior Mediastinotomy: A Feasible Tool for Critical Lymphoma Patients

A significant part of all neoplasms growing in anterior mediastinum are lymphomas (25%). Achieving a correct diagnosis and a clear definition of a lymphoma’s subtype is crucial for beginning chemotherapy as soon as possible. However, most patients present a large mediastinal mass that compresses vessels and airway, with serious cardiorespiratory repercussions. Therefore, having multiple tools available to biopsy the lesion without worsening morbidity becomes fundamental. Patients enrolled in this study were unfit for a surgical biopsy in general anesthesia and the need to begin chemotherapy as fast as possible prompted us to avoid percutaneous fine needle aspiration to prevent diagnostic failures. Our observational study included 13 consecutive patients with radiological findings of anterior mediastinal mass. Ultrasonography was performed directly in the theatre to mark the lesion and to localize vessels and vascularized neoplastic tissue. Open biopsy was carried out in spontaneous breathing with a laryngeal mask and with short-acting medications for a rapid anesthesia, performing an anterior mediastinotomy. The mean operative time was 33.4 ± 6.2 min and spontaneous respiration was maintained throughout the procedure. No complications were reported. All patients were discharged in the first or second postoperative day after a chest X-ray (1.38 ± 0.5 days). The diagnostic yield of this approach was 100%. With the addition of ultrasonography right before the procedure and with spontaneous breathing, anterior mediastinotomy still represents a useful tool in critical patients that could hardly tolerate a general anesthesia. The diagnostic yield is high, and the low postoperative morbidity allows a rapid onset of chemotherapy.

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.

A genome-wide association study identifies susceptibility loci for primary central nervous system lymphoma at 6p25.3 and 3p22.1: a LOC Network study

BACKGROUND

Primary central nervous system lymphoma (PCNSL) is a rare form of extra-nodal non-Hodgkin lymphoma. PCNSL is a distinct subtype of non-Hodgkin lymphoma, with over 95% of tumors belonging to the diffuse large B-cell lymphoma (DLBCL) group. We have conducted a genome-wide association study (GWAS) on immunocompetent patients to address the possibility that common genetic variants influence the risk of developing PCNSL.

METHODS

We performed a meta-analysis of 2 new GWASs of PCNSL totaling 475 cases and 1134 controls of European ancestry. To increase genomic resolution, we imputed >10 million single nucleotide polymorphisms using the 1000 Genomes Project combined with UK10K as reference. In addition we performed a transcription factor binding disruption analysis and investigated the patterns of local chromatin by Capture Hi-C data.

RESULTS

We identified independent risk loci at 3p22.1 (rs41289586, ANO10, P = 2.17 × 10-8) and 6p25.3 near EXOC2 (rs116446171, P = 1.95 x 10-13). In contrast, the lack of an association between rs41289586 and DLBCL suggests distinct germline predisposition to PCNSL and DLBCL. We found looping chromatin interactions between noncoding regions at 6p25.3 (rs11646171) with the IRF4 promoter and at 8q24.21 (rs13254990) with the MYC promoter, both genes with strong relevance to B-cell tumorigenesis.

CONCLUSION

To our knowledge this is the first study providing insight into the genetic predisposition to PCNSL. Our findings represent an important step in defining the contribution of common genetic variation to the risk of developing PCNSL.

No Evidence for a Causal Relationship Between Cancers and Parkinson's Disease

BACKGROUND

Epidemiological data suggest that cancer patients have a reduced risk of subsequent Parkinson's disease (PD) development, but the prevalence of PD in melanoma patients is often reported to be increased. Causal relationships between cancers and PD have not been fully explored.

OBJECTIVE

To study causal relationship between different cancers and PD.

METHODS

We used GWAS summary statistics of 15 different types of cancers and two-sample Mendelian randomization to study the causal relationship with PD.

RESULTS

There was no evidence to support a causal relationship between the studied cancers and PD. We also performed reverse analyses between PD and cancers with available full summary statistics (melanoma, breast, prostate, endometrial and keratinocyte cancers) and did not find evidence of causal relationship.

CONCLUSION

We found no evidence to support a causal relationship between cancers and PD and the previously reported associations could be a result of genetic pleiotropy, shared biology or biases.

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.

The genetic polymorphism down-regulating HLA-DRB1 enhancer activity facilitates HBV persistence, evolution and hepatocarcinogenesis in the Chinese Han population

Genetic predisposition of human leucocyte antigen (HLA)-DR has been linked to nonresponse to hepatitis B virus (HBV) vaccination. We sought to reveal their effects on chronic infection and evolution of HBV and development of hepatocellular carcinoma (HCC). Genetic polymorphisms at HLA-DR enhancer regions were genotyped in 4588 participants using quantitative PCR. HBV mutations were determined by sequencing. A dual-luciferase assay was applied to detect the enhancer activity. Associations between HLA-DR polymorphisms and postoperative prognosis were investigated in another cohort of 397 HBV-infected HCC patients. Variant alleles (rs3135395-T, rs3135338-C and rs477515-T) were significantly associated with a decreased risk of HBV persistence in Chinese patients. rs3135395-T, rs3135338-C, rs477515-T and rs2395178-G also significantly decreased HCC risk. rs3135395-T, rs477515-T and rs2395178-G were inversely associated with the generation of A1762T/G1764A, T1753V and C1653T, the HCC-risk HBV mutations. Multiplicative interactions of the variant genotypes with the HCC-risk HBV mutations were significantly associated with a decreased risk of HCC. In multivariate Cox analysis, rs477515-T independently predicted a favourable prognosis, with a hazard ratio of 0.48 (P = .002). The activity of the HLA-DRB1 enhancer with rs477515-T was significantly higher than that with rs477515-C. The activity of the HLA-DRB1 enhancer with rs477515-T and that with rs477515-C was significantly up-regulated by interferon-γ and interleukin-4, respectively. Interleukin-6 significantly inhibited the HLA-DRB1 enhancer activity, and this effect was more evident in those carrying rs477515-T. Polymorphisms predisposing to down-regulation of HLA-DR facilitate the Th1-to-Th2 transition and promote HCC development, possibly via selecting the HCC-risk HBV mutations. This can be transformed into specific prophylaxis of HCC.

Concordance in survival among first-degree relatives diagnosed with indolent lymphoid malignancies including chronic lymphocytic leukemia

Objectives

To investigate concordance in survival time among first‐degree relatives with lymphoid malignancies.

Methods

By linkage of national Swedish registers, we identified 66 430 patients diagnosed with a lymphoid malignancy 1958‐2016 with information on first‐degree relationships and follow‐up until 2017. Among these, we identified pairs of first‐degree relatives with any (N = 3326) or a similar (N = 690) lymphoid malignancy subtype. We defined survival in the first‐degree relative as good, expected, or poor based on tertiles of deviance residuals from a multivariable Cox regression model. Next, we used Cox regression to estimate hazard ratios (HR) of death with 95% confidence intervals (CI) among patients, using the survival of their first‐degree relative as exposure and adjusting for confounders.

Results

There was no concordance in survival among first‐degree relatives with any lymphoid malignancy (HR_good = 1.00 (reference), HR_Expected = 1.02, 95% CI: 0.89‐1.17, HR_Poor = 1.12, 95% CI: 0.98‐1.27, P_trend = .08). Among first‐degree relatives with indolent lymphoma, including chronic lymphocytic leukemia, those with a first‐degree relative to an expected or poor survival had worse outcome compared to those with a first‐degree relative with good survival (HR_Expected = 1.44, 95% CI: 0.82‐2.53, HR_Poor = 1.79, 95% CI: 1.07‐3.00, P_trend = .03).

Conclusion

Our results support a role of inherited factors in the outcome of indolent lymphoma, including chronic lymphocytic leukemia.

Active monitoring, radical prostatectomy and radical radiotherapy in PSA-detected clinically localised prostate cancer: the ProtecT three-arm RCT

BACKGROUND

Prostate cancer is the most common cancer among men in the UK. Prostate-specific antigen testing followed by biopsy leads to overdetection, overtreatment as well as undertreatment of the disease. Evidence of treatment effectiveness has lacked because of the paucity of randomised controlled trials comparing conventional treatments.

OBJECTIVES

To evaluate the effectiveness of conventional treatments for localised prostate cancer (active monitoring, radical prostatectomy and radical radiotherapy) in men aged 50-69 years.

DESIGN

A prospective, multicentre prostate-specific antigen testing programme followed by a randomised trial of treatment, with a comprehensive cohort follow-up.

SETTING

Prostate-specific antigen testing in primary care and treatment in nine urology departments in the UK.

PARTICIPANTS

Between 2001 and 2009, 228,966 men aged 50-69 years received an invitation to attend an appointment for information about the Prostate testing for cancer and Treatment (ProtecT) study and a prostate-specific antigen test; 82,429 men were tested, 2664 were diagnosed with localised prostate cancer, 1643 agreed to randomisation to active monitoring (n = 545), radical prostatectomy (n = 553) or radical radiotherapy (n = 545) and 997 chose a treatment.

INTERVENTIONS

The interventions were active monitoring, radical prostatectomy and radical radiotherapy.

TRIAL PRIMARY OUTCOME MEASURE

Definite or probable disease-specific mortality at the 10-year median follow-up in randomised participants.

SECONDARY OUTCOME MEASURES

Overall mortality, metastases, disease progression, treatment complications, resource utilisation and patient-reported outcomes.

RESULTS

There were no statistically significant differences between the groups for 17 prostate cancer-specific (p = 0.48) and 169 all-cause (p = 0.87) deaths. Eight men died of prostate cancer in the active monitoring group (1.5 per 1000 person-years, 95% confidence interval 0.7 to 3.0); five died of prostate cancer in the radical prostatectomy group (0.9 per 1000 person-years, 95% confidence interval 0.4 to 2.2 per 1000 person years) and four died of prostate cancer in the radical radiotherapy group (0.7 per 1000 person-years, 95% confidence interval 0.3 to 2.0 per 1000 person years). More men developed metastases in the active monitoring group than in the radical prostatectomy and radical radiotherapy groups: active monitoring, n = 33 (6.3 per 1000 person-years, 95% confidence interval 4.5 to 8.8); radical prostatectomy, n = 13 (2.4 per 1000 person-years, 95% confidence interval 1.4 to 4.2 per 1000 person years); and radical radiotherapy, n = 16 (3.0 per 1000 person-years, 95% confidence interval 1.9 to 4.9 per 1000 person-years; p = 0.004). There were higher rates of disease progression in the active monitoring group than in the radical prostatectomy and radical radiotherapy groups: active monitoring (n = 112; 22.9 per 1000 person-years, 95% confidence interval 19.0 to 27.5 per 1000 person years); radical prostatectomy (n = 46; 8.9 per 1000 person-years, 95% confidence interval 6.7 to 11.9 per 1000 person-years); and radical radiotherapy (n = 46; 9.0 per 1000 person-years, 95% confidence interval 6.7 to 12.0 per 1000 person years; p < 0.001). Radical prostatectomy had the greatest impact on sexual function/urinary continence and remained worse than radical radiotherapy and active monitoring. Radical radiotherapy's impact on sexual function was greatest at 6 months, but recovered somewhat in the majority of participants. Sexual and urinary function gradually declined in the active monitoring group. Bowel function was worse with radical radiotherapy at 6 months, but it recovered with the exception of bloody stools. Urinary voiding and nocturia worsened in the radical radiotherapy group at 6 months but recovered. Condition-specific quality-of-life effects mirrored functional changes. No differences in anxiety/depression or generic or cancer-related quality of life were found. At the National Institute for Health and Care Excellence threshold of £20,000 per quality-adjusted life-year, the probabilities that each arm was the most cost-effective option were 58% (radical radiotherapy), 32% (active monitoring) and 10% (radical prostatectomy).

LIMITATIONS

A single prostate-specific antigen test and transrectal ultrasound biopsies were used. There were very few non-white men in the trial. The majority of men had low- and intermediate-risk disease. Longer follow-up is needed.

CONCLUSIONS

At a median follow-up point of 10 years, prostate cancer-specific mortality was low, irrespective of the assigned treatment. Radical prostatectomy and radical radiotherapy reduced disease progression and metastases, but with side effects. Further work is needed to follow up participants at a median of 15 years.

TRIAL REGISTRATION

Current Controlled Trials ISRCTN20141297.

FUNDING

This project was funded by the National Institute for Health Research Health Technology Assessment programme and will be published in full in Health Technology Assessment; Vol. 24, No. 37. See the National Institute for Health Research Journals Library website for further project information.

GWAS in cancer: progress and challenges

The genome-wide association study (GWAS) is an effective method to detect single-nucleotide polymorphisms (SNPs) of multiple individual genes based on linkage disequilibrium (LD). GWAS examines genotypes and distinguishing gene characteristics that are exhibited in diseases. In the past few decades, more and more literature has reported the results of applying GWAS to study tumors. Although many pleiotropic loci associated with complex phenotypes have been identified by GWAS, the biological functions of many genetic variation loci remain unclear, and the genetic mechanisms of most complex phenotypes cannot be systematically explained. In this article, we will review the new findings of several tumor types, and categorize the new sites and mechanisms that have recently been discovered. We linked the mechanisms of action of various tumors and searched for links to related gene expression pathways. We found that susceptible sites can be divided into hub genes and peripheral genes; the two interact to link gene expression in a variety of diseases.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.

Precision Medicine in Lymphoma by Innovative Instrumental Platforms

In recent years, many efforts have been addressed to the growing field of precision medicine in order to offer individual treatments to every patient on the basis of his/her genetic background. Formerly adopted to achieve new disease classifications as it is still done, innovative platforms, such as microarrays, genome-wide association studies (GWAS), and next generation sequencing (NGS), have made the progress in pharmacogenetics faster and cheaper than previously expected. Several studies in lymphoma patients have demonstrated that these platforms can be used to identify biomarkers predictive of drug efficacy and tolerability, discovering new possible druggable proteins. Indeed, GWAS and NGS allow the investigation of the human genome, finding interesting associations with putative or unexpected targets, which in turns may represent new therapeutic possibilities. Importantly, some objective difficulties have initially hampered the translation of findings in clinical routines, such as the poor quantity/quality of genetic material or the paucity of targets that could be investigated at the same time. At present, some of these technical issues have been partially solved. Furthermore, these analyses are growing in parallel with the development of bioinformatics and its capabilities to manage and analyze big data. Because of pharmacogenetic markers may become important during drug development, regulatory authorities (i.e., EMA, FDA) are preparing ad hoc guidelines and recommendations to include the evaluation of genetic markers in clinical trials. Concerns and difficulties for the adoption of genetic testing in routine are still present, as well as affordability, reliability and the poor confidence of some patients for these tests. However, genetic testing based on predictive markers may offers many advantages to caregivers and patients and their introduction in clinical routine is justified.

Global Analysis of Enhancer Targets Reveals Convergent Enhancer-Driven Regulatory Modules

Single-cell screens enable high-throughput functional assessment of enhancers in their endogenous genomic context. However, the design of current studies limits their application to identifying the primary gene targets of enhancers. Here, we improve the experimental and computational parameters of single-cell enhancer screens to identify the secondary gene targets of enhancers. Our analysis of >500 putative enhancers in K562 cells reveals an interwoven enhancer-driven gene regulatory network. We find that enhancers from distinct genomic loci converge to modulate the expression of common sub-modules, including the α- and β-globin loci, by directly regulating transcription factors. Our analysis suggests that several genetic variants associated with myeloid blood cell traits alter the activity of a distal enhancer of MYB (~140 kb away), with downstream consequences on hemoglobin genes expression and cell state. These data have implications for the understanding of enhancer-associated traits and emphasize the flexibility of controlling transcriptional systems by modifying enhancer activity.

Xie et al. apply improved strategies for single-cell screens to identify an enhancer-driven transcriptional regulatory network in K562 cells. They demonstrate that the same group of genes can be indirectly regulated by enhancers from distinct genomic loci. These data have implications for the understanding of enhancer-associated traits.

Analysis of 153 115 patients with hematological malignancies refines the spectrum of familial risk

Estimating familial cancer risks is clinically important in being able to discriminate between individuals in the population at differing risk for malignancy. To gain insight into the familial risk for the different hematological malignancies and their possible inter-relationship, we analyzed data on more than 16 million individuals from the Swedish Family-Cancer Database. After identifying 153 115 patients diagnosed with a primary hematological malignancy, we quantified familial relative risks (FRRs) by calculating standardized incident ratios (SIRs) in 391 131 of their first-degree relatives. The majority of hematological malignancies showed increased FRRs for the same tumor type, with the highest FRRs being observed for mixed cellularity Hodgkin lymphoma (SIR, 16.7), lymphoplasmacytic lymphoma (SIR, 15.8), and mantle cell lymphoma (SIR, 13.3). There was evidence for pleiotropic relationships; notably, chronic lymphocytic leukemia was associated with an elevated familial risk for other B-cell tumors and myeloproliferative neoplasms. Collectively, these data provide evidence for shared etiological factors for many hematological malignancies and provide information for identifying individuals at increased risk, as well as informing future gene discovery initiatives.

The role of HLA variation in lymphoma aetiology and survival

Epidemiologic and laboratory evidence has consistently supported a strong inflammatory and immune component for lymphoma aetiology. These studies have consistently implicated variation in the immune gene, human leucocyte antigen (HLA), to be associated with lymphoma risk. In this review, we summarize the historical and recent evidence of HLA in both lymphoma aetiology and survival. The recent momentum in uncovering HLA associations has been propelled by the conduct of genome-wide association studies (GWAS), which has permitted the evaluation of imputed HLA alleles in much larger sample sizes than historically feasible with allelotyping studies. Based on the culmination of smaller HLA typing studies and larger GWAS, we now recognize several HLA associations with Hodgkin (HL) and non-Hodgkin lymphomas (NHLs) and their subtypes. Although other genetic variants have also been implicated with lymphoma risk, it is notable that HLA associations have been reported in every NHL and HL subtype evaluated to date. Both HLA class I and class II alleles have been linked with NHL and HL risk. It is notable that the associations identified are largely specific to each lymphoma subtype. However, pleiotropic HLA associations have also been observed. For example, rs10484561, which is in linkage disequilibrium with HLA-DRB1*01:01˜DQA1*01:01˜DQB1*05:01, has been implicated in increased FL and DLBCL risk. Opposing HLA associations across subtypes have also been reported, such as for HLA-A*01:01 which is associated with increased risk of EBV-positive cHL but decreased risk of EBV-negative cHL and chronic lymphocytic leukaemia/small cell lymphoma. Due to extensive linkage disequilibrium and allele/haplotypic variation across race/ethnicities, identification of causal alleles/haplotypes remains challenging. Follow-up functional studies are needed to identify the specific immunological pathways responsible in the multifactorial aetiology of HL and NHL. Correlative studies linking HLA alleles with known molecular subtypes and HLA expression in the tumours are also needed. Finally, additional association studies investigating HLA diversity and lymphoma survival are also required to replicate initial associations reported to date.

Primary Immunodeficiency and Cancer Predisposition Revisited: Embedding Two Closely Related Concepts Into an Integrative Conceptual Framework

Common understanding suggests that the normal function of a "healthy" immune system safe-guards and protects against the development of malignancies, whereas a genetically impaired one might increase the likelihood of their manifestation. This view is primarily based on and apparently supported by an increased incidence of such diseases in patients with specific forms of immunodeficiencies that are caused by high penetrant gene defects. As I will review and discuss herein, such constellations merely represent the tip of an iceberg. The overall situation is by far more varied and complex, especially if one takes into account the growing difficulties to define what actually constitutes an immunodeficiency and what defines a cancer predisposition. The enormous advances in genome sequencing, in bioinformatic analyses and in the functional in vitro and in vivo assessment of novel findings together with the availability of large databases provide us with a wealth of information that steadily increases the number of sequence variants that concur with clinically more or less recognizable immunological problems and their consequences. Since many of the newly identified hard-core defects are exceedingly rare, their tumor predisposing effect is difficult to ascertain. The analyses of large data sets, on the other hand, continuously supply us with low penetrant variants that, at least in statistical terms, are clearly tumor predisposing, although their specific relevance for the respective carriers still needs to be carefully assessed on an individual basis. Finally, defects and variants that affect the same gene families and pathways in both a constitutional and somatic setting underscore the fact that immunodeficiencies and cancer predisposition can be viewed as two closely related errors of development. Depending on the particular genetic and/or environmental context as well as the respective stage of development, the same changes can have either a neutral, predisposing and, in some instances, even a protective effect. To understand the interaction between the immune system, be it "normal" or "deficient" and tumor predisposition and development on a systemic level, one therefore needs to focus on the structure and dynamic functional organization of the entire immune system rather than on its isolated individual components alone.

Epstein–Barr Virus and Cancer

Epstein–Barr virus (EBV) contributes to about 1.5% of all cases of human cancer worldwide, and viral genes are expressed in the malignant cells. EBV also very efficiently causes the proliferation of infected human B lymphocytes. The functions of the viral proteins and small RNAs that may contribute to EBV-associated cancers are becoming increasingly clear, and a broader understanding of the sequence variation of the virus genome has helped to interpret their roles. The improved understanding of the mechanisms of these cancers means that there are great opportunities for the early diagnosis of treatable stages of EBV-associated cancers and the use of immunotherapy to target EBV-infected cells or overcome immune evasion. There is also scope for preventing disease by immunization and for developing therapeutic agents that target the EBV gene products expressed in the cancers.

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.

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.