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.

Abstract 3155: Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse

The purpose of this study is to identify key regulatory pathways that potentially drive abnormal gene expression program in relapsed Acute Myeloid Leukemia (AML) patients, by integrative computational analyses on multi-omics molecular profiles. Relapsed AML remains a clinical challenge. Epigenetic heterogeneity may contribute to transcriptional dysregulation and disease progression in AML. However, what specific transcriptional programs and potential regulatory mechanisms contribute to disease relapse are not yet well understood. To characterize the global transcriptional landscapes in relapsed AML, we integrated genomics data from two cohorts of matched diagnosis and relapse patient specimens. We identified 5,416 differentially expressed genes (DEGs) between diagnosis and relapse in Cohort I. Unsupervised clustering yielded three distinct DEG groups: group A, B and C genes that were predominantly (88%) down-regulated, divergently regulated, or predominantly (65%) up-regulated, respectively, upon relapse. The expression pattern of all DEGs separated the patients into two clusters, most robustly by Group B genes. Interestingly, the majority of DEGs did not associate with changes in gene promoter methylation. Similar patterns were observed in Cohort II. We used Binding Analysis for Regulation of Transcription (BART) to identify transcriptional regulators (TRs) that potentially regulated the DEGs not associated with DNA methylation changes, and assessed the differential expression of identified TRs during disease progression. PU.1 was identified as a potential TR for Group A genes and was down-regulated upon relapse. GATA1 and TAL1 were identified as regulating Group B genes and were up-regulated in patient cluster1 and down-regulated in cluster2, consistent with the expression pattern of Group B genes. RBBP5 was a top predicted TR for Group C genes and was up-regulated upon relapse. We next validated the potential functionality of those predicted factors. In NSG mice transplanted with a human AML specimen, TAL1 and GATA1 were downregulated in AML cells collected four weeks after chemotherapy treatment, and were inferred as TRs for the down-regulated genes, similar to the patient data. PU.1 was inferred as regulating the up-regulated genes. Furthermore, we found that the level of differential expression of TAL1, GATA1, and PU.1 in each patient specimen associated with the correlation of DEG profiles between the patient specimen and TR perturbation in human-derived hematopoietic cell lines. Our results support the possibility that in some AML patients, TRs with roles in hematopoiesis and leukemia might contribute to disease relapse. Further mechanistic studies deciphering the molecular and phenotypic events facilitated by these TRs will yield significant insight into disease biology and possible therapeutic targeting approaches in relapsed AML.

Citation Format: Zhenjia Wang, Yaseswini Neelamraju, Cem Meydan, Nicholas Dunham, Jorge Gandara, Tak Lee, Subhash Prajapati, Franck Rapaport, Caroline Sheridan, Paul Zumbo, Michael Becker, Lars Bullinger, Martin Carroll, Richard D’Andrea, Richard Dillon, Ross Levine, Christopher E. Mason, Ari Melnick, Donna Neuberg, Stefan Bekiranov, Chongzhi Zang, Francine E. Garrett-Bakelman. Gene expression profiles reveal distinct regulatory activities of transcription factors GATA1 and TAL1 upon AML relapse [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3155.

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.

Abstract P5-04-03: Aggressiveness of epithelial cancers is independent of epithelial-to-mesenchymal transition

Background: Epithelial-to-Mesenchymal Transition (EMT)is postulated to be an important step in cancer progression and controlled by multiple mechanisms including EMT transcription factors (EMT-TFs) and splicing factors such as Epithelial Splicing Regulatory Proteins (ESRP1 and ESRP2). We previously have shown that the expression of ESRP1 and ESRP2 have significantly elevated in cases with high Oncotype DX scores and in ERα-positive cells with acquired endocrine resistance (SABCS 2013). This study seeks to identify the role of EMT-TFs and ESRP1s in the determination of outcomes of patients with ER+ breast cancer.

Patients and Methods: The expression of EMT-TFs and ESRP1 was analyzed in the Affymetrix microarray and TCGA BRCA databases. Next, we generated genetically engineered knockdown models of ESRP1 to understand its functional role in endocrine resistance. We performed RNA-seq and MATS analysis to identify alternative splicing events (ASEs) between ESRP1 knockdown and control cell lines [(2C3 vs 2-control (LCC2 set) and 9C2 vs 9-control(LCC9 set)]. Validation of the ASEs was performed using a probe-based platform [Human Transcriptome Array 2.0 (HTA)] and TCGA SpliceSeq from breast tumors.

Results: High levels of ESRP1 mRNA, but not EMT-TFs, are associated with poor prognosis in human ER+ breast tumors (Affymetrix; P=2.8e-07 and TCGA; P=0.00011). Knockdown of ESRP1 in ER+ endocrine resistant breast cancer induced glandular differentiation, rather than mesenchymal features. This was associated with significant reduction in cell and tumor growth in mammary fat pad orthotopic xenograft mice models of LCC2 and LCC9. No alterations in EMT-TFs were observed in these cells. Transcriptome profiling of ESRP1 knockdown cells further revealed altered ASEs in EMT splicing gene signature, but not at the gene level. These alterations (SE-skipped exon) were further validated for ARHGEF11, ENAH, FNIP1, SCRIB, and SLK using probe based HTA platform for ESRP1 knockdown cells and TCGA-SpliceSeq ER+ BRCA tumors in ER+ ESRP1low versus ESRP1high breast tumors.

Conclusions: Our data demonstrates for the first time that high ESRP1 is associated with poor prognosis in ER+ breast cancer. Despite its involvement in regulation of EMT splicing signature, low levels (or knockdown) of ESRP1 were not associated with EMT phenotype in tumors or in endocrine-resistant ER+ cells. Taken together, our findings indicate that EMT is not important in determining prognosis in ER+ breast cancer and that ESRP1 exerts a different role in aggressive ER+ breast cancers.

Citation Format: Badve SS, Neelamraju Y, Goswami CP, Gu X, Nallamothu G, Gu Y, Vieth E, Janga SC, Ryan M, Gokmen-Polar Y. Aggressiveness of epithelial cancers is independent of epithelial-to-mesenchymal transition [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr P5-04-03.

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.

Abstract 1097: Knockdown of splicing factor ESRP1 affects multiple splicing factors

Introduction: Tissue-specific alternative splicing (AS) is an important mechanism for regulating gene expression. Epithelial Splicing Regulatory Proteins 1 and 2 (ESRP1 and ESRP2), two RNA-binding proteins (RBPs) that promote splicing, are potential candidates that contribute to breast cancer recurrence and resistance to therapies. In our prior studies, we have shown that ESRPs are associated with endocrine resistance. In this study, we seek to investigate the impact of the ESRP1 knockdown in endocrine resistant breast cancer cells. Methods: Expression of ESRP1 was analyzed in endocrine resistant-LCC2 and LCC9 cells. Knockdown of ESRP1 expression in LCC2 and LCC9 cells was achieved using lenti-viral based shRNA approach (Mission TRC human shRNA constructs, Sigma). We next performed paired-end RNA-sequencing of LCC2 and LCC9 cells transfected with control and two clones representative of shRNA knockdown (Illumina HiSeq platform, SeqWright Genomic Services) according to the manufacturer's instructions. After processing and quality check of raw FASTQ data (Illumina CASAVA, FASTQ and FASTX toolkits), reads were aligned, and quantified using bioinformatics software tools (Tophat/Cufflinks and Sailfish). Results: Expression of ESRP1 is significantly elevated in ERα-positive endocrine -resistant cells (LCC2 and LCC9 cells). ESRP1 knockdown was confirmed using qRT-PCR analysis in LCC2 cell lines (clone1 - 80% reduction and clone 3- 95% reduction) and in LCC9 cells (clone 2- 90% reduction and clone 3- 50% reduction). RNA-seq analysis identified 5153 differentially regulated genes between ESRP1 knockdown and control LCC2 cell lines. Of the 5153 genes, 1117 and 2109 genes were unique to clone 1 and clone 3, respectively, while 1927 genes were common in both clones. Targeted RNA-seq analysis also revealed that ESRP1 regulates the expression of several genes including other splicing factors. Furthermore, a functional enrichment analysis of these genes revealed a possible role of altering important pathways such as cell-cell adhesion and cell signaling. Additionally, it was observed that neither ESR1 nor any of its downstream genes were affected by the knockdown of ESRP1. This data was similar to that obtained in knockdown of LCC9 cells. Conclusions: ESRP1 regulates splicing of multiple genes in ER+ breast cancer. The process is independent of ER pathway. Targeting splicing could be a novel modality for treating endocrine therapy resistant breast cancer.

Citation Format: Yesim Gökmen-Polar, Yaseswini Neelamraju, Xiaoping Gu, Gouthami Nallamothu, Sarath C. Janga, Sunil Badve. Knockdown of splicing factor ESRP1 affects multiple splicing factors. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1097. doi:10.1158/1538-7445.AM2015-1097

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.

Abstract P3-05-20: ESRP1 adds sp(l)ice to endocrine resistance

Introduction:

De novo or acquired resistance to endocrine therapy limits its utility in a significant number of estrogen receptor (ER) breast cancers. An increasing number of molecular assays predict the likelihood of distant recurrence in tamoxifen-treated patients with node-negative, ER+ breast cancer. However, these do not provide the mechanistic basis for endocrine resistance. It is crucial to identify novel targets and improve the success of endocrine therapies. We previously shown that epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2), RNA binding proteins that promote splicing, are significantly elevated in cases with high Oncotype DX scores (innate resistance) and in ERα-positive cells with acquired tamoxifen resistance (MCF-7 LCC2 cells) and fulvestrant and tamoxifen resistance (MCF-7/LCC9 cells). The aim of this study was to investigate the ESRP1/ESRP2-regulated alternative splicing events leading to innate and acquired endocrine resistance.

Methods:

A combinatorial bioinformatics approach was employed to identify genes altered through alternative splicing by ESRP1/ESRP2 in endocrine resistance. Briefly, genes with ESRP1/ESRP2 motifs were screened in the human genome. This data was integrated with tamoxifen-treated datasets, and filtered using protein-protein interactions (PPI; BIOGRID) and clinical outcome (KM plotter. Potential target gene transcripts were further narrowed down using an algorithm based on motif binding location and the Cancer Genome Atlas (TCGA) breast cancer datasets with low and high ESRP1 cases. The alternative transcripts were further validated using splice variant specific-custom qRT-PCR in low and high RS Oncotype cases as surrogate for endocrine sensitivity and resistance.

Results:

Motif analysis in combination with tamoxifen-treated datasets identified 2212 differentially expressed genes. Further collective analysis of number of PPI (&gt;50) and survival data from KM plotter (P&lt;0.001) narrowed down 79 candidate genes that are associated with tamoxifen resistance. TCGA data confirmed presence of distinct transcripts (splicing variants) based on ESPR1 expression levels. Splice specific RT-PCR confirmed alternative splicing events involving cycle related genes such as AURKA, FZR1, and MDM2 in low ESRP1 (and low RS) and high ESRP1 (high RS) cases.

Conclusion:

ESRP1/ESRP2 by inducing alternative splicing play an important role in tamoxifen resistance and recurrence of ER+ breast cancer. Targeting alternative splicing may offer novel avenues for combating endocrine-resistance in breast cancer.

Citation Format: Yesim Gokmen-Polar, Yaseswini Neelamraju, Chirayu Pankaj Goswami, Harikrishna Nakshatri, Sarath Chandra Janga, Sunil Badve. ESRP1 adds sp(l)ice to endocrine resistance [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P3-05-20.

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.