- Research
- Open access
- Published:
Genomic landscape and distinct molecular subtypes of primary testicular lymphoma
Journal of Translational Medicine volume 22, Article number: 414 (2024)
Abstract
Primary testicular lymphoma (PTL) is a rare lymphoma predominantly occurring in the elderly male population. It is characterized by a limited response to treatment and a heightened tendency towards relapse. Histologically, approximately 90% of PTL cases are classified as diffuse large B-cell lymphomas (DLBCL). Genetic features of PTL were delineated in a limited scope within several independent studies. Some of the articles which analyzed the genetic characterization of DLBCL have incorporated PTL samples, but these have been constrained by small sample sizes. In addition, there have been an absence of independent molecular typing studies of PTL. This report summarizes the common mutational features, copy number variations (CNVs) and molecular typing of PTL patients, based on whole-exome sequencing (WES) conducted on a cohort of 25 PTL patients. Among them, HLA, CDKN2A and MYD88 had a high mutation frequency. In addition, we found two core mutational characteristics in PTL including mutation in genes linked to genomic instability (TP53 and CDKN2A) and mutation in immune-related genes (HLA, MYD88, CD79B). We performed molecular typing of 25 PTL patients into C1 subtype with predominantly TP53 mutations and C2 subtype with predominantly HLA mutations. Notably, mutations in the TP53 gene predicted a poor outcome in most types of lymphomas. However, the C1 subtype, dominated by TP53 mutations, had a better prognosis compared to the C2 subtype in PTL. C2 subtype exhibited a worse prognosis, aligning with our finding that the mechanism of immune escape in PTL was primarily the deletions of HLA rather than PD-L1/PD-L2 alterations, a contrast to other DLBCLs. Moreover, we calculated the tumor mutation burden (TMB) and identified that TMB can predict prognosis and recurrence rate in PTL. Our study underscores the significance of molecular typing in PTL based on mutational characteristics, which plays a crucial role in prognostication and guiding therapeutic strategies for patients.
Introduction
PTL is an extremely rare and aggressive tumor, mainly occurring in older males. Histologically, about 90% of PTL cases are classified as diffuse large B-cell lymphomas (DLBCL). Immunohistochemical analyses and DNA microarray studies have indicated that 60% to 90% of PTL-DLBCL cases belong to the activated B-cell-like (ABC) subtype. Approximately 40% of PTL cases have the propensity to disseminate to extranodal organs, including central nervous system (CNS), the contralateral testis, skin, lung and other soft tissue [1,2,3,4]. The occurrence of PTL in what is considered an immune-privileged site contributes to its poor responsiveness to various treatment modalities and a relatively high rate of recurrence [2, 5, 6].
The negative prognosis and high recurrence rate of PTL are associated with immune escape mechanisms, including activation of immune checkpoints and loss of human leukocyte antigen (HLA). Programmed cell death 1 (PD-1) and cytotoxic T lymphocyte-associated Antigen-4 (CTLA-4) are currently the most investigated immune checkpoints. PD-L1, a ligand protein produced by tumor cells, can bind to PD-1 of T lymphocytes. The combination between PD-L1 and PD-1 induces deactivation, apoptosis, and depletion of T cells, which in turn inhibits the activation, proliferation and anti-tumor function of tumor antigen-specific CD8 + T cells. PD-L1/PD-L2 copy number alterations and additional translocations at relevant sites are common in PTL. These changes are thought to be responsible for immune escape in PTL [7,8,9,10]. CTLA-4 plays a crucial role in modulating T cell activation and tolerance [11, 12]. Another mechanism of tumor escape is the loss of HLA. Deficiencies in both HLA class I (HLA-A, B, C) and class II (HLA-DR, DQ, DP) can lead to inherently weak anti-tumor immune responses and undermine the efficacy of T cell-targeted immunotherapies. In addition, loss of HLA I and II expression significantly contributes to the phenomenon of immune privilege [13, 14]. In a study by Chapuy et al. rearrangements of PD-L1/PD-L2 in PTL were found to be prevalent in over 50% of PTL cases [15]. Contrastingly, Minderman's study revealed expression of PD-L1 or 9p24.1/PD-L1/2 CNA was detected in only a small number of samples. In addition, their team extracted PD-L1 messenger RNA (mRNA) expression data from publicly available databases for validation, and they found that there was no enhancement of PD-L1 expression in PTL compared to primary mediastinal B-cell lymphomas (PMBCLs) [16,17,18,19]. Instead, HLA deletions were common in their PTL cohort. Their article suggests that the mechanism of immune escape in PTL and CNS lymphoma is primarily HLA loss rather than the use of the PD1/PD-L1/2 checkpoint [20].
A limited number of prior studies focusing exclusively on PTL, as well as investigations of DLBCL that include PTL cases, have identified a range of genetic alterations in PTL patients [1, 15, 21,22,23,24,25]. High-frequency of somatic mutations in the MYD88 and CD79B genes (NF-κ-B pathway-associated genes) have been shown to be broadly present in PTL. Additionally, previous studies have also demonstrated the presence of BCL2 and BCL6 translocations in PTL [22, 26, 27].
Currently, both independent studies on PTL and those incorporating PTL within broader DLBCL research have been limited by small sample sizes. And there has been a lack of systematic molecular typing studies in PTL. Here, we performed whole exome sequencing (WES) analysis on 25 PTL patients, enabling us to perform mutation typing in this group. We anticipate that our molecular typing approach will shed light on the molecular heterogeneity of PTL, thereby facilitating more precise, subtype-specific treatments. Moreover, we calculated the tumor mutation burden (TMB) and established its potential as a prognostic and recurrence predictor in PTL.
Methods
Patients and samples
Formalin-fixed, paraffin-embedded specimens of 25 PTL patients (14 cancer samples and matched normal tissues, 11 cancer samples) were obtained and identified at the time of diagnosis. We statistically analyzed the clinical information of 25 PTLs, including age, histologic type, pathology, stage, B-symptoms, ECOG, IPI and recurrence (Table 1). In each case, the diagnosis of PTL was established employing appropriate diagnostic criteria (WHO classification) of lymphoid tumors with combinations of histologic, immunohistochemical, flow cytometric, and genetic evaluation (Additional file 12: Table S1). The study protocols described above comply with the ethical requirements of Peking University Third Hospital.
Whole exome sequencing and quality control of raw sequencing data
We used the traditional CTAB method to extract genomic DNA from tumor samples, followed by quality assessment of the extracted DNA using agarose gel electrophoresis and Nanodrop. The quality-checked genomic DNA was randomly broken into 180–280 bp fragments using a Covaris fragmentation instrument. Subsequently, DNA libraries were prepared and whole exome capture was conducted via hybridization with biotinylated single-strand DNA capture probe. Post-hybridization, the libraries underwent PCR amplification and quality control. Libraries passing quality checks were sequenced on the Novaseq platform. The read length of this this sequencing is 2*150 bp. The highest read depth is 155.24, the lowest read depth is 32.91, the average read depth is 72.98.
In the high-throughput sequencing data, a minor fraction of reads may contain low-quality bases or undetected bases (N). To ensure the quality of information analysis, raw reads were subjected to quality controlled using the Fastp, resulting in clean reads [28]. These clean reads were aligned to the reference genome (version hg19) using the software bwa to generate alignment data [29]. And subsequent analyses were performed based on these clean reads. The quality control process involved trimming low-quality bases from both the 5' and 3' ends of the sequences, discarding sequences with a high proportion of undetected bases (N), and removing adaptor sequences from the reads.
Somatic SNVs/Indels variant detection
Somatic mutation detection was performed using the Mutect2 module in GATK to obtain information on SNVs/Indels variants [30]. The results of Mutect2 were annotated through the simultaneous association with multiple databases (e.g. dbSNP, 1000g, ESP6500, HGMD, OMIM, etc.) using the software ANNOVAR [31]. To ensure high-quality somatic mutation data, further filtration was applied based on specific criteria: (a) Mut_ratio > 5%, Alt >  = 3 and Depth >  = 20; (b) remove the results with SNP effect of unknown and nonframeshift and keep the results with gene region of exonic and splicing; (c) keep the results with mutation frequency below 0.001 in 1000g, ESP6500 and ExAC databases; (d) remove mutation results annotated to interfering genes. Heat maps representing SNV and CNV mutation results were generated using the ComplexHeatmap R package [32]. Additional analyses, including the statistics of SNV mutation patterns, TMB statistics calculation, survival analysis, forest plotting, sample clustering analysis, and CNV mutation frequency statistics on samples were performed using the R language internal script (version 4.1.0).
CNV variant detection
For CNV detection in our samples, we employed CNVkit software, which determines copy number changes by comparing the depth distribution of sequencing reads to that of a reference genome [33]. The occurrence of CNVs in specific genes was subsequently analyzed for each sample using the CNVRanger package in R [34].
Pathway enrichment analysis
Pathway enrichment analysis was performed using Reactome to elucidate the functions of genes undergoing SNV mutations [35]. And the relationship network between pathways and genes was visualized using Cytoscape [36].
Gene expression data and analysis
To further our research, we accessed gene expression arrays for 22 cases of PTL and 232 nodal DLBCLs from the NCBI Gene Expression Omnibus (GEO) database (GSE10524, GSE10846, GSE61578) [1, 15, 37]. Gene expression levels were quantified from these microarray data using the robust multiarray averaging (RMA) algorithm. These expression levels were then transformed using a log2 scale for enhanced analytical clarity.
Result
Genetic characteristics and signaling pathways of 25 PTL patients
Based on the SNV and CNV profiles of 25 PTL patients, we can find that both HLA-A and HLA-C exhibit an overall mutation rate of 68% and the total mutation rate of HLA is as high as 84% (Fig. 1A). Although previous studies have reported CNV deletion in HLA-C. However, it is noteworthy that we also identified somatic mutations in HLA-C in PTL, which had not been reported in previous studies. 4 out of 25 patients had mutations in HLA-C, three of which were stop gain alterations. All three stop-gains are on the second exon of the transcript. The change of its cDNA from C to T at position 232 caused a change in the amino acid of the protein sequence at position 78 from glutamine (Q) to stop-gains. One of the nonsense mutations was in exon 4 of the transcript. Its cDNA was changed from A to G at position 749, causing a change in amino acid from glutamine (Q) to arginine (R) at position 250 of the protein sequence (Fig. 1 and Additional file 11: Table S4). These mutations lead to changes in the encoded proteins and thus loss of function. Moreover, we discovered that somatic mutations in HLA-C, as well as somatic mutations in ASH1L, contributed to worse prognosis for PTL patients (Additional files 2, 3: Figs. S2, S3). We also found that the CDKN2A gene had the third-highest mutation rate at 52% in 25 PTL patients (Fig. 1A and Additional file 1 Fig. S1A). CNV deletions in CDKN2A primarily manifested as homozygous deletions in seven out of the 25 PTL patients, corroborating previous findings by Chapuy et al. [15]. Analysis of the genes on chromosome 9 reveals a very large proportion of deletions in CDKN2A and CDKN2B (Additional file 7: Fig. S7A). Compared with nodal DLBCLs, the expression of CDKN2A on PTL patients is significantly low (Additional file 11: Fig. S11A, Wilcoxon test, P = 5.8e–06). In addition to the deletion of CDKN2A in PTL patients, which has been reported to be associated with genomic instability, we also identified deletions and mutations in TP53, another gene associated with genomic instability in PTL patients. We found that 4 patients with PTL had abnormal CNV profiles in TP53, with three presenting heterozygous deletions and one demonstrating both heterozygous deletions and mutations (Fig. 1A).
Then we performed TMB in 25 PTL patients and subsequently found out that high TMB predicted better progression-free survival (PFS) and overall survival (OS) time (Fig. 1B and Additional file 4: Fig. S4A, PFS, P = 0.032; OS, P = 0.041). Additionally, patients with high TMB had a low relapse rate (Fig. 1B and Additional files 5, 6: Figs. S5, S6, Chi-square test, P = 0.01).
CNVs profiles in PTL
As for CNV profile, we can find a substantial number of deletions on chromosome 6 where HLA locates (Fig. 1C). In addition to the high-frequency deletions of HLA class I genes, including HLA-A, HLA-B, HLA-C, HLA-E, HLA-F and HLA-G, we also found high-frequency deletions of CNV in HLA class II genes (Fig. 1D and Additional file 7: Fig. S7A). Notably, the expression of HLA-A and HLA-C in PTL is lower than that in nodal DLBCLs. However, compared with the expression of PD-L1 in nodal DLBCLs, the expression of PD-L1 did not decrease in PTL patients. And the expression of PD-L2 in PTL patients was higher than that in nodal DLBCLs (Fig. 1D, Wilcoxon test). Compared with HLA genes expression on nodal DLBCLs, the expression of HLA genes was decreased in PTL patients (Additional file 8: Fig. S8A). Furthermore, we identified genes displaying deletions such as P4HTM and WDR6, which were associated with poor prognosis for PTL patients (Fig. 1C and Additional file 9 Fig. S9, log-rank test). Additionally, amplified genes in PTL patients such as TLE1 and ADAMTSL4-AS1 also contributed to the unfavorable prognosis of PTL patients (Additional file 10: Fig. S10, log-rank test). Besides, pathway analysis of the mutant genes in these 25 PTL samples revealed significant enrichment in pathways related to the immune system, antigen processing cross-presentation, and signal transduction (Figs. 1E–F).
Molecular subtyping of 25 PTL patients
We performed a cluster analysis of the 25 PTL patients based on the results of SNV, CNV and TMB profiles, which divided the patients into two distinct subtypes (Fig. 2A). C1 is a genomic-instability subtype characterized primarily by the mutation of TP53. C2 is an immune-escape subtype characterized by the mutation of immune-related genes such as HLA-A, HLA-C, PIM1 and so on, with HLA genes having the highest frequency. Moreover, the total counts of gene with CNV in C1 is higher than that in C2, signifying a significantly greater degree of genomic instability in C1 (Fig. 2B). Interestingly, despite TP53 mutations typically predicting poorer prognosis in lymphomas, C1, predominantly featuring TP53 mutations in PTL, exhibits improved PFS and OS (Fig. 2C, OS, P = 0.01, PFS, P = 0.045). By comparing the clinical information between C1 and C2, we found that there was only a significant difference in stage between the two groups. Patients in C2 were basically at stage III-IV when they were diagnosed (Table 2).
Discussion
Several studies in recent years have delved into the mechanisms underlying immune escape in PTL and to elucidate the factors contributing to its poor response to treatment and prognosis. In addition, there have been a few studies on the mutation spectrum of PTL genes, but there are currently no studies that characterize the mutation profiles of PTL patients in order to perform typing and prognostic analyses. Such a clinically applicable typing approach we have created could stratify PTL, thereby guiding the precise treatment of patients. Here, we performed the WES on 25 PTL patients to delineate the mutation characteristics, chromosomal rearrangement, and CNVs in PTL, ultimately classifying these patients into two molecular subtypes.
By analyzing the mutation spectrum of PTL, we can find that HLA-A, HLA-C, CDKN2A, MYD88, PIM1, SETD1B, ASH1L, IRF4, TP53, DTX1, SMG7, CD79B, ADGRL2 and CREBBP exhibit high mutation rates in PTL. Among them, HLA-A, HLA-C, PIM1 and MYD88 are immune-related genes. In addition to mutations in immune-related genes, we identified a high frequency of mutations in genes affecting genomic stability, such as CDKN2A, TP53, SETD1B and SMG7. Genomic instability is often caused by DNA mutations and aberrant epigenetic modifications (DNA methylation and histone modifications) [38,39,40,41,42,43,44]. Many studies have demonstrated that genetic alterations in TP53 and CDKN2A can induce or tolerate genomic instability [45,46,47,48,49,50,51,52]. In early deep sequencing analyses, somatic mutations in TP53 were seen in only 20% of DLBCL. Somatic mutations in TP53 were even rarer in PTL [15, 21]. However, in our cohort, we found TP53 somatic mutations in 24% of PTLs, including four nonsynonymous as well as one stop-gain. In our subsequent molecular typing of PTL, somatic mutations in TP53 were a major feature present in one of the subtypes. We also found a high frequency of homozygous deletions in the CDKN2A gene. This was a more characteristic genetic manifestation of PTL compared to DLBCL [15].
Currently there is no consensus in the study of PTL regarding its main immune escape mechanism. Some scholars consider the gain of PD-L1/PD-L2 in tumors as the main immune escape mechanism of PTL. Alterations in PD-L1/PD-L2 (9p24.1) copy number were common (> 50%) in PTL in their studies [15, 53, 54]. Other scholars have suggested that HLA deletion is the main immune escape mechanism in PTL. They found a high frequency of HLA I and II deletions in PTL but little expression of PD-L1, resulting in no 9p24.1/PD-L1/2 CNA [1, 55]. Our analysis of PTL by whole exome sequencing suggests that the main tumor immune escape mechanism in PTL is HLA deletion on rather than PD-L1/PD-L2 alteration. We found a large number of deletions at the chromosome 6 locus, with very frequent expression deletions in HLA I and HLA II. Compared to CNV in HLA I, a higher percentage of CNV deletions are found in the HLA-DR genes (HLA-DRA, HLA-DRB5, HLA-DRB6, and HLA-DRB1, etc.) and the HLA-DQ genes (HLA-DQA1, HLA-DQB1, and HLA-DQHLA-DQB2, etc.). CNV deletion of the HLA-DR genes can severely disrupts the immune response to PTL. Then we compared the CNVs of PTL patients with those of nodal DLBCLs patients for validation. Our analysis revealed that the CNVs of HLA-A and HLA-C were significantly lower in PTL patients compared with those of nodal DLBCL patients. Conversely, the CNVs of PD-L2 were significantly higher in PTL patients compared to nodal DLBCL patients. However, there was no significant difference in CNVs of PD-L1 between the two groups. These findings suggest that copy deletion at HLA loci is the main mechanism of immune escape in PTL tumors. In addition, we found somatic mutations in the HLA genes, including HLA-A and HLA-C. Among them, the mutation in HLA-C was a somatic mutation that had not been reported in previous articles. Additionally, by calculating the TMB of the PTL patients in our cohort, we found that PTL patients with higher TMB were less likely to relapse and had longer PFS. PTL patients with high TMB may indirectly reflect the tumor's ability to produce neoantigens, especially considering that most of the genes with high mutation rates in our cohort are immune-related genes. Consequently, such patients may benefit from immunotherapy, making TMB a potential indicator of treatment efficacy.
We successfully categorized PTL patients into C1 and C2 based on their mutation characteristics by cluster analysis. C1 is the subtype characterized by TP53 mutations, while C2 is the subtype characterized by mutations in immune-related genes such as HLA-A, HLA-C and PIM1. A comparison of the two subtypes revealed that C2 had worse PFS and OS compared to C1. It's worth noting that, to date, there have been no studies that have typed patients with PTL. Our mutation typing of PTL patients not only facilitates the prediction of prognosis in patients with primary treatment of PTL but also lays the groundwork for future precision medicine in PTL.
In conclusion, we found two types of core genetic variants in PTL patients including mutations in immune-related genes (HLA, MYD88, CD79B and PIM1) and genomic instability-related genes (TP53 and CDKN2A). Among them, somatic mutations in HLA-C and TP53 have hardly been identified in previous genetic characterization studies of PTL. This group has a better prognosis compared to C2. In addition, we discovered that TMB can predict prognosis and recurrence rate in PTL. Finally, PTL patients can be divided into two groups based on molecular subtyping. The first group, C1, is characterized by mutations in genes that maintain genomic stability, such as TP53. The second group, C2, is characterized by mutations in immune-related genes, including those in the HLA complex. Notably, patients within the C2 subtype exhibit a poorer prognosis, indicating a significant correlation between the prognosis of PTL patients and their immune functionality. This molecular stratification presented in our study holds potential for enhancing the precision in diagnosis, prognostic assessment, tumor staging, treatment guidance, monitoring of recurrence, and the development of targeted therapies for PTL.
Availability of date and materials
The authors thank GSE10524, GSE10846, GSE61578 dataset from GEO database which are all available online. And the authors confirm that the data supporting the fundings of this study are available within our Additional file table.
References
Booman M, et al. Genomic alterations and gene expression in primary diffuse large B-cell lymphomas of immune-privileged sites: the importance of apoptosis and immunomodulatory pathways. J Pathol. 2008;216(2):209–17.
Cheah CY, Wirth A, Seymour JF. Primary testicular lymphoma. Blood. 2014;123(4):486–93.
Shahab N, Doll DC. Testicular lymphoma. Semin Oncol. 1999;26(3):259–69.
Crellin AM, et al. Non-Hodgkin’s lymphoma of the testis. Radiother Oncol. 1993;27(2):99–106.
Vitolo U, Ferreri AJ, Zucca E. Primary testicular lymphoma. Crit Rev Oncol Hematol. 2008;65(2):183–9.
Horne MJ, Adeniran AJ. Primary diffuse large B-cell lymphoma of the testis. Arch Pathol Lab Med. 2011;135(10):1363–7.
Boussiotis VA. Molecular and biochemical aspects of the PD-1 checkpoint pathway. N Engl J Med. 2016;375(18):1767–78.
Gordon SR, et al. PD-1 expression by tumour-associated macrophages inhibits phagocytosis and tumour immunity. Nature. 2017;545(7655):495–9.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64.
Zou W, Wolchok JD, Chen L. PD-L1 (B7–H1) and PD-1 pathway blockade for cancer therapy: mechanisms, response biomarkers, and combinations. Sci Transl Med. 2016;8(328):328rv4.
Leach DR, Krummel MF, Allison JP. Enhancement of antitumor immunity by CTLA-4 blockade. Science. 1996;271(5256):1734–6.
Schwartz RH. Costimulation of T lymphocytes: the role of CD28, CTLA-4, and B7/BB1 in interleukin-2 production and immunotherapy. Cell. 1992;71(7):1065–8.
Ferguson TA, Griffith TS. A vision of cell death: insights into immune privilege. Immunol Rev. 1997;156:167–84.
Riemersma SA, et al. Extensive genetic alterations of the HLA region, including homozygous deletions of HLA class II genes in B-cell lymphomas arising in immune-privileged sites. Blood. 2000;96(10):3569–77.
Chapuy B, et al. Targetable genetic features of primary testicular and primary central nervous system lymphomas. Blood. 2016;127(7):869–81.
Sanjari Nia AH, et al. Non-coding RNAs: a new frontier in benzene-mediated toxicity. Toxicology. 2023;500: 153660.
Pordel S, et al. The role of microRNA-185 in the pathogenesis of human diseases: a focus on cancer. Pathol Res Pract. 2023;249: 154729.
Fattahi M, et al. microRNA-382 as a tumor suppressor? Roles in tumorigenesis and clinical significance. Int J Biol Macromol. 2023;250: 125863.
Huang C, et al. Non-coding RNAs/DNMT3B axis in human cancers: from pathogenesis to clinical significance. J Transl Med. 2023;21(1):621.
Minderman M, et al. Immune evasion in primary testicular and central nervous system lymphomas HLA loss rather than 9p241/PD-L.1/PD-L2 alterations. Blood. 2021;138(13):1194–7.
Monti S, et al. Integrative analysis reveals an outcome-associated and targetable pattern of p53 and cell cycle deregulation in diffuse large B cell lymphoma. Cancer Cell. 2012;22(3):359–72.
Ngo VN, et al. Oncogenically active MYD88 mutations in human lymphoma. Nature. 2011;470(7332):115–9.
Lenz G, et al. Oncogenic CARD11 mutations in human diffuse large B cell lymphoma. Science. 2008;319(5870):1676–9.
Davis RE, et al. Chronic active B-cell-receptor signalling in diffuse large B-cell lymphoma. Nature. 2010;463(7277):88–92.
Menter T, et al. Phenotype profiling of primary testicular diffuse large B-cell lymphomas. Hematol Oncol. 2014;32(2):72–81.
Kraan W, et al. High prevalence of oncogenic MYD88 and CD79B mutations in primary testicular diffuse large B-cell lymphoma. Leukemia. 2014;28(3):719–20.
Oishi N, et al. High prevalence of the MYD88 mutation in testicular lymphoma: Immunohistochemical and genetic analyses. Pathol Int. 2015;65(10):528–35.
Chen S, et al. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.
gatk, . Homepage. https://gatk.broadinstitute.org/hc/en-us/articles/2-Mutect.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16): e164.
Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–9.
Talevich E, et al. CNVkit: genome-wide copy number detection and visualization from targeted DNA sequencing. PLoS Comput Biol. 2016;12(4): e1004873.
da Silva V, et al. CNVRanger: association analysis of CNVs with gene expression and quantitative phenotypes. Bioinformatics. 2020;36(3):972–3.
Fabregat A, et al. Reactome diagram viewer: data structures and strategies to boost performance. Bioinformatics. 2018;34(7):1208–14.
Cytoscape. Homepage. https://cytoscape.org/.
Lenz G, et al. Stromal gene signatures in large-B-cell lymphomas. N Engl J Med. 2008;359(22):2313–23.
Rogakou EP, et al. DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139. J Biol Chem. 1998;273(10):5858–68.
Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100(1):57–70.
Burma S, et al. ATM phosphorylates histone H2AX in response to DNA double-strand breaks. J Biol Chem. 2001;276(45):42462–7.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Chapuy B, et al. Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes. Nat Med. 2018;24(5):679–90.
Vij M, et al. P16 immunohistochemistry is a sensitive and specific surrogate marker for CDKN2A homozygous deletion in gliomas. Acta Neuropathol Commun. 2023;11(1):73.
Clurman BE, Groudine M. The CDKN2A tumor-suppressor locus–a tale of two proteins. N Engl J Med. 1998;338(13):910–2.
Annaratone L, et al. The multifaceted nature of tumor microenvironment in breast carcinomas. Pathobiology. 2020;87(2):125–42.
Fridman JS, Lowe SW. Control of apoptosis by p53. Oncogene. 2003;22(56):9030–40.
Kastan MB, Canman CE, Leonard CJ. P53, cell cycle control and apoptosis: Implications for cancer. Cancer Metastasis Rev. 1995;14(1):3–15.
Vaseva AV, Moll UM. The mitochondrial p53 pathway. Biochim Biophys Acta. 2009;1787(5):414–20.
Speidel D. Transcription-independent p53 apoptosis: an alternative route to death. Trends Cell Biol. 2010;20(1):14–24.
Hassin O, Oren M. Drugging p53 in cancer: one protein, many targets. Nat Rev Drug Discov. 2023;22(2):127–44.
Mortaza Hosainzadegan AE. Rovshan Khalilov, Aygun Nasibova, Amir Hasanzadeh, Parviz Vahedi, Hasan Hosain zadegan, are microbial infections and some antibiotics causes cancer? Adv Biol Earth Sci. 2020;5(1):58.
Rovshan Khalilov, S.A. Mechanisms of insulin action and insulin resistance. Advances in Biology & Earth Sciences, 2023. 8 2
Kiyasu J, et al. Expression of programmed cell death ligand 1 is associated with poor overall survival in patients with diffuse large B-cell lymphoma. Blood. 2015;126(19):2193–201.
Pollari M, et al. PD-L1(+) tumor-associated macrophages and PD-1(+) tumor-infiltrating lymphocytes predict survival in primary testicular lymphoma. Haematologica. 2018;103(11):1908–14.
Chowell D, et al. Patient HLA class I genotype influences cancer response to checkpoint blockade immunotherapy. Science. 2018;359(6375):582–7.
Acknowledgements
The authors thank GSE10524, GSE10846, GSE61578 dataset from GEO database which are all available online.
Funding
Financial support for this research was provided by National Natural Science Foundation of China (81800195), Key Clinical Projects of Peking University Third Hospital (BYSYZD2019026 and BYSYZD2023014), Beijing Xisike Clinical Oncology Research Foundation (Y-NCJH202201-0049), the special fund of the National Clinical Key Specialty Construction Program, P. R. China (2023).
Author information
Authors and Affiliations
Contributions
HMJ conceived the project. WLZ, PY, YRY and SZL analyzed the data. YDX, CLW, JW, CLL, HL, SSL, WH and HMJ contributed towards the interpretation of the data. All authors wrote and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Study protocols described above comply with the ethical requirements of Peking University Third Hospital.
Consent for publication
All authors approved the final manuscript and the submission to this journal.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Mutation spectrum of 25 PTL patients, showing the mutation frequency of each gene (right) and clinical data of each sample (bottom).
Additional file 2: Figure S2.
Effect of mutations in HLA-C on prognosis in PTL patients. A The Kaplan-Meier curves for OS of HLA-C mutation (log-rank test, P=2.9e−03)
Additional file 3: Figure S3.
Effect of mutations in ASH1L on prognosis in PTL patients. A The Kaplan-Meier curves for PFS of ASH1L mutation (log-rank test, P=3e−02)
Additional file 4: Figure S4.
The Kaplan-Meier curves for OS of 25 PTL patients in TMB-low and TMB-high groups (log-rank test, P=0.041)
Additional file 5: Figure S5.
Mutational signature analysis of 25 PTL patients. A–B Proportion of six mutation classes in the exome of 25 PTL patients. C Distribution of the six mutation classes in 25 PTL patients
Additional file 6: Figure S6.
Distribution of five mutation types in 25 PTL patients
Additional file 7: Figure S7.
The percentage of CNV and SNV profiles on chromosome. A The percentage of CNV and SNV profiles for multiple genes on chromosome 6. B The percentage of CNV and SNV profiles for multiple genes on chromosome 9
Additional file 8: Figure S8.
Different expression of each HLA types between PTL patients and nodal DLBCLs (22 PTL patients and 232 nodal DLBCLs patients from GSE10524,GSE10846,GSE61578)
Additional file 9: Figure S9.
Effect of deletion-based genes on prognosis in PTL patients. A Forest plots of deleted genes. The Kaplan-Meier curves for PFS of the CNV in P4HTM (log-rank test, P=7.6e−03) and WDR6 (log-rank test, P=7.6e−03).
Additional file 10: Figure S10.
Effect of amplification-based genes on prognosis in patients with PTL. A Forest plots of amplified genes. B The Kaplan-Meier curves for OS of the CNV in TLE1 (log-rank test, P=9e−04) and ADAMTSL4-AS1 (log-rank test, P=3e−03)
Additional file 11: Figure S11.
Different gene expression between PTL patients and nodal DLBCLs. (A) The expression of CDKN2A between PTL patients and nodal DLBCLs.
Additional file 12: Table S1.
Clinical information of PTL patients in this study. Table S2. Sequencing fastq files quality before and after trim. Table S3. Map statistics for aligning to the genome. Table S4. SNV results after filter. Table S5. CNV results after CNVRanger processing. Table S6. Reactome enrichment results using the gene with SNV and CNV.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhang, W., Yang, P., Yang, Y. et al. Genomic landscape and distinct molecular subtypes of primary testicular lymphoma. J Transl Med 22, 414 (2024). https://doi.org/10.1186/s12967-024-05140-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-024-05140-8