Hypoxia-induced changes in intragenic DNA methylation correlate with alternative splicing in breast cancer

Bayesian Multi-View Clustering given complex inter-view structure

Multi-view datasets are becoming increasingly prevalent. These datasets consist of different modalities that provide complementary characterizations of the same underlying system. They can include heterogeneous types of information with complex relationships, varying degrees of missingness, and assorted sample sizes, as is often the case in multi-omic biological studies. Clustering multi-view data allows us to leverage different modalities to infer underlying systematic structure, but most existing approaches are limited to contexts in which entities are the same across views or have clear one-to-one relationships across data types with a common sample size. Many methods also make strong assumptions about the similarities of clusterings across views. We propose a Bayesian multi-view clustering approach (BMVC) which can handle the realities of multi-view datasets that often have complex relationships and diverse structure. BMVC incorporates known and complex many-to-many relationships between entities via a probabilistic graphical model that enables the joint inference of clusterings specific to each view, but where each view informs the others. Additionally, BMVC estimates the strength of the relationships between each pair of views, thus moderating the degree to which it imposes dependence constraints. We benchmarked BMVC on simulated data to show that it accurately estimates varying degrees of inter-view dependence when inter-view relationships are not limited to one-to-one correspondence. Next, we demonstrated its ability to capture visually interpretable inter-view structure in a public health survey of individuals and households in Puerto Rico following Hurricane Maria. Finally, we showed that BMVC clusters integrate the complex relationships between multi-omic profiles of breast cancer patient data, improving the biological homogeneity of clusters and elucidating hypotheses for functional biological mechanisms. We found that BMVC leverages complex inter-view structure to produce higher quality clusters than those generated by standard approaches. We also showed that BMVC is a valuable tool for real-world discovery and hypothesis generation.

The role of mitochondrial/metabolic axis in development of tamoxifen resistance in breast cancer

Only a few investigations, to our knowledge, have examined the bioenergetics of Tamoxifen (TMX) resistant individuals and reported altered mitochondrial activity and metabolic profile. The primary cause of TMX resistance is firmly suggested to be metabolic changes. Metabolic variations and hypoxia have also been linked in a bidirectional manner. Increased hypoxic levels correlate with early recurrence and proliferation and have a negative therapeutic impact on breast cancer (BC) patients. Hypoxia, carcinogenesis, and patient death are all correlated, resulting in more aggressive traits, a higher chance of metastasis, and TMX resistance. Consequently, we sought to investigate the possible role of the metabolic/hypoxial axis Long non-coding RNA (LncRNA) Taurine up-regulated 1 (TUG-1), Micro-RNA 186-5p (miR-186), Sirtuin-3 (SIRT3), Peroxisome Proliferator Activator Receptor alpha (PPAR-α), and Hypoxia-Inducible Factor-1 (HIF-1) in the development of TMX resistance in BC patients and to correlate this axis with tumor progression. Interestingly, this will be the first time to explore epigenetic regulation of this axis in BC.

Hypoxia-induced loss of SRSF2-dependent DNA methylation promotes CTCF-mediated alternative splicing of VEGFA in breast cancer

Alternative splicing of vascular endothelial growth factor A (VEGFA) generates numerous isoforms with unique roles in tumor angiogenesis, and investigating the underlying mechanism during hypoxia necessitates diligent pursuance. Our research systematically demonstrated that the splicing factor SRSF2 causes the inclusion of exon-8b, leading to the formation of the anti-angiogenic VEGFA-165b isoform under normoxic conditions. Additionally, SRSF2 interacts with DNMT3A and maintains methylation on exon-8a, inhibiting CCCTC-binding factor (CTCF) recruitment and RNA polymerase II (pol II) occupancy, causing exon-8a exclusion and decreased expression of pro-angiogenic VEGFA-165a. Conversely, SRSF2 is downregulated by HIF1α-induced miR-222-3p under hypoxic conditions, which prevents exon-8b inclusion and reduces VEGFA-165b expression. Furthermore, reduced SRSF2 under hypoxia promotes hydroxymethylation on exon-8a, increasing CTCF recruitment, pol II occupancy, exon-8a inclusion, and VEGFA-165a expression. Overall, our findings unveil a specialized dual mechanism of VEGFA-165 alternative splicing, instrumented by the cross-talk between SRSF2 and CTCF, which promotes angiogenesis under hypoxic conditions.

Sex-specific DNA methylation and associations with in utero tobacco smoke exposure at nuclear-encoded mitochondrial genes

Sex-linked differences in mitochondrial ATP production, enzyme activities, and reactive oxygen species generation have been reported in multiple tissue and cell types. While the effects of reproductive hormones underlie many of these differences, regulation of sexually dimorphic mitochondrial function has not been fully characterized. We hypothesized that sex-specific DNA methylation contributes to sex-specific expression of nuclear genes that influence mitochondrial function. Herein, we analysed DNA methylation data specifically focused on nuclear-encoded mitochondrial genes in 191 males and 190 females. We found 596 differentially methylated sites (DMSs) (FDR p < 0.05), corresponding to 324 genes, with at least a 1% difference in methylation between sexes. To investigate the potential functional significance, we utilized gene expression microarray data. Of the 324 genes containing DMSs, 17 showed differences in gene expression by sex. Particularly striking was that ATP5G2, encoding subunit C of ATP synthase, contains seven DMSs and exhibits a sex difference in expression (p = 0.04). Finally, we also found that alterations in DNA methylation associated with in utero tobacco smoke exposure were sex-specific in these nuclear-encoded mitochondrial genes. Interestingly, the level of sex differences in DNA methylation at nuclear-encoded mitochondrial genes and the level of methylation changes associated with smoke exposure were less prominent than that of other genes. This suggests more conservative regulation of DNA methylation at these nuclear-encoded mitochondrial genes as compared to others. Overall, our findings suggest that sex-specific DNA methylation may help establish sex differences in expression and function and that sex-specific alterations in DNA methylation in response to exposures could contribute to sex-variable toxicological responses.

Adverse maternal environment affects hippocampal HTR2c variant expression and epigenetic characteristics in mouse offspring

BACKGROUND

An adverse maternal environment (AME) predisposes progeny towards cognitive impairment in humans and mice. Cognitive impairment associates with hippocampal dysfunction. An important regulator of hippocampal function is the hippocampal serotonergic system. Dysregulation of hippocampal serotonin receptor 2c (HTR2c) expression is linked with cognitive impairment. HTR2c contains multiple mRNA variants and isoforms that are epigenetically regulated including DNA methylation, histone modifications, and small nucleolar RNA MBII-52. We tested the hypotheses that AME increases HTR2c variant expression and alters epigenetic modifications along the HTR2c gene locus.

METHODS

We create an AME through maternal Western diet and prenatal environmental stress in the mouse. We analyzed hippocampal HTR2c and variants' expression, DNA methylation and histone modifications along the gene locus, and MBII-52 levels in postnatal day 21 offspring.

RESULTS

AME significantly increased the expressions of total HTR2c and full-length variants (V201 and V202) concurrently with an altered epigenetic profile along the HTR2c gene locus in male offspring hippocampi. Moreover, increased full-length variants' expression in AME males was in line with increased MBII-52 levels.

CONCLUSIONS

AME affects male offspring hippocampal expression of HTR2c and full-length variants via epigenetic mechanisms. Altered hippocampal HTR2c expression may contribute to cognitive impairment seen in adult males in this model.

IMPACT

The key message of our article is that an adverse maternal environment increases expression of total HTR2c mRNA and protein, alters proportions of HTR2c mRNA variants, and impacts HTR2c epigenetic modifications in male offspring hippocampi relative to controls. Our findings add to the literature by providing the first report of altered HTR2c mRNA variant expression in association with altered epigenetic modifications in the hippocampus of offspring mice exposed to an adverse maternal environment. Our findings suggest that an adverse maternal environment affects the expression of genes previously determined to regulate cognitive function through an epigenetic mechanism in a sex-specific manner.

DEAH-box polypeptide 32 promotes hepatocellular carcinoma progression via activating the β-catenin pathway

PURPOSE

Hepatocellular carcinoma (HCC) is refractory cancer with high morbidity and high mortality. DEAH-box polypeptide 32 (DHX32) was upregulated in several types of malignancies and predicted poor prognosis. Herein, we investigated the role of DHX32 in HCC progression.

METHODS

The expression of DHX32, β-catenin, and epithelial-mesenchymal transition (EMT)-related makers were determined by Western blot and quantitative real-time PCR assays. Cell proliferation was tested by EdU cell proliferation assay. The effect of DHX32 and β-catenin on cell migration and invasion were detected by wound-healing and Traswell invasion assays. Tumour xenografts were performed to determine the effect of DHX32 on HCC tumour growth.

RESULTS

High level of DHX32 expression was associated with reduced overall survival in HCC patients. DHX32 expression was upregulated in human HCC cells and ectopic expression of DHX32 induced EMT, promoted the mobility and proliferation of HCC cells, and enhanced tumour growth in vivo. Silencing DHX32 reversed EMT, inhibited the malignancy behaviors of HCC cells, and suppressed tumour growth. Mechanistically, silencing DHX32 decreased the expression of β-cateninin in nucleus and β-catenin siRNA abrogated DHX32-mediated HCC progression.

CONCLUSION

DHX32 was an attractive regulator of HCC progression and indicated DHX32 canserve as a potential biomarker and therapeutic target for HCC patients.

Development and Verification of the Hypoxia-Related and Immune-Associated Prognosis Signature for Hepatocellular Carcinoma

Background

It has been widely suggested that the association of hypoxia with the immune status within the microenvironment of hepatocellular carcinoma (HCC) is of great clinical significance. The present work was carried out aiming to establish the hypoxia-related and immune-associated gene signature to stratify the risks in HCC.

Patients and Methods

The ssGSEA and t-SNE algorithms were utilized to estimate the immune and hypoxia statuses, respectively, using the TCGA database-derived cohort transcriptome profiles. Different immune groups are distinguished according to the ssGSEA scores, while the hypoxia-high and -low groups are inferred based on the distinct overall survival (OS) of the two groups of patients. Moreover, prognostic genes were identified using the Cox regression model in combination with the LASSO approach, which were later used to establish the hypoxia-related and immune-associated gene signature. At the same time, an ICGC cohort was used for external validation.

Results

A total of 13 genes, namely, HAVCR1, PSRC1, CCNJL, PDSS1, MEX3A, EID3, EPO, PLOD2, KPNA2, CDCA8, ADAMTS5, SLC1A7 and PIGZ, were discovered by the LASSO approach for constructing a gene signature to stratify the risk of HCC. Those low-risk cases showed superior prognosis (OS) to the high-risk counterparts (p<0.05). Moreover, it was suggested by multivariate analysis that our constructed hypoxia-related and immune-associated prognosis signature might be used as the independent factor for prognosis prediction (p<0.001). Patients in high-risk groups had severe hypoxia, higher immune checkpoint expression such as PD-L1, and different immunocyte infiltration states (eg, higher infiltration of regulatory T cells in the high-risk group) compared with those low-risk patients.

Conclusion

Our as-constructed hypoxia-related and immune-associated prognosis signature can be used as an approach to stratify the risk of HCC.

Machine learning analysis of DNA methylation in a hypoxia-immune model of oral squamous cell carcinoma

BACKGROUND

Hypoxia status and immunity are related with the development and prognosis of oral squamous cell carcinoma (OSCC). Here, we constructed a hypoxia-immune model to explore its upstream mechanism and identify potential CpG sites.

METHODS

The hypoxia-immune model was developed and validated by the iCluster algorithm. The LASSO, SVM-RFE and GA-ANN were performed to screen CpG sites correlated to the hypoxia-immune microenvironment.

RESULTS

We found seven hypoxia-immune related CpG sites. Lasso had the best classification performance among three machine learning algorithms.

CONCLUSION

We explored the clinical significance of the hypoxia-immune model and found seven hypoxia-immune related CpG sites by multiple machine learning algorithms. This model and candidate CpG sites may have clinical applications to predict the hypoxia-immune microenvironment.

Hypoxia-induced TGF-β-RBFOX2-ESRP1 axis regulates human MENA alternative splicing and promotes EMT in breast cancer

Hypoxic microenvironment heralds epithelial-mesenchymal transition (EMT), invasion and metastasis in solid tumors. Deregulation of alternative splicing (AS) of several cancer-associated genes has been instrumental in hypoxia-induced EMT. Our study in breast cancer unveils a previously unreported mechanism underlying hypoxia-mediated AS of hMENA, a crucial cytoskeleton remodeler during EMT. We report that the hypoxia-driven depletion of splicing regulator ESRP1 leads to skipping of hMENA exon 11a producing a pro-metastatic isoform, hMENAΔ11a. The transcriptional repression of ESRP1 is mediated by SLUG, which gets stimulated via hypoxia-driven TGF-β signaling. Interestingly, RBFOX2, an otherwise RNA-binding protein, is also found to transcriptionally repress ESRP1 while interacting with SLUG. Similar to SLUG, RBFOX2 gets upregulated under hypoxia via TGF-β signaling. Notably, we found that the exosomal delivery of TGF-β contributes to the elevation of TGF-β signaling under hypoxia. Moreover, our results show that in addition to hMENA, hypoxia-induced TGF-β signaling contributes to global changes in AS of genes associated with EMT. Overall, our findings reveal a new paradigm of hypoxia-driven AS regulation of hMENA and insinuate important implications in therapeutics targeting EMT.

PRKCSH Alternative Splicing Involves in Silica-Induced Expression of Epithelial-Mesenchymal Transition Markers and Cell Proliferation

Background:

Mounting evidence suggests that alternative splicing is one of the ways for cells to adapt to environmental stress insults. The aim of this study was firstly to examine the effect of silica on the alternative splicing of lung fibrosis–associated genes.

Methods:

Microarray analysis was used to construct the alternative splicing profile. Functional experiments were conducted using Cell Counting Kit-8, cell cycle, apoptosis, and epithelial–mesenchymal transition (EMT) analyses. Alternative splicing variants were verified by quantitative real-time polymerase chain reaction (qRT-PCR) polymerase chain reaction method.

Results:

A total of 1850 genes that have alternative splices in response to silica insult were identified. PCDHB11, MALAT1, MT2A, RP11-126D17.1, and RP11-415I12.2 are the top 5 upregulated genes with occurrence of alternative splice, whereas NDE1, RNPEPL1, TREML2, CSF2RB, and PRKCSH are the top 5 downregulated genes with occurrence of alternative splice. Bioinformatic analysis showed these genes with the occurrence of alternative splice mainly are associated with EMT pathway, N-Glycan biosynthesis, and leukocyte transendothelial migration. Further study indicated that PRKCSH-2 knockdown promotes A549 cell proliferation potential by partially promoting EMT signals.

Conclusions:

Significant changes in alternative splicing of silicosis-associated genes occur in patients with silicosis in silica conditions. Our study provides basic founding for further investigation into the detail molecular mechanisms underlying silica-induced silicosis.