Depletion of tet2 results in age-dependent changes in DNA methylation and gene expression in a zebrafish model of myelodysplastic syndrome.

Introduction: Myelodysplastic syndrome (MDS) is a heterogeneous group of clonal hematopoietic disorders characterized by ineffective hematopoiesis, cytopenias, and dysplasia. The gene encoding ten-eleven translocation 2 (tet2), a dioxygenase enzyme that catalyzes the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine, is a recurrently mutated tumor suppressor gene in MDS and other myeloid malignancies. Previously, we reported a stable zebrafish line with a loss-of-function mutation in the tet2 gene. The tet2m/m-mutant zebrafish developed a pre-MDS state with kidney marrow dysplasia, but normal circulating blood counts by 11 months of age and accompanying anemia, signifying the onset of MDS, by 24 months of age. Methods: In the current study, we collected progenitor cells from the kidney marrows of the adult tet2m/m and tet2wt/wt fish at 4 and 15 months of age and conducted enhanced reduced representation of bisulfite sequencing (ERRBS) and bulk RNA-seq to measure changes in DNA methylation and gene expression of hematopoietic stem and progenitor cells (HSPCs). Results and discussion: A global increase in DNA methylation of gene promoter regions and CpG islands was observed in tet2m/m HSPCs at 4 months of age when compared with the wild type. Furthermore, hypermethylated genes were significantly enriched for targets of SUZ12 and the metal-response-element-binding transcription factor 2 (MTF2)-involved in the polycomb repressive complex 2 (PRC2). However, between 4 and 15 months of age, we observed a paradoxical global decrease in DNA methylation in tet2m/m HSPCs. Gene expression analyses identified upregulation of genes associated with mTORC1 signaling and interferon gamma and alpha responses in tet2m/m HSPCs at 4 months of age when compared with the wild type. Downregulated genes in HSPCs of tet2-mutant fish at 4 months of age were enriched for cell cycle regulation, heme metabolism, and interleukin 2 (IL2)/signal transducer and activator of transcription 5 (STAT5) signaling, possibly related to increased self-renewal and clonal advantage in HSPCs with tet2 loss of function. Finally, there was an overall inverse correlation between overall increased promoter methylation and gene expression.

Blocking UBE2N abrogates oncogenic immune signaling in acute myeloid leukemia.

Dysregulation of innate immune signaling pathways is implicated in various hematologic malignancies. However, these pathways have not been systematically examined in acute myeloid leukemia (AML). We report that AML hematopoietic stem and progenitor cells (HSPCs) exhibit a high frequency of dysregulated innate immune-related and inflammatory pathways, referred to as oncogenic immune signaling states. Through gene expression analyses and functional studies in human AML cell lines and patient-derived samples, we found that the ubiquitin-conjugating enzyme UBE2N is required for leukemic cell function in vitro and in vivo by maintaining oncogenic immune signaling states. It is known that the enzyme function of UBE2N can be inhibited by interfering with thioester formation between ubiquitin and the active site. We performed in silico structure-based and cellular-based screens and identified two related small-molecule inhibitors UC-764864/65 that targeted UBE2N at its active site. Using these small-molecule inhibitors as chemical probes, we further revealed the therapeutic efficacy of interfering with UBE2N function. This resulted in the blocking of ubiquitination of innate immune- and inflammatory-related substrates in human AML cell lines. Inhibition of UBE2N function disrupted oncogenic immune signaling by promoting cell death of leukemic HSPCs while sparing normal HSPCs in vitro. Moreover, baseline oncogenic immune signaling states in leukemic cells derived from discrete subsets of patients with AML exhibited a selective dependency on UBE2N function in vitro and in vivo. Our study reveals that interfering with UBE2N abrogates leukemic HSPC function and underscores the dependency of AML cells on UBE2N-dependent oncogenic immune signaling states.

Frequent somatic TET2 mutations in chronic NK-LGL leukemia with distinct patterns of cytopenias.

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher's exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch's t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

Targeted detection and quantitation of histone modifications from 1,000 cells.

Histone post-translational modifications (PTMs) create a powerful regulatory mechanism for maintaining chromosomal integrity in cells. Histone acetylation and methylation, the most widely studied histone PTMs, act in concert with chromatin-associated proteins to control access to genetic information during transcription. Alterations in cellular histone PTMs have been linked to disease states and have crucial biomarker and therapeutic potential. Traditional bottom-up mass spectrometry of histones requires large numbers of cells, typically one million or more. However, for some cell subtype-specific studies, it is difficult or impossible to obtain such large numbers of cells and quantification of rare histone PTMs is often unachievable. An established targeted LC-MS/MS method was used to quantify the abundance of histone PTMs from cell lines and primary human specimens. Sample preparation was modified by omitting nuclear isolation and reducing the rounds of histone derivatization to improve detection of histone peptides down to 1,000 cells. In the current study, we developed and validated a quantitative LC-MS/MS approach tailored for a targeted histone assay of 75 histone peptides with as few as 10,000 cells. Furthermore, we were able to detect and quantify 61 histone peptides from just 1,000 primary human stem cells. Detection of 37 histone peptides was possible from 1,000 acute myeloid leukemia patient cells. We anticipate that this revised method can be used in many applications where achieving large cell numbers is challenging, including rare human cell populations.

Splicing factor ESRP1 controls ER-positive breast cancer by altering metabolic pathways.

The epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) control the epithelial-to-mesenchymal transition (EMT) splicing program in cancer. However, their role in breast cancer recurrence is unclear. In this study, we report that high levels of ESRP1, but not ESRP2, are associated with poor prognosis in estrogen receptor positive (ER+) breast tumors. Knockdown of ESRP1 in endocrine-resistant breast cancer models decreases growth significantly and alters the EMT splicing signature, which we confirm using TCGA SpliceSeq data of ER+ BRCA tumors. However, these changes are not accompanied by the development of a mesenchymal phenotype or a change in key EMT-transcription factors. In tamoxifen-resistant cells, knockdown of ESRP1 affects lipid metabolism and oxidoreductase processes, resulting in the decreased expression of fatty acid synthase (FASN), stearoyl-CoA desaturase 1 (SCD1), and phosphoglycerate dehydrogenase (PHGDH) at both the mRNA and protein levels. Furthermore, ESRP1 knockdown increases the basal respiration and spare respiration capacity. This study reports a novel role for ESRP1 that could form the basis for the prevention of tamoxifen resistance in ER+ breast cancer.

Mutational landscape of RNA-binding proteins in human cancers.

RNA Binding Proteins (RBPs) are a class of post-transcriptional regulatory molecules which are increasingly documented to be dysfunctional in cancer genomes. However, our current understanding of these alterations is limited. Here, we delineate the mutational landscape of ∼1300 RBPs in ∼6000 cancer genomes. Our analysis revealed that RBPs have an average of ∼3 mutations per Mb across 26 cancer types. We identified 281 RBPs to be enriched for mutations (GEMs) in at least one cancer type. GEM RBPs were found to undergo frequent frameshift and inframe deletions as well as missense, nonsense and silent mutations when compared to those that are not enriched for mutations. Functional analysis of these RBPs revealed the enrichment of pathways associated with apoptosis, splicing and translation. Using the OncodriveFM framework, we also identified more than 200 candidate driver RBPs that were found to accumulate functionally impactful mutations in at least one cancer. Expression levels of 15% of these driver RBPs exhibited significant difference, when transcriptome groups with and without deleterious mutations were compared. Functional interaction network of the driver RBPs revealed the enrichment of spliceosomal machinery, suggesting a plausible mechanism for tumorogenesis while network analysis of the protein interactions between RBPs unambiguously revealed the higher degree, betweenness and closeness centrality for driver RBPs compared to non-drivers. Analysis to reveal cancer-specific Ribonucleoprotein (RNP) mutational hotspots showed extensive rewiring even among common drivers between cancer types. Knockdown experiments on pan-cancer drivers such as SF3B1 and PRPF8 in breast cancer cell lines, revealed cancer subtype specific functions like selective stem cell features, indicating a plausible means for RBPs to mediate cancer-specific phenotypes. Hence, this study would form a foundation to uncover the contribution of the mutational spectrum of RBPs in dysregulating the post-transcriptional regulatory networks in different cancer types.

Dissecting the expression relationships between RNA-binding proteins and their cognate targets in eukaryotic post-transcriptional regulatory networks.

RNA-binding proteins (RBPs) are pivotal in orchestrating several steps in the metabolism of RNA in eukaryotes thereby controlling an extensive network of RBP-RNA interactions. Here, we employed CLIP (cross-linking immunoprecipitation)-seq datasets for 60 human RBPs and RIP-ChIP (RNP immunoprecipitation-microarray) data for 69 yeast RBPs to construct a network of genome-wide RBP- target RNA interactions for each RBP. We show in humans that majority (~78%) of the RBPs are strongly associated with their target transcripts at transcript level while ~95% of the studied RBPs were also found to be strongly associated with expression levels of target transcripts when protein expression levels of RBPs were employed. At transcript level, RBP - RNA interaction data for the yeast genome, exhibited a strong association for 63% of the RBPs, confirming the association to be conserved across large phylogenetic distances. Analysis to uncover the features contributing to these associations revealed the number of target transcripts and length of the selected protein-coding transcript of an RBP at the transcript level while intensity of the CLIP signal, number of RNA-Binding domains, location of the binding site on the transcript, to be significant at the protein level. Our analysis will contribute to improved modelling and prediction of post-transcriptional networks.

Database of RNA binding protein expression and disease dynamics (READ DB).

RNA Binding Protein (RBP) Expression and Disease Dynamics database (READ DB) is a non-redundant, curated database of human RBPs. RBPs curated from different experimental studies are reported with their annotation, tissue-wide RNA and protein expression levels, evolutionary conservation, disease associations, protein-protein interactions, microRNA predictions, their known RNA recognition sequence motifs as well as predicted binding targets and associated functional themes, providing a one stop portal for understanding the expression, evolutionary trajectories and disease dynamics of RBPs in the context of post-transcriptional regulatory networks.

The human RBPome: from genes and proteins to human disease.

RNA binding proteins (RBPs) play a central role in mediating post transcriptional regulation of genes. However less is understood about them and their regulatory mechanisms. In this study, we construct a catalogue of 1344 experimentally confirmed RBPs. The domain architecture of RBPs enabled us to classify them into three groups - Classical (29%), Non-classical (19%) and unclassified (52%). A higher percentage of proteins with unclassified domains reveals the presence of various uncharacterised motifs that can potentially bind RNA. RBPs were found to be highly disordered compared to Non-RBPs (p<2.2e-16, Fisher's exact test), suggestive of a dynamic regulatory role of RBPs in cellular signalling and homeostasis. Evolutionary analysis in 62 different species showed that RBPs are highly conserved compared to Non-RBPs (p<2.2e-16, Wilcox-test), reflecting the conservation of various biological processes like mRNA splicing and ribosome biogenesis. The expression patterns of RBPs from human proteome map revealed that ~40% of them are ubiquitously expressed and ~60% are tissue-specific. RBPs were also seen to be highly associated with several neurological disorders, cancer and inflammatory diseases. Anatomical contexts like B cells, T-cells, foetal liver and foetal brain were found to be strongly enriched for RBPs, implying a prominent role of RBPs in immune responses and different developmental stages. The catalogue and meta-analysis presented here should form a foundation for furthering our understanding of RBPs and the cellular networks they control, in years to come. This article is part of a Special Issue entitled: Proteomics in India.

Uncovering RNA binding proteins associated with age and gender during liver maturation.

In the present study, we perform an association analysis focusing on the expression changes of 1344 RNA Binding proteins (RBPs) as a function of age and gender in human liver. We identify 88 and 45 RBPs to be significantly associated with age and gender respectively. Experimental verification of several of the predicted associations in mice confirmed our findings. Our results suggest that a small fraction of the gender-associated RBPs (~40%) are expressed higher in males than females. Altogether, these observations show that several of these RBPs are important and conserved regulators in maintaining liver function. Further analysis of the protein interaction network of RBPs associated with age and gender based on the centrality measures like degree, betweenness and closeness revealed that several of these RBPs might be prominent players in aging liver and impart gender specific alterations in gene expression via the formation of protein complexes. Indeed, both age and gender-associated RBPs in liver were found to show significantly higher clustering coefficients and network centrality measures compared to non-associated RBPs. The compendium of RBPs and this study will help us gain insight into the role of post-transcriptional regulatory molecules in aging and gender specific expression of genes.

Prognostic impact of HOTAIR expression is restricted to ER-negative breast cancers.

Expression of HOX transcript antisense intergenic RNA (HOTAIR), a large intergenic noncoding RNA (lincRNA), has been described as a metastases-associated lincRNA in various cancers including breast, liver and colon cancer cancers. We sought to determine if expression of HOTAIR could be used as a surrogate for assessing nodal metastases and evaluated RNA in situ hybridization (RNA-ISH) assay in a tissue microarray constructed from 133 breast cancer patients. The prognostic value of HOTAIR was further validated in large cohorts using The Cancer Genome Atlas (TCGA) breast cancer subjects. RNA-ISH analysis was successful in 94 cases (17% cases scored 0, 32.9% scored 1, 30.8% scored 2, and 19.1% scored 3). The expression of HOTAIR did not correlate with nodal metastasis regardless of the scoring intensity or with other study parameters (age, tumor size and grade, expression status). Further analysis of TCGA dataset showed that HOTAIR expression was lower in ductal carcinomas but higher in ER-negative tumors. Overexpression of HOTAIR was not associated with nodal metastases or prognosis in ER-positive patients. Its function as a poor prognostic indicator in ER-negative patients was restricted to node-positive patients. HOTAIR appears to be a marker for lymphatic metastases rather than hematogenous metastases in ER-negative patients.

Expression levels of SF3B3 correlate with prognosis and endocrine resistance in estrogen receptor-positive breast cancer.

De novo or acquired resistance to endocrine therapy limits its utility in a significant number of estrogen receptor-positive (ER-positive) breast cancers. It is crucial to identify novel targets for therapeutic intervention and improve the success of endocrine therapies. Splicing factor 3b, subunit 1 (SF3B1) mutations are described in luminal breast cancer albeit in low frequency. In this study, we evaluated the role of SF3B1 and SF3B3, critical parts of the SF3b splicing complex, in ER-positive endocrine resistance. To ascertain the role of SF3B1/SF3B3 in endocrine resistance, their expression levels were evaluated in ER-positive/endocrine-resistant cell lines (MCF-7/LCC2 and MCF-7/LCC9) using a real-time quantitative reverse transcription PCR (qRT-PCR). To further determine their clinical relevance, expression analysis was performed in a cohort of 60 paraffin-embedded ER-positive, node-negative breast carcinomas with low, intermediate, and high Oncotype DX recurrence scores. Expression levels of SF3B1 and SF3B3 and their prognostic value were validated in large cohorts using publicly available gene expression data sets including The Cancer Genome Atlas. SF3B1 and SF3B3 levels were significantly increased in ERα-positive cells with acquired tamoxifen (MCF-7/LCC2; both P<0.0002) and fulvestrant/tamoxifen resistance (MCF-7/LCC9; P=0.008 for SF3B1 and P=0.0006 for SF3B3). Expression levels of both MCF-7/LCC2 and MCF-7/LCC9 were not affected by additional treatments with E2 and/or tamoxifen. Furthermore, qRT-PCR analysis confirmed that SF3B3 expression is significantly upregulated in Oncotype DX high-risk groups when compared with low risk (P=0.019). Similarly, in publicly available breast cancer gene expression data sets, overexpression of SF3B3, but not SF3B1, was significantly correlated with overall survival. Furthermore, the correlation was significant in ER-positive, but not in ER-negative tumors.This is the first study to document the role of SF3B3 in endocrine resistance and prognosis in ER-positive breast cancer. Potential strategies for therapeutic targeting of the splicing mechanism(s) need to be evaluated.