Novel antiapoptotic effect of TBX15: overexpression of TBX15 reduces apoptosis in cancer cells

Epigenomic and transcriptomic profiling of solitary fibrous tumors identifies site-specific patterns and candidate genes regulated by DNA methylation

Solitary fibrous tumor (SFT) is a rare mesenchymal neoplasm which can arise at any anatomic site and is characterized by recurrent NAB2::STAT6 fusions and metastatic progression in 10-30%. The cell of origin has not been identified. Despite some progress in understanding the contribution of heterogeneous fusion types and secondary mutations to SFT biology, epigenetic alterations in extrameningeal SFT remain largely unexplored, and most sarcoma research to date has focused on the use of methylation profiling for tumor classification. We interrogated genome-wide DNA methylation in 79 SFTs to identify informative epigenetic changes. RNA-seq data from targeted panels and data from the Cancer Genome Atlas (TCGA) were used for orthogonal validation of selected findings. In unsupervised clustering analysis, the top 500 most variable CpGs segregated SFTs by primary anatomic site. Differentially methylated genes (DMGs) associated with primary SFT site included EGFR, TBX15, multiple HOX genes and their cofactors EBF1, EBF3, and PBX1, as well as RUNX1 and MEIS1. Of the 20 DMGs that were interrogated on the RNA-seq panel, twelve were significantly differentially expressed according to site. However, with the exception of TBX15, most of these also showed differential expression according to NAB2::STAT6 fusion type, suggesting that the fusion oncogene contributes to transcriptional regulation of these genes. Transcriptomic data confirmed an inverse correlation between gene methylation and the expression of TBX15 in both SFT and TCGA sarcomas. TBX15 also showed differential mRNA expression and 5' UTR methylation between tumors located in different anatomic sites in TCGA data. In all analyses, TBX15 methylation and mRNA expression retained the strongest association with tissue of origin in SFT and other sarcomas, suggesting a possible marker to distinguish metastatic tumors from new primaries without genomic profiling. Epigenetic signatures may further help to identify SFT progenitor cells at different anatomic sites.

The overexpression and clinical significance of TBX15 in human gliomas

T-box transcription factor 15 (TBX15) is upregulated in a variety of tumors and has been reported to promote uncontrolled proliferation of tumor cells and induce tumor cells to avoid apoptosis, thus accelerating the malignant transformation of malignant tumors. However, the prognostic value of TBX15 in glioma and its relationship with immune infiltration remain unknown. In this study, we intended to explore the prognostic value of TBX15 and its link to glioma immune infiltration and examine TBX15 expression in pan-cancer using RNAseq data in TPM format from TCGA and GTEx. TBX15 mRNA and protein expressions in glioma cells and adjacent normal tissue were detected and compared by RT-qPCR and Western blot. The effect of TBX15 on survival was assessed by Kaplan-Meier Method. The correlation between TBX15 upregulation and the clinicopathological characteristics of glioma patients was assessed by using TCGA databases, and the relationship between TBX15 and other genes in glioma was evaluated by using TCGA data. The top 300 genes most significantly associated with TBX15 were selected to establish a PPI network through the STRING database. The relationship between TBX15 mRNA expression and immune cell infiltration was explored by using ssGSEA and the TIMER Database. It was found that TBX15 mRNA expression in glioma tissues was significantly higher than that in the adjacent normal tissues, and this difference was most obvious in high-grade gliomas. TBX15 expression was increased in human gliomas and associated with worse clinicopathological characteristics and poorer survival prognosis in glioma patients. In addition, elevated TBX15 expression was linked to a collection of genes involved in immunosuppression. In conclusion, TBX15 played an important role in immune cell infiltration in glioma and may prove to be a predictor of the prognosis in glioma patients.

miR-4732-3p prevents lung cancer progression via inhibition of the TBX15/TNFSF11 axis

Aim: Possible roles of miRNAs in cancer treatment have been highly studied. This study aimed to elucidate the role of miR-4732-3p in lung cancer. Methods: Bioinformatics analysis was conducted to predict miR-4732-3p-related mRNA targets in lung cancer. Following interaction determination between miR-4732-3p and TBX15 as well as between TBX15 and TNFSF11, their in vitro and in vivo roles were assayed. Results: miR-4732-3p negatively targeted TBX15, which upregulated TNFSF11 by enhancing the activity of the TNFSF11 promoter. Overexpression of miR-4732-3p or silencing of TBX15 or TNFSF11 inhibited the malignant phenotype of lung cancer cells and reduced tumorigenicity in vivo. Conclusion: Overall, this study highlighted the inhibitory role of miR-4732-3p in lung cancer progression through the TBX15/TNFSF11 axis.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

TBX15 rs98422, DNM3 rs1011731, RAD51B rs8017304, and rs2588809 Gene Polymorphisms and Associations With Pituitary Adenoma

BACKGROUND

Pituitary adenoma (PA) is a benign tumor of parenchymal cells in the adenohypophysis, and it's development is strongly associated with genetic factors.This study aim was to find whether TBX15 rs98422, DNM3 rs1011731, RAD51B rs8017304, and rs2588809 single nucleotide polymorphisms can be associated with pituitary adenoma. While the TBX15 gene belongs to the T-box family of genes and is a transcription factor involved in many developmental processes, the DNM3 encodes a protein that is a member of the dynamin family with mechanochemical properties involved in actin-membrane processes, predominantly in membrane budding, and the RAD51B gene plays a significant role in homologous recombination in DNA repair for genome stability.

MATERIALS AND METHODS

The study enrolled 113 patients with pituitary adenoma and 283 healthy control subjects. DNA samples were extracted and purified from peripheral blood leukocytes. Genotyping was carried out using real-time polymerase chain reaction. The results were assessed using binomial logistic regression.

RESULTS

Our study revealed that RAD51B rs2588809 TT genotype could be associated with PA development in the co-dominant (OR=6.833; 95% CI=2.557-18.262; p<0.001) and recessive (OR=7.066; 95% CI=2.667-18.722; p<0.001) models. The same results were observed in females but not in males and PA without recurrence, while in PA with recurrence, no statistically significant results were obtained.

CONCLUSION

RAD51B rs2588809 TT genotype may increase the odds of PA development in women; it may also be associated with non-recurrent PA development.

Establishment of the TALE-code reveals aberrantly activated homeobox gene PBX1 in Hodgkin lymphoma

Homeobox genes encode transcription factors which regulate basic processes in development and cell differentiation and are grouped into classes and subclasses according to sequence similarities. Here, we analyzed the activities of the 20 members strong TALE homeobox gene class in early hematopoiesis and in lymphopoiesis including developing and mature B-cells, T-cells, natural killer (NK)-cells and innate lymphoid cells (ILC). The resultant expression pattern comprised eleven genes and which we termed TALE-code enables discrimination of normal and aberrant activities of TALE homeobox genes in lymphoid malignancies. Subsequent expression analysis of TALE homeobox genes in public datasets of Hodgkin lymphoma (HL) patients revealed overexpression of IRX3, IRX4, MEIS1, MEIS3, PBX1, PBX4 and TGIF1. As paradigm we focused on PBX1 which was deregulated in about 17% HL patients. Normal PBX1 expression was restricted to hematopoietic stem cells and progenitors of T-cells and ILCs but absent in B-cells, reflecting its roles in stemness and early differentiation. HL cell line SUP-HD1 expressed enhanced PBX1 levels and served as an in vitro model to identify upstream regulators and downstream targets in this malignancy. Genomic studies of this cell line therein showed a gain of the PBX1 locus at 1q23 which may underlie its aberrant expression. Comparative expression profiling analyses of HL patients and cell lines followed by knockdown experiments revealed NFIB and TLX2 as target genes activated by PBX1. HOX proteins operate as cofactors of PBX1. Accordingly, our data showed that HOXB9 overexpressed in HL coactivated TLX2 but not NFIB while activating TNFRSF9 without PBX1. Further downstream analyses showed that TLX2 activated TBX15 which operated anti-apoptotically. Taken together, we discovered a lymphoid TALE-code and identified an aberrant network around deregulated TALE homeobox gene PBX1 which may disturb B-cell differentiation in HL by reactivation of progenitor-specific genes. These findings may provide the framework for future studies to exploit possible vulnerabilities of malignant cells in therapeutic scenarios.

Reduction of T-Box 15 gene expression in tumor tissue is a prognostic biomarker for patients with hepatocellular carcinoma

Genome-wide analysis is widely applied to detect molecular alterations during oncogenesis and tumor progression. We analyzed DNA methylation profiles of hepatocellular carcinoma (HCC), and investigated the clinical role of most heypermethylated of tumor, encodes T-box 15 (TBX15), which was originally involved in mesodermal differentiation. We conducted a genome-wide analysis of DNA methylation of tumor and non-tumor tissue of 15 patients with HCC, and revealed TBX15 was the most hypermethylated gene of tumor (Beta-value in tumor tissue = 0.52 compared with non-tumor tissue). Another validation set, which comprised 58 HCC with radical resection, was analyzed to investigate the relationships between tumor phenotype and TBX15 mRNA expression. TBX15 mRNA levels in tumor tissues were significantly lower compared with those of nontumor tissues (p < 0.0001). When we assigned a cutoff value = 0.5-fold, the overall survival 5-year survival rates of the low-expression group (n = 17) were significantly shorter compared with those of the high-expression group (n = 41) (43.3% vs. 86.2%, p = 0.001). Multivariate analysis identified low TBX15 expression as an independent prognostic factor for overall and disease-free survival. Therefore, genome-wide DNA methylation profiling indicates that hypermethylation and reduced expression of TBX15 in tumor tissue represents a potential biomarker for predicting poor survival of patients with HCC.

Identification of contributing genes of Huntington's disease by machine learning

Background

Huntington’s disease (HD) is an inherited disorder caused by the polyglutamine (poly-Q) mutations of the HTT gene results in neurodegeneration characterized by chorea, loss of coordination, cognitive decline. However, HD pathogenesis is still elusive. Despite the availability of a wide range of biological data, a comprehensive understanding of HD’s mechanism from machine learning is so far unrealized, majorly due to the lack of needed data density.

Methods

To harness the knowledge of the HD pathogenesis from the expression profiles of postmortem prefrontal cortex samples of 157 HD and 157 controls, we used gene profiling ranking as the criteria to reduce the dimension to the order of magnitude of the sample size, followed by machine learning using the decision tree, rule induction, random forest, and generalized linear model.

Results

These four Machine learning models identified 66 potential HD-contributing genes, with the cross-validated accuracy of 90.79 ± 4.57%, 89.49 ± 5.20%, 90.45 ± 4.24%, and 97.46 ± 3.26%, respectively. The identified genes enriched the gene ontology of transcriptional regulation, inflammatory response, neuron projection, and the cytoskeleton. Moreover, three genes in the cognitive, sensory, and perceptual systems were also identified.

Conclusions

The mutant HTT may interfere with both the expression and transport of these identified genes to promote the HD pathogenesis.

Osteoporosis- and obesity-risk interrelationships: an epigenetic analysis of GWAS-derived SNPs at the developmental gene TBX15

A major challenge in translating findings from genome-wide association studies (GWAS) to biological mechanisms is pinpointing functional variants because only a very small percentage of variants associated with a given trait actually impact the trait. We used an extensive epigenetics, transcriptomics, and genetics analysis of the TBX15/WARS2 neighbourhood to prioritize this region's best-candidate causal variants for the genetic risk of osteoporosis (estimated bone density, eBMD) and obesity (waist-hip ratio or waist circumference adjusted for body mass index). TBX15 encodes a transcription factor that is important in bone development and adipose biology. Manual curation of 692 GWAS-derived variants gave eight strong candidates for causal SNPs that modulate TBX15 transcription in subcutaneous adipose tissue (SAT) or osteoblasts, which highly and specifically express this gene. None of these SNPs were prioritized by Bayesian fine-mapping. The eight regulatory causal SNPs were in enhancer or promoter chromatin seen preferentially in SAT or osteoblasts at TBX15 intron-1 or upstream. They overlap strongly predicted, allele-specific transcription factor binding sites. Our analysis suggests that these SNPs act independently of two missense SNPs in TBX15. Remarkably, five of the regulatory SNPs were associated with eBMD and obesity and had the same trait-increasing allele for both. We found that WARS2 obesity-related SNPs can be ascribed to high linkage disequilibrium with TBX15 intron-1 SNPs. Our findings from GWAS index, proxy, and imputed SNPs suggest that a few SNPs, including three in a 0.7-kb cluster, act as causal regulatory variants to fine-tune TBX15 expression and, thereby, affect both obesity and osteoporosis risk.

Tumor suppressor miR-33b-5p regulates cellular function and acts a prognostic biomarker in RCC

BACKGROUND

Renal cell carcinoma (RCC) is a renal parenchyma neoplasm with a 30% recurrence rate even when treated properly. MicroRNAs are noncoding small RNAs that are involved in cellular communication and may participate in cancer development. This study aimed to explore the relationship between miR-33b-5p expression and RCC progression and prognosis.

METHOD

RT-qPCR, CCK-8 assay, wound scratch assay, transwell assay and flow cytometry assay were used to evaluate the expression and function of miR-33b-5p in RCC. Additionally, RCC samples and survival data from The Cancer Genome Atlas were used to analyze the prognostic functions of miR-33b-5p.

RESULTS

miR-33b-5p expression in RCC tissues and cell lines (786-O, ACHN) were found to be significantly downregulated, compared with normal tissues and cell lines (P<0.001). The miR-33b-5p mimic transfected cells showed a slower proliferation rate (P<0.01), while its invasion ability decreased by 38.16% (786-O, P<0.001) and 49.19% (ACHN, P<0.05), compared with the negative control (NC). The migration ability of both RCC lines were found to be as follows: miR-33b-5p inhibitor > NC or NC inhibitor > miR-33b-5p mimic. Additionally, TCGA and RCC samples reveal that low miR-33b-5p expression is related to poor survival outcomes (univariate analysis, P=0.029; multivariate analysis, P=0.024; Kaplan-Meier survival curves, P=0.014). Target genes prediction suggests that miR-33b-5p performs its tumor-suppressive effects and prognostic role through targeting TBX15, SLC12A5, and PTGFRN.

CONCLUSIONS

miR-33b-5p may function as a tumor-suppressive regulator and prognostic biomarker in RCC.

A Hybrid Ensemble Approach for Identifying Robust Differentially Methylated Loci in Pan-Cancers

DNA methylation is a widely investigated epigenetic mark that plays a vital role in tumorigenesis. Advancements in high-throughput assays, such as the Infinium 450K platform, provide genome-scale DNA methylation landscapes in single-CpG locus resolution, and the identification of differentially methylated loci has become an insightful approach to deepen our understanding of cancers. However, the situation with extremely unbalanced numbers of samples and loci (approximately 1:1,000) makes it rather difficult to explore differential methylation between the sick and the normal. In this article, a hybrid approach based on ensemble feature selection for identifying differentially methylated loci (HyDML) was proposed by incorporating instance perturbation and multiple function models. Experiments on data from The Cancer Genome Atlas showed that HyDML not only achieved effective DML identification, but also outperformed the single-feature selection approach in terms of classification performance and the robustness of feature selection. The intensive analysis of the DML indicated that different types of cancers have mutual patterns, and the stable DML sharing in pan-cancers is of the great potential to be biomarkers, which may strengthen the confidence of domain experts to implement biological validations.

Regulation of antitumor miR-144-5p targets oncogenes: Direct regulation of syndecan-3 and its clinical significance

In the human genome, miR-451a, miR-144-5p (passenger strand), and miR-144-3p (guide strand) reside in clustered microRNA (miRNA) sequences located within the 17q11.2 region. Low expression of these miRNAs is significantly associated with poor prognosis of patients with renal cell carcinoma (RCC) (miR-451a: P = .00305; miR-144-5p: P = .00128; miR-144-3p: P = 9.45 × 10-5 ). We previously reported that miR-451a acted as an antitumor miRNA in RCC cells. Involvement of the passenger strand of the miR-144 duplex in the pathogenesis of RCC is not well understood. Functional assays showed that miR-144-5p and miR-144-3p significantly reduced cancer cell migration and invasive abilities, suggesting these miRNAs acted as antitumor miRNAs in RCC cells. Analyses of miR-144-5p targets identified a total of 65 putative oncogenic targets in RCC cells. Among them, high expression levels of 9 genes (FAM64A, F2, TRIP13, ANKRD36, CENPF, NCAPG, CLEC2D, SDC3, and SEMA4B) were significantly associated with poor prognosis (P < .001). Among these targets, expression of SDC3 was directly controlled by miR-144-5p, and its expression enhanced cancer cell aggressiveness. We identified genes downstream by SDC3 regulation. Data showed that expression of 10 of the downstream genes (IL18RAP, SDC3, SH2D1A, GZMH, KIF21B, TMC8, GAB3, HLA-DPB2, PLEK, and C1QB) significantly predicted poor prognosis of the patients (P = .0064). These data indicated that the antitumor miR-144-5p/oncogenic SDC3 axis was deeply involved in RCC pathogenesis. Clustered miRNAs (miR-451a, miR-144-5p, and miR-144-3p) acted as antitumor miRNAs, and their targets were intimately involved in RCC pathogenesis.

Unforeseen clonal evolution of tumor cell population in recurrent and metastatic dermatofibrosarcoma protuberans

Dermatofibrosarcoma protuberans (DFSP) is a very rare soft tissue sarcoma, generally of low-grade malignancy. DFSP is locally aggressive with a high recurrence rate, but metastasis occurs rarely. To investigate the mechanism of metastasis in DFSP, we analyzed the whole exome sequencing data of serial tumor samples obtained from a patient who had a 10-year history of recurrent and metastatic DFSP. Tracking various genomic alterations, namely somatic mutations, copy number variations, and chromosomal rearrangements, we observed a dramatic change in tumor cell population during the occurrence of metastasis in this DFSP case. The new subclone that emerged in metastatic DFSP harbored a completely different set of somatic mutations and new focal amplifications, which had not been observed in the primary clone before metastasis. The COL1A1-PDGFB fusion, characteristic of DFSP, was found in all of the serial samples. Moreover, the break position on the fusion gene was identical in all samples. Based on these observations, we suggest a clonal evolution model to explain the mechanism underlying metastasis in DFSP and identified several candidate target genes responsible for metastatic DFSP by utilizing The Cancer Genome Atlas database. This is the first study to observe clonal evolution in metastatic DFSP and provide insight for a possible therapeutic strategy for imatinib-resistant or metastatic DFSP.

NF-κB Mediates the Expression of TBX15 in Cancer Cells

TBX15 is a T-box transcription factor essential for development, also proposed as a marker in prostate cancer; and, recently, its antiapoptotic function indicates a role in carcinogenesis. Regulation of TBX15 is uncovered. In this study, we investigated the regulation of TBX15 expression in human cancer cells, by analyzing the regulatory function of a 5’-distal conserved region of TBX15. Bisulfite sequencing showed high methylation of the CpG island contained in this region that was not correlated with TBX15 mRNA levels, in the cancer cell lines analyzed; however, after 5-aza-dC treatment of TPC-1 cells an increase of TBX15 expression was observed. We also found a significant response of TBX15 to TNF-α activation of the NF-κB pathway using five cancer cell lines, and similar results were obtained when NF-κB was activated with PMA/ionomycin. Next, by luciferase reporter assays, we identified the TBX15 regulatory region containing two functional NF-κB binding sites with response to NF-κBp65, mapping on the -3302 and -3059 positions of the TBX15 gene. Moreover, a direct interaction of NF-κBp65 with one of the two NF-κB binding sites was indicated by ChIP assays. In summary, we provide novel data showing that NF-κB signaling up-regulates TBX15 expression in cancer cells. Furthermore, the link between TBX15 and NF-κB found in this study may be important to understand cancer and development processes.