Weighted gene correlation network analysis identifies RSAD2, HERC5, and CCL8 as prognostic candidates for breast cancer

MicroRNAs Associated with Keloids Identified by Microarray Analysis and In Vitro Experiments

MicroRNAs (miRNAs) play a crucial role in gene regulation and the development of keloid. This research aimed to identify and verify miRNAs associated with keloids by microarray analysis and in vitro experiments, shedding light on seeking for potential therapeutic molecular targets. In this study, the weighted gene co-expression network analysis was performed based on the GSE113620. The key miRNA module most relevant to the keloid was further screened to identify hub miRNAs, and then hub miRNAs was verified by the microarray analysis and qRT-PCR experiments. Additionally, targeted genes of hub miRNAs were predicted and verified. Gene ontology (GO) analysis and KEGG enrichment analysis were also conducted. Five miRNA modules were divided, and the blue module exhibited the highest correlation with keloids. Then, hsa-miR-127-3p, hsa-miR-214-3p, hsa-miR-155-5p, hsa-miR-409-5p, and hsa-miR-542-5p were identified as the hub miRNAs. Subsequently, the microarray analysis and qRT-PCR results demonstrated that the expression of five miRNAs were upregulated in keloid tissues. The GO analysis revealed that the target genes of these miRNAs were mainly enriched in biological processes including gene transcription, protein phosphorylation and the MAPK (mitogen-activated protein kinase) cascade, and the KEGG pathway enrichment analysis showed that the PI3K-AKT signaling pathway were significantly enriched. In conclusion, these five miRNAs (hsa-miR-127-3p, hsa-miR-155-5p, hsa-miR-214-3p, hsa-miR-409-5p, and hsa-miR-542-5p) play vital roles in the pathogenesis of keloid and might be potential therapeutic targets. These miRNAs might regulate genes enriched in gene transcription, protein phosphorylation, the MAPK cascade, and the PI3K-Akt signaling pathway.

Multi-omics and in-silico approach reveals AURKA, AURKB, and RSAD2 as therapeutic biomarkers in OSCC progression

Oral squamous cell carcinoma (OSCC), a prevalent form of head and neck cancer, poses a significant health challenge with limited improvements in patient outcomes over the years. Its development is influenced by a complex interplay of genetic alterations and environmental factors. While progress has been made in understanding the molecular mechanisms driving OSCC, pinpointing critical molecular markers and potential drug candidates has proven elusive. This study uniquely endeavors to conduct a meta-analysis to unveil therapeutic genes responsible for OSCC tumorigenesis. A multi-omics approach identified 951 differentially expressed genes (DEGs) associated with OSCC by analyzing microarray data from the NCBI GEO database. Weighted gene co-expression network analysis (WGCNA) detected a significant hub gene module comprising 805 genes, followed by the construction of protein-protein interaction network resulting in two small clusters of 7 gene-encoded proteins each. These clusters were filtered out based on top 10 significant pathways and gene ontology terms to identify six key target proteins with elevated expression levels, acting as potential therapeutic biomarkers for OSCC. Notably, RSAD2 emerged as a novel biomarker linked to OSCC progression. Furthermore, we identified potential inhibitors targeting AURKA, AURKB, and RSAD2 proteins and validated their interactions through molecular dynamics simulation studies. The simulations confirmed the stability of receptor-ligand complexes, suggesting ZINC03839281, ZINC04026167, and ZINC00718292 compounds hold promise as potential inhibitors for therapeutically targeting AURKA, AURKB, and RSAD2 as significant OSCC biomarkers. We recommend further comprehensive studies, including experimental and preclinical investigations, to validate the effectiveness of these lead compounds for OSCC treatment.

CircJPH1 regulates the NF-κB/HERC5 axis to promote the malignant progression of esophageal squamous cell carcinoma through binding to XRCC6

Esophageal squamous cell carcinoma (ESCC) is a prevalent and malignant cancer with an unknown pathogenesis and a poor prognosis; therefore, the identification of effective biomarkers and targets is crucial for its diagnosis and treatment. Circular (circ)RNAs are prominent functional biomarkers and therapeutic targets in various diseases, particularly cancer, due to their widespread expression and regulatory mechanisms. Our study aimed to investigate the therapeutic potential of circRNA for ESCC. We identified Hsa_circ_0137111 for the first time as one of the most significantly up-regulated genes in ESCC sequencing and named it circJPH1. The results of the present study demonstrated an enhanced expression of circJPH1 in ESCC tissues. Moreover, circJPH1-knockdown could significantly inhibit the proliferation, migration, and invasion of ESCC cells, while its overexpression promoted these characteristics. In addition, circJPH1 promoted ESCC cell tumor growth in vivo. For the first time, mass spectrometry and RNA pull-down analysis revealed the interaction of X-ray repair cross-complementary 6 (XRCC6) protein with circJPH1, thereby promoting its nuclear translocation. Consequently, the nuclear factor kappa-B (NF-κB) signaling pathway was activated, leading to an up-regulation of HECT and RLD domain containing E3 ubiquitin protein ligase 5 (HERC5), thereby promoting ESCC progression. In summary, the present study elucidated the regulatory impact of circJPH1 on ESCC progression in vitro and in vivo, thereby indicating its potential role in ESCC treatment.

Identification of novel protein biomarkers and therapeutic targets for ankylosing spondylitis using human circulating plasma proteomics and genome analysis

The proteome serves as the primary basis for identifying targets for treatment. This study conducted proteomic range two-sample Mendelian randomization (MR) analysis to pinpoint potential protein markers and treatment targets for ankylosing spondylitis (AS). A total of 4907 data points on circulating protein expression were collected from a large-scale protein quantitative trait locus investigation involving 35,559 individuals. Using data from a Finnish study on AS as the outcome, the dataset comprised 166,144 individuals of European ancestry (1462 cases and 164,682 controls), and causal relationships were determined through bidirectional Mendelian randomization of two samples. Proteins were further validated and identified through single-cell expression analysis, certain cells showing enriched expression levels were detected, and possible treatment targets were optimized. Increased HERC5 expression predicted by genes was related to increased AS risk, whereas the expression of the remaining five circulating proteins, AIF1, CREB3L4, MLN, MRPL55, and SPAG11B, was negatively correlated with AS risk. For each increase in gene-predicted protein levels, the ORs of AS were 2.11 (95% CI 1.44-3.09) for HERC5, 0.14 (95% CI 0.05-0.41) for AIF1, 0.48 (95% CI 0.34-0.68) for CREB3L4, 0.54 (95% CI 0.42-0.68) for MLN, 0.23 (95% CI 0.13-0.38) for MRPL55, and 0.26 (95% CI 0.17-0.39) for SPAG11B. The hypothesis of a reverse causal relationship between these six circulating proteins and AS is not supported. Three of the six protein-coding genes were expressed in both the AS and healthy control groups, while CREB3L4, MLN, and SPAG11B were not detected. Increased levels of HERC5 predicted by genes are related to increased AS risk, whereas the levels of the remaining five circulating proteins, AIF1, CREB3L4, MLN, MRPL55, and SPAG11B, negatively correlate with AS risk. HERC5, AIF1, and MRPL55 are potential therapeutic targets for AS. This study advanced the field by employing a novel combination of proteomic range two-sample MR analysis and single-cell expression analysis to identify potential protein markers and therapeutic targets for AS. This approach enabled a comprehensive understanding of the causal relationships between circulating proteins and AS, which has not been extensively explored in previous studies.

The Effects of Bacillus amyloliquefaciens SC06 on Behavior and Brain Function in Broilers Infected by Clostridium perfringens

Poultry studies conducted on Clostridium perfringens (CP) mainly focus on the effects of intestinal health and productive performance. Notably, the probiotic Bacillus amyloliquefaciens SC06 (BaSC06) is known to play a role in preventing bacterial infection. However, whether CP could induce the changes in brain function and behaviors and whether BaSC06 could play roles in these parameters is yet to be reported. The aim of this study was to evaluate the effects of BaSC06 on stress-related behaviors and gene expression, as well as the brain morphology and mRNA sequence of the hypothalamus in broiler chickens. A total of 288 one-day-old chicks were randomly divided into four groups: (1) a control group with no treatment administered or infection; (2) birds treated with the BaSC06 group; (3) a CP group; and (4) a BaSC06 plus CP (Ba_CP) group. The results showed that stress and fear-related behaviors were significantly induced by a CP infection and decreased due to the treatment of BaSC06. CP infection caused pathological damage to the pia and cortex of the brain, while BaSC06 showed a protective effect. CP significantly inhibited hypothalamic GABA and promoted HTR1A gene expression, while BaSC06 promoted GABA and decreased HTR1A gene expression. The different genes were nearly found between the comparisons of control vs. Ba group and Ba vs. CP group, while there were a great number of different genes between the comparisons of control vs. Ba_CP as well as CP vs. Ba_CP. Several different gene expression pathways were found that were related to disease, energy metabolism, and nervous system development. Our results will help to promote poultry welfare and health, as well as provide insights into probiotics to replace antibiotics and reduce resistance in the chicken industry.

HERC5 downregulation in non-small cell lung cancer is associated with altered energy metabolism and metastasis

BACKGROUND

Metastasis is the leading cause of cancer-related death in non-small cell lung cancer (NSCLC) patients. We previously showed that low HERC5 expression predicts early tumor dissemination and a dismal prognosis in NSCLC patients. Here, we performed functional studies to unravel the mechanism underlying the "metastasis-suppressor" effect of HERC5, with a focus on mitochondrial metabolism pathways.

METHODS

We assessed cell proliferation, colony formation potential, anchorage-independent growth, migration, and wound healing in NSCLC cell line models with HERC5 overexpression (OE) or knockout (KO). To study early tumor cell dissemination, we used these cell line models in zebrafish experiments and performed intracardial injections in nude mice. Mass spectrometry (MS) was used to analyze protein changes in whole-cell extracts. Furthermore, electron microscopy (EM) imaging, cellular respiration, glycolytic activity, and lactate production were used to investigate the relationships with mitochondrial energy metabolism pathways.

RESULTS

Using different in vitro NSCLC cell line models, we showed that NSCLC cells with low HERC5 expression had increased malignant and invasive properties. Furthermore, two different in vivo models in zebrafish and a xenograft mouse model showed increased dissemination and metastasis formation (in particular in the brain). Functional enrichment clustering of MS data revealed an increase in mitochondrial proteins in vitro when HERC5 levels were high. Loss of HERC5 leads to an increased Warburg effect, leading to improved adaptation and survival under prolonged inhibition of oxidative phosphorylation.

CONCLUSIONS

Taken together, these results indicate that low HERC5 expression increases the metastatic potential of NSCLC in vitro and in vivo. Furthermore, HERC5-induced proteomic changes influence mitochondrial pathways, ultimately leading to alterations in energy metabolism and demonstrating its role as a new potential metastasis suppressor gene.

CXADR promote epithelial-mesenchymal transition in endometriosis by modulating AKT/GSK-3β signaling

Endometriosis is a benign high prevalent disease exhibiting malignant features. However, the underlying pathogenesis and key molecules of endometriosis remain unclear. By integrating and analysis of existing expression profile datasets, we identified coxsackie and adenovirus receptor (CXADR), as a novel key gene in endometriosis. Based on the results of immunohistochemistry (IHC), we confirmed significant down-regulation of CXADR in ectopic endometrial tissues obtained from women with endometriosis compared with healthy controls. Further in vitro investigation indicated that CXADR regulated the stability and function of the phosphatases and AKT inhibitors PHLPP2 (pleckstrin homology domain and leucine-rich repeat protein phosphatase 2) and PTEN (phosphatase and tensin homolog). Loss of CXADR led to phosphorylation of AKT and glycogen synthase kinase-3β (GSK-3β), which resulted in stabilization of an epithelial-mesenchymal transition (EMT) factor, SNAIL1 (snail family transcriptional repressor 1). Therefore, EMT processs was induced, and the proliferation, migration and invasion of Ishikawa cells were enhanced. Over-expression of CXADR showed opposite effects. These findings suggest a previously undefined role of AKT/GSK-3β signaling axis in regulating EMT and reveal the involvement of a CXADR-induced EMT, in pathogenic progression of endometriosis.

Metastasis-Related Signature for Clinically Predicting Prognosis and Tumor Immune Microenvironment of Osteosarcoma Patients

Osteosarcoma is the most prevalent clinical malignant bone tumor in adolescents. The prognosis of metastatic osteosarcoma is still very poor. The aim of our study was to investigate the clinical diagnosis and prognostic significance of metastasis related genes (MRGs) in patients with osteosarcoma. Clinical information and RNA sequencing data with osteosarcoma patients were obtained and set as the training set from UCSC databases. GSE21257 were downloaded and chosen as the verification cohort. An eight gene metastasis related risk signature including MYC, TAC4, ABCA4, GADD45GIP1, TNFRSF21, HERC5, MAGEA11, and PDE1B was built to predict the overall survival of osteosarcoma patients. Based on risk assessments, patients were classified into high- and low-risk groups. The high-risk patients had higher risk score and shorter survival time. ROC curves revealed that this risk signature can accurately predict survival times of osteosarcoma patients at the 1-, 2-, 3-, 4- and 5- year. GSEA revealed that MYC targets, E2F targets, mTORC1 signaling, Wnt /β-catenin signaling and cell cycle were upregulated, and cell adhesion molecules, and primary immunodeficiency were decreased in high-risk group. MRGs were highly linked with the tumor immune microenvironment and ICB response. These results identified that MRGs as a novel prognostic and diagnostic biomarker in osteosarcoma.

The E3 ligase HERC5 promotes antimycobacterial responses in macrophages by ISGylating the phosphatase PTEN

Innate immune signaling in macrophages during viral infection is regulated by ISGylation, the covalent attachment of the ubiquitin-like protein interferon-stimulated gene 15 (ISG15) to protein targets. Here, we explored the role of ISGylation in the macrophage response to infection with Mycobacterium tuberculosis. In human and mouse macrophages, the E3 ubiquitin ligases HERC5 and mHERC6, respectively, mediated the ISGylation of the phosphatase PTEN, which promoted its degradation. The decreased abundance of PTEN led to an increase in the activity of the PI3K-AKT signaling pathway, which stimulated the synthesis of proinflammatory cytokines. Bacterial growth was increased in culture and in vivo when human or mouse macrophages were deficient in the major E3 ISG15 ligase. The findings expand the role of ISGylation in macrophages to antibacterial immunity and suggest that HERC5 signaling may be a candidate target for adjunct host-directed therapy in patients with tuberculosis.

The Functional Roles of ISG15/ISGylation in Cancer

The protein ISG15 encoded by interferon-stimulated gene (ISG) 15 is the first identified member of the ubiquitin-like protein family and exists in the form of monomers and conjugated complexes. Like ubiquitin, ISG15 can mediate an ubiquitin-like modification by covalently modifying other proteins, known as ISGylation. There is growing evidence showing that both the free and conjugated ISG15 are involved in multiple key cellular processes, including autophagy, exosome secretion, DNA repair, immune regulation, and cancer occurrence and progression. In this review, we aim to further clarify the function of ISG15 and ISGylation in cancer, demonstrate the important relationship between ISG15/ISGylation and cancer, and emphasize new insights into the different roles of ISG15/ISGylation in cancer progression. This review may contribute to therapeutic intervention in cancer. However, due to the limitations of current research, the regulation of ISG15/ISGylation on cancer progression is not completely clear, thus further comprehensive and sufficient correlation studies are still needed.

Transcriptomic and Functional Evidence That miRNA193a-3p Inhibits Lymphatic Endothelial Cell (LEC) and LEC + MCF-7 Spheroid Growth Directly and by Altering MCF-7 Secretome

MicroRNA 193a-3p (miR193a-3p) is a short non-coding RNA with tumor suppressor properties. Breast cancer (BC) progression is governed by active interaction between breast cancer cells, vascular (V)/lymphatic (L) endothelial cells (ECs), and BC secretome. We have recently shown that miR193a-3p, a tumor suppressor miRNA, inhibits MCF-7 BC cell-driven growth of VECs via direct antimitogenic actions and alters MCF-7 secretome. Since LEC-BC cross-talk plays a key role in BC progression, we investigated the effects of miR193a-3p on MCF-7 secretome and estradiol-mediated growth effects in LECs and LEC + MCF-7 spheroids, and delineated the underlying mechanisms. Transfection of LECs with miR193a-3p, as well as secretome from MCF-7 transfected cells, inhibited LEC growth, and these effects were mimicked in LEC + MCF-7 spheroids. Moreover, miR193a-3p inhibited ERK1/2 and Akt phosphorylation in LECs and LEC + MCF-7 spheroids, which are importantly involved in promoting cancer development and metastasis. Treatment of LECs and LEC + MCF-7 spheroids with estradiol (E2)-induced growth, as well as ERK1/2 and Akt phosphorylation, and was abrogated by miR193a-3p and secretome from MCF-7 transfected cells. Gene expression analysis (GEA) in LEC + MCF-7 spheroids transfected with miR193a-3p showed significant upregulation of 54 genes and downregulation of 73 genes. Pathway enrichment analysis of regulated genes showed significant modulation of several pathways, including interferon, interleukin/cytokine-mediated signaling, innate immune system, ERK1/2 cascade, apoptosis, and estrogen receptor signaling. Transcriptomic analysis showed downregulation in interferon and anti-apoptotic and pro-growth molecules, such as IFI6, IFIT1, OSA1/2, IFITM1, HLA-A/B, PSMB8/9, and PARP9, which are known to regulate BC progression. The cytokine proteome array of miR193a-3p transfected MCF secretome and confirmed the upregulation of several growth inhibitory cytokines, including IFNγ, Il-1a, IL-1ra, IL-32, IL-33, IL-24, IL-27, cystatin, C-reactive protein, Fas ligand, MIG, and sTIM3. Moreover, miR193a-3p alters factors in MCF-7 secretome, which represses ERK1/2 and Akt phosphorylation, induces pro-apoptotic protein and apoptosis in LECs, and downregulates interferon-associated proteins known to promote cancer growth and metastasis. In conclusion, miR193a-3p can potentially modify the tumor microenvironment by altering pro-growth BC secretome and inhibiting LEC growth, and may represent a therapeutic molecule to target breast tumors/cancer.

Identification of Key Genes and FUNCTIONAL Pathway in Radioresistance of Non-Small Cell Lung Cancer

Purpose

For better understanding of radiotherapy resistance and its potential mechanism.

Methods

We established radioresistance cell lines of non-small cell lung cancer (NSCLC) followed by microarray analysis. 529 differentially expressed genes (DEGs) were then screened between radiation resistant cell lines compared with the sensitive cell lines. The biological functions and enrichment pathways of the above DEGs were identified using Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analyses. Gene Set Enrichment Analysis (GSEA) revealed that the radiation resistance group had the most gene sets enriched in altered immune response, such as TNF signaling pathway, when compared to the radiation sensitive group. Protein-protein interaction (PPI) network was carried out through the STRING database, and then five hub genes (CXCL10, IFIH1, DDX58, CXCL11, RSAD2) were screened by Cytoscape software. RT-PCR confirmed the expression of the above hub genes. ChIP-X Enrichment Analysis showed that STAT1 might be the transcription factor of the above hub genes. Considering that PD-L1 could be activated by STAT1 in a variety of tumors and ultimately lead to immune exhaustion, RT-PCR and Western blot verified the expression level of PD-L1.

Results

Five hub genes (CXCL10, IFIH1, DDX58, CXCL11, RSAD2) were screened and verified to be highly expressed in radioresistance group, STAT1 might be the transcription factor of the above hub genes. Our study found that the expression level of PD-L1 was increased after radiotherapy resistance.

Conclusion

Although immune system activation occurs followed by radiation resistance, we hypothesized that the upregulation of PD-L1 expression caused by STAT1 activation might be one of the mechanisms of radiotherapy resistance.

Micro-RNA193a-3p Inhibits Breast Cancer Cell Driven Growth of Vascular Endothelial Cells by Altering Secretome and Inhibiting Mitogenesis: Transcriptomic and Functional Evidence

Breast cancer (BC) cell secretome in the tumor microenvironment (TME) facilitates neo-angiogenesis by promoting vascular endothelial cell (VEC) growth. Drugs that block BC cell growth or angiogenesis can restrict tumor growth and are of clinical relevance. Molecules that can target both BC cell and VEC growth as well as BC secretome may be more effective in treating BC. Since small non-coding microRNAs (miRs) regulate cell growth and miR193a-3p has onco-suppressor activity, we investigated whether miR193a-3p inhibits MCF-7-driven growth (proliferation, migration, capillary formation, signal transduction) of VECs. Using BC cells and VECs grown in monolayers or 3D spheroids and gene microarrays, we demonstrate that: pro-growth effects of MCF-7 and MDA-MB231 conditioned medium (CM) are lost in CM collected from MCF-7/MDA-MB231 cells pre-transfected with miR193a-3p (miR193a-CM). Moreover, miR193a-CM inhibited MAPK and Akt phosphorylation in VECs. In microarray gene expression studies, miR193a-CM upregulated 553 genes and downregulated 543 genes in VECs. Transcriptomic and pathway enrichment analysis of differentially regulated genes revealed downregulation of interferon-associated genes and pathways that induce angiogenesis and BC/tumor growth. An angiogenesis proteome array confirmed the downregulation of 20 pro-angiogenesis proteins by miR193a-CM in VECs. Additionally, in MCF-7 cells and VECs, estradiol (E2) downregulated miR193a-3p expression and induced growth. Ectopic expression of miR193a-3p abrogated the growth stimulatory effects of estradiol E2 and serum in MCF-7 cells and VECs, as well as in MCF-7 and MCF-7+VEC 3D spheroids. Immunostaining of MCF-7+VEC spheroid sections with ki67 showed miR193a-3p inhibits cell proliferation. Taken together, our findings provide first evidence that miR193a-3p abrogates MCF-7-driven growth of VECs by altering MCF-7 secretome and downregulating pro-growth interferon signals and proangiogenic proteins. Additionally, miR193a-3p inhibits serum and E2-induced growth of MCF-7, VECs, and MCF-7+VEC spheroids. In conclusion, miRNA193a-3p can potentially target/inhibit BC tumor angiogenesis via a dual mechanism: (1) altering proangiogenic BC secretome/TME and (2) inhibiting VEC growth. It may represent a therapeutic molecule to target breast tumor growth.

Systematic Analysis of the Grafting-Related Glucanase-Encoding GH9 Family Genes in Pepper, Tomato and Tobacco

Grafting is an important agricultural practice to control soil-borne diseases, alleviate continuous cropping problems and improve stress tolerance in vegetable industry, but it is relatively less applied in pepper production. A recent study has revealed the key roles of β-1, 4-glucanase in graft survival. We speculated that the GH9 family gene encoding glucanase may be involved in the obstacles of pepper grafting. Therefore, we performed a systematic analysis of the GH9 family in pepper, tomato and tobacco. A total of 25, 24 and 42 GH9 genes were identified from these three species. Compared with the orthologues of other solanaceous crops, the deduced pepper GH9B3 protein lacks a conserved motif (Motif 5). Promoter cis-element analysis revealed that a wound-responsive element exists in the promoter of tobacco NbGH9B3, but it is absent in the GH9B3 promoter of most solanaceous crops. The auxin-responsive related element is absent in CaGH9B3 promoter, but it presents in the promoter of tobacco, tomato, potato and petunia GH9B3. Tissue and induction expression profiles indicated that GH9 family genes are functionally differentiated. Nine GH9 genes, including CaGH9B3, were detected expressing in pepper stem. The expression patterns of NbGH9B3 and CaGH9B3 in grafting were different in our test condition, with obvious induction in tobacco but repression in pepper. Furthermore, weighted correlation network analysis (WGCNA) revealed 58 transcription factor genes highly co-expressed with NbGH9B3. Eight WRKY binding sites were detected in the promoter of NbGH9B3, and several NbWRKYs were highly co-expressed with NbGH9B3. In conclusion, the missing of Motif 5 in CaGH9B3, and lacking of wound- and auxin-responsive elements in the gene promoter are the potential causes of grafting-related problems in pepper. WRKY family transcription factors could be important regulator of NbGH9B3 in tobacco grafting. Our analysis points out the putative regulators of NbGH9B3, which would be helpful to the functional validation and the study of signal pathways related to grafting in the future.

Construction of miRNA-lncRNA-mRNA co-expression network affecting EMT-mediated cisplatin resistance in ovarian cancer

Platinum resistance is one of the major concerns in ovarian cancer treatment. Recent evidence shows the critical role of epithelial-mesenchymal transition (EMT) in this resistance. Epithelial-like ovarian cancer cells show decreased sensitivity to cisplatin after cisplatin treatment. Our study prospected the association between epithelial phenotype and response to cisplatin in ovarian cancer. Microarray dataset GSE47856 was acquired from the GEO database. After identifying differentially expressed genes (DEGs) between epithelial-like and mesenchymal-like cells, the module identification analysis was performed using weighted gene co-expression network analysis (WGCNA). The gene ontology (GO) and pathway analyses of the most considerable modules were performed. The protein-protein interaction network was also constructed. The hub genes were specified using Cytoscape plugins MCODE and cytoHubba, followed by the survival analysis and data validation. Finally, the co-expression of miRNA-lncRNA-TF with the hub genes was reconstructed. The co-expression network analysis suggests 20 modules relating to the Epithelial phenotype. The antiquewhite4, brown and darkmagenta modules are the most significant non-preserved modules in the Epithelial phenotype and contain the most differentially expressed genes. GO, and KEGG pathway enrichment analyses on these modules divulge that these genes were primarily enriched in the focal adhesion, DNA replication pathways and stress response processes. ROC curve and overall survival rate analysis show that the co-expression pattern of the brown module's hub genes could be a potential prognostic biomarker for ovarian cancer cisplatin resistance.

A Novel ERK2 Degrader Z734 Induces Apoptosis of MCF–7 Cells via the HERC3/p53 Signaling Pathway

Breast cancer is one of the leading causes of death worldwide, and synthetic chemicals targeting specific proteins or various molecular pathways for tumor suppression, such as ERK inhibitors and degraders, have been intensively investigated. The targets of ERK participate in the regulation of critical cellular mechanisms and underpin the progression of anticancer therapy. In this study, we identified a novel small molecule, which we named Z734, as a new mitogen–activated protein kinase 1 (ERK2) degrader and demonstrated that Z734 inhibits cell growth by inducing p53–mediated apoptotic pathways in human breast cancer cells. Treatment with Z734 resulted in the inhibition of cancer cell proliferation, colony formation and migration invasion, as well as cancer cell death via apoptosis. In addition, the Co–IP and GST pulldown assays indicated that the HECT and RLD domains containing E3 ubiquitin protein ligase 3 (HERC3) could directly interact with ERK2 through the HECT domain, promoting ERK2 ubiquitination. We also observed a strong link between HERC3 and p53 for the modulation of apoptosis. HERC3 can increase the protein and phosphorylation levels of p53, which further promotes apoptotic activity. In a xenograft mouse model, the effect was obtained in a treatment group that combined Z734 with lapatinib compared with that of the single–treatment groups. In summary, our results indicated that Z734 actively controls the development of breast cancer through apoptosis, and HERC3 may mediate ERK2 and p53 signaling, which offers new potential targets for clinical therapy.

HERC5/IFI16/p53 signaling mediates breast cancer cell proliferation and migration

AIMS

This study aims to find differentially expressed ubiquitination-related gene(s) and elucidates their biological significance in breast cancer.

MAIN METHODS

Differentially expressed genes were profiled in MCF-7 and MDA-MB-231 cells by using PCR array method. Abnormal expression of HERC5 was studied in the cells and in breast cancer specimens via Quantitative Real-time PCR and western blot. Cell proliferation and cell migration abilities were evaluated by using cell counting kits, or through colony formation, wound healing and trans-well assays. HERC5 target proteins were investigated via proteomic, co-immunoprecipitation and western blot methods. Down-stream signaling pathways were investigated through gene expression/knockdown methods.

KEY FINDINGS

Huge increase of HERC5 expression was found in MCF-7 and MDA-MB-231 cells, knockdown of which repressed the cell proliferation and migration. HERC5 interacted with IFI16, mediated IFI16 ISGylation at K274 and facilitated IFI16 proteasomal degradation. IFI16 acted as a tumor suppressor and to some extent mediated the HERC5 function in the breast cancer (BC) cells. HERC5 was negatively correlated with IFI16 protein, while IFI16 was positively correlated to p53 expression at mRNA and protein levels, which indicates a novel signaling pathway - HERC5/IFI16/p53. HERC5 expression was increased in glucose-starved BC cells and in human breast cancer tissues, accompanied with the decrease of IFI16 and P53.

SIGNIFICANCE

Our work reveals the abnormal expression of HERC5 and its carcinogenic role in breast cancer cells, which is probably mediated by an HERC5/IFI16/p53 signaling pathway. This work also provides potential diagnostic/therapeutic biomarkers for breast cancer diagnosis and treatment.

Using Weighted Gene Co-Expression Network Analysis to Identify Increased MND1 Expression as a Predictor of Poor Breast Cancer Survival

Objective

We used bioinformatics analysis to identify potential biomarker genes and their relationship with breast cancer (BC).

Materials and Methods

We used a weighted gene co-expression network analysis (WGCNA) to create a co-expression network based on the top 25% genes in the GSE24124, GSE33926, and GSE86166 datasets obtained from the Gene Expression Omnibus. We used the DAVID online platform to perform GO and KEGG pathway enrichment analyses and the Cytoscape CytoHubba plug-in to screen the potential genes. Then, we related the genes to prognostic values in BC using the Oncomine, GEPIA, and Kaplan–Meier Plotter databases. Findings were validated by immunohistochemical (IHC) staining in the Human Protein Atlas and the TCGA-BRCA cohort. LinkedOmics identified the interactive expressions of hub genes. We used UALCAN to evaluate the methylation levels of these hub genes. MethSurv and SurvivalMeth were used to assess the multilevel prognostic value. Finally, we assessed hub gene association with immune cell infiltration using TIMER.

Results

The mRNA levels of MKI67, UBE2C, GTSE1, CCNA2, and MND1 were significantly upregulated in BC, whereas ESR1, THSD4, TFF1, AGR2, and FOXA1 were significantly downregulated. The DNA methylation signature analysis showed a better prognosis in the low-risk group. Further subgroup analyses revealed that MND1 might serve as an independent risk factor for unfavorable BC prognosis. Additionally, MND1 expression levels positively correlate with the immune infiltration statuses of CD4+ T cells, CD8+ T cells, B cells, neutrophils, dendritic cells, and macrophages.

Conclusion

Our results indicate that the ten hub genes may be involved in BC’s carcinogenesis, development, or metastasis, and MND1 may be a potential biomarker and therapeutic target for BC.

HERC3 regulates epithelial-mesenchymal transition by directly ubiquitination degradation EIF5A2 and inhibits metastasis of colorectal cancer

E3 ligase is widely reported to exert fundamental functions in cancers. Through rigorous bioinformatic analysis concentrating E3 ligases based on data from Genotype-Tissue Expression (GTEx) and data from The Cancer Genome Atlas (TCGA), HERC3 was indicated to be downregulated in colorectal cancer (CRC) and HERC3 downregulation showed poor overall survival (OS) and disease-free survival (DFS). Through qRT-PCR, western blotting and Immunohistochemistry (IHC), analytical results were validated based on tissues in Zhongshan hospital. Functionally, HERC3 was indicated to inhibit the migration, invasion and metastasis in vitro and in vivo through transwell assays, wound healing assays and vivo experiments. And HERC3 could regulate epithelial-mesenchymal transition (EMT) in CRC. Furthermore, immunoprecipitation (IP), coimmunoprecipitation (co-IP) and GST-pulldown assays indicated that HERC3 could directly interact with EIF5A2 in vitro and in vivo through the RCC1 domain in HERC3. And HERC3 could function as an E3 to promote the K27 and K48-linked ubiquitination degradation of EIF5A2 via the HECT domain in HERC3, besides, K47, K67, K85, and K121 in EIF5A2 were identified as ubiquitination sites. In addition, HERC3 was indicated to affect the migration, invasion and metastasis and further regulatE EMT via EIF5A2/TGF-/Smad2/3 signal. The present study may provide insight into the mechanism of EMT in CRC.

The Highly Expressed IFIT1 in Nasopharyngeal Carcinoma Enhances Proliferation, Migration, and Invasion of Nasopharyngeal Carcinoma Cells

In this study, we aimed to identify potential targets modulating the progression of nasopharyngeal carcinoma (NPC) using integrated bioinformatics analysis and functional assays. Differentially expressed genes (DEGs) between NPC and normal tissues samples were obtained from publicly availably microarray datasets (GSE68799, GSE34573, and GSE53819) in the Gene Expression Omnibus (GEO) database. The bioinformatics analysis identified 49 common DEGs from three GEO datasets, which were mainly enriched in cytokine/chemokine pathways and extracellular matrix organization pathway. Further protein-protein interaction network analysis identified 11 hub genes from the 49 DEGs. The 11 hub genes were significantly up-regulated in the NPC tissues when compared to normal tissues by analyzing the Oncomine database. The 8 hub genes including COL5A1, COL7A1, COL22A1, CXCL11, IFI44L, IFIT1, RSAD2, and USP18 were significantly up-regulated in the NPC tissues when compared to normal tissues by using the Oncomine database. Further validation studies showed that IFIT1 was up-regulated in the NPC cells. Knockdown of IFI1T1 suppressed the proliferation, migration, and invasion of NPC cells; while IFIT1 overexpression promoted the proliferation, migration, and invasion of NPC cells. In conclusion, a total of 49 DEGs and 11 hub genes in NPC using the integrated bioinformatics analysis. IFIT1 was up-regulated in the NPC cells lines, and IFIT1 may act as an oncogene by promoting NPC cell proliferation, migration, and invasion.

Identification of potential prognostic markers associated with lung metastasis in breast cancer by weighted gene co-expression network analysis

Breast cancer (BC) is an aggressive cancer with a high percentage recurrence and metastasis. As one of the most common distant metastasis organ in BC, lung metastasis has a worse prognosis than that of liver and bone. Therefore, it's important to explore some potential prognostic markers associated with the lung metastasis in BC for preventive treatment. In this study, transcriptomic data and clinical information of BC patients were downloaded from The Cancer Genome Atlas (TCGA) database. Co-expression modules constructed by weighted gene co-expression network analysis (WGCNA) found the royal blue module was significantly associated with lung metastasis in BC. Then, co-expression genes of this module were analyzed for functional enrichment. Furthermore, the prognostic value of these genes was assessed by GEPIA Database and Kaplan-Meier Plotter. Results showed that the hub genes, LMNB and CDC20, were up-regulated in BC and had a worse survival of the patients. Therefore, we speculate that these two genes play crucial roles in the process of lung metastasis in BC, which can be used as potential prognostic markers in lung metastasis of BC. Collectively, our study identified two potential key genes in the lung metastasis of BC, which might be applied as the prognostic markers of the precise treatment in breast cancer with lung metastasis.

Potential Novel Modules and Hub Genes as Prognostic Candidates of Thyroid Cancer by Weighted Gene Co-Expression Network Analysis

Background

Although thyroid cancer (THCA) is one of the most common type of endocrine malignancy, its highly complex molecular mechanisms of carcinogenesis are not completely known.

Materials and Methods

In this study, weighted gene co-expression network analysis (WGCNA) was utilized to construct gene co-expression networks and evaluate the relations between modules and clinical traits to identify potential prognostic biomarkers for THCA patients. RNA-seq data and clinical data were downloaded from The Cancer Genome Atlas (TCGA). Other independent datasets from the Gene Expression Omnibus (GEO) database and the Human Protein Atlas database were performed to validate findings.

Results

Finally, 11 co-expression modules were constructed and four hub genes, CCDC146, SLC4A4, TDRD9 and MUM1L1, were identified and validated statistically, which were considerably interrelated to worse survival of THCA patients.

Conclusion

This research study revealed four hub genes may be considered candidate prognostic biomarkers and potential therapeutic targets for THCA patients in the future.

CD200 Blockade Modulates Tumor Immune Microenvironment but Fails to Show Efficacy in Inhibiting Tumor Growth in a Murine Model of Melanoma

CD200-CD200R pathway regulates immune responses and has been implicated in the pathogenesis of a number of cancer types. CD200 blockade is considered a strategy for immunotherapy of CD200-positive cancers such as melanoma. Thus, it is critical to understand the potential impacts of CD200 blockade in a more human relevant tumor model. In this study, we evaluated these issues using the CD200⁺ Yumm1.7 mouse melanoma model. Yumm1.7 cells bear Braf/Pten mutations resembling human melanoma. We found that Yumm1.7 tumors grow significantly faster in CD200R^–/– mice compared to wild type mice. Analysis of tumor immune microenvironment (TIME) revealed that tumors from CD200R^–/– or anti-CD200 treated mice had downregulated immune cell contents and reduced TCR clonality compared to tumors from untreated wild type mice. T cells also showed impaired effector functions, as reflected by reduced numbers of IFN-γ⁺ and TNF-α⁺ T cells. Mechanistically, we found upregulation of the CCL8 gene in CD200R^–/– tumors. In vitro co-culture experiments using Yumm1.7 tumor cells with bone marrow derived macrophages (BMDM) from WT and CD200R^–/– mice confirmed upregulation of macrophage CCL8 in the absence of CD200-CD200R interaction. Finally, we found that anti-CD200 therapy failed to show efficacy either alone or in combination with checkpoint inhibitors such as anti-PD-1 or anti-CTLA4 in inhibiting Yumm1.7 tumor growth. Given that CD200R-deficiency or anti-CD200 treatment leads to reduced T cell responses in TME, using blockade of CD200 as an immunotherapy for cancers such as melanoma should be practiced with caution.

Gene expression network related to DNA methylation and miRNA regulation during the process of aflatoxin B1-induced malignant transformation of L02 cells

Aflatoxin is a secondary metabolite secreted by Aspergillus flavus, parasitic Aspergillus, and other fungi through the polyketone pathway, and it can be detected in many foods. Aflatoxin has strong toxicity and carcinogenicity, and many studies have shown that aflatoxin is highly associated with liver cancer. In the present study, malignant transformation of L02 cells was induced by aflatoxin B1 (AFB1), and the gene expression, miRNA expression, and methylation level were detected by high-throughput sequencing. The gene and miRNA expression results showed that 2547 genes and 315 miRNAs were changed in the AFB1-treated group compared with the DMSO group. Among them, RSAD2 and SCIN were significantly upregulated, whereas TRAPPC3L and UBE2L6 were significantly downregulated. Has-miR-33b-3p was significantly upregulated, whereas Has-miR-3613-5p was significantly downregulated. The methylation results showed that 2832 CpG sites were methylated on the promoter or coding DNA sequence (CDS) of the gene, whereas the expression of DNMT3a and DNMT3b was significantly upregulated. Moreover, hypermethylation occurred in TRAPPC3L, CDH13, and SPINK13. The results of GO and KEGG pathway analyses showed that significantly changed genes and miRNAs were mainly involved in tumor formation, proliferation, invasion, and migration. The results of network map analysis showed that Hsa-miR-3613-5p, Hsa-miR-615-5p, Hsa-miR-615-3p, and Hsa-miR-3158-3p were the key miRNAs for malignant transformation of L02 cells induced by AFB1. In addition, the expression of ONECUT2, RAP1GAP2, and FSTL4 was regulated by DNA methylation and miRNAs. These results suggested that the gene expression network regulated by DNA methylation and miRNAs may play a vital role in AFB1-induced hepatocellular carcinoma.

Identification of latent biomarkers in connection with progression and prognosis in oral cancer by comprehensive bioinformatics analysis

BACKGROUND

Oral cancer (OC) is a common and dangerous malignant tumor with a low survival rate. However, the micro level mechanism has not been explained in detail.

METHODS

Gene and miRNA expression micro array data were extracted from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) and miRNAs (DE miRNAs) were identified by R software. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of genes and genomes (KEGG) pathway analysis were used to assess the potential molecular mechanisms of DEGs. Cytoscape software was utilized to construct protein-protein interaction (PPI) network and miRNA-gene network. Central genes were screened out with the participation of gene degree, molecular complex detection (MCODE) plugin, and miRNA-gene network. Then, the identified genes were checked by The Cancer Genome Atlas (TCGA) gene expression profile, Kaplan-Meier data, Oncomine, and the Human Protein Atlas database. Receiver operating characteristic (ROC) curve was drawn to predict the diagnostic efficiency of crucial gene level in normal and tumor tissues. Univariate and multivariate Cox regression were used to analyze the effect of dominant genes and clinical characteristics on the overall survival rate of OC patients.

RESULTS

Gene expression data of gene expression profiling chip(GSE9844, GSE30784, and GSE74530) were obtained from GEO database, including 199 tumor and 63 non-tumor samples. We identified 298 gene mutations, including 200 upregulated and 98 downregulated genes. GO functional annotation analysis showed that DEGs were enriched in extracellular structure and extracellular matrix containing collagen. In addition, KEGG pathway enrichment analysis demonstrated that the DEGs were significantly enriched in IL-17 signaling pathway and PI3K-Akt signaling pathway. Then, we detected three most relevant modules in PPI network. Central genes (CXCL8, DDX60, EIF2AK2, GBP1, IFI44, IFI44L, IFIT1, IL6, MMP9,CXCL1, CCL20, RSAD2, and RTP4) were screened out with the participation of MCODE plugin, gene degree, and miRNA-gene network. TCGA gene expression profile and Kaplan-Meier analysis showed that high expression of CXCL8, DDX60, IL6, and RTP4 was associated with poor prognosis in OC patients, while patients with high expression of IFI44L and RSAD2 had a better prognosis. The elevated expression of CXCL8, DDX60, IFI44L, RSAD2, and RTP44 in OC was verified by using Oncomine database. ROC curve showed that the mRNA levels of these five genes had a helpful diagnostic effect on tumor tissue. The Human Protein Atlas database showed that the protein expressions of DDX60, IFI44L, RSAD2, and RTP44 in tumor tissues were higher than those in normal tissues. Finally, univariate and multivariate Cox regression showed that DDX60, IFI44L, RSAD2, and RTP44 were independent prognostic indicators of OC.

CONCLUSION

This study revealed the potential biomarkers and relevant pathways of OC from publicly available GEO database, and provided a theoretical basis for elucidating the diagnosis, treatment, and prognosis of OC.

Survival-Based Biomarker Module Identification Associated with Oral Squamous Cell Carcinoma (OSCC)

Simple Summary

In this study, four OSCC-specific hub genes were identified using high-throughput RNA-Seq data from TCGA cohort. The significant genes within tumor and normal samples were used for weighted PPI network construction based on survival of patients along with their expression profiles. The analysis revealed the most significant module in the training and test datasets. The genes from this module were used for pathway enrichment analysis followed by hub gene selection. These novel biomarkers might have clinical utility for diagnosis and prognosis prediction in OSCC, providing diagnosis at a very early stage. Moreover, a combination of all these biomarkers might distinguish the OSCC patients with low risk and high risk for cancer progression and recurrence, which will provide useful guidance for personalized and precision therapy. However, the results in the present study were obtained by integrative theoretical analysis, and the findings remain to be confirmed by further experimental validations.

Abstract

Head and neck squamous cell carcinoma (HNSC) is one of the most common malignant tumors worldwide with a high rate of morbidity and mortality, with 90% of predilections occurring for oral squamous cell carcinoma (OSCC). Cancers of the mouth account for 40% of head and neck cancers, including squamous cell carcinomas of the tongue, floor of the mouth, buccal mucosa, lips, hard and soft palate, and gingival. OSCC is the most devastating and commonly occurring oral malignancy, with a mortality rate of 500,000 deaths per year. This has imposed a strong necessity to discover driver genes responsible for its progression and malignancy. In the present study we filtered oral squamous cell carcinoma tissue samples from TCGA-HNSC cohort, which we followed by constructing a weighted PPI network based on the survival of patients and the expression profiles of samples collected from them. We found a total of 46 modules, with 18 modules having more than five edges. The KM and ME analyses revealed a single module (with 12 genes) as significant in the training and test datasets. The genes from this significant module were subjected to pathway enrichment analysis for identification of significant pathways and involved genes. Finally, the overlapping genes between gene sets ranked on the basis of weighted PPI module centralities (i.e., degree and eigenvector), significant pathway genes, and DEGs from a microarray OSCC dataset were considered as OSCC-specific hub genes. These hub genes were clinically validated using the IHC images available from the Human Protein Atlas (HPA) database.

Weighted Gene Co-Expression Network Analysis Identifies Five Hub Genes Associated with Metastasis in Synovial Sarcoma

BACKGROUND

Synovial sarcoma (SS) refers to a malignant soft tissue sarcoma (STS) which often occurs in children and adults and has a poor prognosis in elderly patients. Patients with local lesions can be treated with extensive surgical resection combined with adjuvant or radiotherapy, whereas about half of the cases have recurrent diseases and metastatic lesions, and five-year survival ratio is assessed within the range of 27% - 55% only.

METHOD

We downloaded a set of expression profile data (GSE40021) related to SS metastasis based on the Gene Expression Omnibus (GEO) database, and selected distinctly represented genes (DEGs) related to tumor metastasis. WGCNA was used to emphasize the DEGs related to tumor metastasis and obtain co-expression modules. Then, the module most related to SS metastasis was screened out. The genes enriched in this module were analyzed by Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway improvement analysis. Cytoscape software was used for constructing protein-protein interaction (PPI) networks, and hub genes were screened in Oncomine analysis.

RESULT

We selected 514 DEGs, consisting of 210 up-regulated genes and 304 down-regulated genes. Through WGCAN, we got seven co-expression modules and the module most related to SS metastasis was the turquoise module, which contained 66 genes. Finally, we screened out five hub genes (HJURP, NCAPG, TPX2, CENPA, NDC80) through CytoHubba and Oncomine analysis.

CONCLUSION

In this study, we screened five hub genes that may help in clinical diagnosis and serve as the latent purpose of SS treatment.

Predicting the Key Genes Involved in Aortic Valve Calcification Through Integrated Bioinformatics Analysis

Background: Valvular heart disease is obtaining growing attention in the cardiovascular field and it is believed that calcific aortic valve disease (CAVD) is the most common valvular heart disease (VHD) in the world. CAVD does not have a fully effective treatment to delay its progression and the specific molecular mechanism of aortic valve calcification remains unclear.

Materials and Methods: We obtained the gene expression datasets GSE12644 and GSE51472 from the public comprehensive free database GEO. Then, a series of bioinformatics methods, such as GO and KEGG analysis, STING online tool, Cytoscape software, were used to identify differentially expressed genes in CAVD and healthy controls, construct a PPI network, and then identify key genes. In addition, immune infiltration analysis was used via CIBERSORT to observe the expression of various immune cells in CAVD.

Results: A total of 144 differential expression genes were identified in the CAVD samples in comparison with the control samples, including 49 up-regulated genes and 95 down-regulated genes. GO analysis of DEGs were most observably enriched in the immune response, signal transduction, inflammatory response, proteolysis, innate immune response, and apoptotic process. The KEGG analysis revealed that the enrichment of DEGs in CAVD were remarkably observed in the chemokine signaling pathway, cytokine-cytokine receptor interaction, and PI3K-Akt signaling pathway. Chemokines CXCL13, CCL19, CCL8, CXCL8, CXCL16, MMP9, CCL18, CXCL5, VCAM1, and PPBP were identified as the hub genes of CAVD. It was macrophages that accounted for the maximal proportion among these immune cells. The expression of macrophages M0, B cells memory, and Plasma cells were higher in the CAVD valves than in healthy valves, however, the expression of B cells naïve, NK cells activated, and macrophages M2 were lower.

Conclusion: We detected that chemokines CXCL13, CXCL8, CXCL16, and CXCL5, and CCL19, CCL8, and CCL18 are the most important markers of aortic valve disease. The regulatory macrophages M0, plasma cells, B cells memory, B cells naïve, NK cells activated, and macrophages M2 are probably related to the occurrence and the advancement of aortic valve stenosis. These identified chemokines and these immune cells may interact with a subtle adjustment relationship in the development of calcification in CAVD.

Osteomodulin is a Potential Genetic Target for Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is one of the most common genetic heart diseases. Its features include abnormal cardiomyocyte hypertrophy, microvascular dysfunction, and increased accumulation of intercellular matrix. We aim to unravel genes associated with the pathogenesis of HCM and provide a potential target for diagnosis and treatment. Key modules were identified by weighted gene co-expression network analysis (WGCNA). A miRNA-mRNA network was constructed with the predicted miRNA and the most likely hub gene was screened out for gene set enrichment analysis (GSEA). The diagnostic capacity of hub gene was verified by receiver operating characteristic (ROC) curves. Single-cell sequencing (sc-RNA seq) data of normal adult hearts were used to further predict the specific cell types expressing the hub gene. WGCNA assigned genes into different modules and found that the genes contained in the red module had the strongest positive correlation with HCM disease. 2.5% of the genes were common between DEG and hub genes. With the miRNA-mRNA network, osteomodulin (OMD) was identified as the most potential hub gene. GSEA showed that OMD was mainly involved in the synthesis of extracellular matrix and had a certain inhibitory effect on the immune system. The expression of OMD in HCM was validated and ROC curve analysis showed that OMD could distinguish HCM from controls with the area under the curve (AUC) > 0.7. The sc-RNA seq revealed that OMD was mainly expressed in the later stages of cardiac fibroblasts, suggesting that OMD may have an effect on fibroblasts, participating in the pathogenesis of HCM. OMD may serve as a biomarker and therapeutic target for HCM in the future.

Overexpression of TICRR and PPIF confer poor prognosis in endometrial cancer identified by gene co-expression network analysis

The incidence of endometrial cancer (EC) is intensively increasing. However, due to the complexity and heterogeneity of EC, the molecular targeted therapy is still limited. The reliable and accurate biomarkers for tumor progression are urgently demanded. After normalizing the data from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA), we utilized limma and WGCNA packages to identify differentially expressed genes (DEGs). The copy number variations of candidate genes were investigated by cBioPortal. Enrichment pathways analysis was performed by ClueGO and CluePedia. The methylation status was explored by UALCAN. ROC curve and survival analysis were conducted by SPSS and Kaplan–Meier. Infiltration immune cells in microenvironment were analyzed by TISIDB. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA) were applied to explore potential biological pathways. Immunohistochemistry staining (IHC), cell proliferation, cell apoptosis, colony formation, migration, invasion and scratch-wound assays were performed to investigate the function of key genes in vitro. In this study, four expression profile datasets were integrated to identify candidate genes. Combined with WGCNA analysis, the top ten candidates were screened out, whose abnormal methylation patterns were extremely correlated with their expression level and they were associated with tumor grades and predicted poor survival. GSEA and GSVA demonstrated they were involved in DNA replication and cell cycle transition in EC. Gene silencing of TICRR and PPIF dramatically inhibited cell growth, migration and epithelial-mesenchymal transition (EMT) and enhanced progesterone sensitivity. Additionally, from DrugBank database, cyclosporine may be effective for PPIF targeted therapy. By integrative bioinformatics analysis and in vitro experiments, our study shed novel light on the molecular mechanisms of EC. TICRR and PPIF may promise to be potential therapeutic targets for endometrial cancer.

Coexpression of PBX1 and EMP2 as Prognostic Biomarkers in Estrogen Receptor-Negative Breast Cancer via Data Mining

Previous studies revealed that PBX1 ranked the third in the differentially expressed genes about development and progression of breast cancer (BC). Nevertheless, the role of PBX1 contributing to progression of BC has been unevaluated. Here, on the basis of ONCOMINE and GOBO databases, we compared BC samples with normal controls about the expression of PBX1 in various types of cancers, as well as their related expression levels in cancer cell lines by Cancer Cell Line Encyclopedia (CCLE) analysis. It was also found that, when compared with normal controls, PBX1 was markedly higher expressed not only in BC samples but also in BC cell lines, and coexpressed with EMP2 by ONCOMINE and CCLE coexpression analysis, which was also expressed higher in BC samples and BC cell lines similarly. According to Kaplan-Meier plotter, we further explored the prognostic functions of PBX1 and EMP2 in different molecular subtypes of BC, respectively. We demonstrated that overexpression of PBX1 mRNA was correlated with worse survival in luminal B subtype BC, whereas increased EMP2 expression was associated with shorter relapse-free survival in estrogen receptor (ER)-negative patients. Combining with previous studies, we could make a conclusion that coexpression of PBX1 and EMP2 predicts poor prognosis in ER-negative BC, which could be effective biomarkers for BC.

Suppression of CXCL-1 Could Restore Necroptotic Pathway in Chronic Lymphocytic Leukemia

Purpose

To clarify the role of different cytokines and selenite in the defective necroptotic pathway of chronic lymphocytic leukemia (CLL).

Patients and Methods

We randomly collected the peripheral blood samples of 11 untreated CLL patients and 10 healthy volunteers, and then separated B lymphocytes from peripheral blood. Then, real-time polymerase chain reaction (PCR), enzyme-linked immunosorbent assay (ELISA) and Western Blot were performed to detect the expression of different cytokines, including CXC-motif chemokine ligand 1 (CXCL-1). Finally, we used flow cytometry to analyze the percentage of surviving cells to figure out whether CLL cells or normal B lymphocytes underwent necroptosis.

Results

1) The high expression of CXCL-1 was seen in CLL cells compared with normal B lymphocytes (p = 0.0001, adjusted p =0.0012); 2) The downregulation of CXCL-1 was shown in normal B lymphocytes after induction by TNF-α and z-VAD; 3) CLL cells could restore necroptosis induced by TNF-α and z-VAD after knockdown of CXCL-1; 4) The transcriptional and translational expression of LEF-1 were downregulated after the knockdown of CXCL-1 in CLL cells; 5. 3.2μM selenite could help CLL cells restore necroptosis (p = 0.0102) and inhibit the transcriptional and translational expression of CXCL-1.

Conclusion

CXCL-1 played an important role in the defective necroptosis of CLL cells and regulated the expression of LEF-1. Selenite could inhibit the expression of CXCL-1 and help CLL cells restore necroptosis together with TNF-α and z-VAD. Selenite might be the potential medication of CLL in the future.

HERC Ubiquitin Ligases in Cancer

HERC proteins are ubiquitin E3 ligases of the HECT family. The HERC subfamily is composed of six members classified by size into large (HERC1 and HERC2) and small (HERC3-HERC6). HERC family ubiquitin ligases regulate important cellular processes, such as neurodevelopment, DNA damage response, cell proliferation, cell migration, and immune responses. Accumulating evidence also shows that this family plays critical roles in cancer. In this review, we provide an integrated view of the role of these ligases in cancer, highlighting their bivalent functions as either oncogenes or tumor suppressors, depending on the tumor type. We include a discussion of both the molecular mechanisms involved and the potential therapeutic strategies.

Identification of Hub Genes as Biomarkers Correlated with the Proliferation and Prognosis in Lung Cancer: A Weighted Gene Co-Expression Network Analysis

Lung cancer is one of the most malignant tumors in the world. Early diagnosis and treatment of lung cancer are vitally important to reduce the mortality of lung cancer patients. In the present study, we attempt to identify the candidate biomarkers for lung cancer by weighted gene co-expression network analysis (WGCNA). Gene expression profile of GSE30219 was downloaded from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) were analyzed by the limma package, and the co-expression modules of genes were built by WGCNA. UALCAN was used to analyze the relative expression of normal group and tumor subgroups based on tumor individual cancer stages. Survival analysis for the hub genes was performed by Kaplan–Meier plotter analysis with the TCGA database. A total of 2176 genes (745 upregulated and 1431 downregulated genes) were obtained from the GSE30219 database. Seven gene co-expression modules were conducted by WGCNA and the blue module might be inferred as the most crucial module in the pathogenesis of lung cancer. In the pathway enrichment analysis of KEGG, the candidate genes were enriched in the “DNA replication,” “Cell cycle,” and “P53 signaling pathway” pathways. Among these, the cell cycle pathway was the most significant pathway in the blue module with four hub genes CCNB1, CCNE2, MCM7, and PCNA which were selected in our study. Kaplan–Meier plotter analysis indicated that the high expressions of four hub genes were correlated with a worse overall survival (OS) and advanced tumors. qRT-PCR showed that mRNA expression levels of MCM7 (p = 0.038) and CCNE2 (0.003) were significantly higher in patients with the TNM stage. In summary, the high expression of the MCM7 and CCNE2 were significantly related with advanced tumors and worse OS in lung cancer. Thus, the MCM7 and CCNE2 genes can be good indicators for cellular proliferation and prognosis in lung cancer.

Identification of hub methylated-CpG sites and associated genes in oral squamous cell carcinoma

To improve personalized diagnosis and prognosis for oral squamous cell carcinoma (OSCC) by identification of hub methylated‐CpG sites and associated genes, weighted gene comethylation network analysis (WGCNA) was performed to examine and identify hub modules and CpG sites correlated with OSCC. Here, WGCNA modeling yielded blue and brown comethylation modules that were significantly associated with OSCC status. Following screening of the differentially expressed genes (DEGs) from gene expression microarrays and differentially methylated‐CpG sites (DCGs), integrated multiomics analysis of the DEGs, DCGs, and hub CpG sites from the modules was performed to investigate their correlations. Expression levels of 16 CpG sites‐associated genes were negatively correlated with methylation patterns of promoter. Moreover, Kaplan‐Meier survival analysis of the hub CpG sites and associated genes was carried out using 2 public databases, MethSurv and GEPIA. Only 5 genes, ACTA1, ACTN2, OSR1, SYNGR1, and ZNF677, had significant overall survival using GEPIA. Hypermethylated‐CpG sites ACTN2‐cg21376883 and OSR1‐cg06509239 were found to be associated with poor survival by MethSurv. Methylation status of specific site and expression levels of associated genes were determined using clinical samples by quantitative methylation‐specific PCR and real‐time PCR. Pearson's correlation analysis showed that methylation levels of cg06509239 and cg18335068 were negatively related to OSR1 and ZNF677 expression levels, respectively. Our classification schema using multiomics analysis represents a screening framework for identification of hub CpG sites and associated genes.

Immune-related key gene CLDN10 correlates with lymph node metastasis but predicts favorable prognosis in papillary thyroid carcinoma

The functions of immune cells in lymph node metastasis (LNM) have attracted considerable attention. This study aimed to screen the key immune-related and LNM-related genes in PTC. In the discovery phase, the immune-related genes in LNM were screened by using bioinformatics methods. In the validation phases, the association of the genes with LNM was first confirmed in a cohort from The Cancer Genome Atlas and a cohort based on a tissue chip. Then, the relationship of the genes with immune cell infiltration was further explored. Consequently, CLDN10 was identified, and its high expression was correlated with the presence of LNM in PTC but predicted a favorable prognosis. High CLDN10 expression was positively correlated with the infiltration of several immune cells, such as B cells, CD8+T cells, and macrophages. High CLDN10 expression may improve the outcomes of patients with PTC by increasing immune cell infiltration, although it might be associated with LNM. In conclusion, although CLDN1 might be correlated with LNM, it may also increase the infiltration of immune cells, including CD8+T cells and macrophages, and improve the clinical outcomes of patients with PTC. The effects of tumor purity and immune cell infiltration need to be considered in prognosis evaluation.

Weighted gene co-expression network analysis identifies CCNA2 as a treatment target of prostate cancer through inhibiting cell cycle

Prostate cancer is a malignant tumor disease that seriously harms the lives of middle-aged and elderly men. Weighted gene co-expression analysis can be used to construct gene co-expression networks to explore gene sets and genes that are significantly correlated with clinical features. In this study, the transcriptome data of prostate cancer on TCGA was analyzed by weighted gene co-expression network, and the gene with a significant correlation with disease Gleason stage and tumor T stage was identified: CCNA2. CCNA2 was significantly associated with biochemical recurrence, disease-free survival and overall survival rate of prostate cancer. The ability of cancer cell proliferation, invasion and metastasis was decreased after down-regulated expression of CCNA2 in prostate cancer cell lines. Flow cytometry revealed that tumor cells were arrested in the S phase after down-regulated the expression of CCNA2. Taken together, we used WGCNA and obtain a gene CCNA2 which is significantly associated with the prognosis of prostate cancer, which may be an indicator of the prognosis of prostate cancer and a new therapeutic target.