Dynamics in protein translation sustaining T cell preparedness.

In response to pathogenic threats, naive T cells rapidly transition from a quiescent to an activated state, yet the underlying mechanisms are incompletely understood. Using a pulsed SILAC approach, we investigated the dynamics of mRNA translation kinetics and protein turnover in human naive and activated T cells. Our datasets uncovered that transcription factors maintaining T cell quiescence had constitutively high turnover, which facilitated their depletion following activation. Furthermore, naive T cells maintained a surprisingly large number of idling ribosomes as well as 242 repressed mRNA species and a reservoir of glycolytic enzymes. These components were rapidly engaged following stimulation, promoting an immediate translational and glycolytic switch to ramp up the T cell activation program. Our data elucidate new insights into how T cells maintain a prepared state to mount a rapid immune response, and provide a resource of protein turnover, absolute translation kinetics and protein synthesis rates in T cells ( https://www.immunomics.ch ).

Systems analyses of the Fabry kidney transcriptome and its response to enzyme replacement therapy identified and cross-validated enzyme replacement therapy-resistant targets amenable to drug repurposing.

Fabry disease is a rare disorder caused by variations in the alpha-galactosidase gene. To a degree, Fabry disease is manageable via enzyme replacement therapy (ERT). By understanding the molecular basis of Fabry nephropathy (FN) and ERT's long-term impact, here we aimed to provide a framework for selection of potential disease biomarkers and drug targets. We obtained biopsies from eight control individuals and two independent FN cohorts comprising 16 individuals taken prior to and after up to ten years of ERT, and performed RNAseq analysis. Combining pathway-centered analyses with network-science allowed computation of transcriptional landscapes from four nephron compartments and their integration with existing proteome and drug-target interactome data. Comparing these transcriptional landscapes revealed high inter-cohort heterogeneity. Kidney compartment transcriptional landscapes comprehensively reflected differences in FN cohort characteristics. With exception of a few aspects, in particular arteries, early ERT in patients with classical Fabry could lastingly revert FN gene expression patterns to closely match that of control individuals. Pathways nonetheless consistently altered in both FN cohorts pre-ERT were mostly in glomeruli and arteries and related to the same biological themes. While keratinization-related processes in glomeruli were sensitive to ERT, a majority of alterations, such as transporter activity and responses to stimuli, remained dysregulated or reemerged despite ERT. Inferring an ERT-resistant genetic module of expressed genes identified 69 drugs for potential repurposing matching the proteins encoded by 12 genes. Thus, we identified and cross-validated ERT-resistant gene product modules that, when leveraged with external data, allowed estimating their suitability as biomarkers to potentially track disease course or treatment efficacy and potential targets for adjunct pharmaceutical treatment.

Protein Expression of AEBP1, MCM4, and FABP4 Differentiate Osteogenic, Adipogenic, and Mesenchymal Stromal Stem Cells.

Mesenchymal stem cells (MSCs) gain an increasing focus in the field of regenerative medicine due to their differentiation abilities into chondrocytes, adipocytes, and osteoblastic cells. However, it is apparent that the transformation processes are extremely complex and cause cellular heterogeneity. The study aimed to characterize differences between MSCs and cells after adipogenic (AD) or osteoblastic (OB) differentiation at the proteome level. Comparative proteomic profiling was performed using tandem mass spectrometry in data-independent acquisition mode. Proteins were quantified by deep neural networks in library-free mode and correlated to the Molecular Signature Database (MSigDB) hallmark gene set collections for functional annotation. We analyzed 4108 proteins across all samples, which revealed a distinct clustering between MSCs and cell differentiation states. Protein expression profiling identified activation of the Peroxisome proliferator-activated receptors (PPARs) signaling pathway after AD. In addition, two distinct protein marker panels could be defined for osteoblastic and adipocytic cell lineages. Hereby, overexpression of AEBP1 and MCM4 for OB as well as of FABP4 for AD was detected as the most promising molecular markers. Combination of deep neural network and machine-learning algorithms with data-independent mass spectrometry distinguish MSCs and cell lineages after adipogenic or osteoblastic differentiation. We identified specific proteins as the molecular basis for bone formation, which could be used for regenerative medicine in the future.

The contribution of active and passive mechanisms of 5mC and 5hmC removal in human T lymphocytes is differentiation- and activation-dependent.

In mammals, the 5'-methylcytosine (5mC) modification in the genomic DNA contributes to the dynamic control of gene expression. 5mC erasure is required for the activation of developmental programs and occurs either by passive dilution through DNA replication, or by enzymatic oxidation of the methyl mark to 5-hydroxymethylcytosine (5hmC), which can persist as such or undergo further oxidation and enzymatic removal. The relative contribution of each mechanism to epigenetic control in dynamic biological systems still remains a compelling question. To explore this critical issue, we used primary human T lymphocytes, in which two cellular states can be clearly identified, namely quiescent naïve T cells, which are slowly or rarely proliferating, and rapidly proliferating activated T cells. We found that active mechanisms of methylation removal were selectively at work in naïve T cells, while memory T lymphocytes entirely relied on passive, replication-dependent dilution, suggesting that proliferative capacity influences the choice of the preferential demethylation mechanism. Active processes of demethylation appear to be critical in quiescent naïve T lymphocytes for the maintenance of regulatory regions poised for rapid responses to physiological stimuli.

KMT2D mutations and TP53 disruptions are poor prognostic biomarkers in mantle cell lymphoma receiving high-dose therapy: a FIL study.

In recent years, the outcome of mantle cell lymphoma (MCL) has improved, especially in younger patients, receiving cytarabine-containing chemoimmunotherapy and autologous stem cell transplantation. Nevertheless, a proportion of MCL patients still experience early failure. To identify biomarkers anticipating failure of intensive chemotherapy in MCL, we performed target resequencing and DNA profiling of purified tumor samples collected from patients enrolled in the prospective FIL-MCL0208 phase 3 trial (high-dose chemoimmunotherapy followed by autologous transplantation and randomized lenalidomide maintenance). Mutations of KMT2D and disruption of TP53 by deletion or mutation associated with an increased risk of progression and death, both in univariate and multivariate analysis. By adding KMT2D mutations and TP53 disruption to the MIPI-c backbone, we derived a new prognostic index, the "MIPI-genetic" ("MIPI- g"). The "MIPI-g" improved the model discrimination ability compared to the MIPI-c alone, defining three risk groups: i) low-risk patients (4-year progression free survival and overall survival of 72.0% and 94.5%); ii) inter-mediate-risk patients (4-year progression free survival and overall survival of 42.2% and 65.8%) and iii) high-risk patients (4-year progression free survival and overall survival of 11.5% and 44.9%). Our results: i) confirm that TP53 disruption identifies a high-risk population characterized by poor sensitivity to conventional or intensified chemotherapy; ii) provide the pivotal evidence that patients harboring KMT2D mutations share the same poor outcome as patients harboring TP53 disruption; and iii) allow to develop a tool for the identification of high-risk MCL patients for whom novel therapeutic strategies need to be investigated. (Trial registered at clinicaltrials.gov identifier: NCT02354313).

cuRnet: an R package for graph traversing on GPU.

BACKGROUND: R has become the de-facto reference analysis environment in Bioinformatics. Plenty of tools are available as packages that extend the R functionality, and many of them target the analysis of biological networks. Several algorithms for graphs, which are the most adopted mathematical representation of networks, are well-known examples of applications that require high-performance computing, and for which classic sequential implementations are becoming inappropriate. In this context, parallel approaches targeting GPU architectures are becoming pervasive to deal with the execution time constraints. Although R packages for parallel execution on GPUs are already available, none of them provides graph algorithms. RESULTS: This work presents cuRnet, a R package that provides a parallel implementation for GPUs of the breath-first search (BFS), the single-source shortest paths (SSSP), and the strongly connected components (SCC) algorithms. The package allows offloading computing intensive applications to GPU devices for massively parallel computation and to speed up the runtime up to one order of magnitude with respect to the standard sequential computations on CPU. We have tested cuRnet on a benchmark of large protein interaction networks and for the interpretation of high-throughput omics data thought network analysis. CONCLUSIONS: cuRnet is a R package to speed up graph traversal and analysis through parallel computation on GPUs. We show the efficiency of cuRnet applied both to biological network analysis, which requires basic graph algorithms, and to complex existing procedures built upon such algorithms.

Correction to: cuRnet: an R package for graph traversing on GPU.

After publication of this supplement article [1], it was brought to our attention that reference 10 and reference 12 in the article are incorrect.

Identification of a new subclass of ALK-negative ALCL expressing aberrant levels of ERBB4 transcripts.

Anaplastic large-cell lymphoma (ALCL) is a clinical and biological heterogeneous disease that includes systemic anaplastic lymphoma kinase (ALK)-positive and ALK-negative entities. To discover biomarkers and/or genes involved in ALK-negative ALCL pathogenesis, we applied the cancer outlier profile analysis algorithm to a gene expression profiling data set including 249 cases of T-cell non-Hodgkin lymphoma and normal T cells. Ectopic coexpression of ERBB4 and COL29A1 genes was detected in 24% of ALK-negative ALCL patients. RNA sequencing and 5' RNA ligase-mediated rapid amplification of complementary DNA ends identified 2 novel ERBB4-truncated transcripts displaying intronic transcription start sites. By luciferase assays, we defined that the expression of ERBB4-aberrant transcripts is promoted by endogenous intronic long terminal repeats. ERBB4 expression was confirmed at the protein level by western blot analysis and immunohistochemistry. Lastly, we demonstrated that ERBB4-truncated forms show oncogenic potentials and that ERBB4 pharmacologic inhibition partially controls ALCL cell growth and disease progression in an ERBB4-positive patient-derived tumorgraft model. In conclusion, we identified a new subclass of ALK-negative ALCL characterized by aberrant expression of ERBB4-truncated transcripts carrying intronic 5' untranslated regions.

Bromodomain and extra-terminal domain inhibition modulates the expression of pathologically relevant microRNAs in diffuse large B-cell lymphoma.

Aberrant changes in microRNA expression contribute to lymphomagenesis. Bromodomain and extra-terminal domain inhibitors such as OTX015 (MK-8628, birabresib) have demonstrated preclinical and clinical activity in hematologic tumors. MicroRNA profiling of diffuse large B-cell lymphoma cells treated with OTX015 revealed changes in the expression levels of a limited number of microRNAs, including miR-92a-1-5p, miR-21-3p, miR-155-5p and miR-96-5p. Analysis of publicly available chromatin immunoprecipitation sequencing data of diffuse large B-cell lymphoma cells treated with bromodomain and extra-terminal domain (BET) inhibitors showed that the BET family member BRD4 bound to the upstream regulatory regions of multiple microRNA genes and that this binding decreased following BET inhibition. Alignment of our microRNA profiling data with the BRD4 chromatin immunoprecipitation sequencing data revealed that microRNAs downregulated by OTX015 also exhibited reduced BRD4 binding in their promoter regions following treatment with another bromodomain and extra-terminal domain inhibitor, JQ1, indicating that BRD4 contributes directly to microRNA expression in lymphoma. Treatment with bromodomain and extra-terminal domain inhibitors also decreased the expression of the arginine methyltransferase PRMT5, which plays a crucial role in B-cell transformation and negatively modulates the transcription of miR-96-5p. The data presented here indicate that in addition to previously observed effects on the expression of coding genes, bromodomain and extra-terminal domain inhibitors also modulate the expression of microRNAs involved in lymphomagenesis.

Population-based outcome analysis of diffuse large B-cell lymphoma in people living with HIV infection and competent individuals.

The prognostic factors and outcome of 58 acquired immunodeficiency syndrome-related diffuse large B-cell lymphoma (AR-DLBCL) patients from the Swiss HIV Cohort Study, diagnosed from 2004 to 2011, were compared with those of 326 immunocompetent (IC)-DLBCL from the Hematology Division of the Amedeo Avogadro University (Italy) and the Oncology Institute of Southern Switzerland. Median follow-up was 6 years; 5-year overall survival (OS) was 68% (95% CI: 63%-73%) in IC-DLBCL and 63% (95% CI: 49%-75%) in AR-DLBCL (P = .220). The acquired immunodeficiency syndrome-related lymphoma international prognostic index predicted OS in AR-DLBCL. Among 148 patients younger than 61 years (40 AR-DLBCL and 108 IC-DLBCL) treated with RCHOP/RCHOP-like regimens, 20 IC-DLBCL and 9 AR-DLBCL patients died and OS was not significantly different. A higher proportion of early deaths occurred in the AR-DLBCL: indeed, 1-year OS was 94% (95% CI: 87%-97%) in IC-DLBCL and 82% (95% CI: 66%-91%) in AR-DLBCL patients. After rituximab and active antiretroviral therapy introduction, AR-DLBCL and IC-DLBCL patients treated with curative intent have similar long-term survival.

PCSF: An R-package for network-based interpretation of high-throughput data.

With the recent technological developments a vast amount of high-throughput data has been profiled to understand the mechanism of complex diseases. The current bioinformatics challenge is to interpret the data and underlying biology, where efficient algorithms for analyzing heterogeneous high-throughput data using biological networks are becoming increasingly valuable. In this paper, we propose a software package based on the Prize-collecting Steiner Forest graph optimization approach. The PCSF package performs fast and user-friendly network analysis of high-throughput data by mapping the data onto a biological networks such as protein-protein interaction, gene-gene interaction or any other correlation or coexpression based networks. Using the interaction networks as a template, it determines high-confidence subnetworks relevant to the data, which potentially leads to predictions of functional units. It also interactively visualizes the resulting subnetwork with functional enrichment analysis.

Preclinical evaluation of the BET bromodomain inhibitor BAY 1238097 for the treatment of lymphoma.

The epigenome is often deregulated in cancer and treatment with inhibitors of bromodomain and extra-terminal proteins, the readers of epigenetic acetylation marks, represents a novel therapeutic approach. Here, we have characterized the anti-tumour activity of the novel bromodomain and extra-terminal (BET) inhibitor BAY 1238097 in preclinical lymphoma models. BAY 1238097 showed anti-proliferative activity in a large panel of lymphoma-derived cell lines, with a median 50% inhibitory concentration between 70 and 208 nmol/l. The compound showed strong anti-tumour efficacy in vivo as a single agent in two diffuse large B cell lymphoma models. Gene expression profiling showed BAY 1238097 targeted the NFKB/TLR/JAK/STAT signalling pathways, MYC and E2F1-regulated genes, cell cycle regulation and chromatin structure. The gene expression profiling signatures also highly overlapped with the signatures obtained with other BET Bromodomain inhibitors and partially overlapped with HDAC-inhibitors, mTOR inhibitors and demethylating agents. Notably, BAY 1238097 presented in vitro synergism with EZH2, mTOR and BTK inhibitors. In conclusion, the BET inhibitor BAY 1238097 presented promising anti-lymphoma preclinical activity in vitro and in vivo, mediated by the interference with biological processes driving the lymphoma cells. Our data also indicate the use of combination schemes targeting EZH2, mTOR and BTK alongside BET bromodomains.

DNA methylation profiling identifies two splenic marginal zone lymphoma subgroups with different clinical and genetic features.

Splenic marginal zone lymphoma is a rare lymphoma. Loss of 7q31 and somatic mutations affecting the NOTCH2 and KLF2 genes are the commonest genomic aberrations. Epigenetic changes can be pharmacologically reverted; therefore, identification of groups of patients with specific epigenomic alterations might have therapeutic relevance. Here we integrated genome-wide DNA-promoter methylation profiling with gene expression profiling, and clinical and biological variables. An unsupervised clustering analysis of a test series of 98 samples identified 2 clusters with different degrees of promoter methylation. The cluster comprising samples with higher-promoter methylation (High-M) had a poorer overall survival compared with the lower (Low-M) cluster. The prognostic relevance of the High-M phenotype was confirmed in an independent validation set of 36 patients. In the whole series, the High-M phenotype was associated with IGHV1-02 usage, mutations of NOTCH2 gene, 7q31-32 loss, and histologic transformation. In the High-M set, a number of tumor-suppressor genes were methylated and repressed. PRC2 subunit genes and several prosurvival lymphoma genes were unmethylated and overexpressed. A model based on the methylation of 3 genes (CACNB2, HTRA1, KLF4) identified a poorer-outcome patient subset. Exposure of splenic marginal zone lymphoma cell lines to a demethylating agent caused partial reversion of the High-M phenotype and inhibition of proliferation.

The novel atypical retinoid ST5589 down-regulates Aurora Kinase A and has anti-tumour activity in lymphoma pre-clinical models.

Despite the marked improvements in the treatment of lymphomas, there is still a need for new therapeutic agents. Synthetic retinoids represent a class of compounds with anti-cancer activity. Here, we report the preclinical activity of a new member of this class, the ST1926-derivative ST5589, in lymphomas. ST5589 presented a dose-dependent anti-proliferative activity in almost all of the 25 lymphoma cell lines analysed, with a median 50% inhibitory concentration of 433 nM. Apoptosis was observed in 8/11 cell lines. ST5589 induced changes in the gene expression profiles of the cell lines, including the down-regulation of Aurora Kinase A (AURKA). Specific gene expression signatures were associated with a higher sensitivity to the compound and combination of ST5589 with carfilzomib revealed the importance of proteasome activity in mediating the anti-tumour activity of ST5589. In conclusion, we have identified a new mechanism of action of atypical retinoids as anti-cancer compounds, and the encouraging results obtained with the new ST1926-derivative ST5589 provide the basis for further developments of the compound.

Deregulation of ETS1 and FLI1 contributes to the pathogenesis of diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) is the most common form of human lymphoma. DLBCL is a heterogeneous disease characterized by different genetic lesions. We herein report the functional characterization of a recurrent gain mapping on chromosome 11q24.3, found in 23% of 166 DLBCL cases analyzed. The transcription factors ETS1 and FLI1, located within the 11q24.3 region, had significantly higher expression in clinical samples carrying the gain. Functional studies on cell lines showed that ETS1 and FLI1 cooperate in sustaining DLBCL proliferation and viability and regulate genes involved in germinal center differentiation. Taken together, these data identify the 11q24.3 gain as a recurrent lesion in DLBCL leading to ETS1 and FLI1 deregulated expression, which can contribute to the pathogenesis of this disease.

PRDM1/BLIMP1 is commonly inactivated in anaplastic large T-cell lymphoma.

Anaplastic large cell lymphoma (ALCL) is a mature T-cell lymphoma that can present as a systemic or primary cutaneous disease. Systemic ALCL represents 2% to 5% of adult lymphoma but up to 30% of all pediatric cases. Two subtypes of systemic ALCL are currently recognized on the basis of the presence of a translocation involving the anaplastic lymphoma kinase ALK gene. Despite considerable progress, several questions remain open regarding the pathogenesis of both ALCL subtypes. To investigate the molecular pathogenesis and to assess the relationship between the ALK(+) and ALK(-) ALCL subtypes, we performed a genome-wide DNA profiling using high-density, single nucleotide polymorphism arrays on a series of 64 cases and 7 cell lines. The commonest lesions were losses at 17p13 and at 6q21, encompassing the TP53 and PRDM1 genes, respectively. The latter gene, coding for BLIMP1, was inactivated by multiple mechanisms, more frequently, but not exclusively, in ALK(-)ALCL. In vitro and in vivo experiments showed that that PRDM1 is a tumor suppressor gene in ALCL models, likely acting as an antiapoptotic agent. Losses of TP53 and/or PRDM1 were present in 52% of ALK(-)ALCL, and in 29% of all ALCL cases with a clinical implication.

Two main genetic pathways lead to the transformation of chronic lymphocytic leukemia to Richter syndrome.

Richter syndrome (RS) occurs in up to 15% of patients with chronic lymphocytic leukemia (CLL). Although RS, usually represented by the histologic transformation to a diffuse large B-cell lymphoma (DLBCL), is associated with a very poor outcome, especially when clonally related to the preexisting CLL, the mechanisms leading to RS have not been clarified. To better understand the pathogenesis of RS, we analyzed a series of cases including 59 RS, 28 CLL phase of RS, 315 CLL, and 127 de novo DLBCL. RS demonstrated a genomic complexity intermediate between CLL and DLBCL. Cell-cycle deregulation via inactivation of TP53 and of CDKN2A was a main mechanism in the histologic transformation from CLL phase, being present in approximately one half of the cases, and affected the outcome of the RS patients. A second major subgroup was characterized by the presence of trisomy 12 and comprised one third of the cases. Although RS shared some of the lesions seen in de novo DLBCL, its genomic profile was clearly separate. The CLL phase preceding RS had not a generalized increase in genomic complexity compared with untransformed CLL, but it presented clear differences in the frequency of specific genetic lesions.

Distinct patterns of global promoter methylation in early stage chronic lymphocytic leukemia.

Genomic and epigenomic studies of chronic lymphocytic leukemia (CLL) are reshaping our understanding of the disease and have provided new perspectives for a more individualized diagnosis and new potential therapeutic targets. In this study, the global promoter methylation profile was determined in highly purified B-cells from 37 (Binet stage A) CLL patients, using high-resolution methylation microarrays (27,578 CpG). Overall, the methylation pattern correlated with the major biological (ZAP-70 and CD38), and molecular (IGHV mutation) markers, distinguishing CLL cases according to IGHV mutational status. Cell adhesion molecules were enriched in the signature of unmutated (UM) versus mutated (M-) CLL. Moreover, in M-CLL CpG hyper-methylation in three genes, including SPG20, was significantly anti-correlated with the corresponding gene expression level. Finally, the correlation between the methylation pattern and clinical parameters was investigated. Notably, out of 42 methyl-probes that were significantly associated with progression free survival (PFS), hyper-methylation of SPG20 was also positively associated with PFS. These data support the notion that epigenetic changes have clinical impact in CLL and may contribute to the identification of novel candidate disease-associated genes potentially useful to predict the clinical outcome of early stage CLL patients.

Identification of a 3-gene model as a powerful diagnostic tool for the recognition of ALK-negative anaplastic large-cell lymphoma.

Anaplastic large-cell lymphomas (ALCLs) are a group of clinically and biologically heterogeneous diseases including the ALK(+) and ALK(-) systemic forms. Whereas ALK(+) ALCLs are molecularly characterized and can be readily diagnosed, specific immunophenotypic or genetic features to define ALK(-) ALCL are missing, and their distinction from other T-cell non-Hodgkin lymphomas (T-NHLs) remains controversial. In the present study, we undertook a transcriptional profiling meta-analysis of 309 cases, including ALCL and other primary T-NHL samples. Pathway discovery and prediction analyses defined a minimum set of genes capable of recognizing ALK(-) ALCL. Application of quantitative RT-PCR in independent datasets from cryopreserved and formalin-fixed paraffin-embedded samples validated a 3-gene model (TNFRSF8, BATF3, and TMOD1) able to successfully separate ALK(-) ALCL from peripheral T-cell lymphoma not otherwise specified, with overall accuracy near 97%. In conclusion, our data justify the possibility of translating quantitative RT-PCR protocols to routine clinical settings as a new approach to objectively dissect T-NHL and to select more appropriate therapeutic protocols.