Co-expression analysis and ceRNA network reveal eight novel potential lncRNA biomarkers in hepatocellular carcinoma

PAX6/CXCL14 regulatory axis promotes the repair of corneal injury by enhancing corneal epithelial cell proliferation

BACKGROUND

Corneal injuries, often leading to severe vision loss or blindness, have traditionally been treated with the belief that limbal stem cells (LSCs) are essential for repair and homeostasis, while central corneal epithelial cells (CCECs) were thought incapable of such repair. However, our research reveals that CCECs can fully heal and maintain the homeostasis of injured corneas in rats, even without LSCs. We discovered that CXCL14, under PAX6's influence, significantly boosts the stemness, proliferation, and migration of CCECs, facilitating corneal wound healing and homeostasis. This finding introduces CXCL14 as a promising new drug target for corneal injury treatment.

METHODS

To investigate the PAX6/CXCL14 regulatory axis's role in CCECs wound healing, we cultured human corneal epithelial cell lines with either increased or decreased expression of PAX6 and CXCL14 using adenovirus transfection in vitro. Techniques such as coimmunoprecipitation, chromatin immunoprecipitation, immunofluorescence staining, western blot, real-time PCR, cell colony formation, and cell cycle analysis were employed to validate the axis's function. In vivo, a rat corneal epithelial injury model was developed to further confirm the PAX6/CXCL14 axis's mechanism in repairing corneal damage and maintaining corneal homeostasis, as well as to assess the potential of CXCL14 protein as a therapeutic agent for corneal injuries.

RESULTS

Our study reveals that CCECs naturally express high levels of CXCL14, which is significantly upregulated by PAX6 following corneal damage. We identified SDC1 as CXCL14's receptor, whose engagement activates the NF-κB pathway to stimulate corneal repair by enhancing the stemness, proliferative, and migratory capacities of CCECs. Moreover, our research underscores CXCL14's therapeutic promise for corneal injuries, showing that recombinant CXCL14 effectively accelerates corneal healing in rat models.

CONCLUSION

CCECs play a critical and independent role in the repair of corneal injuries and the maintenance of corneal homeostasis, distinct from that of LSCs. The PAX6/CXCL14 regulatory axis is pivotal in this process. Additionally, our research demonstrates that the important function of CXCL14 in corneal repair endows it with the potential to be developed into a novel therapeutic agent for treating corneal injuries.

Identification of lineage-specific epigenetic regulators FOXA1 and GRHL2 through chromatin accessibility profiling in breast cancer cell lines

Breast cancer is a heterogeneous disease, and breast cancer cell lines are invaluable for studying this heterogeneity. However, the epigenetic diversity across these cell lines remains poorly understood. In this study, we performed genome-wide chromatin accessibility analysis on 23 breast cancer cell lines, including 2 estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative (ER+/HER2-), 3 ER+/HER2+, 3 HER2+, and 15 triple-negative breast cancer (TNBC) lines. These cell lines were classified into three groups based on their chromatin accessibility: the receptor-positive group (Group-P), TNBC basal group (Group-B), and TNBC mesenchymal group (Group-M). Motif enrichment analysis revealed that only Group-P exhibited coenrichment of forkhead box A1 (FOXA1) and grainyhead-like 2 (GRHL2) motifs, whereas Group-B was characterized by the presence of the GRHL2 motif without FOXA1. Notably, Group-M did not show enrichment of either FOXA1 or GRHL2 motifs. Furthermore, gene ontology analysis suggested that group-specific accessible regions were associated with their unique lineage characteristics. To investigate the epigenetic landscape regulatory roles of FOXA1 and GRHL2, we performed knockdown experiments targeting FOXA1 and GRHL2, followed by assay for transposase-accessible chromatin sequencing analysis. The findings revealed that FOXA1 maintains Group-P-specific regions while suppressing Group-B-specific regions in Group-P cells. In contrast, GRHL2 preserves commonly accessible regions shared between Group-P and Group-B in Group-B cells, suggesting that FOXA1 and GRHL2 play a pivotal role in preserving distinct chromatin accessibility patterns for each group. Specifically, FOXA1 distinguishes between receptor-positive and TNBC cell lines, whereas GRHL2 distinguishes between basal-like and mesenchymal subtypes in TNBC lines.

M6A-related bioinformatics analysis indicates that LRPPRC is an immune marker for ischemic stroke

Ischemic stroke (IS) is a common cerebrovascular disease whose pathogenesis involves a variety of immune molecules, immune channels and immune processes. 6-methyladenosine (m6A) modification regulates a variety of immune metabolic and immunopathological processes, but the role of m6A in IS is not yet understood. We downloaded the data set GSE58294 from the GEO database and screened for m6A-regulated differential expression genes. The RF algorithm was selected to screen the m6A key regulatory genes. Clinical prediction models were constructed and validated based on m6A key regulatory genes. IS patients were grouped according to the expression of m6A key regulatory genes, and immune markers of IS were identified based on immune infiltration characteristics and correlation. Finally, we performed functional enrichment, protein interaction network analysis and molecular prediction of the immune biomarkers. We identified a total of 7 differentially expressed genes in the dataset, namely METTL3, WTAP, YWHAG, TRA2A, YTHDF3, LRPPRC and HNRNPA2B1. The random forest algorithm indicated that all 7 genes were m6A key regulatory genes of IS, and the credibility of the above key regulatory genes was verified by constructing a clinical prediction model. Based on the expression of key regulatory genes, we divided IS patients into 2 groups. Based on the expression of the gene LRPPRC and the correlation of immune infiltration under different subgroups, LRPPRC was identified as an immune biomarker for IS. GO enrichment analyses indicate that LRPPRC is associated with a variety of cellular functions. Protein interaction network analysis and molecular prediction indicated that LRPPRC correlates with a variety of immune proteins, and LRPPRC may serve as a target for IS drug therapy. Our findings suggest that LRPPRC is an immune marker for IS. Further analysis based on LRPPRC could elucidate its role in the immune microenvironment of IS.

Alternative Polyadenylation Regulatory Factors Signature for Survival Prediction in Kidney Renal Cell Carcinoma

Background

Alternative polyadenylation (APA) plays a vital regulatory role in various diseases. It is widely accepted that APA is regulated by APA regulatory factors.

Objective

Whether APA regulatory factors affect the prognosis of renal cell carcinoma remains unclear, and this is the main topic of this study.

Methods

We downloaded the transcriptome and clinical data from The Cancer Genome Atlas (TCGA) database. We used the Lasso regression system to construct an APA model for analyzing the relationship between common APA regulatory factors and renal cell carcinoma. We also validated our APA model using independent GEO datasets (GSE29609, GSE76207).

Results

It was found that the expression levels of 5 APA regulatory factors (CPSF1, CPSF2, CSTF2, PABPC1, and PABPC4) were significantly associated with tumor gene mutation burden (TMB) score in renal clear cell carcinoma, and the risk score constructed using the expression level of 5 key APA regulatory factors could be used to predict the outcome of renal clear cell carcinoma. The TMB score is associated with the remodeling of the immune microenvironment.

Conclusions

By identifying key APA regulatory factors in renal cell carcinoma and constructing risk scores for key APA regulatory factors, we showed that key APA regulators affect prognosis of renal clear cell carcinoma patients. In addition, the risk score level is associated with TMB, indicating that APA may affect the efficacy of immunotherapy through immune microenvironment-related genes. This helps us better understand the mRNA processing mechanism of renal clear cell carcinoma.

Multi-omics approaches for drug-response characterization in primary biliary cholangitis and autoimmune hepatitis variant syndrome

BACKGROUND

Primary biliary cholangitis (PBC) and autoimmune hepatitis (AIH) variant syndrome (VS) exhibit a complex overlap of AIH features with PBC, leading to poorer prognoses than those with PBC or AIH alone. The biomarkers associated with drug response and potential molecular mechanisms in this syndrome have not been fully elucidated.

METHODS

Whole-transcriptome sequencing was employed to discern differentially expressed (DE) RNAs within good responders (GR) and poor responders (PR) among patients with PBC/AIH VS. Subsequent gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted for the identified DE RNAs. Plasma metabolomics was employed to delineate the metabolic profiles distinguishing PR and GR groups. The quantification of immune cell profiles and associated cytokines was achieved through flow cytometry and immunoassay technology. Uni- and multivariable logistic regression analyses were conducted to construct a predictive model for insufficient biochemical response. The performance of the model was assessed by computing the area under the receiver operating characteristic (AUC) curve, sensitivity, and specificity.

FINDINGS

The analysis identified 224 differentially expressed (DE) mRNAs, 189 DE long non-coding RNAs, 39 DE circular RNAs, and 63 DE microRNAs. Functional pathway analysis revealed enrichment in lipid metabolic pathways and immune response. Metabolomics disclosed dysregulated lipid metabolism and identified PC (18:2/18:2) and PC (16:0/20:3) as predictors. CD4+ T helper (Th) cells, including Th2 cells and regulatory T cells (Tregs), were upregulated in the GR group. Pro-inflammatory cytokines (IFN-γ, TNF-α, IL-9, and IL-17) were downregulated in the GR group, while anti-inflammatory cytokines (IL-10, IL-4, IL-5, and IL-22) were elevated. Regulatory networks were constructed, identifying CACNA1H and ACAA1 as target genes. A predictive model based on these indicators demonstrated an AUC of 0.986 in the primary cohort and an AUC of 0.940 in the validation cohort for predicting complete biochemical response.

CONCLUSION

A combined model integrating genomic, metabolic, and cytokinomic features demonstrated high accuracy in predicting insufficient biochemical response in patients with PBC/AIH VS. Early recognition of individuals at elevated risk for insufficient response allows for the prompt initiation of additional treatments.

ACACA reduces lipid accumulation through dual regulation of lipid metabolism and mitochondrial function via AMPK- PPARα- CPT1A axis

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a multifaceted metabolic disorder, whose global prevalence is rapidly increasing. Acetyl CoA carboxylases 1 (ACACA) is the key enzyme that controls the rate of fatty acid synthesis. Hence, it is crucial to investigate the function of ACACA in regulating lipid metabolism during the progress of NAFLD.

METHODS

Firstly, a fatty liver mouse model was established by high-fat diet at 2nd, 12th, and 20th week, respectively. Then, transcriptome analysis was performed on liver samples to investigate the underlying mechanisms and identify the target gene of the occurrence and development of NAFLD. Afterwards, lipid accumulation cell model was induced by palmitic acid and oleic acid (PA ∶ OA molar ratio = 1∶2). Next, we silenced the target gene ACACA using small interfering RNAs (siRNAs) or the CMS-121 inhibitor. Subsequently, experiments were performed comprehensively the effects of inhibiting ACACA on mitochondrial function and lipid metabolism, as well as on AMPK- PPARα- CPT1A pathway.

RESULTS

This data indicated that the pathways significantly affected by high-fat diet include lipid metabolism and mitochondrial function. Then, we focus on the target gene ACACA. In addition, the in vitro results suggested that inhibiting of ACACA in vitro reduces intracellular lipid accumulation, specifically the content of TG and TC. Furthermore, ACACA ameliorated mitochondrial dysfunction and alleviate oxidative stress, including MMP complete, ATP and ROS production, as well as the expression of mitochondria respiratory chain complex (MRC) and AMPK proteins. Meanwhile, ACACA inhibition enhances lipid metabolism through activation of PPARα/CPT1A, leading to a decrease in intracellular lipid accumulation.

CONCLUSION

Targeting ACACA can reduce lipid accumulation by mediating the AMPK- PPARα- CPT1A pathway, which regulates lipid metabolism and alleviates mitochondrial dysfunction.

Enhancing nasopharyngeal carcinoma cell radiosensitivity by suppressing AKT/mTOR via CENP-N knockdown

OBJECTIVE

Investigating the impact of centromere protein N (CENP-N) on radiosensitivity of nasopharyngeal carcinoma (NPC) cells.

METHODS

Using immunohistochemistry and immunofluorescence to detect CENP-N expression in tissues from 35 patients with radiosensitive or radioresistant NPC. Assessing the effect of combined CENP-N knockdown and radiotherapy on various cellular processes by CCK-8, colony formation, flow cytometry, and Western blotting. Establishing a NPC xenograft model. When the tumor volume reached 100 mm3, a irradiation dose of 6 Gy was given, and the effects of the combined treatment were evaluated in vivo using immunofluorescence and Western blotting techniques.

RESULTS

The level of CENP-N was significantly reduced in radiosensitive tissues of NPC (p < 0.05). Knockdown of CENP-N enhanced NPC radiosensitivity, resulting in sensitizing enhancement ratios (SER) of 1.44 (5-8 F) and 1.16 (CNE-2Z). The combined treatment showed significantly higher levels of proliferation suppression, apoptosis, and G2/M phase arrest (p < 0.01) compared to either CENP-N knockdown alone or radiotherapy alone. The combined treatment group showed the highest increase in Bax and γH2AX protein levels, whereas the protein Cyclin D1 exhibited the greatest decrease (p < 0.01). However, the above changes were reversed after treatment with AKT activator SC79. In vivo, the mean volume and weight of tumors in the radiotherapy group were 182 ± 54 mm3 and 0.16 ± 0.03 g. The mean tumor volume and weight in the combined treatment group were 84 ± 42 mm3 and 0.04 ± 0.01 g.

CONCLUSION

Knockdown of CENP-N can enhance NPC radiosensitivity by inhibiting AKT/mTOR.

Weighted gene co-expression network analysis revealed T cell differentiation associated with the age-related phenotypes in COVID-19 patients

The risk of severe condition caused by Corona Virus Disease 2019 (COVID-19) increases with age. However, the underlying mechanisms have not been clearly understood. The dataset GSE157103 was used to perform weighted gene co-expression network analysis on 100 COVID-19 patients in our analysis. Through weighted gene co-expression network analysis, we identified a key module which was significantly related with age. This age-related module could predict Intensive Care Unit status and mechanical-ventilation usage, and enriched with positive regulation of T cell receptor signaling pathway biological progress. Moreover, 10 hub genes were identified as crucial gene of the age-related module. Protein-protein interaction network and transcription factors-gene interactions were established. Lastly, independent data sets and RT-qPCR were used to validate the key module and hub genes. Our conclusion revealed that key genes were associated with the age-related phenotypes in COVID-19 patients, and it would be beneficial for clinical doctors to develop reasonable therapeutic strategies in elderly COVID-19 patients.

Bioinformatics Analysis of the Inflammation-Associated lncRNA-mRNA Coexpression Network in Type 2 Diabetes

Introduction

Diabetes is a chronic inflammatory state, and a key role of lncRNAs in diabetes complications is a new area of research.

Methods

In this study, key lncRNAs related to diabetes inflammation were identified by RNA-chip mining and lncRNA-mRNA coexpression network construction and finally verified by RT-qPCR.

Results

We ultimately obtained 12 genes, including A1BG-AS1, AC084125.4, RAMP2-AS1, FTX, DBH-AS1, LOXL1-AS1, LINC00893, LINC00894, PVT1, RUSC1-AS1, HCG25, and ATP1B3-AS1. RT-qPCR assays verified that LOXL1-AS1, A1BG-AS1, FTX, PVT1, and HCG25 were upregulated in the HG+LPS-induced THP-1 cells, and LINC00893, LINC00894, RUSC1-AS1, DBH-AS1, and RAMP2-AS1 were downregulated in the HG+LPS-induced THP-1 cells.

Conclusions

lncRNAs and mRNAs are extensively linked and form a coexpression network, and lncRNAs may influence the development of type 2 diabetes by regulating the corresponding mRNAs. The ten key genes obtained may become biomarkers of inflammation in type 2 diabetes in the future.

A universal co-expression gene network and prognostic model for hepatic-biliary-pancreatic cancers identified by integrative analyses

Hepatic, biliary and pancreatic cancers are a diverse set of malignancies with poor prognoses. It is possible that common molecular mechanisms are involved in the carcinogenesis of these cancers. Here, we identified LINC01537 and seven protein-coding genes by integrative analysis of transcriptomes of mRNAs, microRNAs and long non-coding RNAs from cholangiocarcinoma, hepatocellular carcinoma and pancreatic adenocarcinoma cohorts in TCGA. A predictive model constructed from seven biomarkers was established to successfully predict the survival rate of patients, which was then further verified in external cohorts. Additionally, patients with high-risk scores in our model were prone to epithelial-mesenchymal transition. Finally, activation of the biomarker PDE2A significantly attenuated migration and epithelial-mesenchymal transition in the HepG2 liver cancer cell line.

Identifying Potential Tumor Antigens and Antigens-Related Subtypes in Hepatocellular Carcinoma for mRNA Vaccine Development

Background

The mRNA vaccine has become a promising platform for cancer therapy. Lots of studies have been focusing on discovering novel potent cancer-associated antigens to develop mRNA vaccines against cancers. Besides, immunotyping shows the immune status, and immune microenvironment of immunotyping is related with therapeutic reaction. However, potential antigens for mRNA vaccines and immunotyping of liver hepatocellular carcinoma (LIHC) remain far from being understood.

Methods

In this study, we collected gene expression data and clinical information data from ICGC and TCGA databases. Using GEPIA2, we calculated differential expression genes and prognostic indices. We applied TIMER to calculate the correlation coefficient between immune infiltrating cells and each gene. Consensus cluster was used for immunotyping of LIHC.

Results

We uncovered four most potential candidates including PES1, MCM3, PPM1G, and KPNA2, which were all related with antigen-presenting cell (APC) infiltration and poor survival in LIHC in two independent datasets. Furthermore, three immune-related subtypes (IS1-IS3) of LIHC were identified. All these results were validated in two independent datasets. Furthermore, we validated our results in vitro.

Conclusions

The above candidates will be expected to be potential antigen genes for developing anti-LIHC mRNA vaccine, and furthermore, patients with IS2 and IS3 tumors are supposed to be appropriate for mRNA vaccine in LIHC.

Long noncoding RNA NR2F1-AS1 promotes the malignancy of non-small cell lung cancer via sponging microRNA-493-5p and thereby increasing ITGB1 expression

Several studies have reported that the long noncoding ribonucleic acid (lncRNA) NR2F1 antisense RNA 1 (NR2F1-AS1) affects multiple cellular pathways that are involved in tumorigenesis and tumor progression. The present study aimed to detect NR2F1-AS1 expression in non-small cell lung cancer (NSCLC), investigate the role of NR2F1-AS1 in promoting the tumorigenic behavior of NSCLC cells, and elucidate the mechanism underlying the effect of NR2F1-AS1 on NSCLC progression. Our results showed that NR2F1-AS1 expression was upregulated in NSCLC cells, and notably, its upregulation was correlated with adverse clinical characteristics and shorter overall survival in patients with NSCLC. The absence of NR2F1-AS1 functionally decreased NSCLC cell proliferation, migration, and invasion and promoted tumor cell apoptosis. In addition, the tumor growth of NSCLC cells in vivo was inhibited after NR2F1-AS1 silencing. Mechanistically, NR2F1-AS1 functioned as a competing endogenous RNA for miR-493-5p and consequently increased ITGB1 expression. Rescue assays further validated that an increased output of the miR-493-5p/ITGB1 axis could neutralize the regulatory impact of NR2F1-AS1 knockdown on the malignant phenotype of NSCLC cells. In summary, the NR2F1-AS1/miR-493-5p/ITGB1 pathway initiates pro-oncogenic behavior in NSCLC tumor progression, and the NR2F1-AS1/miR-493-5p/ITGB1 axis may provide new molecular targets for anticancer therapy against NSCLC.