Genome-wide association study identifies multiple susceptibility loci for diffuse large B cell lymphoma

Exosome secretion related gene signature predicts chemoresistance in patients with colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a highly heterogeneous malignancy, and patients often have different responses to treatment. In this study, the genetic characteristics related to exosome formation and secretion procedure were used to predict chemoresistance and guide the individualized treatment of patients.

METHODS

Firstly, seven microarray datasets in Gene Expression Omnibus (GEO) and RNA-Seq dataset from the Cancer Genome Atlas (TCGA) were used to analysis the transcriptome profiles and associated characteristics of CRC patients. Then, a predictive model based on gene features linked to exosome formation and secretion was created and validated using Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis and Support Vector Machine-Recursive Feature Elimination (SVM-RFE) machine learning. Finally, we evaluated the model using chemoresistant/chemosensitive cells and tissues by immunofluorescence (IF), western blot (WB), quantitative real-time PCR (qRT-PCR) and immunocytochemistry (IHC) experiments, and the predictive value of integrated model in the clinical validation cohort were performed by Receiver Operating Characteristic (ROC) and Kaplan-Meier (K-M) curves analyses.

RESULTS

We established a risk score signature based on three genes related to exosome secretion in CRC. Better Overall Survival (OS) and greater chemosensitivity were seen in the low-risk group, whereas the high-risk group exhibited chemoresistance and a subpar response to immune checkpoint blockade (ICB) therapy. Higher expression of the model genes EXOC2, EXOC3 and STX4 were observed in chemoresistant cells and specimens. The AUC of 5-year disease-free survival (DFS) was 0.804. Compared with that in the low-risk group, patients' DFS was found to be significantly worse in the high-risk group.

CONCLUSIONS

In summary, the gene signature related to exosome formation and secretion could reliably predict patients' chemosensitivity and ICB treatment response, which providing new independent biomarkers for the treatment of CRC.

HLA variations in patients with diffuse large B-cell lymphoma and association with disease risk and prognosis: a case-control study

Introduction

Human leukocyte antigen (HLA) polymorphisms have been associated with the development of various autoimmune diseases, as well as malignant neoplasms. Non-Hodgkin lymphomas (NHLs) are a heterogenous group of lymphoid malignancies in which a genetic substrate has been established and is deemed to play a crucial role in disease pathogenesis. This study aimed to identify whether variations in the HLA gene region were associated with diffuse large B-cell lymphoma (DLBCL) risk and prognosis.

Methods

We defined HLA class I (HLA-A, HLA-B, HLA-C) and class II (HLA-DRB1, HLA-DQB1) alleles in 60 patients with DLBCL and compared the results to those found by 236 healthy adult donors from the bone marrow bank of Northern Greece. HLA typing was performed by two molecular methods, Sequence - Specific Oligonucleotide HLA typing (SSO) and Sequence - Specific Primer HLA typing (SSP), from white blood cells recovered from peripheral blood. The phenotypic frequencies of HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 between patients and controls were compared with the 2-sided Fisher's exact test. Results with p-value <0.05 were considered statistically significant. Odds Ratios with 95% Confidence Intervals were calculated to further strengthen the results. The 2-sided Fisher's exact test was also applied to alleles found only in one of the two groups, while the odds ratios together with the confidence intervals were corrected with Haldane-Anscombe method.

Results

Among the studied HLA polymorphisms, the frequency HLA-C*12 allele was significantly lower in patients with DLBCL compared with control subjects (6.7% vs. 34.7%, OR = 0.16, 95% CI: 0.04-0.44). Frequency of HLA-B*39 was significantly lower in patients with DLBCL compared with controls, but due to the low frequency of this polymorphism in the studied population and small sample size, determinations regarding the significance of this findings were limited. Survival analysis revealed that the presence of HLA-C*12 was not associated with improved or worsened overall and progression-free survival. No statistically significant associations were observed in the phenotypic frequencies of HLA-A, HLA-DQB1, HLA-DRB1 and the rest of HLA-B alleles between the control and DLBCL groups.

Discussion

Collectively, our results provide valuable insight regarding the role of HLA variations on DLBCL risk. Further studies are required to consolidate our findings and ascertain the clinical implications of these genetic variations on DLBCL management and prognosis.

A genome-wide association study on hematopoietic stem cell transplantation reveals novel genomic loci associated with transplant outcomes

Introduction

Data on genomic susceptibility for adverse outcomes after hematopoietic stem cell transplantation (HSCT) for recipients are scarce.

Methods

We performed a genome wide association study (GWAS) to identify genes associated with survival/mortality, relapse, and severe graft-versus-host disease (sGvHD), fitting proportional hazard and subdistributional models to data of n=1,392 recipients of European ancestry from three centres.

Results

The single nucleotide polymorphism (SNP) rs17154454, intronic to the neuronal growth guidant semaphorin 3C gene (SEMA3C), was genome-wide significantly associated with event-free survival (p=7.0x10-8) and sGvHD (p=7.5x10-8). Further associations were detected for SNPs in the Paxillin gene (PXN) with death without prior relapse or sGvHD, as well as for SNPs of the Plasmacytoma Variant Translocation 1 gene (PVT1, a long non-coding RNA gene), the Melanocortin 5 Receptor (MC5R) gene and the WW Domain Containing Oxidoreductase gene (WWOX), all associated with the occurrence of sGvHD. Functional considerations support the observed associations.

Discussion

Thus, new genes were identified, potentially influencing the outcome of HSCT.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomization analysis

BACKGROUND

People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer.

METHODS

We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers.

RESULTS

We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08-1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship.

CONCLUSIONS

These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Genomic Analysis of Lymphoma Risk in Bullmastiff Dogs

Lymphoma is the most common haematological malignancy affecting dogs and has a high incidence in the Bullmastiff breed. The aim of this study was to identify risk loci predisposing this breed to the disease. The average age of lymphoma diagnosis in 55 cases was less than 6 years, similar to the median age of 64 cases from our clinical and pathology databases. When fine-scale population structure was explored using NETVIEW, cases were distributed throughout an extended pedigree. When genotyped cases (n = 49) and dogs from the control group (n = 281) were compared in a genome-wide association analysis of lymphoma risk, the most prominent associated regions were detected on CFA13 and CFA33. The top SNPs in a 5.4 Mb region on CFA13 were significant at a chromosome-wide level, and the region was fine-mapped to ~1.2 Mb (CFA13: 25.2-26.4 Mb; CanFam3.1) with four potential functional candidates, including the MYC proto-oncogene bHLH transcription factor (MYC) and a region syntenic with the human and mouse lncRNA Pvt1 oncogene (PVT1). A 380 Kb associated region at CFA33: 7.7-8.1 Mb contained the coding sequence for SUMO specific peptidase7 (SENP7) and NFK inhibitor zeta (NFKBIZ) genes. These genes have annotations related to cancer, amongst others, and both have functional links to MYC regulation. Genomic signatures identified in lymphoma cases suggest that increased risk contributed by the regions identified by GWAS may complement a complex predisposing genetic background.

Causal relationships between circulating inflammatory cytokines and diffuse large B cell lymphoma: a bidirectional Mendelian randomization study

Diffuse large B cell lymphoma (DLBCL) is the most common type of non-Hodgkin lymphoma. Studies indicated that inflammatory cytokines involved in the occurrence and progression of DLBCL and it is challenging to discern causality from the effects due to the presence of feedback loops. We conducted a bidirectional Mendelian randomization (MR) study to investigate the potential causal relationship between DLBCL and inflammatory cytokines. The genetic variants associated with inflammatory cytokines were obtained from a genome-wide association study (GWAS) involving 8293 European participants, and the data on 1010 individuals with DLBCL were sourced from the FinnGen consortium. The primary method employed in this study was the inverse-variance weighted (IVW) method, with supplementary analyses conducted using the MR-Egger, weighted median, and MR-PRESSO approaches. Based on the IVW method, genetically predicted that increasing level of Monokine induced by interferon gamma (MIG/CXC chemokine ligand 9, CXCL9) [OR: 1.31; 95% CI: 1.05-1.62; P = 0.01] and interferon gamma-induced protein 10(IP-10/CXC chemokine ligand 10, CXCL10) [OR: 1.30; 95% CI: 1.02-1.66; P = 0.03] showed suggestive associations with DLBCL risk. DLBCL may increase the level of macrophage colony-stimulating factor (M-CSF) [OR: 1.12; 95% CI: 1.01-1.2; P = 0.03], tumor necrosis factor beta (TNF-β) [OR: 1.16; 95% CI: 1.02-1.31; P = 0.02] and TNF-related apoptosis-inducing ligand (TRAIL) [OR: 1.07; 95% CI: 1.01-1.13; P = 0.02]. This study presents evidence supporting a causal relationship between inflammation cytokines and DLBCL. Specifically, MIG/CXCL9 and IP-10/CXCL10 were identified as indicators of upstream causes of DLBCL; while, DLBCL itself was found to elevate the levels of M-CSF, TNF-β, and TRAIL. These findings suggest that targeting specific inflammatory factors through regulation and intervention could serve as a potential approach for the treatment and prevention of DLBCL.

Epidemiology and etiology of diffuse large B-cell lymphoma

As the most common non-Hodgkin lymphoma subtype, diffuse large B-cell lymphoma (DLBCL) incidence patterns generally parallel that for NHL overall. Globally, DLBCL accounts for a third of all NHLs, ranging between 20% and 50% by country. Based on United States (U.S.) cancer registry data, age-standardized incidence rate for DLBCL was 7.2 per 100,000. DLBCL incidence rises with age and is generally higher in males than females; in the U.S., incidence is highest among non-Hispanic whites (9.2/100,000). Like NHL incidence, DLBCL incidence rose in the first half of the 20th century but has largely plateaued. However, there is some evidence that incidence rates are rising in areas of historically low rates, such as Asia; there are also estimates for rising DLBCL incidence in the near future due to the changing demographics in developed countries whose aging population is growing. Established risk factors for DLBCL include those that result in severe immune deficiency such as HIV/AIDS, inherited immunodeficiency syndromes, and organ transplant recipients. Factors that lead to chronic immune dysregulations are also established risk factors, and include a number of autoimmune conditions (eg, Sjögren syndrome, systemic lupus erythematosus, rheumatoid arthritis), viral infections (eg, HIV, KSHV/HHV8, HCV, EBV), and obesity. Family history of NHL/DLBCL, personal history of cancer, and multiple genetic susceptibility loci are also well-established risk factors for DLBCL. There is strong evidence for multiple environmental exposures in DLBCL etiology, including exposure to trichloroethylene, benzene, and pesticides and herbicides, with recent associations noted with glyphosate. There is also strong evidence for associations with other viruses, such as HBV. Recent estimates suggest that obesity accounts for nearly a quarter of DLBCLs that develop, but despite recent gains in the understanding of DLBCL etiology, the majority of disease remain unexplained. An understanding of the host and environmental contributions to disease etiology, and concerted efforts to expand our understanding to multiple race/ethnic groups, will be essential for constructing clinically relevant risk prediction models and develop effective strategies for disease prevention.

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).

Polygenic Risk and Chemotherapy-Related Subsequent Malignancies in Childhood Cancer Survivors: A Childhood Cancer Survivor Study and St Jude Lifetime Cohort Study Report

PURPOSE

Chemotherapeutic exposures are associated with subsequent malignant neoplasm (SMN) risk. The role of genetic susceptibility in chemotherapy-related SMNs should be defined as use of radiation therapy (RT) decreases.

PATIENTS AND METHODS

SMNs among long-term childhood cancer survivors of European (EUR; N = 9,895) and African (AFR; N = 718) genetic ancestry from the Childhood Cancer Survivor Study and St Jude Lifetime Cohort Study were evaluated. An externally validated 179-variant polygenic risk score (PRS) associated with pleiotropic adult cancer risk from the UK Biobank Study (N > 400,000) was computed for each survivor. SMN cumulative incidence comparing top and bottom PRS quintiles was estimated, along with hazard ratios (HRs) from proportional hazards models.

RESULTS

A total of 1,594 survivors developed SMNs, with basal cell carcinomas (n = 822), breast cancers (n = 235), and thyroid cancers (n = 221) being the most frequent. Although SMN risk associations with the PRS were extremely modest in RT-exposed EUR survivors (HR, 1.22; P = .048; n = 4,630), the increase in 30-year SMN cumulative incidence and HRs comparing top and bottom PRS quintiles was statistically significant among nonirradiated EUR survivors (n = 4,322) treated with alkylating agents (17% v 6%; HR, 2.46; P < .01), anthracyclines (20% v 8%; HR, 2.86; P < .001), epipodophyllotoxins (23% v 1%; HR, 12.20; P < .001), or platinums (46% v 7%; HR, 8.58; P < .01). This PRS also significantly modified epipodophyllotoxin-related SMN risk among nonirradiated AFR survivors (n = 414; P < .01). Improvements in prediction attributable to the PRS were greatest for epipodophyllotoxin-exposed (AUC, 0.71 v 0.63) and platinum-exposed (AUC,0.68 v 0.58) survivors.

CONCLUSION

A pleiotropic cancer PRS has strong potential for improving SMN clinical risk stratification among nonirradiated survivors treated with specific chemotherapies. A polygenic risk screening approach may be a valuable complement to an early screening strategy on the basis of treatments and rare cancer-susceptibility mutations.

The causal effect of inflammatory bowel disease on diffuse large B-cell lymphoma: two-sample Mendelian randomization study

Background

It has been reported that inflammatory bowel disease (IBD) is associated with an increased risk of malignancies, including lymphoma. A number of large observational studies have been devoted to exploring the causal link between IBD and malignant lymphoma. However, no consensus exists on whether there is a causal relationship between IBD and malignant lymphoma.

Methods

The summary dataset of the IBD and lymphoma genome-wide association studies (GWAS) was obtained from the OPEN GWAS website. Single-nucleotide polymorphisms (SNPs) were selected as genetic instrumental variants (IVs) for fulling P < 5 × 10-8 and linkage disequilibrium (LD) of r2 = 0.001 in the IBD GWAS. The proxy SNPs with LD of r2 > 0.8 were identified. Palindromic SNPs and outlier SNPs were excluded. The assessments of sensitivity employed the Cochran's Q test, Mendelian randomization (MR)-Egger intercept test, and leave-one-out analysis.

Results

The MR analysis results proved the causality of IBD on diffuse large B-cell lymphoma (DLBCL). The risk of developing DLBCL is increased by 28.6% in patients with IBD [odds ratio (OR)IVW = 1.286, 95% confidence interval (CI) 1.066-1.552, P = 0.009]. The results of the subgroup analysis showed that Crohn's disease (ORIVW = 1.218, 95% CI 1.030-1.441, P = 0.021) rather than ulcerative colitis (ORIVW = 1.206, 95% CI 0.984-1.478, P = 0.072) had a causal effect on DLBCL. No horizontal and directional pleiotropy was observed in the MR studies.

Conclusions

The above MR study concluded that IBD itself is causally responsible for DLBCL, especially Crohn's disease. Further investigations are needed to elucidate the mechanism underlying this direct causal link.

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.

Gene and pathway based burden analyses in familial lymphoid cancer cases: Rare variants in immune pathway genes

Genome-wide association studies have revealed common genetic variants with small effect sizes associated with diverse lymphoid cancers. Family studies have uncovered rare variants with high effect sizes. However, these variants explain only a portion of the heritability of these cancers. Some of the missing heritability may be attributable to rare variants with small effect sizes. We aim to identify rare germline variants associated with familial lymphoid cancers using exome sequencing. One case per family was selected from 39 lymphoid cancer families based on early onset of disease or rarity of subtype. Control data was from Non-Finnish Europeans in gnomAD exomes (N = 56,885) or ExAC (N = 33,370). Gene and pathway-based burden tests for rare variants were performed using TRAPD. Five putatively pathogenic germline variants were found in four genes: INTU, PEX7, EHHADH, and ASXL1. Pathway-based association tests identified the innate and adaptive immune systems, peroxisomal pathway and olfactory receptor pathway as associated with lymphoid cancers in familial cases. Our results suggest that rare inherited defects in the genes involved in immune system and peroxisomal pathway may predispose individuals to lymphoid cancers.

Causal relationships between risk of venous thromboembolism and 18 cancers: a bidirectional Mendelian randomisation analysis

Background

People with cancer experience high rates of venous thromboembolism (VTE). Additionally, risk of subsequent cancer is increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer.

Methods

We used data from large genome-wide association study meta-analyses to perform bi-directional Mendelian randomisation analyses to estimate causal associations between genetically-proxied lifetime risk of VTE and risk of 18 different cancers.

Results

We found no conclusive evidence that genetically-proxied lifetime risk of VTE was causally associated with an increased incidence of cancer, or vice-versa. We observed an association between VTE and pancreatic cancer risk (odds ratio for pancreatic cancer 1.23 (95% confidence interval 1.08 - 1.40) per log-odds increase in risk of VTE, P = 0.002). However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence from Mendelian randomisation to suggest a causal relationship.

Conclusions

These findings do not support the hypothesis that genetically-proxied lifetime risk of VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesise evidence for these mechanisms.

Selective Enhancer Dependencies in MYC -Intact and MYC -Rearranged Germinal Center B-cell Diffuse Large B-cell Lymphoma

High expression of MYC and its target genes define a subset of germinal center B-cell diffuse large B-cell lymphoma (GCB-DLBCL) associated with poor outcomes. Half of these high-grade cases show chromosomal rearrangements between the MYC locus and heterologous enhancer-bearing loci, while focal deletions of the adjacent non-coding gene PVT1 are enriched in MYC -intact cases. To identify genomic drivers of MYC activation, we used high-throughput CRISPR-interference (CRISPRi) profiling of candidate enhancers in the MYC locus and rearrangement partner loci in GCB-DLBCL cell lines and mantle cell lymphoma (MCL) comparators that lacked common rearrangements between MYC and immunoglobulin (Ig) loci. Rearrangements between MYC and non-Ig loci were associated with unique dependencies on specific enhancer subunits within those partner loci. Notably, fitness dependency on enhancer modules within the BCL6 super-enhancer ( BCL6 -SE) cluster regulated by a transcription factor complex of MEF2B, POU2F2, and POU2AF1 was higher in cell lines bearing a recurrent MYC::BCL6 -SE rearrangement. In contrast, GCB-DLBCL cell lines without MYC rearrangement were highly dependent on a previously uncharacterized 3' enhancer within the MYC locus itself (GCBME-1), that is regulated in part by the same triad of factors. GCBME-1 is evolutionarily conserved and active in normal germinal center B cells in humans and mice, suggesting a key role in normal germinal center B cell biology. Finally, we show that the PVT1 promoter limits MYC activation by either native or heterologous enhancers and demonstrate that this limitation is bypassed by 3' rearrangements that remove PVT1 from its position in cis with the rearranged MYC gene.

Key points

CRISPR-interference screens identify a conserved germinal center B cell MYC enhancer that is essential for GCB-DLBCL lacking MYC rearrangements. Functional profiling of MYC partner loci reveals principles of MYC enhancer-hijacking activation by non-immunoglobulin rearrangements.

Circulating immune markers and risks of non-Hodgkin lymphoma subtypes: A pooled analysis

Although prediagnostic circulating concentrations of the immune activation markers soluble CD27 (sCD27), sCD30 and chemokine ligand-13 (CXCL13) have been associated with non-Hodgkin lymphoma (NHL) risk, studies have been limited by sample size in associations with NHL subtypes. We pooled data from eight nested case-control studies to investigate subtype-specific relationships for these analytes. Using polytomous regression, we calculated odds ratios (ORs) with 95% confidence intervals (CIs) relating study-specific analyte tertiles to selected subtypes vs controls (n = 3310): chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL; n = 623), diffuse large B cell lymphoma (DLBCL; n = 621), follicular lymphoma (FL; n = 398), marginal zone lymphoma (MZL; n = 138), mantle cell lymphoma (MCL; n = 82) and T cell lymphoma (TCL; n = 92). We observed associations with DLBCL for elevated sCD27 [OR for third vs first tertile (ORT3 ) = 2.2, 95% CI = 1.6-3.1], sCD30 (ORT3  = 2.0, 95% CI = 1.6-2.5) and CXCL13 (ORT3  = 2.3, 95% CI = 1.8-3.0). We also observed associations with sCD27 for CLL/SLL (ORT3  = 3.3, 95% CI = 2.4-4.6), MZL (ORT3  = 7.7, 95% CI = 3.0-20.1) and TCL (ORT3  = 3.4, 95% CI = 1.5-7.7), and between sCD30 and FL (ORT3  = 2.7, 95% CI = 2.0-3.5). In analyses stratified by time from phlebotomy to case diagnosis, the sCD27-TCL and all three DLBCL associations were equivalent across both follow-up periods (<7.5, ≥7.5 years). For other analyte-subtype comparisons, associations were stronger for the follow-up period closer to phlebotomy, particularly for indolent subtypes. In conclusion, we found robust evidence of an association between these immune markers and DLBCL, consistent with hypotheses that mechanisms related to immune activation are important in its pathogenesis. Our other findings, particularly for the rarer subtypes MZL and TCL, require further investigation.

The epidemiology of Waldenström macroglobulinemia

Waldenström macroglobulinemia (WM) is a rare subtype of non-Hodgkin lymphoma characterized by the presence of lymphoplasmacytic lymphoma (LPL) in the bone marrow accompanied by a monoclonal immunoglobulin type M (IgM) in the serum. WM was first described only 80 years ago and became reportable in the US as a malignancy in 1988. Very little systematic research was conducted prior to 2000 to characterize incidence, clinical characteristics, risk factors or diagnostic and prognostic criteria, and there were essentially no WM-specific clinical interventional trials. Since the inaugural meeting of the International Workshop in Waldenström's Macroglobulinemia (IWWM) in 2000, WM has become the focus of a steadily increasing and productive body of research, engaging a growing number of investigators throughout the world. This introductory overview provides summary of the current understanding of the epidemiology of WM/LPL as a backdrop for a series of consensus panel recommendations arising from research presented at the 11th IWWM.

Functional impact of multi-omic interactions in breast cancer subtypes

Multi-omic approaches are expected to deliver a broader molecular view of cancer. However, the promised mechanistic explanations have not quite settled yet. Here, we propose a theoretical and computational analysis framework to semi-automatically produce network models of the regulatory constraints influencing a biological function. This way, we identified functions significantly enriched on the analyzed omics and described associated features, for each of the four breast cancer molecular subtypes. For instance, we identified functions sustaining over-representation of invasion-related processes in the basal subtype and DNA modification processes in the normal tissue. We found limited overlap on the omics-associated functions between subtypes; however, a startling feature intersection within subtype functions also emerged. The examples presented highlight new, potentially regulatory features, with sound biological reasons to expect a connection with the functions. Multi-omic regulatory networks thus constitute reliable models of the way omics are connected, demonstrating a capability for systematic generation of mechanistic hypothesis.

Selecting the appropriate hurdles and endpoints for pentilludin, a novel antiaddiction pharmacotherapeutic targeting the receptor type protein tyrosine phosphatase D

Substance use disorders provide challenges for development of effective medications. Use of abused substances is likely initiated, sustained and "quit" by complex brain and pharmacological mechanisms that have both genetic and environmental determinants. Medical utilities of prescribed stimulants and opioids provide complex challenges for prevention: how can we minimize their contribution to substance use disorders while retaining medical benefits for pain, restless leg syndrome, attention deficit hyperactivity disorder, narcolepsy and other indications. Data required to support assessments of reduced abuse liability and resulting regulatory scheduling differs from information required to support licensing of novel prophylactic or therapeutic anti-addiction medications, adding further complexity and challenges. I describe some of these challenges in the context of our current efforts to develop pentilludin as a novel anti-addiction therapeutic for a target that is strongly supported by human and mouse genetic and pharmacologic studies, the receptor type protein tyrosine phosphatase D (PTPRD).

COVID-19 and non-Hodgkin's lymphoma: A common susceptibility pattern?

OBJECTIVE

To explore the link between COVID-19 incidence, socio-economic covariates, and NHL incidence.

DESIGN

Ecological study design.

SETTING

Sardinia, Italy.

PARTICIPANTS

We used official reports on the total cases of COVID-19 in 2020, published data on NHL incidence, and socio-economic indicators by administrative unit, covering the whole regional population.

MAIN OUTCOMES AND MEASURES

We used multivariable regression analysis to explore the association between the natural logarithm (ln) of the 2020 cumulative incidence of COVID-19 and the ln-transformed NHL incidence in 1974-2003, weighing by population size and adjusting by socioeconomic deprivation and other covariates.

RESULTS

The cumulative incidence of COVID-19 increased in relation to past incidence of NHL (p < 0.001), socioeconomic deprivation (p = 0.006), and proportion of elderly residents (p < 0.001) and decreased with urban residency (p = 0.001). Several sensitivity analyses confirmed the finding of an association between COVID-19 and NHL.

CONCLUSION

This ecological study found an ecological association between NHL and COVID-19. If further investigation would confirm our findings, shared susceptibility factors should be investigated among the plausible underlying mechanisms.

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

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

Evolutionary Origins of Metabolic Reprogramming in Cancer

Metabolic changes that facilitate tumor growth are one of the hallmarks of cancer. These changes are not specific to tumors but also take place during the physiological growth of tissues. Indeed, the cellular and tissue mechanisms present in the tumor have their physiological counterpart in the repair of tissue lesions and wound healing. These molecular mechanisms have been acquired during metazoan evolution, first to eliminate the infection of the tissue injury, then to enter an effective regenerative phase. Cancer itself could be considered a phenomenon of antagonistic pleiotropy of the genes involved in effective tissue repair. Cancer and tissue repair are complex traits that share many intermediate phenotypes at the molecular, cellular, and tissue levels, and all of these are integrated within a Systems Biology structure. Complex traits are influenced by a multitude of common genes, each with a weak effect. This polygenic component of complex traits is mainly unknown and so makes up part of the missing heritability. Here, we try to integrate these different perspectives from the point of view of the metabolic changes observed in cancer.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

Background

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

Case presentation

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

Conclusions

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

Typical Enhancers, Super-Enhancers, and Cancers

Simple Summary

The cancer genome has been exhaustively studied upon the advent of Next-Generation Sequencing technologies. Coding and non-coding sequences have been defined as hotspots of genomic variations that affect the naïve gene expression programs established in normal cells, thus working as endogenous drivers of carcinogenesis. In this review, we comprehensively summarize fundamental aspects of gene expression regulation, with emphasis on the impact of sequence and structural variations mapped across non-coding cis-acting elements of genes encoding for tumor-related transcription factors. Chromatin architecture, epigenome reprogramming, transcriptional enhancers and Super-enhancers, oncogene regulation, cutting-edge technologies, and pharmacological treatment are substantially highlighted.

Abstract

Non-coding segments of the human genome are enriched in cis-regulatory modules that constitute functional elements, such as transcriptional enhancers and Super-enhancers. A hallmark of cancer pathogenesis is the dramatic dysregulation of the “archetype” gene expression profiles of normal human cells. Genomic variations can promote such deficiencies when occurring across enhancers and Super-enhancers, since they affect their mechanistic principles, their functional capacity and specificity, and the epigenomic features of the chromatin microenvironment across which these regulatory elements reside. Here, we comprehensively describe: fundamental mechanisms of gene expression dysregulation in cancers that involve genomic abnormalities within enhancers’ and Super-enhancers’ (SEs) sequences, which alter the expression of oncogenic transcription factors (TFs); cutting-edge technologies applied for the analysis of variation-enriched hotspots of the cancer genome; and pharmacological approaches for the treatment of Super-enhancers’ aberrant function. Finally, we provide an intratumor meta-analysis, which highlights that genomic variations in transcription-factor-driven tumors are accompanied overexpression of genes, a portion of which encodes for additional cancer-related transcription factors.

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.

Genomic characterization of lymphomas in patients with inborn errors of immunity

Inborn errors of immunity-associated lymphomas are characterized by distinct clinical features and genetic signatures.

Both germline and somatic alterations contribute to lymphomagenesis in patients with inborn errors of immunity.

Patients with inborn errors of immunity (IEI) have a higher risk of developing cancer, especially lymphoma. However, the molecular basis for IEI-related lymphoma is complex and remains elusive. Here, we perform an in-depth analysis of lymphoma genomes derived from 23 IEI patients. We identified and validated disease-causing or -associated germline mutations in 14 of 23 patients involving ATM, BACH2, BLM, CD70, G6PD, NBN, PIK3CD, PTEN, and TNFRSF13B. Furthermore, we profiled somatic mutations in the lymphoma genome and identified 8 genes that were mutated at a significantly higher level in IEI-associated diffuse large B-cell lymphomas (DLBCLs) than in non-IEI DLBCLs, such as BRCA2, NCOR1, KLF2, FAS, CCND3, and BRWD3. The latter, BRWD3, is furthermore preferentially mutated in tumors of a subgroup of activated phosphoinositide 3-kinase δ syndrome patients. We also identified 5 genomic mutational signatures, including 2 DNA repair deficiency-related signatures, in IEI-associated lymphomas and a strikingly high number of inter- and intrachromosomal structural variants in the tumor genome of a Bloom syndrome patient. In summary, our comprehensive genomic characterization of lymphomas derived from patients with rare genetic disorders expands our understanding of lymphomagenesis and provides new insights for targeted therapy.

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

BACKGROUND

A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non-Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes.

METHODS

In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell-mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression.

RESULTS

We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08-1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59-2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell-mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction.

CONCLUSIONS

Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions.

IMPACT

Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.

The association between polymorphism of the long noncoding RNA, Plasmacytoma variant translocation 1, and the risk of gastric cancer

Genetic polymorphisms of plasmacytoma variant translocation 1 can affect various tumors including gastro-intestinal, sexual hormone sensitive cancers and lymphoma. Accumulated evidence have shown that plasmacytoma variant translocation 1 acts as an oncogene and tumor suppressor in various cancers. In fact, the rs13255292 and rs2608053 single nucleotide polymorphisms of plasmacytoma variant translocation 1are known to affect lymphoma; however, their effects on gastric cancer are primarily unknown. In this study, we evaluated the association between these plasmacytoma variant translocation 1 polymorphisms and the risk of gastric cancer.In the present study, 462 patients diagnosed with gastric cancer and 377 cancer-free controls were enrolled. The TaqMan genotyping assay was used to analyze the association between rs13255292 and rs2608053 single nucleotide polymorphisms and the risk of gastric cancer.The rs2608053 dominant model (CT + TT) was associated with a decreased risk of gastric cancer in T3 + T4 (odds ratio [OR] = 0.61, confidence interval (CI) = 0.41 - 0.92, P = .019), and stage III Gastric cancer subgroups (OR = 0.59, 95% CI = 0.38 - 0.91, P = .017) compared to the CC genotype. When stratified analysis by sex was carried out, the rs13255292 dominant model (CT + TT) had a significant association with an increased risk of gastric cancer in the female negative lymph node metastasis gastric cancer subgroup, compared to the CC genotype (OR = 1.96, 95% CI = 1.16 - 3.30, P = .012). The recessive model (TT) of rs13255292 was associated with an increased risk of gastric cancer in the male T3 + T4 gastric cancer subgroups compared to the CC + CT genotype (OR = 3.82, 95% CI = 1.02 - 14.33, P = .047). The dominant model (CT + TT) of rs2608053 was related to a decreased risk of gastric cancer in male T3 + T4 (OR = 0.57, 95% CI = 0.33 - 0.98, P = .042) and stage III gastric cancer subgroups (OR = 0.49, 95% CI = 0.27 - 0.89, P = .020) compared to the CC genotype.The rs13255292 and rs2608053 single nucleotide polymorphisms in plasmacytoma variant translocation 1 may contribute to susceptibility of gastric cancer. Further studies with more subjects and different ethnic groups are needed to validate our results.

Role of PVT1 polymorphisms in the glioma susceptibility and prognosis

BACKGROUND

Genetic factors play a crucial role in the glioma risk and prognosis of glioma patients. To explore the role of plasmacytoma variant translocation 1 (PVT1) polymorphism in the susceptibility and survival of glioma in the Chinese Han population, we conducted a case-control study.

METHODS

The three single-nucleotide polymorphisms (SNPs) in PVT1 were genotyped using Agena MassARRAY from 575 patients with glioma and 500 healthy controls. We used the χ2 test to analyze the differences in distribution of allele and genotype between the cases and controls. Odds ratio and 95% confidence interval (CI) were calculated by logistic regression analysis to evaluate the association SNPs with glioma risk. The effects of polymorphisms and clinical features on survival of glioma patients were evaluated using the log-rank test, Kaplan-Meier and Cox regression analysis.

RESULTS

We found that rs13255292 was associated with a decreased risk of glioma in the recessive model in overall or male; and rs4410871 was significantly associated with an increased the risk of glioma in age ≤40 years old or female. Moreover, the extent of resection and chemotherapy were found to be key prognostic factors in survival of glioma patients. However, the gender, age, tumor grade, radiotherapy and PVT1 polymorphisms have no effect on prognosis of glioma patients.

CONCLUSIONS

Our results indicated that PVT1 polymorphisms (rs13255292 and rs4410871) were associated with glioma susceptibility, but have no effect on prognosis of glioma patients. Further studies with large samples are required to confirm the results.

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.

Genome-wide homozygosity and risk of four non-Hodgkin lymphoma subtypes

AIM

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

METHODS

We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.

RESULTS

We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10-6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10-5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.

CONCLUSION

Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

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.

Genome instability and lymphoma

Lymphoma is one of the most common malignant tumor of the hematologic system. The genome instability is not only an important molecular basis for the development of lymphoma, but also has important value in the diagnosis and prognosis of lymphoma. There are 2 types of genome instability: Microsatellite instability (MSI/MIN) at gene level and chromosomal instability at chromosome level. Through the study on genes associated with lymphoma, the unstable genes associated with lymphoma could be found, meanwhile the mechanism of its occurrence and development of lymphoma could be explored, and the important basis of molecular biology could also be provided in the field of current hot lymphoma precision medical research.

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.

Diffuse Large B-Cell Lymphoma

Large B-cell lymphomas, with an estimated 150,000 new cases annually worldwide, represent almost 30% of all cases of non-Hodgkin’s lymphoma. Patients typically present with progressive lymphadenopathy, extranodal disease, or both and require therapy. Despite the advanced stage at presentation in the majority of patients, more than 60% can be cured with R-CHOP (rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone) immunochemotherapy (Fig. 1A). Patients with treatment failure after R-CHOP often have a poor outcome — in particular, those with disease that is refractory to frontline or subsequent therapies — although some patients can have a durable remission and be cured after secondary therapies. Over the past two decades, improved insights into large B-cell lymphomas, in terms of epidemiology, prognostic factors, and biologic heterogeneity, have led to a refinement of disease classification and the development of new therapeutic approaches.

Cross-cancer evaluation of polygenic risk scores for 16 cancer types in two large cohorts

Even distinct cancer types share biological hallmarks. Here, we investigate polygenic risk score (PRS)-specific pleiotropy across 16 cancers in European ancestry individuals from the Genetic Epidemiology Research on Adult Health and Aging cohort (16,012 cases, 50,552 controls) and UK Biobank (48,969 cases, 359,802 controls). Within cohorts, each PRS is evaluated in multivariable logistic regression models against all other cancer types. Results are then meta-analyzed across cohorts. Ten positive and one inverse cross-cancer associations are found after multiple testing correction. Two pairs show bidirectional associations; the melanoma PRS is positively associated with oral cavity/pharyngeal cancer and vice versa, whereas the lung cancer PRS is positively associated with oral cavity/pharyngeal cancer, and the oral cavity/pharyngeal cancer PRS is inversely associated with lung cancer. Overall, we validate known, and uncover previously unreported, patterns of pleiotropy that have the potential to inform investigations of risk prediction, shared etiology, and precision cancer prevention strategies.

Tracking the Genetic Susceptibility Background of B-Cell Non-Hodgkin's Lymphomas from Genome-Wide Association Studies

B-cell non-Hodgkin’s lymphoma (NHL) risk associations had been mainly attributed to family history of the disease, inflammation, and immune components including human leukocyte antigen (HLA) genetic variations. Nevertheless, a broad range of genome-wide association studies (GWAS) have shed light into the identification of several genetic variants presumptively associated with B-cell NHL etiologies, survival or shared genetic risk with other diseases. The present review aims to overview HLA structure and diversity and summarize the evidence of genetic variations, by GWAS, on five NHL subtypes (diffuse large B-cell lymphoma DLBCL, follicular lymphoma FL, chronic lymphocytic leukemia CLL, marginal zone lymphoma MZL, and primary central nervous system lymphoma PCNSL). Evidence indicates that the HLA zygosity status in B-cell NHL might promote immune escape and that genome-wide significance variants can give biological insight but also potential therapeutic markers such as WEE1 in DLBCL. However, additional studies are needed, especially for non-DLBCL, to replicate the associations found to date.

Pigmentary traits, sun exposure, and risk of non-Hodgkin's lymphoma/chronic lymphocytic leukemia: A study within the French E3N prospective cohort

To investigate whether risk factors for keratinocyte carcinomas (KCs), namely pigmentary traits and sun exposure, are associated with risk of non‐Hodgkin's lymphomas (NHL) and chronic lymphocytic leukemia (CLL). E3N is a prospective cohort of French women aged 40–65 years at inclusion in 1990. Cancer data were collected at baseline and updated every 2–3 years. Hazard Ratios (HRs) and 95% confidence intervals (CIs) for associations between pigmentary traits and sun exposure, and risk of CLL/NHL were estimated using Cox models. With a median follow‐up of 24 years, 622 incident cases of CLL/NHL were ascertained among the 92,097 included women. The presence of nevi was associated with CLL/NHL risk: HR for “many or very many nevi” relative to “no nevi”: 1.56 [1.15; 2.11]. Such association with number of nevi appears to be mostly limited to risk of CLL: HR for “many or very many nevi”: 3.00 [1.38; 6.52]; versus 1.32 [0.94; 1.84] for NHL. Women whose skin was highly sensitive to sunburn also had a higher risk of CLL: HR = 1.96 [1.21; 3.18], while no increase in risk of NHL was observed. Skin or hair color, number of freckles, and average daily ultraviolet (UV) dose during spring and summer in location of residence at birth or at inclusion (kJ/m²) were not associated with CLL/NHL risk. Some pigmentary traits (presence of nevi and skin sensitivity), but not sun exposure, were associated with CLL/NHL. These observations suggest that CLL may share some constitutional risk factors with keratinocyte cancers.

Mutations in the ARAP3 Gene in Three Families with Primary Lymphedema Negative for Mutations in Known Lymphedema-Associated Genes

Background

ARAP3 is a small GTPase-activating protein regulator, which has important functions in lymphatic vessel organogenesis and modulation of cell adhesion and migration. Mutations in the ARAP3 gene are associated with impaired lymphatic vessel formation.

Objective

The aim of our study was to determine the genotypes of lymphedema patients in relation to variants in the ARAP3 gene in order to explore its role in the development of lymphedema.

Methods and Results

We applied next-generation sequencing to DNA samples of a cohort of 246 Italian patients with lymphatic malformations. When we tested probands for known lymphedema genes, 235 out of 246 were negative. Retrospectively, we tested the DNA of these 235 patients for new candidate lymphedema-associated genes, including ARAP3. Three out of 235 probands proved to carry rare missense heterozygous variants in ARAP3. In the case of two families, other family members were also tested and proved negative for the ARAP3 variant, besides being unaffected by lymphedema. According to in silico analysis, alterations due to these variants have a significant impact on the overall structure and stability of the resulting proteins.

Conclusions

Based on our results, we propose that variants in ARAP3 could be included in genetic testing for lymphedema.

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.

Insecticide use and risk of non-Hodgkin lymphoma subtypes: A subset meta-analysis of the North American Pooled Project

Insecticide use has been linked to increased risk of non-Hodgkin lymphoma (NHL), however, findings of epidemiologic studies have been inconsistent, particularly for NHL subtypes. We analyzed 1690 NHL cases and 5131 controls in the North American Pooled Project (NAPP) to investigate self-reported insecticide use and risk of NHL overall and by subtypes: follicular lymphoma (FL), diffuse large B-cell lymphoma (DLBCL) and small lymphocytic lymphoma (SLL). Odds ratios (OR) and 95% confidence intervals for each insecticide were estimated using logistic regression. Subtype-specific associations were evaluated using ASSET (Association analysis for SubSETs). Increased risks of multiple NHL subtypes were observed for lindane (OR = 1.60, 1.20-2.10: FL, DLCBL, SLL), chlordane (OR = 1.59, 1.17-2.16: FL, SLL) and DDT (OR = 1.36, 1.06-1.73: DLBCL, SLL). Positive trends were observed, within the subsets with identified associations, for increasing categories of exposure duration for lindane (P_trend = 1.7 × 10^-4 ), chlordane (P_trend = 1.0 × 10^-3 ) and DDT (P_trend = 4.2 × 10^-3 ), however, the exposure-response relationship was nonlinear. Ever use of pyrethrum was associated with an increased risk of FL (OR = 3.65, 1.45-9.15), and the relationship with duration of use appeared monotonic (OR for >10 years: OR = 5.38, 1.75-16.53; P_trend = 3.6 × 10^-3 ). Our analysis identified several novel associations between insecticide use and specific NHL subtypes, suggesting possible etiologic heterogeneity in the context of pesticide exposure.

Genome-defined African ancestry is associated with distinct mutations and worse survival in patients with diffuse large B-cell lymphoma

BACKGROUND

Significant racial differences have been observed in the incidence and clinical outcomes of diffuse large B-cell lymphoma (DLBCL) in the United States, but to the authors' knowledge it remains unclear whether genomic differences contribute to these disparities.

METHODS

To understand the influences of genetic ancestry on tumor genomic alterations, the authors estimated the genetic ancestry of 1001 previously described patients with DLBCL using unsupervised model-based Admixture global ancestry analysis applied to exome sequencing data and examined the mutational profile of 150 DLBCL driver genes in tumors obtained from this cohort.

RESULTS

Global ancestry prediction identified 619 patients with >90% European ancestry, 81 patients with >90% African ancestry, and 50 patients with >90% Asian ancestry. Compared with patients with DLBCL with European ancestry, patients with African ancestry were aged >10 years younger at the time of diagnosis and were more likely to present with B symptoms, elevated serum lactate dehydrogenase, extranodal disease, and advanced stage disease. Patients with African ancestry demonstrated worse overall survival compared with patients with European ancestry (median, 4.9 years vs 8.8 years; P = .04). Recurrent mutations of MLL2 (KMT2D), HIST1H1E, MYD88, BCL2, and PIM1 were found across all ancestry groups, suggesting shared mechanisms underlying tumor biology. The authors also identified 6 DLBCL driver genes that were more commonly mutated in patients with African ancestry compared with patients with European ancestry: ATM (21.0% vs 7.75%; P < .001), MGA (19.7% vs 5.33%; P < .001), SETD2 (17.3% vs 5.17%; P < .001), TET2 (12.3% vs 5.82%; P = .029), MLL3 (KMT2C) (11.1% vs 4.36%; P = .013), and DNMT3A (11.1% vs 4.52%; P = .016).

CONCLUSIONS

Distinct prevalence and patterns of mutation highlight an important difference in the mutational landscapes of DLBCL arising in different ancestry groups. To the authors' knowledge, the results of the current study provide the first-ever characterization of genetic alterations among patients with African descent who are diagnosed with DLBCL.

Significance of Single-Nucleotide Variants in Long Intergenic Non-protein Coding RNAs

Single-nucleotide variants (SNVs) are the most common genetic variants and universally present in the human genome. Genome-wide association studies (GWASs) have identified a great number of disease or trait-associated variants, many of which are located in non-coding regions. Long intergenic non-protein coding RNAs (lincRNAs) are the major subtype of long non-coding RNAs; lincRNAs play crucial roles in various disorders and cellular models via multiple mechanisms. With rapid growth in the number of the identified lincRNAs and genetic variants, there is great demand for an investigation of SNVs in lincRNAs. Hence, in this article, we mainly summarize the significant role of SNVs within human lincRNA regions. Some pivotal variants may serve as risk factors for the development of various disorders, especially cancer. They may also act as important regulatory signatures involved in the modulation of lincRNAs in a tissue- or disorder-specific manner. An increasing number of researches indicate that lincRNA variants would potentially provide additional options for genetic testing and disease risk assessment in the personalized medicine era.

Lipid Trait Variants and the Risk of Non-Hodgkin Lymphoma Subtypes: A Mendelian Randomization Study

BACKGROUND

Lipid traits have been inconsistently linked to risk of non-Hodgkin lymphoma (NHL). We examined the association of genetically predicted lipid traits with risk of diffuse large B-cell lymphoma (DLBCL), chronic lymphocytic leukemia (CLL), follicular lymphoma (FL), and marginal zone lymphoma (MZL) using Mendelian randomization (MR) analysis.

METHODS

Genome-wide association study data from the InterLymph Consortium were available for 2,661 DLBCLs, 2,179 CLLs, 2,142 FLs, 824 MZLs, and 6,221 controls. SNPs associated (P < 5 × 10-8) with high-density lipoprotein (HDL, n = 164), low-density lipoprotein (LDL, n = 137), total cholesterol (TC, n = 161), and triglycerides (TG, n = 123) were used as instrumental variables (IV), explaining 14.6%, 27.7%, 16.8%, and 12.8% of phenotypic variation, respectively. Associations between each lipid trait and NHL subtype were calculated using the MR inverse variance-weighted method, estimating odds ratios (OR) per standard deviation and 95% confidence intervals (CI).

RESULTS

HDL was positively associated with DLBCL (OR = 1.14; 95% CI, 1.00-1.30) and MZL (OR = 1.09; 95% CI, 1.01-1.18), while TG was inversely associated with MZL risk (OR = 0.90; 95% CI, 0.83-0.99), all at nominal significance (P < 0.05). A positive trend was observed for HDL with FL risk (OR = 1.08; 95% CI, 0.99-1.19; P = 0.087). No associations were noteworthy after adjusting for multiple testing.

CONCLUSIONS

We did not find evidence of a clear or strong association of these lipid traits with the most common NHL subtypes. While these IVs have been previously linked to other cancers, our findings do not support any causal associations with these NHL subtypes.

IMPACT

Our results suggest that prior reported inverse associations of lipid traits are not likely to be causal and could represent reverse causality or confounding.

Linear and circular PVT1 in hematological malignancies and immune response: two faces of the same coin

Non coding RNAs (ncRNAs) have emerged as regulators of human carcinogenesis by affecting the expression of key tumor suppressor genes and oncogenes. They are divided into short and long ncRNAs, according to their length. Circular RNAs (circRNAs) are included in the second group and were recently discovered as being originated by back-splicing, joining either single or multiple exons, or exons with retained introns. The human Plasmacytoma Variant Translocation 1 (PVT1) gene maps on the long arm of chromosome 8 (8q24) and encodes for 52 ncRNAs variants, including 26 linear and 26 circular isoforms, and 6 microRNAs. PVT1 genomic locus is 54 Kb downstream to MYC and several interactions have been described among these two genes, including a feedback regulatory mechanism. MYC-independent functions of PVT1/circPVT1 have been also reported, especially in the regulation of immune responses. We here review and discuss the role of both PVT1 and circPVT1 in the hematopoietic system. No information is currently available concerning their transforming ability in hematopoietic cells. However, present literature supports their cooperation with a more aggressive and/or undifferentiated cell phenotype, thus contributing to cancer progression. PVT1/circPVT1 upregulation through genomic amplification or rearrangements and/or increased transcription, provides a proliferative advantage to malignant cells in acute myeloid leukemia, acute promyelocytic leukemia, Burkitt lymphoma, multiple myeloma (linear PVT1) and acute lymphoblastic leukemia (circPVT1). In addition, PVT1 and circPVT1 regulate immune responses: the overexpression of the linear form in myeloid derived suppressor cells induced immune tolerance in preclinical tumor models and circPVT1 showed immunosuppressive properties in myeloid and lymphoid cell subsets. Overall, these recent data on PVT1 and circPVT1 functions in hematological malignancies and immune responses reflect two faces of the same coin: involvement in cancer progression by promoting a more aggressive phenotype of malignant cells and negative regulation of the immune system as a novel potential therapy-resistance mechanism.

Integrating transcriptome-wide association study and copy number variation study identifies candidate genes and pathways for diffuse non-Hodgkin's lymphoma

BACKGROUND

The genetic basis of diffuse non-Hodgkin's lymphoma (DNHL) is largely unknown now. We conducted a large-scale transcriptome-wide association study (TWAS) of DNHL to identify novel candidates for DNHL.

METHODS

The GWAS summary data of DNHL was obtained from the UKBiobank, involving 685 cases and 451,579 controls. TWAS of DNHL was performed using tissue-specific gene expression weights generated from the Genotype-Tissue Expression (GTEx) data. The DNHLTWAS results were further validated by a previous published copy number alterations (CNA) study of DNHL. Gene ontology (GO) and pathway enrichment analysis of identified candidate genes were conducted by the DAVID 6.8.

RESULTS

We identified 214 genes with TWAS P value < 0.05 for DNHL, such as MRPL19 (P_TWAS = 0.0010), CRCP (P_TWAS = 0.0010) and SEMA3C (P_TWAS = 0.0010). After further comparing the 214 genes with copy number variations of DNHL patients, we found 1 overlapped gene, BCL10 (P_TWAS = 0.0100). We also detected 6 common GO terms shared between gene set enrichment analysis results of TWAS and CNAs, such as cytosol (P_TWAS = 0.0003, PCNAs = 4.99 × 10^-7) and membrane (P_TWAS = 0.0048, P_CNAs = 0.0046). The pathway enrichment analysis of TWAS and CNAs detected 3 common pathways, including HIF-1 signaling pathway (P_TWAS = 0.0195, P_CNAs = 1.96 × 10^-5), mTOR signaling pathway (P_TWAS = 0.0242, P_CNAs = 6.75 × 10^-5) and adipocytokine signaling pathway (P_TWAS = 0.0392, P_CNAs = 0.0103).

CONCLUSIONS

Our study identified multiple DNHL associated genes and pathways, providing novel useful information for the pathogenetic studies of DNHL.

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.