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.

A time-resolved meta-analysis of consensus gene expression profiles during human T-cell activation

BACKGROUND

The coordinated transcriptional regulation of activated T-cells is based on a complex dynamic behavior of signaling networks. Given an external stimulus, T-cell gene expression is characterized by impulse and sustained patterns over the course. Here, we analyze the temporal pattern of activation across different T-cell populations to develop consensus gene signatures for T-cell activation.

RESULTS

Here, we identify and verify general biomarker signatures robustly evaluating T-cell activation in a time-resolved manner. We identify time-resolved gene expression profiles comprising 521 genes of up to 10 disjunct time points during activation and different polarization conditions. The gene signatures include central transcriptional regulators of T-cell activation, representing successive waves as well as sustained patterns of induction. They cover sustained repressed, intermediate, and late response expression rates across multiple T-cell populations, thus defining consensus biomarker signatures for T-cell activation. In addition, intermediate and late response activation signatures in CAR T-cell infusion products are correlated to immune effector cell-associated neurotoxicity syndrome.

CONCLUSION

This study is the first to describe temporally resolved gene expression patterns across T-cell populations. These biomarker signatures are a valuable source, e.g., monitoring transcriptional changes during T-cell activation with a reasonable number of genes, annotating T-cell states in single-cell transcriptome studies, or assessing dysregulated functions of human T-cell immunity.

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.

Transcriptomic signatures reveal a shift towards an anti-inflammatory gene expression profile but also the induction of type I and type II interferon signaling networks through aryl hydrocarbon receptor activation in murine macrophages

Introduction

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that regulates a broad range of target genes involved in the xenobiotic response, cell cycle control and circadian rhythm. AhR is constitutively expressed in macrophages (Mϕ), acting as key regulator of cytokine production. While proinflammatory cytokines, i.e., IL-1β, IL-6, IL-12, are suppressed through AhR activation, anti-inflammatory IL-10 is induced. However, the underlying mechanisms of those effects and the importance of the specific ligand structure are not yet completely understood.

Methods

Therefore, we have compared the global gene expression pattern in activated murine bone marrow-derived macrophages (BMMs) subsequently to exposure with either benzo[a]pyrene (BaP) or indole-3-carbinol (I3C), representing high-affinity vs. low-affinity AhR ligands, respectively, by means of mRNA sequencing. AhR dependency of observed effects was proved using BMMs from AhR-knockout (Ahr^-/-) mice.

Results and discussion

In total, more than 1,000 differentially expressed genes (DEGs) could be mapped, covering a plethora of AhR-modulated effects on basal cellular processes, i.e., transcription and translation, but also immune functions, i.e., antigen presentation, cytokine production, and phagocytosis. Among DEGs were genes that are already known to be regulated by AhR, i.e., Irf1, Ido2, and Cd84. However, we identified DEGs not yet described to be AhR-regulated in Mϕ so far, i.e., Slpi, Il12rb1, and Il21r. All six genes likely contribute to shifting the Mϕ phenotype from proinflammatory to anti-inflammatory. The majority of DEGs induced through BaP were not affected through I3C exposure, probably due to higher AhR affinity of BaP in comparison to I3C. Mapping of known aryl hydrocarbon response element (AHRE) sequence motifs in identified DEGs revealed more than 200 genes not possessing any AHRE, and therefore being not eligible for canonical regulation. Bioinformatic approaches modeled a central role of type I and type II interferons in the regulation of those genes. Additionally, RT-qPCR and ELISA confirmed a AhR-dependent expressional induction and AhR-dependent secretion of IFN-γ in response to BaP exposure, suggesting an auto- or paracrine activation pathway of Mϕ.

Atypical fibroxanthoma and pleomorphic dermal sarcoma - gene expression analysis compared with undifferentiated cutaneous squamous cell carcinoma

BACKGROUND

The histogenetic origin of atypical fibroxanthoma (AFX) and pleomorphic dermal sarcoma (PDS) has not been definitively elucidated. In addition to a fibroblastic origin, a keratinocytic differentiation is discussed due to strong clinical, histomorphological and molecular genetic similarities with undifferentiated cutaneous squamous cell carcinoma (cSCC).

PATIENTS AND METHODS

56 cases (36 AFXs, 8 PDSs, 12 undifferentiated cSCCs) were evaluated for their clinical, histomorphological, and immunohistochemical characteristics. RNA transcriptome analysis was performed on 18 cases (6 AFXs/PDSs, 6 undifferentiated cSCCs, 6 differentiated cSCCs).

RESULTS

Clinically, the strong similarities in age, gender and tumor location were confirmed. Without further immunohistochemical staining, histomorphological differentiation between AFX/PDS and undifferentiated cSCC is often impossible. Principal component analysis of the RNA transcriptome analysis showed that AFX/PDS and differentiated cSCC each formed their own cluster, while the undifferentiated cSCCs fall in between these two groups, but without forming a cluster of their own. When examining differentially expressed genes (DEGs), the heat maps showed that there were cases within the undifferentiated cSCC that were more likely to be AFX/PDS than differentiated cSCC based on their expression profile.

CONCLUSIONS

The results provide evidence of molecular similarities between AFX/PDS and undifferentiated cSCC and suggest a common histogenetic origin.

RNA sequencing of glioblastoma tissue slice cultures reveals the effects of treatment at the transcriptional level

One of the major challenges in cancer research is finding models that closely resemble tumors within patients. Human tissue slice cultures are a promising approach to provide a model of the patient's tumor biology ex vivo. Recently, it was shown that these slices can be successfully analyzed by whole transcriptome sequencing as well as automated histochemistry, increasing their usability as preclinical model. Glioblastoma multiforme (GBM) is a highly malignant brain tumor with poor prognosis and little is known about its genetic background and heterogeneity regarding therapy success. In this study, tissue from the tumors of 25 patients with primary GBM was processed into slice cultures and treated with standard therapy (irradiation and temozolomide). Total RNA sequencing and automated histochemistry were performed to enable analysis of treatment effects at a transcriptional and histological level. Slice cultures from long‐term survivors (overall survival [OS] > 24 months) exhibited more apoptosis than cultures from patients with shorter OS. Proliferation within these slices was slightly increased in contrast to other groups, but not significantly. Among all samples, 58 protein‐coding genes were upregulated and 32 downregulated in treated vs. untreated slice cultures. In general, an upregulation of DNA damage‐related and cell cycle checkpoint genes as well as enrichment of genotoxicity pathways and p53‐dependent signaling was found after treatment. Overall, the current study reproduces knowledge from former studies regarding the feasibility of transcriptomic analyses and automated histology in tissue slice cultures. We further demonstrate that the experimental data merge with the clinical follow‐up of the patients, which improves the applicability of our model system.

Low Energy Electron Irradiation Is a Potent Alternative to Gamma Irradiation for the Inactivation of (CAR-)NK-92 Cells in ATMP Manufacturing

Background

With increasing clinical use of NK-92 cells and their CAR-modified derivatives in cancer immunotherapy, there is a growing demand for efficient production processes of these “off-the-shelf” therapeutics. In order to ensure safety and prevent the occurrence of secondary tumors, (CAR-)NK-92 cell proliferation has to be inactivated before transfusion. This is commonly achieved by gamma irradiation. Recently, we showed proof of concept that low energy electron irradiation (LEEI) is a new method for NK-92 inactivation. LEEI has several advantages over gamma irradiation, including a faster reaction time, a more reproducible dose rate and much less requirements on radiation shielding. Here, LEEI was further evaluated as a promising alternative to gamma irradiation yielding cells with highly maintained cytotoxic effector function.

Methods

Effectiveness and efficiency of LEEI and gamma irradiation were analyzed using NK-92 and CD123-directed CAR-NK-92 cells. LEE-irradiated cells were extensively characterized and compared to gamma-irradiated cells via flow cytometry, cytotoxicity assays, and comet assays, amongst others.

Results

Our results show that both irradiation methods caused a progressive decrease in cell viability and are, therefore, suitable for inhibition of cell proliferation. Notably, the NK-mediated specific lysis of tumor cells was maintained at stable levels for three days post-irradiation, with a trend towards higher activities after LEEI treatment as compared to gamma irradiation. Both gamma irradiation as well as LEEI led to substantial DNA damage and an accumulation of irradiated cells in the G2/M cell cycle phases. In addition, transcriptomic analysis of irradiated cells revealed approximately 12-fold more differentially expressed genes two hours after gamma irradiation, compared to LEEI. Analysis of surface molecules revealed an irradiation-induced decrease in surface expression of CD56, but no changes in the levels of the activating receptors NKp46, NKG2D, or NKp30.

Conclusions

The presented data show that LEEI inactivates (CAR-)NK-92 cells as efficiently as gamma irradiation, but with less impact on the overall gene expression. Due to logistic advantages, LEEI might provide a superior alternative for the manufacture of (CAR-)NK-92 cells for clinical application.

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.

New insights into valve-related intramural and intracellular bacterial diversity in infective endocarditis

Aims

In infective endocarditis (IE), a severe inflammatory disease of the endocardium with an unchanged incidence and mortality rate over the past decades, only 1% of the cases have been described as polymicrobial infections based on microbiological approaches. The aim of this study was to identify potential biodiversity of bacterial species from infected native and prosthetic valves. Furthermore, we compared the ultrastructural micro-environments to detect the localization and distribution patterns of pathogens in IE.

Material and methods

Using next-generation sequencing (NGS) of 16S rDNA, which allows analysis of the entire bacterial community within a single sample, we investigated the biodiversity of infectious bacterial species from resected native and prosthetic valves in a clinical cohort of 8 IE patients. Furthermore, we investigated the ultrastructural infected valve micro-environment by focused ion beam scanning electron microscopy (FIB-SEM).

Results

Biodiversity was detected in 7 of 8 resected heart valves. This comprised 13 bacterial genera and 16 species. In addition to 11 pathogens already described as being IE related, 5 bacterial species were identified as having a novel association. In contrast, valve and blood culture-based diagnosis revealed only 4 species from 3 bacterial genera and did not show any relevant antibiotic resistance. The antibiotics chosen on this basis for treatment, however, did not cover the bacterial spectra identified by our amplicon sequencing analysis in 4 of 8 cases. In addition to intramural distribution patterns of infective bacteria, intracellular localization with evidence of bacterial immune escape mechanisms was identified.

Conclusion

The high frequency of polymicrobial infections, pathogen diversity, and intracellular persistence of common IE-causing bacteria may provide clues to help explain the persistent and devastating mortality rate observed for IE. Improved bacterial diagnosis by 16S rDNA NGS that increases the ability to tailor antibiotic therapy may result in improved outcomes.

Cell cycle, oncogenic and tumor suppressor pathways regulate numerous long and macro non-protein-coding RNAs

Background

The genome is pervasively transcribed but most transcripts do not code for proteins, constituting non-protein-coding RNAs. Despite increasing numbers of functional reports of individual long non-coding RNAs (lncRNAs), assessing the extent of functionality among the non-coding transcriptional output of mammalian cells remains intricate. In the protein-coding world, transcripts differentially expressed in the context of processes essential for the survival of multicellular organisms have been instrumental in the discovery of functionally relevant proteins and their deregulation is frequently associated with diseases. We therefore systematically identified lncRNAs expressed differentially in response to oncologically relevant processes and cell-cycle, p53 and STAT3 pathways, using tiling arrays.

Results

We found that up to 80% of the pathway-triggered transcriptional responses are non-coding. Among these we identified very large macroRNAs with pathway-specific expression patterns and demonstrated that these are likely continuous transcripts. MacroRNAs contain elements conserved in mammals and sauropsids, which in part exhibit conserved RNA secondary structure. Comparing evolutionary rates of a macroRNA to adjacent protein-coding genes suggests a local action of the transcript. Finally, in different grades of astrocytoma, a tumor disease unrelated to the initially used cell lines, macroRNAs are differentially expressed.

Conclusions

It has been shown previously that the majority of expressed non-ribosomal transcripts are non-coding. We now conclude that differential expression triggered by signaling pathways gives rise to a similar abundance of non-coding content. It is thus unlikely that the prevalence of non-coding transcripts in the cell is a trivial consequence of leaky or random transcription events.

Analysis of the STAT3 interactome using in-situ biotinylation and SILAC

Signal transducer and activator of transcription 3 (STAT3) is activated by a variety of cytokines and growth factors. To generate a comprehensive data set of proteins interacting specifically with STAT3, we applied stable isotope labeling with amino acids in cell culture (SILAC). For high-affinity pull-down using streptavidin, we fused STAT3 with a short peptide tag allowing biotinylation in situ (bio-tag), which did not affect STAT3 functions. By this approach, 3642 coprecipitated proteins were detected in human embryonic kidney-293 cells. Filtering using statistical and functional criteria finally extracted 136 proteins as putative interaction partners of STAT3. Both, a physical interaction network analysis and the enrichment of known and predicted interaction partners suggested that our filtering criteria successfully enriched true STAT3 interactors. Our approach identified numerous novel interactors, including ones previously predicted to associate with STAT3. By reciprocal coprecipitation, we were able to verify the physical association between STAT3 and selected interactors, including the novel interaction with TOX4, a member of the TOX high mobility group box family. Applying the same method, we next investigated the activation-dependency of the STAT3 interactome. Again, we identified both known and novel interactions. Thus, our approach allows to study protein-protein interaction effectively and comprehensively.

BIOLOGICAL SIGNIFICANCE

The location, activity, function, degradation, and synthesis of proteins are significantly regulated by interactions of proteins with other proteins, biopolymers and small molecules. Thus, the comprehensive characterization of interactions of proteins in a given proteome is the next milestone on the path to understanding the biochemistry of the cell. In order to generate a comprehensive interactome dataset of proteins specifically interacting with a selected bait protein, we fused our bait protein STAT3 with a short peptide tag allowing biotinylation in situ (bio-tag). This bio-tag allows an affinity pull-down using streptavidin but affected neither the activation of STAT3 by tyrosine phosphorylation nor its transactivating potential. We combined SILAC for accurate relative protein quantification, subcellular fractionation to increase the coverage of interacting proteins, high-affinity pull-down and a stringent filtering method to successfully analyze the interactome of STAT3. With our approach we confirmed several already known and identified numerous novel STAT3 interactors. The approach applied provides a rapid and effective method, which is broadly applicable for studying protein-protein interactions and their dependency on post-translational modifications.

Cyclophilins contribute to Stat3 signaling and survival of multiple myeloma cells

Signal transducer and activator of transcription 3 (Stat3) is the major mediator of interleukin-6 (IL-6) family cytokines. In addition, Stat3 is known to be involved in the pathophysiology of many malignancies. Here, we show that the cis-trans peptidyl-prolyl isomerase cyclophilin (Cyp) B specifically interacts with Stat3, whereas the highly related CypA does not. CypB knockdown inhibited the IL-6-induced transactivation potential but not the tyrosine phosphorylation of Stat3. Binding of CypB to Stat3 target promoters and alteration of the intranuclear localization of Stat3 on CypB depletion suggested a nuclear function of Stat3/CypB interaction. By contrast, CypA knockdown inhibited Stat3 IL-6-induced tyrosine phosphorylation and nuclear translocation. The Cyp inhibitor cyclosporine A (CsA) caused similar effects. However, Stat1 activation in response to IL-6 or interferon-gamma was not affected by Cyp silencing or CsA treatment. As a result, Cyp knockdown shifted IL-6 signaling to a Stat1-dominated pathway. Furthermore, Cyp depletion or treatment with CsA induced apoptosis in IL-6-dependent multiple myeloma cells, whereas an IL-6-independent line was not affected. Thus, Cyps support the anti-apoptotic action of Stat3. Taken together, CypA and CypB both play pivotal roles, yet at different signaling levels, for Stat3 activation and function. These data also suggest a novel mechanism of CsA action.

Interleukin-6 dependent survival of multiple myeloma cells involves the Stat3-mediated induction of microRNA-21 through a highly conserved enhancer

Signal transducer and activator of transcription 3 (Stat3) is implicated in the pathogenesis of many malignancies and essential for IL-6-dependent survival and growth of multiple myeloma cells. Here, we demonstrate that the gene encoding oncogenic microRNA-21 (miR-21) is controlled by an upstream enhancer containing 2 Stat3 binding sites strictly conserved since the first observed evolutionary appearance of miR-21 and Stat3. MiR-21 induction by IL-6 was strictly Stat3 dependent. Ectopically raising miR-21 expression in myeloma cells in the absence of IL-6 significantly reduced their apoptosis levels. These data provide strong evidence that miR-21 induction contributes to the oncogenic potential of Stat3.