ST13, a proliferation regulator, inhibits growth and migration of colorectal cancer cell lines

Transcriptomics analysis reveals distinct mechanism of breast cancer stem cells regulation in mammospheres from MCF-7 and T47D cells

Luminal A breast cancer, constituting 70 % of breast cancer cases, presents a challenge due to the development of resistance and recurrence caused by breast cancer stem cells (BCSC). Luminal breast tumors are characterized by TP53 expression, a tumor suppressor gene involved in maintaining stem cell attributes in cancer. Although a previous study successfully developed mammospheres (MS) from MCF-7 (with wild-type TP53) and T47D (with mutant TP53) luminal breast cancer cells for BCSC enrichment, their transcriptomic profiles remain unclear. We aimed to elucidate the transcriptomic disparities between MS of MCF-7 and T47D cells using bioinformatics analyses of differentially expressed genes (DEGs), including the KEGG pathway, Gene Ontology (GO), drug-gene association, disease-gene association, Gene Set Enrichment Analysis (GSEA), DNA methylation analysis, correlation analysis of DEGs with immune cell infiltration, and association analysis of genes and small-molecule compounds via the Connectivity Map (CMap). Upregulated DEGs were enriched in metabolism-related KEGG pathways, whereas downregulated DEGs were enriched in the MAPK signaling pathway. Drug-gene association analysis revealed that both upregulated and downregulated DEGs were associated with fostamatinib. The KEGG pathway GSEA results indicated that the DEGs were enriched for oxidative phosphorylation, whereas the downregulated DEGs were negatively enriched for the p53 signaling pathway. Examination of DNA methylation revealed a noticeable disparity in the expression patterns of the PKM2, ERO1L, SLC6A6, EPAS1, APLP2, RPL10L, and NEDD4 genes when comparing cohorts with low- and high-risk breast cancer. Furthermore, a significant positive correlation was identified between SLC6A6 expression and macrophage presence, as well as MSN, and AKR1B1 expression and neutrophil and dentritic cell infiltration. CMap analysis unveiled SA-83851 as a potential candidate to counteract the effects of DEGs, specifically in cells harbouring mutant TP53. Further research, including in vitro and in vivo validations, is warranted to develop drugs targeting BCSCs.

T cell-related prognostic risk model and tumor immune environment modulation in lung adenocarcinoma based on single-cell and bulk RNA sequencing

BACKGROUND

T cells are present in all stages of tumor formation and play an important role in the tumor microenvironment. We aimed to explore the expression profile of T cell marker genes, constructed a prognostic risk model based on these genes in Lung adenocarcinoma (LUAD), and investigated the link between this risk model and the immunotherapy response.

METHODS

We obtained the single-cell sequencing data of LUAD from the literature, and screened out 6 tissue biopsy samples, including 32,108 cells from patients with non-small cell lung cancer, to identify T cell marker genes in LUAD. Combined with TCGA database, a prognostic risk model based on T-cell marker gene was constructed, and the data from GEO database was used for verification. We also investigated the association between this risk model and immunotherapy response.

RESULTS

Based on scRNA-seq data 1839 T-cell marker genes were identified, after which a risk model consisting of 9 gene signatures for prognosis was constructed in combination with the TCGA dataset. This risk model divided patients into high-risk and low-risk groups based on overall survival. The multivariate analysis demonstrated that the risk model was an independent prognostic factor. Analysis of immune profiles showed that high-risk groups presented discriminative immune-cell infiltrations and immune-suppressive states. Risk scores of the model were closely correlated with Linoleic acid metabolism, intestinal immune network for IgA production and drug metabolism cytochrome P450.

CONCLUSION

Our study proposed a novel prognostic risk model based on T cell marker genes for LUAD patients. The survival of LUAD patients as well as treatment outcomes may be accurately predicted by the prognostic risk model, and make the high-risk population present different immune cell infiltration and immunosuppression state.

Identification of blood-based key biomarker and immune infiltration in Immunoglobulin A nephropathy by comprehensive bioinformatics analysis and a cohort validation

Background

To identify the critical genes in the onset and progression of Immunoglobulin A nephropathy (IgAN) and to explore its immune cell infiltration feature.

Methods

Differentially expressed genes (DEGs) were firstly screened from 1 blood-derived dataset GSE73953 and a glomerulus derived dataset GSE93798 through limma analysis, overlap genes omitting and weighted gene correlation network analysis (WGCNA) and further reduced according to expression pattern and correlation with the clinical features: eGFR and proteinuria, followed by external validation using the GSE37460 dataset and an IgAN cohort. In addition, the CIBERSORT tool for immune cell infiltration analysis, ceRNA network construction and Connectivity Map (CMAP) were also performed.

Results

A total of 195 DEGs were found, and among them, 3 upregulated (ORMDL2, NRP1, and COL4A1) and 3 downregulated genes (ST13, HSPA8 and PKP4) are verified to correlate clinically, and finally ORMDL2, NRP1 and COL4A1 were validated in patient cohort and with the ability of IgAN discrimination (highest AUC was COL4A1: 97.14%). The immune cell infiltration results revealed that significant differences could be found on resting memory CD4 T cells, activated NK cells, and M2 macrophages between control and IgAN.

Conclusions

Our results demonstrated here that significantly upregulated DEGs: ORMDL2, NRP1 and COL4A1, could be served as the diagnostic marker for IgAN, and dysregulated immune cell infiltration hinted possible the immune system intervention point in the setting of IgAN.

Oncolytic Adenovirus Expressing ST13 Increases Antitumor Effect of Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand Against Pancreatic Ductal Adenocarcinoma

Oncolytic adenoviruses (OAds) are promising agents for cancer therapy, representing a novel therapeutic strategy for pancreatic ductal adenocarcinoma (PDAC). However, there are challenges associated with the successful use of an OAd alone, involving the security of the viral vector and screening of an effective antitumor gene. In the present study, a novel OAd CD55-ST13-tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) was constructed in which the dual therapeutic genes ST13 and TRAIL were inserted, featuring the carcinoembryonic antigen (CEA) as a promoter to control E1A and deletion of the 55 kDa E1B gene. ST13, known as a colorectal cancer suppressor gene, exhibited lower expression in PDAC than in tumor-adjacent tissues and was associated with poor prognosis in PDAC patients. In vitro studies demonstrated that CD55-ST13-TRAIL was effective in promoting the expression of ST13 and TRAIL in CEA-positive pancreatic cancer cells. Moreover, CD55-ST13-TRAIL exhibited a synergistic effect toward tumor cell death compared with CD55-ST13 alone or CD55-TRAIL alone, and inhibited tumor cell proliferation and induced cell apoptosis dependent on caspase pathways in PDAC cells. Furthermore, xenograft experiments in a mouse model indicated that CD55-ST13-TRAIL significantly inhibited tumor growth and improved the survival of animals with xenografts. The findings demonstrate that oncolytic virotherapy under the control of the promoter CEA enables safe and efficient treatment of PDAC, and suggest that it represents a promising candidate for the treatment of metastatic diseases.

Differential proteomic analysis of fetal and geriatric lumbar nucleus pulposus: immunoinflammation and age-related intervertebral disc degeneration

BACKGROUND

Intervertebral disc degeneration (IVDD) is a major cause of low back pain. Although the mechanism of degeneration remains unclear, aging has been recognized as a key risk factor for IVDD. Most studies seeking to identify IVDD-associated molecular alterations in the context of human age-related IVDD have focused only on a limited number of proteins. Differential proteomic analysis is an ideal method for comprehensively screening altered protein profiles and identifying the potential pathways related to pathological processes such as disc degeneration.

METHODS

In this study, tandem mass tag (TMT) labeling was combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) for differential proteomic analysis of human fetal and geriatric lumbar disc nucleus pulposus (NP) tissue. Parallel reaction monitoring (PRM) and Western blotting (WB) techniques were used to identify target proteins. Bioinformatic analyses, including Gene Ontology (GO) annotation, domain annotation, pathway annotation, subcellular localization and functional enrichment analyses, were used to interpret the potential significance of the protein alterations in the mechanism of IVDD. Student's t-tests and two-tailed Fisher's exact tests were used for statistical analysis.

RESULTS

Six hundred forty five proteins were significantly upregulated and 748 proteins were downregulated in the geriatric group compared with the fetal group. Twelve proteins were verified to have significant differences in abundance between geriatric and fetal NP tissue; most of these have not been previously identified as being associated with human IVDD. The potential significance of the differentially expressed proteins in age-related IVDD was analyzed from multiple perspectives, especially with regard to the association of the immunoinflammatory response with IVDD.

CONCLUSIONS

Differential proteomic analysis was used as a comprehensive strategy for elucidating the protein alterations associated with age-related IVDD. The findings of this study will aid in the screening of new biomarkers and molecular targets for the diagnosis and therapy of IVDD. The results may also significantly enhance our understanding of the pathophysiological process and mechanism of age-related IVDD.

Construction and Validation of an Autophagy-Related Prognostic Risk Signature for Survival Predicting in Clear Cell Renal Cell Carcinoma Patients

Background: Clear cell renal cell carcinoma (ccRCC) is a common type of malignant tumors in urinary system. Evaluating the prognostic outcome at the time of initial diagnosis is essential for patients. Autophagy is known to play a significant role in tumors. Here, we attempted to construct an autophagy-related prognostic risk signature based on the expression profile of autophagy-related genes (ARGs) for predicting the long-term outcome and effect of precise treatments for ccRCC patients.

Methods: We obtained the expression profile of ccRCC from the cancer genome atlas (TCGA) database and extract the portion of ARGs. We conducted differentially expressed analysis on ARGs and then performed enrichment analyses to confirm the anomalous autophagy-related biological functions. Then, we performed univariate Cox regression to screen out overall survival (OS)-related ARGs. With these genes, we established an autophagy-related risk signature by least absolute shrinkage and selection operator (LASSO) Cox regression. We validated the reliability of the risk signature with receiver operating characteristic (ROC) analysis, survival analysis, clinic correlation analysis, and Cox regression. Then we analyzed the function of each gene in the signature by single-gene gene set enrichment analysis (GSEA). Finally, we analyzed the correlation between our risk score and expression level of several targets of immunotherapy and targeted therapy.

Results: We established a seven-gene prognostic risk signature, according to which we could divide patients into high or low risk groups and predict their outcomes. ROC analysis and survival analysis validated the reliability of the signature. Clinic correlation analysis found that the risk group is significantly correlated with severity of ccRCC. Multivariate Cox regression revealed that the risk score could act as an independent predictor for the prognosis of ccRCC patients. Correlation analysis between risk score and targets of precise treatments showed that our risk signature could predict the effects of precise treatment powerfully.

Conclusion: Our study provided a brand new autophagy-related seven-gene prognostic risk signature, which could perform as a prognostic indicator for ccRCC. Meanwhile, our study provides a novel sight to understand the role of autophagy and suggest therapeutic strategies in the category of precise treatment in ccRCC.

Quinolinate as a Marker for Kynurenine Metabolite Formation and the Unresolved Question of NAD+ Synthesis During Inflammation and Infection

Quinolinate (Quin) is a classic example of a biochemical double-edged sword, acting as both essential metabolite and potent neurotoxin. Quin is an important metabolite in the kynurenine pathway of tryptophan catabolism leading to the de novo synthesis of nicotinamide adenine dinucleotide (NAD+). As a precursor for NAD+, Quin can direct a portion of tryptophan catabolism toward replenishing cellular NAD+ levels in response to inflammation and infection. Intracellular Quin levels increase dramatically in response to immune stimulation [e.g., lipopolysaccharide (LPS) or pokeweed mitogen (PWM)] in macrophages, microglia, dendritic cells, and other cells of the immune system. NAD+ serves numerous functions including energy production, the poly ADP ribose polymerization (PARP) reaction involved in DNA repair, and the activity of various enzymes such as the NAD+-dependent deacetylases known as sirtuins. We used highly specific antibodies to protein-coupled Quin to delineate cells that accumulate Quin as a key aspect of the response to immune stimulation and infection. Here, we describe Quin staining in the brain, spleen, and liver after LPS administration to the brain or systemic PWM administration. Quin expression was strong in immune cells in the periphery after both treatments, whereas very limited Quin expression was observed in the brain even after direct LPS injection. Immunoreactive cells exhibited diverse morphology ranging from foam cells to cells with membrane extensions related to cell motility. We also examined protein expression changes in the spleen after kynurenine administration. Acute (8 h) and prolonged (48 h) kynurenine administration led to significant changes in protein expression in the spleen, including multiple changes involved with cytoskeletal rearrangements associated with cell motility. Kynurenine administration resulted in several expression level changes in proteins associated with heat shock protein 90 (HSP90), a chaperone for the aryl-hydrocarbon receptor (AHR), which is the primary kynurenine metabolite receptor. We propose that cells with high levels of Quin are those that are currently releasing kynurenine pathway metabolites as well as accumulating Quin for sustained NAD+ synthesis from tryptophan. Further, we propose that the kynurenine pathway may be linked to the regulation of cell motility in immune and cancer cells.

A Framework for Identifying Genotypic Information from Clinical Records: Exploiting Integrated Ontology Structures to Transfer Annotations between ICD Codes and Gene Ontologies

Although some methods are proposed for automatic ontology generation, none of them address the issue of integrating large-scale heterogeneous biomedical ontologies. We propose a novel approach for integrating various types of ontologies efficiently and apply it to integrate International Classification of Diseases, Ninth Revision, Clinical Modification (ICD9CM), and Gene Ontologies. This approach is one of the early attempts to quantify the associations among clinical terms (e.g., ICD9 codes) based on their corresponding genomic relationships. We reconstructed a merged tree for a partial set of GO and ICD9 codes and measured the performance of this tree in terms of associations' relevance by comparing them with two well-known disease-gene datasets (i.e., MalaCards and Disease Ontology). Furthermore, we compared the genomic-based ICD9 associations to temporal relationships between them from electronic health records. Our analysis shows promising associations supported by both comparisons suggesting a high reliability. We also manually analyzed several significant associations and found promising support from literature.

New Phage Display-Isolated Heptapeptide Recognizing the Regulatory Carboxy-Terminal Domain of Human Tumour Protein p53

The transcription factor tumor protein p53 (P53) controls a variety of genes most involved in cell cycle and is at the origin of apoptosis when DNA is irreparably damaged. We planned to select novel tumor protein p53-interacting peptides through the screening of hepta-peptide phage-display libraries. For this aim, human tumor suppressor protein p53 was expressed in Escherichia coli as Glutathione S-transferase fusion and purified by affinity chromatography. The phage library was then screened on this immobilized protein target. After three rounds of panning, phages were sequenced and shown to contain a consensus sequence NPNSAQG. Thereafter, either free p53 liberated from the fusion protein through thrombin treatment or Histidine-tagged p53 were recognized efficiently by the selected phage. To locate the p53-binding epitope of the selected hepta-peptide, three long peptides parts of the three known domains of the protein were synthesized and screened by the selected phage/peptide. Thus, the Carboxy-terminal p53 region was shown to be the target of the isolated phage as well as by its derived Fluorescein isothiocyanate-peptide. Molecular docking showed Lysine 386 as an important residue potentially engaged in this interaction. The selected hepta-peptide is a novel p53-interacting peptide, not described by other studies, and could be used as therapeutic tool in the future.

Exosomal microRNA in plasma as a non-invasive biomarker for the recurrence of non-small cell lung cancer

Predictive biomarkers for the recurrence of non-small cell lung cancer (NSCLC) in patients who have received curative resection are important for cancer treatment. The functional microRNAs (miRNAs/miRs) in the exosomes of plasma and serum samples are of interest as stable and non-invasive biomarkers for recurrence in cancer patients. The aim of the present study was to clarify the usefulness of plasma exosomal miRNAs as biomarkers for the prediction of recurrence in NSCLC following curative resection. First, microarray-based expression profiling of miRNAs derived from exosomes in the plasma of 6 patients was employed to identify a biomarker that distinguishes between patients with and without NSCLC recurrence. In the miRNA microarray analyses, the exosomal miR-21 and miR-4257 levels of the NSCLC patients showed marked upregulation in those individuals with recurrence compared with those without recurrence and healthy individuals. These two miRNAs were thus selected as recurrence-specific biomarkers and their potential was evaluated in a separate cohort of 195 NSCLC patients. In comparison to the levels in 30 healthy individuals, exosomal miR-21 and miR-4257 levels showed a significant increase in the NSCLC patients (P<0.01). When evaluating the clinicopathological significance of these miRNAs, exosomal miR-21 showed a significant association with tumor size and tumor-node-metastasis (TNM) stage (P<0.05). Exosomal miR-4257 showed a significant association with histological type, lymphatic invasion and TNM stage (P<0.05). The disease-free survival (DFS) rates of high exosomal miR-21 patients were significantly worse than those of low exosomal miR-21 patients (P<0.05), and the DFS rates of patients with high exosomal miR-4257 levels were significantly worse than those with low exosomal miR-4257 levels (P<0.01). In the Cox multivariate analysis, plasma exosomal miR-21 and miR-4257 expression showed a significance association with DFS (P<0.05). These results suggest that plasma exosomal miR-21 and mir-4257 expression has potential as a predictive biomarker for recurrence in NSCLC patients who have received curative resection.

Molecular Analysis of Pediatric Oligodendrogliomas Highlights Genetic Differences with Adult Counterparts and Other Pediatric Gliomas

Oligodendroglioma represents a distinctive neoplasm in adults but similar neoplasms occur rarely in children. We studied 20 cases of pediatric oligodendroglioma by SNP array (median age 9 years, range 1-19; 15 grade II and 5 grade III). Cytogenetic abnormalities were present in 8 (53%) grade II and all five anaplastic oligodendrogliomas. Most changes were in the form of deletion and copy neutral loss of heterozygosity (LOH). The most common abnormality was 1p deletion (n = 5). Whole arm 1p19q co-deletion was present in three cases from adolescent patients and 9p loss in 3, including one low-grade oligodendroglioma with CDKN2A homozygous deletion. Common losses were largely limited to the anaplastic subset (n = 5) and included 3q29 (n = 3), 11p (n = 3), 17q (n = 3), 4q (n = 2), 6p (n = 2), 13q (n = 2), 14q (n = 2), 17p (n = 2) and whole Ch 18 loss (n = 2). Gains were non-recurrent except for whole Ch 7 (n = 2) and gain on 12q (n = 2) including the MDM2 locus. Possible germ line LOH (or uniparental disomy) was present in seven cases (35%), with one focal abnormality (22q13.1-13.2) in two. BRAF-KIAA1549 fusions and BRAF p.V600E mutations were absent (n = 13 and 8). In summary, cytogenetic alterations in pediatric oligodendrogliomas are characterized mostly by genomic losses, particularly in anaplastic tumors.

A Lotus japonicus Cochaperone Protein Interacts With the Ubiquitin-Like Domain Protein CIP73 and Plays a Negative Regulatory Role in Nodulation

The calcium/calmodulin-dependent protein kinase CCaMK forms a complex with its phosphorylation target CIP73 (CCaMK-interacting protein of 73 kDa). In this work, a homolog of the animal HSC/HSP70 interacting protein (HIP) was identified as an interacting partner of CIP73 in Lotus japonicus. L. japonicus HIP contains all functional domains characteristic of animal HIP proteins. The C-terminal STI1-like domain of L. japonicus HIP was found to be necessary and sufficient for interaction with CIP73. The interaction between CIP73 and HIP occurred in both the nuclei and cytoplasm in Nicotiana benthamiana leaf cells. The interactions between CIP73 and HIP and between CIP73 and CCaMK could take place simultaneously in the same nuclei. HIP transcripts were detected in all plant tissues tested. As nodule primordia developed into young nodules, the expression of HIP was down-regulated and the HIP transcript level became very low in mature nodules. More nodules were formed in transgenic hairy roots of L. japonicus expressing HIP RNA interference at 16 days postinoculation as compared with the control hairy roots expressing the empty vector. It appears that HIP may play a role as a negative regulator for nodulation.