A reliable transcriptomic risk-score applicable to formalin-fixed paraffin-embedded biopsies improves outcome prediction in localized prostate cancer

BACKGROUND

Clinical manifestation of prostate cancer (PCa) is highly variable. Aggressive tumors require radical treatment while clinically non-significant ones may be suitable for active surveillance. We previously developed the prognostic ProstaTrend RNA signature based on transcriptome-wide microarray and RNA-sequencing (RNA-Seq) analyses, primarily of prostatectomy specimens. An RNA-Seq study of formalin-fixed paraffin-embedded (FFPE) tumor biopsies has now allowed us to use this test as a basis for the development of a novel test that is applicable to FFPE biopsies as a tool for early routine PCa diagnostics.

METHODS

All patients of the FFPE biopsy cohort were treated by radical prostatectomy and median follow-up for biochemical recurrence (BCR) was 9 years. Based on the transcriptome data of 176 FFPE biopsies, we filtered ProstaTrend for genes susceptible to FFPE-associated degradation via regression analysis. ProstaTrend was additionally restricted to genes with concordant prognostic effects in the RNA-Seq TCGA prostate adenocarcinoma (PRAD) cohort to ensure robust and broad applicability. The prognostic relevance of the refined Transcriptomic Risk Score (TRS) was analyzed by Kaplan-Meier curves and Cox-regression models in our FFPE-biopsy cohort and 9 other public datasets from PCa patients with BCR as primary endpoint. In addition, we developed a prostate single-cell atlas of 41 PCa patients from 5 publicly available studies to analyze gene expression of ProstaTrend genes in different cell compartments.

RESULTS

Validation of the TRS using the original ProstaTrend signature in the cohort of FFPE biopsies revealed a relevant impact of FFPE-associated degradation on gene expression and consequently no significant association with prognosis (Cox-regression, p-value > 0.05) in FFPE tissue. However, the TRS based on the new version of the ProstaTrend-ffpe signature, which included 204 genes (of originally 1396 genes), was significantly associated with BCR in the FFPE biopsy cohort (Cox-regression p-value < 0.001) and retained prognostic relevance when adjusted for Gleason Grade Groups. We confirmed a significant association with BCR in 9 independent cohorts including 1109 patients. Comparison of the prognostic performance of the TRS with 17 other prognostically relevant PCa panels revealed that ProstaTrend-ffpe was among the best-ranked panels. We generated a PCa cell atlas to associate ProstaTrend genes with cell lineages or cell types. Tumor-specific luminal cells have a significantly higher TRS than normal luminal cells in all analyzed datasets. In addition, TRS of epithelial and luminal cells was correlated with increased Gleason score in 3 studies.

CONCLUSIONS

We developed a prognostic gene-expression signature for PCa that can be applied to FFPE biopsies and may be suitable to support clinical decision-making.

Characterization and evaluation of gene fusions as a measure of genetic instability and disease prognosis in prostate cancer

BACKGROUND

Prostate cancer (PCa) is one of the most prevalent cancers worldwide. The clinical manifestations and molecular characteristics of PCa are highly variable. Aggressive types require radical treatment, whereas indolent ones may be suitable for active surveillance or organ-preserving focal therapies. Patient stratification by clinical or pathological risk categories still lacks sufficient precision. Incorporating molecular biomarkers, such as transcriptome-wide expression signatures, improves patient stratification but so far excludes chromosomal rearrangements. In this study, we investigated gene fusions in PCa, characterized potential novel candidates, and explored their role as prognostic markers for PCa progression.

METHODS

We analyzed 630 patients in four cohorts with varying traits regarding sequencing protocols, sample conservation, and PCa risk group. The datasets included transcriptome-wide expression and matched clinical follow-up data to detect and characterize gene fusions in PCa. With the fusion calling software Arriba, we computationally predicted gene fusions. Following detection, we annotated the gene fusions using published databases for gene fusions in cancer. To relate the occurrence of gene fusions to Gleason Grading Groups and disease prognosis, we performed survival analyses using the Kaplan-Meier estimator, log-rank test, and Cox regression.

RESULTS

Our analyses identified two potential novel gene fusions, MBTTPS2,L0XNC01::SMS and AMACR::AMACR. These fusions were detected in all four studied cohorts, providing compelling evidence for the validity of these fusions and their relevance in PCa. We also found that the number of gene fusions detected in a patient sample was significantly associated with the time to biochemical recurrence in two of the four cohorts (log-rank test, p-value < 0.05 for both cohorts). This was also confirmed after adjusting the prognostic model for Gleason Grading Groups (Cox regression, p-values < 0.05).

CONCLUSIONS

Our gene fusion characterization workflow revealed two potential novel fusions specific for PCa. We found evidence that the number of gene fusions was associated with the prognosis of PCa. However, as the quantitative correlations were only moderately strong, further validation and assessment of clinical value is required before potential application.

Transcriptional states of CAR-T infusion relate to neurotoxicity – lessons from high-resolution single-cell SOM expression portraying

Anti-CD19 CAR-T cell immunotherapy is a hopeful treatment option for patients with B cell lymphomas, however it copes with partly severe adverse effects like neurotoxicity. Single-cell resolved molecular data sets in combination with clinical parametrization allow for comprehensive characterization of cellular subpopulations, their transcriptomic states, and their relation to the adverse effects. We here present a re-analysis of single-cell RNA sequencing data of 24 patients comprising more than 130,000 cells with focus on cellular states and their association to immune cell related neurotoxicity. For this, we developed a single-cell data portraying workflow to disentangle the transcriptional state space with single-cell resolution and its analysis in terms of modularly-composed cellular programs. We demonstrated capabilities of single-cell data portraying to disentangle transcriptional states using intuitive visualization, functional mining, molecular cell stratification, and variability analyses. Our analysis revealed that the T cell composition of the patient’s infusion product as well as the spectrum of their transcriptional states of cells derived from patients with low ICANS grade do not markedly differ from those of cells from high ICANS patients, while the relative abundancies, particularly that of cycling cells, of LAG3-mediated exhaustion and of CAR positive cells, vary. Our study provides molecular details of the transcriptomic landscape with possible impact to overcome neurotoxicity.

Classifying Germinal Center Derived Lymphomas-Navigate a Complex Transcriptional Landscape

Classification of lymphoid neoplasms is based mainly on histologic, immunologic, and (rarer) genetic features. It has been supplemented by gene expression profiling (GEP) in the last decade. Despite the considerable success, particularly in associating lymphoma subtypes with specific transcriptional programs and classifier signatures of up- or downregulated genes, competing molecular classifiers were often proposed in the literature by different groups for the same classification tasks to distinguish, e.g., BL versus DLBCL or different DLBCL subtypes. Moreover, rarer sub-entities such as MYC and BCL2 "double hit lymphomas" (DHL), IRF4-rearranged large cell lymphoma (IRF4-LCL), and Burkitt-like lymphomas with 11q aberration pattern (mnBLL-11q) attracted interest while their relatedness regarding the major classes is still unclear in many respects. We explored the transcriptional landscape of 873 lymphomas referring to a wide spectrum of subtypes by applying self-organizing maps (SOM) machine learning. The landscape reveals a continuum of transcriptional states activated in the different subtypes without clear-cut borderlines between them and preventing their unambiguous classification. These states show striking parallels with single cell gene expression of the active germinal center (GC), which is characterized by the cyclic progression of B-cells. The expression patterns along the GC trajectory are discriminative for distinguishing different lymphoma subtypes. We show that the rare subtypes take intermediate positions between BL, DLBCL, and FL as considered by the 5th edition of the WHO classification of haemato-lymphoid tumors in 2022. Classifier gene signatures extracted from these states as modules of coregulated genes are competitive with literature classifiers. They provide functional-defined classifiers with the option of consenting redundant classifiers from the literature. We discuss alternative classification schemes of different granularity and functional impact as possible avenues toward personalization and improved diagnostics of GC-derived lymphomas.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

Lupus erythematosus: correlation of clinical and histological findings and proposal for a modified disease classification

BACKGROUND

Classification of lupus erythematosus (LE) is conflicting as it is carried out from different starting points. Whereas dermatological classifications categorize LE morphologically based on specific cutaneous lesions, rheumatologic classifications are based on symptomatic aspects. Indeed, LE is a systemic autoimmune disease with variable acuity and organ involvement. All cutaneous disease patterns may occur in both limited-cutaneous and systemic LE.

PATIENTS AND METHODS

76 LE-patients with complete clinical data, clinical photographs and biopsy of cutaneous manifestations as well as paraclinical findings were retrospectively analyzed. Based on a published two-dimensional classification system that considers disease-specific skin manifestations and final disease diagnosis separately, patients' diagnoses were revised and compared with those in medical records. In addition, the extent to which patients could be clustered by diagnosis based on their LE-specific skin manifestations, corresponding histopathological changes, and paraclinical data was investigated.

RESULTS

After re-evaluation, the proportion of patients with limited-cutaneous LE decreased from 82% previously to 24%. More than two-thirds of patients indeed showed intermediate or systemic LE. Disease-specific skin manifestations, histologic characteristics and paraclinical data did not cluster with final diagnoses.

CONCLUSIONS

First, the work underlines the systemic character of the disease. Second, a two-dimensional approach can help overcome classification difficulties in LE, as skin-morphologic and symptomatic aspects can be considered separately.

Distinct Molecular Mechanisms of Altered HLA Class II Expression in Malignant Melanoma

BACKGROUND

The human leukocyte antigen (HLA) class II molecules are constitutively expressed in some melanoma, but the underlying molecular mechanisms have not yet been characterized.

METHODS

The expression of HLA class II antigen processing machinery (APM) components was determined in melanoma samples by qPCR, Western blot, flow cytometry and immunohistochemistry. Immunohistochemical and TCGA datasets were used for correlation of HLA class II expression to tumor grading, T-cell infiltration and patients' survival.

RESULTS

The heterogeneous HLA class II expression in melanoma samples allowed us to characterize four distinct phenotypes. Phenotype I totally lacks constitutive HLA class II surface expression, which is inducible by interferon-gamma (IFN-γ); phenotype II expresses low basal surface HLA class II that is further upregulated by IFN-γ; phenotype III lacks constitutive and IFN-γ controlled HLA class II expression, but could be induced by epigenetic drugs; and in phenotype IV, lack of HLA class II expression is not recovered by any drug tested. High levels of HLA class II APM component expression were associated with an increased intra-tumoral CD4+ T-cell density and increased patients' survival.

CONCLUSIONS

The heterogeneous basal expression of HLA class II antigens and/or APM components in melanoma cells is caused by distinct molecular mechanisms and has clinical relevance.

High-Resolution Cartography of the Transcriptome and Methylome Landscapes of Diffuse Gliomas

Molecular mechanisms of lower-grade (II-III) diffuse gliomas (LGG) are still poorly understood, mainly because of their heterogeneity. They split into astrocytoma- (IDH-A) and oligodendroglioma-like (IDH-O) tumors both carrying mutations(s) at the isocitrate dehydrogenase (IDH) gene and into IDH wild type (IDH-wt) gliomas of glioblastoma resemblance. We generated detailed maps of the transcriptomes and DNA methylomes, revealing that cell functions divided into three major archetypic hallmarks: (i) increased proliferation in IDH-wt and, to a lesser degree, IDH-O; (ii) increased inflammation in IDH-A and IDH-wt; and (iii) the loss of synaptic transmission in all subtypes. Immunogenic properties of IDH-A are diverse, partly resembling signatures observed in grade IV mesenchymal glioblastomas or in grade I pilocytic astrocytomas. We analyzed details of coregulation between gene expression and DNA methylation and of the immunogenic micro-environment presumably driving tumor development and treatment resistance. Our transcriptome and methylome maps support personalized, case-by-case views to decipher the heterogeneity of glioma states in terms of data portraits. Thereby, molecular cartography provides a graphical coordinate system that links gene-level information with glioma subtypes, their phenotypes, and clinical context.

ProstaTrend-A Multivariable Prognostic RNA Expression Score for Aggressive Prostate Cancer

BACKGROUND

Prostate cancer (PCa) is the most prevalent solid cancer among men in Western Countries. The clinical behavior of localized PCa is highly variable. Some cancers are aggressive leading to death, while others can even be monitored safely. Hence, there is a high clinical need for precise biomarkers for identification of aggressive disease in addition to established clinical parameters.

OBJECTIVE

To develop an RNA expression-based score for the prediction of PCa prognosis that facilitates clinical decision making.

DESIGN, SETTING, AND PARTICIPANTS

We assessed 233 tissue specimens of PCa patients with long-term follow-up data from fresh-frozen radical prostatectomies (RPs), from formalin-fixed and paraffin-embedded RP specimens and biopsies by transcriptome-wide next-generation sequencing and customized expression microarrays.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

We applied Cox proportional hazard models to the cohorts from different platforms and specimen types. Evidence from these models was combined by fixed-effect meta-analysis to identify genes predictive of the time to death of disease (DoD). Genes were combined by a weighted median approach into a prognostic score called ProstaTrend and transferred for the prediction of biochemical recurrence (BCR) after RP in an independent cohort of The Cancer Genome Atlas (TCGA).

RESULTS AND LIMITATIONS

ProstaTrend comprising ∼1400 genes was significantly associated with DoD in the training cohort of PCa patients treated by RP (leave-one-out cross-validation, Cox regression: p=2e-09) and with BCR in the TCGA validation cohort (Cox regression: p=3e-06). The prognostic impact persisted after multivariable Cox regression analysis adjusting for Gleason grading group (GG) ≥3 and resection status (p=0.001; DoD, training cohort) and for GG≥3, pathological stage ≥T3, and resection state (p=0.037; BCR, validation cohort).

CONCLUSIONS

ProstaTrend is a transcriptome-based score that predicts DoD and BCR in cohorts of PCa patients treated with RP.

PATIENT SUMMARY

ProstaTrend provides molecular patient risk stratification after radical prostatectomy.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

The impact of SOCS1 mutations in diffuse large B-cell lymphoma

Mutations in SOCS1 are frequent in primary mediastinal B‐cell lymphoma and classical Hodgkin lymphoma. In the latter, SOCS1 mutations affect the length of the encoded protein (major mutations) and are associated with shorter patient survival. Two independent studies examined the prognostic impact of SOCS1 mutations in diffuse large B‐cell lymphoma (DLBCL) and showed differing results. This may be due to the small number of included patients, the heterogeneity of patients’ demographics and the distinct treatment schemes in these studies. To overcome the size limitations of these previous studies, we assessed SOCS1 mutations in the RICOVER‐60 cohort. The cohort uniformly consists of elderly patients (aged 61–80 years) treated with the CHOP‐14 scheme (cyclophosphamide, hydroxydaunorubicin, vincristine, prednisolone at 14‐day intervals) with or without an additional rituximab treatment. Patient outcomes were analysed with regard to overall SOCS1 mutation frequency, major and minor mutations and a novel impact‐based classifier – against the treatment modalities. Patients harbouring putative pathogenic SOCS1 mutations showed significant reduced overall survival within the CHOP plus rituximab group. Hence, putative pathogenic SOCS1 mutations seem to efface the beneficial effect of the therapeutic CD20 antibody. Comparing published data of whole exome and transcriptome sequencing of a large DLBCL cohort confirmed that predicted deleterious SOCS1 mutations forecast pre‐eminent survival in early onset DLBCL.

A modular transcriptome map of mature B cell lymphomas

BACKGROUND

Germinal center-derived B cell lymphomas are tumors of the lymphoid tissues representing one of the most heterogeneous malignancies. Here we characterize the variety of transcriptomic phenotypes of this disease based on 873 biopsy specimens collected in the German Cancer Aid MMML (Molecular Mechanisms in Malignant Lymphoma) consortium. They include diffuse large B cell lymphoma (DLBCL), follicular lymphoma (FL), Burkitt's lymphoma, mixed FL/DLBCL lymphomas, primary mediastinal large B cell lymphoma, multiple myeloma, IRF4-rearranged large cell lymphoma, MYC-negative Burkitt-like lymphoma with chr. 11q aberration and mantle cell lymphoma.

METHODS

We apply self-organizing map (SOM) machine learning to microarray-derived expression data to generate a holistic view on the transcriptome landscape of lymphomas, to describe the multidimensional nature of gene regulation and to pursue a modular view on co-expression. Expression data were complemented by pathological, genetic and clinical characteristics.

RESULTS

We present a transcriptome map of B cell lymphomas that allows visual comparison between the SOM portraits of different lymphoma strata and individual cases. It decomposes into one dozen modules of co-expressed genes related to different functional categories, to genetic defects and to the pathogenesis of lymphomas. On a molecular level, this disease rather forms a continuum of expression states than clearly separated phenotypes. We introduced the concept of combinatorial pattern types (PATs) that stratifies the lymphomas into nine PAT groups and, on a coarser level, into five prominent cancer hallmark types with proliferation, inflammation and stroma signatures. Inflammation signatures in combination with healthy B cell and tonsil characteristics associate with better overall survival rates, while proliferation in combination with inflammation and plasma cell characteristics worsens it. A phenotypic similarity tree is presented that reveals possible progression paths along the transcriptional dimensions. Our analysis provided a novel look on the transition range between FL and DLBCL, on DLBCL with poor prognosis showing expression patterns resembling that of Burkitt's lymphoma and particularly on 'double-hit' MYC and BCL2 transformed lymphomas.

CONCLUSIONS

The transcriptome map provides a tool that aggregates, refines and visualizes the data collected in the MMML study and interprets them in the light of previous knowledge to provide orientation and support in current and future studies on lymphomas and on other cancer entities.

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

Burkitt lymphoma (BL) is the most common pediatric B-cell lymphoma. Here, within the International Cancer Genome Consortium, the authors performed whole genome and transcriptome sequencing of 39 sporadic BL, describing the landscape of mutations, structural variants, and mutational processes that underpin this disease how alterations on different cellular levels cooperate in deregulating key pathways and complexes.

Stress-related hormones in association with periodontal condition in adolescents-results of the epidemiologic LIFE Child study

OBJECTIVES

The aim of this study was to investigate the associations between blood levels of stress-related hormones and early signs of periodontal disease in children and adolescents.

MATERIALS AND METHODS

Within the LIFE (Leipzig research center for civilization diseases) Child study, 498 adolescents (10 to 18 years) were included. Early signs of periodontal inflammation were measured by probing depth (PD) at six index teeth (16, 11, 26, 36, 31, 46). Blood levels of stress-related hormones (cortisol, dehydroepiandosterone-sulfate [DHEA-S]) and, additionally interleukine-6 (IL-6) were measured. Socioeconomic status, oral hygiene, orthodontic appliances, and nutritional status, recorded by body-mass-index-standard-deviation-score (BMI-SDS), were considered as confounding factors. Additionally, in 98 participants, an oral chairside active matrix metalloproteinase-8 (aMMP-8) test was performed. Statistical tests are the Mann-Whitney U tests, chi-squared tests and multivariate logistic regression model.

RESULTS

IL-6, BMI-SDS as well as positive aMMP-8 test result were significantly associated with maximum PD > 3 mm (p < 0.05). However, no statistically significant associations between stress-related hormones (cortisol and DHEA-S) and presence of maximum PD > 3 mm were found (p > 0.05). Higher DHEA-S and BMI were associated with positive aMMP-8 result, even after adjusting for age and gender (p = 0.027, p_adj = 0.026).

CONCLUSION

The results reveal no associations between PD and stress-related hormones cortisol and DHEA-S. aMMP-8 test result might be associated with DHEA-S level. Nutritional status seems to influence periodontal disease in adolescents.

CLINICAL RELEVANCE

DHEA-S and BMI-SDS show associations with early signs of periodontal disease in adolescents aged 10 to 18 years. This association should be confirmed by the investigation of high-risk groups.

Gene expression profiling reveals a close relationship between follicular lymphoma grade 3A and 3B, but distinct profiles of follicular lymphoma grade 1 and 2

A linear progression model of follicular lymphomas (FL) FL1, FL2 and FL3A has been favored, since FL3A often co-exist with an FL1/2 component. FL3B, in contrast, is thought to be more closely related to diffuse large B-cell lymphoma (DLBCL), and both are often simultaneously present in one tumor (DLBCL/FL3B). To obtain more detailed insights into follicular lymphoma progression, a comprehensive analysis of a well-defined set of FL1/2 (n=22), FL3A (n=16), FL3B (n=6), DLBCL/FL3B (n=9), and germinal center B-cell-type diffuse large B-cell lymphoma (n=45) was undertaken using gene expression profiling, immunohistochemical stainings and genetic analyses by fluorescence in situ hybridization. While immunohistochemical (CD10, IRF4/MUM1, Ki67, BCL2, BCL6) and genetic profiles (translocations of BCL2, BCL6 and MYC) delineate FL1-3A from FL3B and DLBCL/FL3B, significant differences were observed between FL1/2 and FL3A upon gene expression profiling. Interestingly, FL3B turned out to be closely related to FL3A, not categorizing within a separate gene expression cluster, and both FL3A and FL3B showed overlapping profiles in between FL1/2 and diffuse large B-cell lymphoma. Finally, based upon their gene expression pattern, DLBCL/FL3B represent a composite form of FL3B and DLBCL, with the majority of samples more closely resembling the latter. The fact that gene expression profiling clearly separated FL1/2 from both FL3A and FL3B suggests a closer biological relationship between the latter. This notion, however, is in contrast to immunohistochemical and genetic profiles of the different histological FL subtypes that point to a closer relationship between FL1/2 and FL3A, and separates them from FL3B.

Development and Implementation of an Analysis Tool for Array-based Comparative Genomic Hybridization

Objectives: Array-comparative genomic hybridization (aCGH) is a high-throughput method to detect and map copy number aberrations in the genome. Multi-step analysis of high-dimensional data requires an integrated suite of bioinformatic tools. In this paperwe detail an analysis pipeline for array CGH data.

Methods: We developed an analysis tool for array CGH data which supports single and multi-chip analyses as well as combined analyses with paired mRNA gene expression data. The functions supporting relevant steps of analysis were implemented using the open source software R and combined as package aCGHPipeline. Analysis methods were illustrated using 189 CGH arrays of aggressive B-cell lymphomas.

Results: The package covers data input, quality control, normalization, segmentation and classification. For multi-chip analysis aCGHPipeline offers an algorithm for automatic delineation of recurrent regions. This task was performed manuallyup to now. The package also supports combined analysis with mRNA gene expression data. Outputs consist of HTML documents to facilitate communication with clinical partners.

Conclusions: The R package aCGHPipeline supports basic tasks of single and multi-chip analysis of array CGH data.

Advanced patient age at diagnosis of diffuse large B-cell lymphoma is associated with molecular characteristics including ABC-subtype and high expression of MYC

The incidence of diffuse large B-cell lymphoma (DLBCL) increases with age being patient age at diagnosis an adverse prognostic factor. However, elderly patients are often underrepresented in common studies. To investigate the effect between age and biological characteristics in DLBCL, we analyzed data of 1534 patients encompassing all adult age groups, enriched for the age ≥75 years. Follicular lymphoma (FL) grade 3B with histopathological characteristics of DLBCLs were included. Gender, centroblastic cytology, FL grade 3B morphology, CD10 expression, and ABC/non-GCB-subtype were significantly associated with age after correction for multiple testing and after adjusting for cohorts. Analysis of a subgroup points towards an association of MYC expression with age. Our data indicate that biological features of DLBCL and FL grade 3B are associated with increasing age among adult patients. The prevalence of the ABC/non-GCB-subtype in elderly patients suggests that therapies targeting this molecular subtype should be specifically explored in this subgroup.

TP53 mutation and survival in aggressive B cell lymphoma

TP53 is mutated in 20-25% of aggressive B-cell lymphoma (B-NHL). To date, no studies have addressed the impact of TP53 mutations in prospective clinical trial cohorts. To evaluate the impact of TP53 mutation to current risk models in aggressive B-NHL, we investigated TP53 gene mutations within the RICOVER-60 trial. Of 1,222 elderly patients (aged 61-80 years) enrolled in the study and randomized to six or eight cycles of CHOP-14 with or without Rituximab (NCT00052936), 265 patients were analyzed for TP53 mutations. TP53 mutations were demonstrated in 63 of 265 patients (23.8%). TP53 mutation was associated with higher LDH (65% vs. 37%; p < 0.001), higher international prognostic index-Scores (IPI 4/5 27% vs. 12%; p = 0.025) and B-symptoms (41% vs. 24%; p = 0.011). Patients with TP53 mutation were less likely to obtain a complete remission CR/CRu (CR unconfirmed) 61.9% (mut) vs. 79.7% (wt) (p = 0.007). TP53 mutations were associated with decreased event-free (EFS), progression-free (PFS) and overall survival (OS) (median observation time of 40.2 months): the 3 year EFS, PFS and OS were 42% (vs. 60%; p = 0.012), 42% (vs. 67.5%; p < 0.001) and 50% (vs. 76%; p < 0.001) for the TP53 mutation group. In a Cox proportional hazard analysis adjusting for IPI-factors and treatment arms, TP53 mutation was shown to be an independent predictor of EFS (HR 1.5), PFS (HR 2.0) and OS (HR 2.3; p < 0.001). TP53 mutations are independent predictors of survival in untreated patients with aggressive CD20+ lymphoma. TP53 mutations should be considered for risk models in DLBCL and strategies to improve outcome for patients with mutant TP53 must be developed.

aMMP-8 in correlation to caries and periodontal condition in adolescents-results of the epidemiologic LIFE child study

OBJECTIVES

The suitability of a chairside aMMP-8 test in determination of periodontal inflammation and caries in adolescents was assessed. Secondly, the influence of orthodontic treatment on aMMP-8 test result was analyzed.

MATERIALS AND METHODS

Within the LIFE Child study, 434 adolescents (10 to 18 years) were included. Clinical dental examinations comprised caries experience (DMF/T-Index), signs of periodontal inflammation (probing pocket depth, PPD; community periodontal index of treatment needs; CPITN) at six index teeth and oral hygiene (OH). Information about orthodontic appliances (OA) and socioeconomic status (SES) were obtained by validated questionnaires. Test's sensitivity and specificity to detect periodontal inflammation and carious lesions were evaluated. The influence of OA on the test result was analyzed (multivariate model).

RESULTS

No associations between age, gender, SES or OH, and test outcome were found (p > 0.05). Positive test results correlated to periodontal findings (CPITN, mean PPD; p < 0.001). However, for the detection of ≥ 1 site(s) with PPD ≥ 4 mm, the test's sensitivity and specificity were found to be 61 and 69%, respectively. Multivariate analysis revealed a higher probability for a positive test result in cases of fixed OA (odds ratio 5.02, 95% confidence interval 1.90-13.19). The test had no diagnostic value considering carious lesions.

CONCLUSIONS

The chairside aMMP-8 test does not reliably identify adolescents with periodontal inflammation. Positive test results were more frequent in case of OA.

CLINICAL RELEVANCE

The chairside aMMP-8 test is no appropriate tool to screen children and adolescents neither for periodontal inflammation nor for carious lesions.

In Europe expression of EBNA2 is associated with poor survival in EBV-positive diffuse large B-cell lymphoma of the elderly

Epstein-Barr virus (EBV)-positive diffuse large B-cell lymphoma (DLBCL) of the elderly occurs by definition in patients above the age of 50 years without any known underlying immunodeficiency. We investigated the incidence and clinical relevance of this subtype in Europe with special attention to the EBV-latency type. Among the 598 DLBCL, 15 EBV-positive lymphomas fulfilling the criteria of EBV-positive DLBCL of the elderly were identified (2.5%). Patients with EBV-positive DLBCL expressing EBNA2 showed a significantly poorer overall survival than patients with EBNA2-negative EBV-positive DLBCL (p = 0.0156). The incidence of EBV-positive DLBCL of the elderly in Europe is much lower than in Asian countries (2.5% of all cases of DLBCL). Interestingly, the likelihood of EBV positivity did not increase with age in patient above 50 years. Among EBV-positive DLBCL of the elderly a subgroup with EBV-latency type III expressing EBNA2 can be identified, which shows a poor outcome.

[Gene expression analyses and their possible clinical benefit in head and neck cancer]

Systems biology approaches for mutational (exome analysis and targeted sequencing) and gene expression analysis (transcriptome-wide gene expression profiling) represent a new and growing scientific field in head and neck oncology. In addition to medical biological expertise, bioinformatic assistance is increasingly required. For squamous cell head and neck cancer (HNSCC), the recent molecular genetic single-gene and signal pathway observations represent basic research. Important aspects of this have now been significantly enhanced by systems biology approaches, which have grown into relevant areas of translational clinical research. It is now known that HPV16 is associated with genetic alterations at various locations, but also that it functionally affects genes not altered in their base sequence at the level of methylation. In transcriptome analyses, various consortia found matching clusters of gene expression and HPV16 association with the spectrum of somatic mutations. The differential methylation of gene promoters discovered in HPV16-driven HNSCC proved predictive for survival—even in HNSCC patients without HPV detection. The authors present an overview of some translationally relevant findings and venture an outlook on possible future clinical developments.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project "Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing", we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

The role of HPV RNA transcription, immune response-related gene expression and disruptive TP53 mutations in diagnostic and prognostic profiling of head and neck cancer

Stratification of head and neck squamous cell carcinomas (HNSCC) based on HPV16 DNA and RNA status, gene expression patterns, and mutated candidate genes may facilitate patient treatment decision. We characterize head and neck squamous cell carcinomas (HNSCC) with different HPV16 DNA and RNA (E6*I) status from 290 consecutively recruited patients by gene expression profiling and targeted sequencing of 50 genes. We show that tumors with transcriptionally inactive HPV16 (DNA+ RNA-) are similar to HPV-negative (DNA-) tumors regarding gene expression and frequency of TP53 mutations (47%, 8/17 and 43%, 72/167, respectively). We also find that an immune response-related gene expression cluster is associated with lymph node metastasis, independent of HPV16 status and that disruptive TP53 mutations are associated with lymph node metastasis in HPV16 DNA- tumors. We validate each of these associations in another large data set. Four gene expression clusters which we identify differ moderately but significantly in overall survival. Our findings underscore the importance of measuring the HPV16 RNA (E6*I) and TP53-mutation status for patient stratification and identify associations of an immune response-related gene expression cluster and TP53 mutations with lymph node metastasis in HNSCC.

Histopathology of rhinophyma - a clinical-histopathologic correlation

BACKGROUND

Whereas early rhinophyma is histopathologically thought to resemble fully developed rosacea, a fibromatous variant has previously been described for severe rhinophyma. In terms of clinical characteristics, recently a new Rhinophyma Severity Index (RHISI) was introduced.

METHODS

We studied 24 patients who had been treated with wide shave excisions for rhinophyma. Specimens were stained with hematoxylin-eosin, periodic acid-Schiff reaction and a panel of immunohistochemical stains and observed for any correlation between clinical severity and histopathologic features as well as for predictive markers of clinical recurrence.

RESULTS

There were no significant histopathologic differences between the groups reflecting the different clinical expressions. From a histopathologic perspective, clinically severe forms did not show with exclusive fibrotic changes. Further, there was no histopathologic marker predicting the clinical course or possible recurrence of the disease after surgical treatment. Only the clinical pre-operative RHISI score correlated with the postoperative outcome, with a high pre-operative RHISI being a risk factor for recurrence.

CONCLUSION

Histopathologic features do not correlate with the clinical expression of rhinophyma. An exclusively 'fibrotic' rhinophyma form does not appear to exist and could possibly be the result of sampling error based on small biopsies studied.

Association of NADPH oxidase polymorphisms with anthracycline-induced cardiotoxicity in the RICOVER-60 trial of patients with aggressive CD20(+) B-cell lymphoma

AIM

To identify gene variants responsible for anthracycline-induced cardiotoxicity.

PATIENTS & METHODS

Polymorphisms of the NADPH oxidase subunits and of the anthracycline transporters ABCC1, ABCC2 and SLC28A3 were genotyped in elderly patients (61-80 years) treated for aggressive CD20(+) B-cell lymphomas with CHOP-14 with or without rituximab and followed up for 3 years.

RESULTS

The accumulation of RAC2 subunit genotypes TA/AA among cases was statistically significant upon adjustment for gender, age and doxorubicin dose in a multivariate logistic regression analysis (OR: 2.3, p = 0.028; univariate: OR: 1.8, p = 0.077). RAC2 and CYBA genotypes were significantly associated with anthracycline-induced cardiotoxicity in a meta-analysis of this and a similar previous study.

CONCLUSION

Our results support the theory that NADPH oxidase is involved in anthracycline-induced cardiotoxicity. Original submitted 9 July 2014; Revision submitted 19 December 2014.

Biological characterization of adult MYC-translocation-positive mature B-cell lymphomas other than molecular Burkitt lymphoma

Chromosomal translocations affecting the MYC oncogene are the biological hallmark of Burkitt lymphomas but also occur in a subset of other mature B-cell lymphomas. If accompanied by a chromosomal break targeting the BCL2 and/or BCL6 oncogene these MYC translocation-positive (MYC(+)) lymphomas are called double-hit lymphomas, otherwise the term single-hit lymphomas is applied. In order to characterize the biological features of these MYC(+) lymphomas other than Burkitt lymphoma we explored, after exclusion of molecular Burkitt lymphoma as defined by gene expression profiling, the molecular, pathological and clinical aspects of 80 MYC-translocation-positive lymphomas (31 single-hit, 46 double-hit and 3 MYC(+)-lymphomas with unknown BCL6 status). Comparison of single-hit and double-hit lymphomas revealed no difference in MYC partner (IG/non-IG), genomic complexity, MYC expression or gene expression profile. Double-hit lymphomas more frequently showed a germinal center B-cell-like gene expression profile and had higher IGH and MYC mutation frequencies. Gene expression profiling revealed 130 differentially expressed genes between BCL6(+)/MYC(+) and BCL2(+)/MYC(+) double-hit lymphomas. BCL2(+)/MYC(+) double-hit lymphomas more frequently showed a germinal center B-like gene expression profile. Analysis of all lymphomas according to MYC partner (IG/non-IG) revealed no substantial differences. In this series of lymphomas, in which immunochemotherapy was administered in only a minority of cases, single-hit and double-hit lymphomas had a similar poor outcome in contrast to the outcome of molecular Burkitt lymphoma and lymphomas without the MYC break. Our data suggest that, after excluding molecular Burkitt lymphoma and pediatric cases, MYC(+) lymphomas are biologically quite homogeneous with single-hit and double-hit lymphomas as well as IG-MYC and non-IG-MYC(+) lymphomas sharing various molecular characteristics.

SOCS1 mutation subtypes predict divergent outcomes in diffuse large B-Cell lymphoma (DLBCL) patients

Suppressor of cytokine signaling 1 (SOCS1) is frequently mutated in primary mediastinal and diffuse large B-cell lymphomas (DLBCL). Currently, the prognostic relevance of these mutations in DLBCL is unknown. To evaluate the value of the SOCS1 mutation status as a prognostic biomarker in DLBCL patients, we performed full-length SOCS1 sequencing in tumors of 154 comprehensively characterized DLBCL patients. We identified 90 SOCS1 mutations in 16% of lymphomas. With respect to molecular consequences of mutations, we defined two distinct subtypes: those with truncating (major) and those with non-truncating mutations (minor), respectively. The SOCS1 mutated subgroup or the minor/major subtypes cannot be predicted on clinical grounds; however, assignment of four established gene-expression profile-based classifiers revealed significant associations of SOCS1 major cases with germinal center and specific pathway activation pattern signatures. Above all, SOCS1 major cases have an excellent overall survival, even better than the GCB-like subgroup. SOCS1 minor cases had a dismal survival, even worse than the ABC gene signature group. The SOCS1 mutation subsets retained prognostic significance in uni- and multivariate analyses. Together our data indicate that assessment of the SOCS1 mutation status is a single gene prognostic biomarker in DLBCL.

Circulating levels of TNF receptor II are prognostic for patients with peripheral T-cell non-Hodgkin lymphoma

PURPOSE

Peripheral T-cell non-Hodgkin lymphomas (T-NHL) represent a small but heterogeneous and clinically aggressive subset of NHLs with a poor outcome. Cytokines or their receptors might be associated with the clinical outcome of these lymphomas. Therefore, we tested whether gene variations and serum levels of soluble TNF receptor (TNFR)I (sTNFRI), sTNFRII, interleukin (IL)-10, or sIL-4R are predictive for treatment response in T-NHLs.

EXPERIMENTAL DESIGN

Peripheral blood DNA from 117 patients with T-NHL treated in prospective clinical trials was subjected to genotyping analysis. Whenever possible, pretreatment sera were obtained, and circulating levels of sTNFRI, sTNFRII, IL-10, and sIL-4R were determined with a specific capture enzyme-linked immunoassay.

RESULTS

Patients characterized by TNFRI-609GG (rs4149570) showed a trend toward better event free survival [EFS; univariate: P = 0.041; multivariate: HR, 1.76; confidence interval (CI), 0.99-3.14 with P = 0.056]. A protective role of IL-10-1087A, -824T, and -597A reported in another study was not confirmed in our cohort. Patients with circulating levels of soluble TNFRII ≥2.16 ng/mL had a 2.07-fold increased relative risk for shorter overall survival (OS; univariate: P = 0.0034; multivariate: HR, 2.07; CI, 0.92-4.70 with P = 0.081) and a 2.49-fold higher risk for shorter EFS (univariate: P = 0.00068; multivariate: HR, 2.49; CI, 1.22-5.08 with P = 0.012). Elevations of circulating levels of sTNFRI, IL-10, and sIL-4R are frequent, but the clinical response in these patients is not significantly different.

CONCLUSIONS

Our findings suggest a critical role for TNF-TNFR signaling for the clinical outcome of patients with peripheral T-NHLs.

Identification of genomic alterations associated with metastasis and cancer specific survival in clear cell renal cell carcinoma

PURPOSE

We identified regions of DNA copy number changes that are significantly associated with metastasis and clinical outcome in patients with clear cell renal cell carcinoma.

MATERIALS AND METHODS

We analyzed 53 primary clear cell renal cell carcinomas, including 31 metastasized and 22 nonmetastasized tumors, by array comparative genomic hybridization with a median resolution of 1 to 1.5 Mbp. To validate copy number aberrations with potential prognostic value we performed fluorescence in situ hybridization analysis using commercially available fluorescent probes.

RESULTS

We identified 5 recurrent chromosomal aberrations that were significantly associated with metastasis, including gains of 1q21.3, 12q13.12, 12q13.3q14.1 and 20q11.21q13.2, and loss of 9p21.3p24.1. The most prominent of them with the highest OR for metastatic risk were loss of 9p21.3p24.1, and gains of 1q21.3 and 20q11.21q13.32. Eight alterations involving chromosomes 7, 9, 12, 16 and 20 significantly correlated with shortened cancer specific survival. The lowest p values on Kaplan-Meier analysis showed losses of 9p21.3p24.1 and 9q32q33.1, and gains of 7q36.3 and 20q11.21q13.32. Fluorescence in situ hybridization done in the same cohort for the 4 select regions 1q21.3, 7q36.3, 9p21.3p24.1 and 20q11.21q13.32 clearly confirmed the results of array comparative genomic hybridization.

CONCLUSIONS

Data suggest that specific chromosomal alterations in clear cell renal cell carcinoma can be used to predict metastasis and cancer specific survival in patients with clear cell renal cell carcinoma. It seems possible to design a combined fluorescence in situ hybridization assay based on these genetic targets for outcome prediction, which can be used for routine diagnostics.

Detection of genomic aberrations in molecularly defined Burkitt's lymphoma by array-based, high resolution, single nucleotide polymorphism analysis

BACKGROUND

Knowledge about the genetic lesions that occur in Burkitt's lymphoma, besides the pathognomonic IG-MYC translocations, is limited.

DESIGN AND METHODS

Thirty-nine molecularly-defined Burkitt's lymphomas were analyzed with high-resolution single-nucleotide polymorphism chips for genomic imbalances and uniparental disomy. Imbalances were correlated to expression profiles and selected micro-RNA analysis. Translocations affecting the MYC locus were studied by fluorescence in situ hybridization.

RESULTS

We detected 528 copy number changes, defining 29 recurrently imbalanced regions. Five hundred and eighteen regions of uniparental disomy were found, but these were rarely recurrent. Combined imbalance mapping and expression profiling revealed a strong correlation between copy number and expression. Several recurrent imbalances affected the MYC pathway: the micro-RNA-supercluster 17-92 was frequently gained and the transcription factor E2F2 was recurrently deleted. Molecular Burkitt's lymphoma lacking MYC translocations showed MYC gains. Amplifications of the polymerase iota gene were associated with increased frequency of positions scored as aberrant.

CONCLUSIONS

The present findings suggest that uniparental disomies do not play a major role in the pathogenesis of Burkitt's lymphoma, whereas some genes may contribute to the development of this lymphoma through gene dosage effects. Amplifications of the polymerase iota gene may be functionally linked with increased genomic alterations in Burkitt's lymphoma. The pattern and rarity of chromosomal changes detectable, even at the high resolution employed here, together with aberrations of genes regulating MYC activity, support the hypothesis that deregulation of the MYC pathway is the major force driving the pathogenesis of Burkitt's lymphoma, but show that this deregulation is more complex than previously known.