Functional dissection of HOXD cluster genes in regulation of neuroblastoma cell proliferation and differentiation

Transcriptome Profiling of Mouse Embryonic Fibroblast Spontaneous Immortalization: A Comparative Analysis

Cell immortalization, a hallmark of cancer development, is a process that cells can undergo on their path to carcinogenesis. Spontaneously immortalized mouse embryonic fibroblasts (MEFs) have been used for decades; however, changes in the global transcriptome during this process have been poorly described. In our research, we characterized the poly-A RNA transcriptome changes after spontaneous immortalization. To this end, differentially expressed genes (DEGs) were screened using DESeq2 and characterized by gene ontology enrichment analysis and protein-protein interaction (PPI) network analysis to identify the potential hub genes. In our study, we identified changes in the expression of genes involved in proliferation regulation, cell adhesion, immune response and transcriptional regulation in immortalized MEFs. In addition, we performed a comparative analysis with previously reported MEF immortalization data, where we propose a predicted gene regulatory network model in immortalized MEFs based on the altered expression of Mapk11, Cdh1, Chl1, Zic1, Hoxd10 and the novel hub genes Il6 and Itgb2.

Dynamic regulation of BDNF gene expression by estradiol and lncRNA HOTAIR

Brain derived neurotrophic factor (BDNF) is a major neurotransmitter that controls growth and maintenance of neurons and its misregulation is linked to neurodegeneration and human diseases. Estradiol (E2) is well-known to regulate the process of differentiation and plasticity of hippocampal neurons. Here we examined the mechanisms of BDNF gene regulation under basal conditions and under stimuli such as E2. Our results demonstrated that BDNF expression is induced by E2 in vitro in HT22 cells (hippocampal neuronal cells) and in vivo (in ovariectomized mouse brain under E2-treatment). Using chromatin immunoprecipitation assay, we demonstrated that estrogen receptors (ERα, ERβ) were enriched at the BDNF promoter in presence of E2. Additionally, ER-coregulators (e.g., CBP/p300, MLL3), histone acetylation, H3K4-trimethylation, and RNA polymerase II levels were also elevated at the BDNF promoter in an E2-dependent manner. Additionally, under the basal conditions (in the absence of E2), the long noncoding RNA HOTAIR and its interacting partners PRC2 and LSD1 complexes binds to the promoter of BDNF and represses its expression. HOTAIR knockdown -relieves the repression resulting in elevation of BDNF expression. Further, levels of HOTAIR-interacting partners, EZH2 and LSD1 were reduced at the BDNF promoter upon HOTAIR-knockdown revealing that HOTAIR plays a regulatory role in BDNF gene expression by modulating promoter histone modifications. Additionally, we showed that E2 induced-BDNF expression is mediated by the displacement of silencing factors, EZH2 and LSD1 at BDNF promoter and subsequent recruitment of active transcription machinery. These results reveal the mechanisms of BDNF gene regulation under the basal condition and in presence of a positive regulator such as E2 in neuronal cells.

A pan-cancer analysis of the role of HOXD1, HOXD3, and HOXD4 and validation in renal cell carcinoma

HOXD1, HOXD3, and HOXD4 are members of the HOXD genes family and are related to tumorigenesis of the tumor. However, whether HOXDs (1, 3, 4) have a crucial role across pan-cancer is still unknown. HOXD1, HOXD3, and HOXD4 expressions were analyzed using public databases in 33 types of tumors. The UCSC Xena website was carried out to investigate the relationship between the expression of genes and the progress of cancers. The biological functions of HOXD3 were tested by colony forming, transwell, wound healing, and xenograft assay in vitro and in vivo. GSEA was used to identify the associated cancer hallmarks with HOXDs expression. Immune cell infiltration analysis was applied to verify the immune cell infiltrations related to genes. The results showed HOXD1, HOXD3, and HOXD4 co-low expressed in BRCA, COAD, KICH, KIRC, KIRP, READ, and TGCT. In the KIRC, all of HOXDs expression was connected with tumor stage and histological grade. Upregulation of HOXDs was associated with improved OS, DSS, and PFI. Down-expression of HOXD3 induced cell proliferation, migration, and invasion in vivo and in vitro. In addition, HOXDs were connected with immune-activated hallmarks and cancer immune cell infiltrations. These findings demonstrated that HOXDs may be indicative biomarkers for the prognosis and immunotherapy in pan-cancer.

Chromosome 1p status in neuroblastoma correlates with higher expression levels of miRNAs targeting neuronal differentiation pathway

Neuroblastoma (NB) is characterized by acquired segmental and numerical chromosome aberrations. Although deletions of distal 1p and 11q are frequent alterations, no candidate tumor suppressor gene residing in these chromosomal sites could be identified so far. In the present study, we detected the genomic imbalances of six neuroblastoma cell lines using the multiplex ligation-dependent probe amplification (MLPA) technique and the microRNA (miRNA) expression profiles of the cell lines by a microarray study. According to MLPA results, we aimed to assess the miRNA expression profiles of the cell lines harboring 11q and 1p deletions. The cell lines with 1p deletions revealed statistically significant higher levels of expression for 29 miRNAs in contrast to the cell lines without 1p deletion in microarray study. We also performed GO enrichment analysis for predicted targets of the differentially expressed miRNAs. According to GO enrichment analysis, miRNAs that showed the high change in expression was associated with neuronal differentiation. We showed that hsa-miR-494, hsa-miR-495, and hsa-miR-543 target most of mRNAs in neuronal differentiation pathway. Although limited to the cell lines, our results highly suggest that NBs with different segmental chromosome abnormalities may have different dysregulated miRNA expression signatures that target the genes involved in neuronal differentiation.

Significance of HOXD transcription factors family in progression, migration and angiogenesis of cancer

The transcription factors (TFs) of the HOX family play significant roles during early embryonic development and cellular processes. They also play a key role in tumorigenesis as tumor oncogenes or suppressors. Furthermore, TFs of the HOXD geFIne cluster affect proliferation, migration, and invasion of tumors. Consequently, dysregulated activity of HOXD TFs has been linked to clinicopathological characteristics of cancer. HOXD TFs are regulated by non-coding RNAs and methylation of DNA on promoter and enhancer regions. In addition, HOXD genes modulate the biological function of cancer cells via the MEK and AKT signaling pathways, thus, making HOXD TFs, a suitable molecular marker for cancer prognosis and therapy. In this review, we summarized the roles of HOXD TFs in different cancers and highlighted its potential as a diagnostic and therapeutic target.

TOP2B Is Required to Maintain the Adrenergic Neural Phenotype and for ATRA-Induced Differentiation of SH-SY5Y Neuroblastoma Cells

The neuroblastoma cell line SH-SY5Y is widely used to study retinoic acid (RA)-induced gene expression and differentiation and as a tool to study neurodegenerative disorders. SH-SY5Y cells predominantly exhibit adrenergic neuronal properties, but they can also exist in an epigenetically interconvertible alternative state with more mesenchymal characteristics; as a result, these cells can be used to study gene regulation circuitry controlling neuroblastoma phenotype. Using a combination of pharmacological inhibition and targeted gene inactivation, we have probed the requirement for DNA topoisomerase IIB (TOP2B) in RA-induced gene expression and differentiation and in the balance between adrenergic neuronal versus mesenchymal transcription programmes. We found that expression of many, but not all genes that are rapidly induced by ATRA in SH-SY5Y cells was significantly reduced in the TOP2B null cells; these genes include BCL2, CYP26A1, CRABP2, and NTRK2. Comparing gene expression profiles in wild-type versus TOP2B null cells, we found that long genes and genes expressed at a high level in WT SH-SY5Y cells were disproportionately dependent on TOP2B. Notably, TOP2B null SH-SY5Y cells upregulated mesenchymal markers vimentin (VIM) and fibronectin (FN1) and components of the NOTCH signalling pathway. Enrichment analysis and comparison with the transcription profiles of other neuroblastoma-derived cell lines supported the conclusion that TOP2B is required to fully maintain the adrenergic neural-like transcriptional signature of SH-SY5Y cells and to suppress the alternative mesenchymal epithelial-like epigenetic state.

Prognostic significance of HOXD4 protein expression in human ovarian cancers

Objectives

Ovarian cancer is the most common gynecological malignancy, ranking as the fifth leading cause of cancer-related deaths among females in the United States. Homeobox D4 (HOXD4) is a transcription factor belonging to the homeobox protein family, which plays a critical role in morphogenesis during embryo development. Here we aimed to study the HOXD4 expression in ovarian serous carcinoma (OSC) and determine its clinical significance.

Materials and Methods

Real-time quantitative PCR and immunohistochemistry targeting human OSC tissues and adjacent ovarian tissues were performed to correlate the patterns of HOXD4 expression with clinical characteristics and survival outcomes. Cell lines and nude mice were used for verifying the role of HOXD4 in OSC.

Results

HOXD4 protein was predominantly expressed in OSC tissues compared with nontumorous tissues. The correlation test demonstrated a significant correlation between HOXD4 with tumor FIGO stage. Univariate and multivariate analyses found that HOXD4 expression was associated with poorer overall survival. Furthermore, high expression of HOXD4 protein was observed in OSC cell lines in vitro. Finally, the oncogenic effect of HOXD4 was confirmed by cellular and xenograft experiments.

Conclusion

HOXD4 protein expression may be associated with a poorer prognosis in OSC. The unfavorable prognostic value of HOXD4 in malignancies and its underlying mechanism are worthy of further investigation.

Integrated analysis identifies a long non-coding RNAs-messenger RNAs signature for prediction of prognosis in hepatitis B virus-hepatocellular carcinoma patients

Hepatitis B virus (HBV) infection is a leading cause of hepatocellular carcinoma (HCC), but HBV-HCC related prognosis signature remains rarely investigated. This study was to identify an integrated long non-coding RNAs-messenger RNAs (lncRNA-mRNA) signature for prediction of overall survival (OS) and explore their underlying functions.One RNA-sequencing dataset (training set, n = 95) and one microarray dataset E-TABM-36 (validation set, n = 44) were collected. Least absolute shrinkage and selection operator analysis was performed to identify an lncRNA-mRNA prognosis signature. The OS difference of patients in the high-risk and low-risk risk groups was evaluated by Kaplan-Meier curve. Area under the receiver operating characteristic curve (AUC), Harrell concordance index (C-index) calculation, and multivariate analyses with clinical characteristics were used to determine the prognostic ability. Furthermore, a coexpression network was constructed to interpret the functions.Nine signature genes (3 lncRNAs and 6 mRNAs) were selected to generate the risk score model. Patients belonging to the high-risk group showed a significantly shorter survival than those of the low-risk group. The prediction accuracy of the risk score for 5-year OS was 0.936 and 0.905 for the training set and validation set, respectively. Also, this risk score was independent of various clinical variables for the prognosis prediction. Incorporation of the risk score remarkably increased the predictive power of the routine clinical prognostic factors (vascular invasion status, tumor recurrence status) (AUC = 0.942 vs 0.628; C-index = 0.7997 vs 0.6908). Furthermore, LncRNA insulin-like growth factor 2 antisense RNA (IGF2-AS) and long intergenic non-protein coding RNA 342 (LINC00342) were predicted to exert tumor suppression effects by regulating homeobox D1 (HOXD1) and secreted frizzled related protein 5 (SFRP5), respectively; while lncRNA rhophilin Rho GTPase binding protein 1 antisense RNA 1 (RHPN1-AS1) may possess carcinogenic potential by promoting the transcription of chromobox 2 (CBX2), cell division cycle 20 (CDC20), matrix metallopeptidase 12 (MMP12), stratifin (SFN), tripartite motif containing 16 (TRIM16), and uroplakin 3A (UPK3A). These mRNAs may be associated with cell proliferation or apoptosis related pathways.This study may provide a novel, effective prognostic biomarker, and some therapeutic targets for HBV-HCC patients.

The rs12826786 in HOTAIR lncRNA Is Associated with Risk of Autism Spectrum Disorder

Autism spectrum disorder (ASD) encompasses a group of neurodevelopmental disorders with complex pattern of inheritance. Several single nucleotide polymorphisms (SNPs) within coding or non-coding regions of genome have been associated with risk of this disorder. In the current study, we genotyped rs12826786, rs1899663, and rs4759314 SNPs within HOX transcript antisense RNA (HOTAIR) in 427 ASD cases and 430 normally developed children. The rs12826786 was associated with ASD in allelic (T vs. C: OR (95% CI) = 1.29 (1.07-1.57), adjusted P value = 0.03) and recessive (TT vs. TC + CC: OR (95% CI) = 1.60 (1.10-2.32), adjusted P value = 0.04) models. However, the other SNPs were not associated with ASD in any inheritance model. No estimated haplotype within HOTAIR was associated with risk of ASD in the assessed population. Based on the results of the current investigation, the rs12826786 can be regarded as a risk locus for ASD in Iranian population.

Histone demethylase KDM6B has an anti-tumorigenic function in neuroblastoma by promoting differentiation

Induction of differentiation is a therapeutic strategy in high-risk neuroblastoma, a childhood cancer of the sympathetic nervous system. Neuroblastoma differentiation requires transcriptional upregulation of neuronal genes. How this process is regulated at epigenetic levels is not well understood. Here we report that the histone H3 lysine 27 demethylase KDM6B is an epigenetic activator of neuroblastoma cell differentiation. KDM6B mRNA expression is downregulated in poorly differentiated high-risk neuroblastomas and upregulated in differentiated tumors, and high KDM6B expression is prognostic for better survival in neuroblastoma patients. In neuroblastoma cell lines, KDM6B depletion promotes cell proliferation, whereas KDM6B overexpression induces neuronal differentiation and inhibits cell proliferation and tumorgenicity. Mechanistically, KDM6B epigenetically activates the transcription of neuronal genes by removing the repressive chromatin marker histone H3 lysine 27 trimethylation. In addition, we show that KDM6B functions downstream of the retinoic acid-HOXC9 axis in inducing neuroblastoma cell differentiation: KDM6B expression is upregulated by retinoic acid via HOXC9, and KDM6B is required for HOXC9-induced neuroblastoma cell differentiation. Finally, we present evidence that KDM6B interacts with HOXC9 to target neuronal genes for epigenetic activation. These findings identify a KDM6B-dependent epigenetic mechanism in the control of neuroblastoma cell differentiation, providing a rationale for reducing histone H3 lysine 27 trimethylation as a strategy for enhancing differentiation-based therapy in high-risk neuroblastoma.

High HOXD4 protein expression in gastric adenocarcinoma tissues indicates unfavorable clinical outcomes

BACKGROUND/AIM

Homeobox D4 (HOXD4) belongs to the homeobox (HOX) family, which plays a crucial role in the early embryo development and cell differentiation. The role of HOXD4 in human gastric adenocarcinoma has not been elucidated. In the present study, we aimed to examine the expression levels of HOXD4 and dissect whether the HOXD4 expression is associated with aggressive clinicopathological outcomes of patients with gastric adenocarcinoma.

PATIENTS AND METHODS

Clinicopathological analyses were performed in 127 patients with gastric adenocarcinoma. Expression of HOXD4 was tested by immunohistochemistry staining and quantitative RT-PCR. Clinical outcomes were evaluated by the Kaplan-Meier method and log-rank test. The prognostic role of HOXD4 in gastric adenocarcinoma patients was assessed by univariate and multivariate analyses. The effects and mechanisms of HOXD4 on cell proliferation, migration and invasion were explored through cellular experiments.

RESULTS

HOXD4 expression was elevated in gastric adenocarcinoma tissues compared to non-tumorous gastric tissues (P = 0.018). High expression of HOXD4 was significantly associated with larger tumor size (P = 0.008), advanced tumor invasion depth (P = 0.014), and positive lymph node metastasis (P < 0.001). Moreover, patients with high HOXD4 expression had poorer overall survival (P = 0.001), and HOXD4 was identified as an independent prognosis factor according to multivariate analysis [hazard ratio (HR) =2.253, 95% confident interval (CI) 1.028-4.979, P = 0.044]. Cellular results revealed that HOXD4 can promote tumor cell proliferation by upregulating c-Myc and cyclin D1.

CONCLUSIONS

Our study demonstrated that overexpression of HOXD4 was significantly correlated with poorer prognosis of gastric adenocarcinoma patients, indicating the potential of HOXD4 as a novel clinical predictive biomarker and drug target.

Identification of a Robust Methylation Classifier for Cutaneous Melanoma Diagnosis

Early diagnosis improves melanoma survival, yet the histopathological diagnosis of cutaneous primary melanoma can be challenging, even for expert dermatopathologists. Analysis of epigenetic alterations, such as DNA methylation, that occur in melanoma can aid in its early diagnosis. Using a genome-wide methylation screening, we assessed CpG methylation in a diverse set of 89 primary invasive melanomas, 73 nevi, and 41 melanocytic proliferations of uncertain malignant potential, classified based on interobserver review by dermatopathologists. Melanomas and nevi were split into training and validation sets. Predictive modeling in the training set using ElasticNet identified a 40-CpG classifier distinguishing 60 melanomas from 48 nevi. High diagnostic accuracy (area under the receiver operator characteristic curve = 0.996, sensitivity = 96.6%, and specificity = 100.0%) was independently confirmed in the validation set (29 melanomas, 25 nevi) and other published sample sets. The 40-CpG melanoma classifier included homeobox transcription factors and genes with roles in stem cell pluripotency or the nervous system. Application of the 40-CpG melanoma classifier to the diagnostically uncertain samples assigned melanoma or nevus status, potentially offering a diagnostic tool to assist dermatopathologists. In summary, the robust, accurate 40-CpG melanoma classifier offers a promising assay for improving primary melanoma diagnosis.

Role of HOX Genes in Stem Cell Differentiation and Cancer

HOX genes encode an evolutionarily conserved set of transcription factors that control how the phenotype of an organism becomes organized during development based on its genetic makeup. For example, in bilaterian-type animals, HOX genes are organized in gene clusters that encode anatomic segment identity, that is, whether the embryo will form with bilateral symmetry with a head (anterior), tail (posterior), back (dorsal), and belly (ventral). Although HOX genes are known to regulate stem cell (SC) differentiation and HOX genes are dysregulated in cancer, the mechanisms by which dysregulation of HOX genes in SCs causes cancer development is not fully understood. Therefore, the purpose of this manuscript was (i) to review the role of HOX genes in SC differentiation, particularly in embryonic, adult tissue-specific, and induced pluripotent SC, and (ii) to investigate how dysregulated HOX genes in SCs are responsible for the development of colorectal cancer (CRC) and acute myeloid leukemia (AML). We analyzed HOX gene expression in CRC and AML using information from The Cancer Genome Atlas study. Finally, we reviewed the literature on HOX genes and related therapeutics that might help us understand ways to develop SC-specific therapies that target aberrant HOX gene expression that contributes to cancer development.

HOTAIR gene polymorphisms contribute to increased neuroblastoma susceptibility in Chinese children

BACKGROUND

Neuroblastoma is the most frequently diagnosed extracranial solid tumor in children. Previous studies have shown that single-nucleotide polymorphisms in some genes are associated with the risk of multiple cancers, including neuroblastoma. Although Hox transcript antisense intergenic RNA (HOTAIR) gene polymorphisms have been investigated in a variety of cancers, to the authors' knowledge the relationships between HOTAIR gene polymorphisms and neuroblastoma susceptibility have not been reported to date. The objective of the current study was to evaluate the correlation between HOTAIR gene polymorphisms and neuroblastoma risk in Chinese children.

METHODS

The authors genotyped 6 polymorphisms (rs920778 A>G, rs12826786 C>T, rs4759314 A>G, rs7958904 G>C, rs874945 C>T, and rs1899663 C>A) of the HOTAIR gene in 2 Chinese populations including 393 neuroblastoma cases and 812 healthy controls. The strength of the associations was evaluated using odds ratios and 95% confidence intervals. Further stratification analyses were conducted to explore the association between the HOTAIR gene polymorphisms rs12826786 C>T, rs874945 C>T, and rs1899663 C>A with neuroblastoma susceptibility in terms of age, sex, clinical stage of disease, and sites of origin.

RESULTS

The authors found that the rs12826786 C>T (P =.013), rs874945 C>T (P =.020), and rs1899663 C>A (P =.029) polymorphisms were significantly associated with increased neuroblastoma risk. In stratification analyses, these associations were more predominant in females and among patients with tumor in the retroperitoneal region or mediastinum. The remaining 3 polymorphisms were not found to be related to neuroblastoma susceptibility.

CONCLUSIONS

The results of the current study verified that HOTAIR gene polymorphisms are associated with increased neuroblastoma risk and suggest that HOTAIR gene polymorphisms might be a potential biomarker for neuroblastoma susceptibility. Cancer 2018;124:2599-606. © 2018 American Cancer Society.

Clinicopathological analysis of HOXD4 expression in diffuse gliomas and its correlation with IDH mutations and 1p/19q co-deletion

Backgrounds

HOX (homologous box) is known as the dominant gene of vertebrate growth and cell differentiation. Abnormal expression of HOX gene in various tumors has attracted the attention of scholars. As a component of HOX clusters, HOXD4 plays a controversial role in the tumorigenesis of central nervous system.

Results

The data demonstrated that and the results demonstrated that HOXD4 was overexpressed in glioma tissues compared to that of normal brain tissues. patients with high HOXD4 expression had a significant shorter survival than those with low HOXD4 expression in total glioma cohort (p<0.001), WHO Grade II cohort (p=0.003) and Grade III cohort (p<0.001), but not in Grade IV cohort when OS (overall survival) was analyzed (p=0.216). The findings were confirmed by the large-scale omics data analysis including lower-grade glioma (LGG) and glioblastoma multiforme (GBM) in TCGA (the cancer genome atlas) and CGGA (Chinese glioma genome atlas). Moreover, it was revealed that the expression of HOXD4 have a significant impact on the OS of Grade IV glioma with IDH wild-type and 1p/19q intact according to TCGA data.

Methods

Clinicopathological analysis of HOXD4 expression in 453 glioma patients was performed in the current study. Expression of HOXD4 was evaluated by qPCR and immunohistochemical (IHC) staining. Univariate and multivariate analysis were conducted to investigate the prognostic role of HOXD4 in glioma patients.

Conclusions

Expression of HOXD4 was closely related to the clinical outcomes of patients with gliomas, and HOXD4 may be a potential prognostic biomarker of gliomas.

Alternative SET/TAFI Promoters Regulate Embryonic Stem Cell Differentiation

Embryonic stem cells (ESCs) are regulated by pluripotency-related transcription factors in concert with chromatin regulators. To identify additional stem cell regulators, we screened a library of endogenously labeled fluorescent fusion proteins in mouse ESCs for fluorescence loss during differentiation. We identified SET, which displayed a rapid isoform shift during early differentiation from the predominant isoform in ESCs, SETα, to the primary isoform in differentiated cells, SETβ, through alternative promoters. SETα is selectively bound and regulated by pluripotency factors. SET depletion causes proliferation slowdown and perturbed neuronal differentiation in vitro and developmental arrest in vivo, and photobleaching methods demonstrate SET's role in maintaining a dynamic chromatin state in ESCs. This work identifies an important regulator of pluripotency and early differentiation, which is controlled by alternative promoter usage.

•

We identify SETα to be rapidly downregulated during ESC differentiation

•

SETα is regulated by pluripotency factors and replaced by SETβ during differentiation

•

SETα/SETβ switch is crucial for ESC differentiation

•

SETα regulates chromatin plasticity in ESCs

HOXD8 exerts a tumor-suppressing role in colorectal cancer as an apoptotic inducer

Homeobox (HOX) genes are conserved transcription factors which determine the anterior-posterior body axis patterning. HOXD8 is a member of HOX genes deregulated in several tumors such as lung carcinoma, neuroblastoma, glioma and colorectal cancer (CRC) in a context-dependent manner. In CRC, HOXD8 is downregulated in cancer tissues and metastatic foci as compared to normal tissues. Whether HOXD8 acts as a tumor suppressor of malignant progression and metastasis is still unclear. Also, the underlying mechanism of its function including the downstream targets is totally unknown. Here, we clarified the lower expression of HOXD8 in clinical colorectal cancer vs. normal colon tissues. Also, we showed that stable expression of HOXD8 in colorectal cancer cells significantly reduced the cell proliferation, anchorage-independent growth and invasion. Further, using The Cancer Genome Atlas (TCGA), we identified the genes associated with HOXD8 in order to demonstrate its function as a suppressor or a promoter of colorectal carcinoma. Among inversely related genes, apoptotic inhibitors like STK38 kinase and MYC were shown to be negatively associated with HOXD8. We demonstrated the ability of HOXD8 to upregulate executioner caspases 6 & 7 and cleaved PARP, thus inducing the apoptotic events in colorectal cancer cells.

DeepCAGE transcriptomics identify HOXD10 as a transcription factor regulating lymphatic endothelial responses to VEGF-C

Lymphangiogenesis plays a crucial role during development, in cancer metastasis and in inflammation. Activation of VEGFR-3 (also known as FLT4) by VEGF-C is one of the main drivers of lymphangiogenesis, but the transcriptional events downstream of VEGFR-3 activation are largely unknown. Recently, we identified a wave of immediate early transcription factors that are upregulated in human lymphatic endothelial cells (LECs) within the first 30 to 80 min after VEGFR-3 activation. Expression of these transcription factors must be regulated by additional pre-existing transcription factors that are rapidly activated by VEGFR-3 signaling. Using transcription factor activity analysis, we identified the homeobox transcription factor HOXD10 to be specifically activated at early time points after VEGFR-3 stimulation, and to regulate expression of immediate early transcription factors, including NR4A1. Gain- and loss-of-function studies revealed that HOXD10 is involved in LECs migration and formation of cord-like structures. Furthermore, HOXD10 regulates expression of VE-cadherin, claudin-5 and NOS3 (also known as e-NOS), and promotes lymphatic endothelial permeability. Taken together, these results reveal an important and unanticipated role of HOXD10 in the regulation of VEGFR-3 signaling in lymphatic endothelial cells, and in the control of lymphangiogenesis and permeability.

A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

FOXP1 inhibits cell growth and attenuates tumorigenicity of neuroblastoma

BACKGROUND

Segmental genomic copy number alterations, such as loss of 11q or 3p and gain of 17q, are well established markers of poor outcome in neuroblastoma, and have been suggested to comprise tumor suppressor genes or oncogenes, respectively. The gene forkhead box P1 (FOXP1) maps to chromosome 3p14.1, a tumor suppressor locus deleted in many human cancers including neuroblastoma. FoxP1 belongs to a family of winged-helix transcription factors that are involved in processes of cellular proliferation, differentiation and neoplastic transformation.

METHODS

Microarray expression profiles of 476 neuroblastoma specimens were generated and genes differentially expressed between favorable and unfavorable neuroblastoma were identified. FOXP1 expression was correlated to clinical markers and patient outcome. To determine whether hypermethylation is involved in silencing of FOXP1, methylation analysis of the 5' region of FOXP1 in 47 neuroblastomas was performed. Furthermore, FOXP1 was re-expressed in three neuroblastoma cell lines to study the effect of FOXP1 on growth characteristics of neuroblastoma cells.

RESULTS

Low expression of FOXP1 is associated with markers of unfavorable prognosis like stage 4, age >18 months and MYCN amplification and unfavorable gene expression-based classification (P < 0.001 each). Moreover, FOXP1 expression predicts patient outcome accurately and independently from well-established prognostic markers. Array-based CGH analysis of 159 neuroblastomas revealed that heterozygous loss of the FOXP1 locus was a rare event (n = 4), but if present, was associated with low FOXP1 expression. By contrast, DNA methylation analysis in 47 neuroblastomas indicated that hypermethylation is not regularly involved in FOXP1 gene silencing. Re-expression of FoxP1 significantly impaired cell proliferation, viability and colony formation in soft agar. Furthermore, induction of FOXP1 expression led to cell cycle arrest and apoptotic cell death of neuroblastoma cells.

CONCLUSIONS

Our results suggest that down-regulation of FOXP1 expression is a common event in high-risk neuroblastoma pathogenesis and may contribute to tumor progression and unfavorable patient outcome.

3'LIFE: a functional assay to detect miRNA targets in high-throughput

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene output at the post-transcriptional level by targeting degenerate elements primarily in 3'untranslated regions (3'UTRs) of mRNAs. Individual miRNAs can regulate networks of hundreds of genes, yet for the majority of miRNAs few, if any, targets are known. Misexpression of miRNAs is also a major contributor to cancer progression, thus there is a critical need to validate miRNA targets in high-throughput to understand miRNAs' contribution to tumorigenesis. Here we introduce a novel high-throughput assay to detect miRNA targets in 3'UTRs, called Luminescent Identification of Functional Elements in 3'UTRs (3'LIFE). We demonstrate the feasibility of 3'LIFE using a data set of 275 human 3'UTRs and two cancer-relevant miRNAs, let-7c and miR-10b, and compare our results to alternative methods to detect miRNA targets throughout the genome. We identify a large number of novel gene targets for these miRNAs, with only 32% of hits being bioinformatically predicted and 27% directed by non-canonical interactions. Functional analysis of target genes reveals consistent roles for each miRNA as either a tumor suppressor (let-7c) or oncogenic miRNA (miR-10b), and preferentially target multiple genes within regulatory networks, suggesting 3'LIFE is a rapid and sensitive method to detect miRNA targets in high-throughput.

The atypical cell cycle regulator Spy1 suppresses differentiation of the neuroblastoma stem cell population

Neuroblastoma is an aggressive pediatric cancer originating embryonically from the neural crest. The heterogeneity of the disease, as most solid tumors, complicates diagnosis and treatment. In neuroblastoma this heterogeneity is well represented in both primary tumours and derived cell lines and has been shown to be driven by a population of stem-like tumour initiating cells. Resolving the molecular mediators driving the division of this population of cells may indicate effective therapeutic options for neuroblastoma patients. This study has determined that the atypical cyclin-like protein Spy1, recently indicated in driving symmetric division of glioma stem cells, is a critical factor in the stem-like properties of neuroblastoma tumor initiating cell populations. Spy1 activates Cyclin Dependent Kinases (CDK) in a manner that is unique from classical cyclins. Hence this discovery may represent an important opportunity to design CDK inhibitor drugs to uniquely target subpopulations of cells within these aggressive neural tumours.

miR-224 promotion of cell migration and invasion by targeting Homeobox D 10 gene in human hepatocellular carcinoma

BACKGROUND AND AIM

MicroRNAs (miRNAs) are small noncoding RNA molecules that control target gene expression and are implicated in the regulation of diverse cellular pathways. In our previous research, we have demonstrated that miR-224 was overexpressed in liver cancer cells and tissues, which was an important factor in the regulation of cell migration and invasion. This study aimed to further explore the regulatory mechanism of miR-224 in the migration and invasion in liver cancer cells.

METHODS

A luciferase reporter assay was used to confirm that the HOXD10 gene was a direct target of miR-224. Quantitative reverse transcriptase-polymerase chain reaction, Western blotting, Transwell migration, and Matrigel invasion assays were performed to clarify the molecular mechanism of miR-224 in the regulation of cell migration and invasion in human hepatocellular carcinoma (HCC).

RESULTS

(i) The expression of miR-224 was strongly upregulated in MHHC97H and MHCC97L cells, and its expression level was significantly associated with cell invasive potential. (ii) The HOXD10 gene was confirmed to be a direct target of miR-224. Compared with normal liver tissues and cells, HOXD10 had lower expression in HCC tissues and cells and inversely regulated HCC cell invasion. (iii) miR-224 promoted expression of the tumor invasion-associated proteins p-PAK4 and MMP-9 by directly targeting HOXD10.

CONCLUSION

Our findings suggest a previously undescribed regulatory pathway in which the miR-224/HOXD10/p-PAK4/MMP-9 signaling pathway contributes to the regulation of cell migration and invasion and provides a new biotarget for HCC treatment.

A 92-gene cancer classifier predicts the site of origin for neuroendocrine tumors

A diagnosis of neuroendocrine carcinoma is often morphologically straight-forward; however, the tumor site of origin may remain elusive in a metastatic presentation. Neuroendocrine tumor subtyping has important implications for staging and patient management. In this study, the novel use and performance of a 92-gene molecular cancer classifier for determination of the site of tumor origin are described in a series of 75 neuroendocrine tumors (44 metastatic, 31 primary; gastrointestinal (n=12), pulmonary (n=22), Merkel cell (n=10), pancreatic (n=10), pheochromocytoma (n=10), and medullary thyroid carcinoma (n=11)). Formalin-fixed, paraffin-embedded samples passing multicenter pathologist adjudication were blinded and tested by a 92-gene molecular assay that predicts tumor type/subtype based upon relative quantitative PCR expression measurements for 87 tumor-related and 5 reference genes. The 92-gene assay demonstrated 99% (74/75; 95% confidence interval (CI) 0.93-0.99) accuracy for classification of neuroendocrine carcinomas and correctly subtyped the tumor site of origin in 95% (71/75; 95% CI 0.87-0.98) of cases. Analysis of gene expression subsignatures within the 92-gene assay panel showed 4 genes with promising discriminatory value for tumor typing and 15 genes for tumor subtyping. The 92-gene classifier demonstrated excellent accuracy for classifying and determining the site of origin in tumors with neuroendocrine differentiation. These results show promise for use of this test to aid in classifying neuroendocrine tumors of indeterminate primary site, particularly in the metastatic setting.

HOXC9 directly regulates distinct sets of genes to coordinate diverse cellular processes during neuronal differentiation

Background

Cellular differentiation is characterized by the acquisition of specialized structures and functions, cell cycle exit, and global attenuation of the DNA damage response. It is largely unknown how these diverse cellular events are coordinated at the molecular level during differentiation. We addressed this question in a model system of neuroblastoma cell differentiation induced by HOXC9.

Results

We conducted a genome-wide analysis of the HOXC9-induced neuronal differentiation program. Microarray gene expression profiling revealed that HOXC9-induced differentiation was associated with transcriptional regulation of 2,370 genes, characterized by global upregulation of neuronal genes and downregulation of cell cycle and DNA repair genes. Remarkably, genome-wide mapping by ChIP-seq demonstrated that HOXC9 bound to 40% of these genes, including a large number of genes involved in neuronal differentiation, cell cycle progression and the DNA damage response. Moreover, we showed that HOXC9 interacted with the transcriptional repressor E2F6 and recruited it to the promoters of cell cycle genes for repressing their expression.

Conclusions

Our results demonstrate that HOXC9 coordinates diverse cellular processes associated with differentiation by directly activating and repressing the transcription of distinct sets of genes.

Role of genomic architecture in the expression dynamics of long noncoding RNAs during differentiation of human neuroblastoma cells

BACKGROUND

Mammalian genomes are extensively transcribed producing thousands of long non-protein-coding RNAs (lncRNAs). The biological significance and function of the vast majority of lncRNAs remain unclear. Recent studies have implicated several lncRNAs as playing important roles in embryonic development and cancer progression. LncRNAs are characterized with different genomic architectures in relationship with their associated protein-coding genes. Our study aimed at bridging lncRNA architecture with dynamical patterns of their expression using differentiating human neuroblastoma cells model.

RESULTS

LncRNA expression was studied in a 120-hours timecourse of differentiation of human neuroblastoma SH-SY5Y cells into neurons upon treatment with retinoic acid (RA), the compound used for the treatment of neuroblastoma. A custom microarray chip was utilized to interrogate expression levels of 9,267 lncRNAs in the course of differentiation. We categorized lncRNAs into 19 architecture classes according to their position relatively to protein-coding genes. For each architecture class, dynamics of expression of lncRNAs was studied in association with their protein-coding partners. It allowed us to demonstrate positive correlation of lncRNAs with their associated protein-coding genes at bidirectional promoters and for sense-antisense transcript pairs. In contrast, lncRNAs located in the introns and downstream of the protein-coding genes were characterized with negative correlation modes. We further classified the lncRNAs by the temporal patterns of their expression dynamics. We found that intronic and bidirectional promoter architectures are associated with rapid RA-dependent induction or repression of the corresponding lncRNAs, followed by their constant expression. At the same time, lncRNAs expressed downstream of protein-coding genes are characterized by rapid induction, followed by transcriptional repression. Quantitative RT-PCR analysis confirmed the discovered functional modes for several selected lncRNAs associated with proteins involved in cancer and embryonic development.

CONCLUSIONS

This is the first report detailing dynamical changes of multiple lncRNAs during RA-induced neuroblastoma differentiation. Integration of genomic and transcriptomic levels of information allowed us to demonstrate specific behavior of lncRNAs organized in different genomic architectures. This study also provides a list of lncRNAs with possible roles in neuroblastoma.

HMCan: a method for detecting chromatin modifications in cancer samples using ChIP-seq data

MOTIVATION

Cancer cells are often characterized by epigenetic changes, which include aberrant histone modifications. In particular, local or regional epigenetic silencing is a common mechanism in cancer for silencing expression of tumor suppressor genes. Though several tools have been created to enable detection of histone marks in ChIP-seq data from normal samples, it is unclear whether these tools can be efficiently applied to ChIP-seq data generated from cancer samples. Indeed, cancer genomes are often characterized by frequent copy number alterations: gains and losses of large regions of chromosomal material. Copy number alterations may create a substantial statistical bias in the evaluation of histone mark signal enrichment and result in underdetection of the signal in the regions of loss and overdetection of the signal in the regions of gain.

RESULTS

We present HMCan (Histone modifications in cancer), a tool specially designed to analyze histone modification ChIP-seq data produced from cancer genomes. HMCan corrects for the GC-content and copy number bias and then applies Hidden Markov Models to detect the signal from the corrected data. On simulated data, HMCan outperformed several commonly used tools developed to analyze histone modification data produced from genomes without copy number alterations. HMCan also showed superior results on a ChIP-seq dataset generated for the repressive histone mark H3K27me3 in a bladder cancer cell line. HMCan predictions matched well with experimental data (qPCR validated regions) and included, for example, the previously detected H3K27me3 mark in the promoter of the DLEC1 gene, missed by other tools we tested.

Regulation of URG4/URGCP and PPARα gene expressions after retinoic acid treatment in neuroblastoma cells

Neuroblastoma (NB), originating from neural crest cells, is the most common extracranial tumor of childhood. Retinoic acid (RA) which is the biological active form of vitamin A regulates differentiation of NB cells, and RA derivatives have been used for NB treatment. PPARα (peroxisome proliferator-activated receptor) plays an important role in the oxidation of fatty acids, carcinogenesis, and differentiation. URG4/URGCP gene is a proto-oncogene and that overexpression of URG4/URGCP is associated with metastasis and tumor recurrence in osteosarcoma. It has been known that URG4/URGCP gene is an overexpressed gene in hepatocellular carcinoma and gastric cancers. This study aims to detect gene expression patterns of PPARα and URG4/URGCP genes in SH-SY5Y NB cell line after RA treatment. Expressions levels of PPARα and URG4/URGCP genes were analyzed after RA treatment for reducing differentiation in SH-SY5Y NB cell line. To induce differentiation, the cells were treated with 10 μM RA in the dark for 3-10 days. Gene expression of URG4/URGCP and PPARα genes were presented as the yield of polymerase chain reaction (PCR) products from target genes compared with the yield of PCR products from the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene. SH-SY5Y cells possess small processes in an undifferentiated state, and after treatment with RA, the cells developed long neurites, resembling a neuronal phenotype. PPARα gene expression increased in RA-treated groups; URG4/URGCP gene expression decreased in SH-SY5Y cells after RA treatment compared with that in the control cells. NB cell differentiation might associate with PPARα and URG4/URGCP gene expression profile after RA treatment.