ZC3HAV1 promotes the proliferation and metastasis via regulating KRAS in pancreatic cancer

Heyingwuzi formulation alleviates diabetic retinopathy by promoting mitophagy via the HIF-1α/BNIP3/NIX axis

BACKGROUND

Diabetic retinopathy (DR) is the primary cause of visual problems in patients with diabetes. The Heyingwuzi formulation (HYWZF) is effective against DR.

AIM

To determine the HYWZF prevention mechanisms, especially those underlying mitophagy.

METHODS

Human retinal capillary endothelial cells (HRCECs) were treated with high glucose (hg), HYWZF serum, PX-478, or Mdivi-1 in vitro. Then, cell counting kit-8, transwell, and tube formation assays were used to evaluate HRCEC proliferation, invasion, and tube formation, respectively. Transmission electron microscopy was used to assess mitochondrial morphology, and Western blotting was used to determine the protein levels. Flow cytometry was used to assess cell apoptosis, reactive oxygen species (ROS) production, and mitochondrial membrane potential. Moreover, C57BL/6 mice were established in vivo using streptozotocin and treated with HYWZF for four weeks. Blood glucose levels and body weight were monitored continuously. Changes in retinal characteristics were evaluated using hematoxylin and eosin, tar violet, and periodic acid-Schiff staining. Protein levels in retinal tissues were determined via Western blotting, immunohistochemistry, and immunostaining.

RESULTS

HYWZF inhibited excessive ROS production, apoptosis, tube formation, and invasion in hg-induced HRCECs via mitochondrial autophagy in vitro. It increased the mRNA expression levels of BCL2-interacting protein 3 (BNIP3), FUN14 domain-containing 1, BNIP3-like (BNIP3L, also known as NIX), PARKIN, PTEN-induced kinase 1, and hypoxia-inducible factor (HIF)-1α. Moreover, it downregulated the protein levels of vascular endothelial cell growth factor and increased the light chain 3-II/I ratio. However, PX-478 and Mdivi-1 reversed these effects. Additionally, PX-478 and Mdivi-1 rescued the effects of HYWZF by decreasing oxidative stress and apoptosis and increasing mitophagy. HYWZF intervention improved the symptoms of diabetes, tissue damage, number of acellular capillaries, and oxidative stress in vivo. Furthermore, in vivo experiments confirmed the results of in vitro experiments.

CONCLUSION

HYWZF alleviated DR and associated damage by promoting mitophagy via the HIF-1α/BNIP3/NIX axis.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

BACKGROUND

Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors.

METHODS

We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers.

RESULTS

We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors.

CONCLUSIONS

This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Multi-System-Level Analysis with RNA-Seq on Pterygium Inflammation Discovers Association between Inflammatory Responses, Oxidative Stress, and Oxidative Phosphorylation

A pterygium is a common conjunctival degeneration and inflammatory condition. It grows onto the corneal surface or limbus, causing blurred vision and cosmetic issues. Ultraviolet is a well-known risk factor for the development of a pterygium, although its pathogenesis remains unclear, with only limited understanding of its hereditary basis. In this study, we collected RNA-seq from both pterygial tissues and conjunctival tissues (as controls) from six patients (a total of twelve biological samples) and retrieved publicly available data, including eight pterygium samples and eight controls. We investigated the intrinsic gene regulatory mechanisms closely linked to the inflammatory reactions of pterygiums and compared Asian (Korea) and the European (Germany) pterygiums using multiple analysis approaches from different perspectives. The increased expression of antioxidant genes in response to oxidative stress and DNA damage implies an association between these factors and pterygium development. Also, our comparative analysis revealed both similarities and differences between Asian and European pterygiums. The decrease in gene expressions involved in the three primary inflammatory signaling pathways-JAK/STAT, MAPK, and NF-kappa B signaling-suggests a connection between pathway dysfunction and pterygium development. We also observed relatively higher activity of autophagy and antioxidants in the Asian group, while the European group exhibited more pronounced stress responses against oxidative stress. These differences could potentially be necessitated by energy-associated pathways, specifically oxidative phosphorylation.

Proximal tubule cells in blood and urine as potential biomarkers for kidney disease biopsy

Early diagnosis and treatment are crucial for managing kidney disease, yet there remains a need to further explore pathological mechanisms and develop minimally invasive diagnostic methods. In this study, we employed single-cell RNA sequencing (scRNA-seq) to assess the cellular heterogeneity of kidney diseases. We analyzed gene expression profiles from renal tissue, peripheral blood mononuclear cells (PBMCs), and urine of four patients with nephritis. Our findings identified 12 distinct cell subsets in renal tissues and leukocytes. These subsets encompassed fibroblast cells, mesangial cells, epithelial cells, proximal tubule cells (PTCs), and six immune cell types: CD8+ T cells, macrophages, natural killer cells, dendritic cells, B cells, and neutrophils. Interestingly, PTCs were present in both PBMCs and urine samples but absent in healthy blood samples. Furthermore, several populations of fibroblast cells, mesangial cells, and PTCs exhibited pro-inflammatory or pro-apoptotic behaviors. Our gene expression analysis highlighted the critical role of inflammatory PTCs and fibroblasts in nephritis development and progression. These cells showed high expression of pro-inflammatory genes, which could have chemotactic and activating effect on neutrophils. This was substantiated by the widespread in these cells. Notably, the gene expression profiles of inflammatory PTCs in PBMCs, urine, and kidney tissues had high similarity. This suggests that PTCs in urine and PBMCs hold significant potential as alternative markers to invasive kidney biopsies.

Identification of CCCH-type zinc finger antiviral protein 1 (ZAP) gene from Pacific white shrimp (Penaeus vannamei): Characterization and expression analysis in response to viral infection

Zinc-finger proteins (ZFPs) are a huge family that exert multiple roles in the cells. ZFPs could be divided into nine types based on the numbers and positions of conserved Cys and His residues, in which CCCH-type ZFP was one of the most widely studied types. CCCH-type zinc finger antiviral protein 1 (ZAP), a CCCH-type ZFP that can inhibit the replication of certain RNA viruses and DNA viruses by mediating degradation of viral RNA and repressing mRNA translation, plays significant roles in the host innate immune defenses against viral infections. Presently, there have been numerous reports investigating the antiviral ability of ZAP, while no data is available about ZAP gene in the species of shrimps or even crustaceans. In this study, a novel protein containing CCCH-type zinc finger motifs (ZnF-CCCH), CCCH-type zinc finger antiviral protein 1 (ZAP) gene, was identified from Pacific white shrimp (Penaeus vannamei) and its role in antiviral immunity was further investigated. Similar to mammalian ZAPs, in addition to ZnF-CCCH, PvZAP also possesses central WWE domains and C-terminal PARP domain. Phylogenetic analysis showed that PvZAP was close to that of the crustacean Pacific oyster, separating from the cluster of vertebrate ZAP proteins. Upon in vivo infection by IHHNV, gene expression of PvZAP was strongly up-regulated in the hepatopancreas and gills of both adult and juvenile shrimps, where adult individuals showed higher fold changes of up-regulation than in juvenile individuals. These results suggested that PvZAP might play an important role in the innate immune defense of Pacific white shrimp against IHHNV infection. This allows us to gain new insights into the immunological function of ZAP in the innate immunity of shrimp species and even crustaceans.

Genomic Association Analysis of Growth and Backfat Traits in Large White Pigs

The pig industry is significantly influenced by complex traits such as growth rate and fat deposition, which have substantial implications for economic returns. Over the years, remarkable genetic advancements have been achieved through intense artificial selection to enhance these traits in pigs. In this study, we aimed to investigate the genetic factors that contribute to growth efficiency and lean meat percentages in Large White pigs. Specifically, we focused on analyzing two key traits: age at 100 kg live weight (AGE100) and backfat thickness at 100 kg (BF100), in three distinct Large White pig populations-500 Canadian, 295 Danish, and 1500 American Large White pigs. By employing population genomic techniques, we observed significant population stratification among these pig populations. Utilizing imputed whole-genome sequencing data, we conducted single population genome-wide association studies (GWAS) as well as a combined meta-analysis across the three populations to identify genetic markers associated with the aforementioned traits. Our analyses highlighted several candidate genes, such as CNTN1-which has been linked to weight loss in mice and is potentially influential for AGE100-and MC4R, which is associated with obesity and appetite and may impact both traits. Additionally, we identified other genes-namely, PDZRN4, LIPM, and ANKRD22-which play a partial role in fat growth. Our findings provide valuable insights into the genetic basis of these important traits in Large White pigs, which may inform breeding strategies for improved production efficiency and meat quality.

STING pathway contributes to the prognosis of hepatocellular carcinoma and identification of prognostic gene signatures correlated to tumor microenvironment

Background

Hepatocellular carcinoma (HCC) is one of the most malignant solid tumors worldwide. Recent evidence shows that the stimulator of interferon genes (STING) pathway is essential for anti-tumor immunity via inducing the production of downstream inflammatory cytokines. However, its impact on the prognosis and tumor microenvironment of HCC was still limited known.

Methods

We obtained gene expression profiles of HCC from GEO, TCGA, and ICGC databases, and immune-related genes (IRGs) from the ImmPort database. Multivariate Cox regression was performed to identify independent prognostic factors. Nomogram was established to predict survival probability for individual patients. Kaplan–Meier curve was used to evaluate the survival difference. Afterward, ESTIMATE, TISCH, and TIMER databases were combined to assess the immune cell infiltration. Furthermore, the qPCR, western blotting, and immunohistochemistry were done to evaluate gene expression, and in vitro cell models were built to determine cell migratory ability.

Results

We found that gene markers of NLRC3, STING1, TBK1, TRIM21, and XRCC6 within STING pathway were independent prognostic factors in HCC patients. Underlying the finding, a predictive nomogram was constructed in TCGA-training cohort and further validated in TCGA-all and ICGC datasets, showing credible performance. Experimentally, up-regulated TBK1 promotes the ability of HCC cell migration. Next, the survival-related immune-related co-expressed gene signatures (IRCGS) (VAV1, RHOA, and ZC3HAV1) were determined in HCC cohorts and their expression was verified in human HCC cells and clinical samples. Furthermore, survival-related IRCGS was associated with the infiltration of various immune cell subtypes in HCC, the transcriptional expression of prominent immune checkpoints, and immunotherapeutic response.

Conclusion

Collectively, we constructed a novel prognostic nomogram model for predicting the survival probability of individual HCC patients. Moreover, an immune-related prognostic gene signature was determined. Both might function as potential therapeutic targets for HCC treatment in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-022-02734-4.

Construction of a breast cancer prognosis model based on alternative splicing and immune infiltration

BACKGROUND

Breast cancer (BC) is the most common malignancy among women in the world. Alternative splicing (AS) is an important mechanism for regulating gene expression and producing proteome diversity, which is closely related to tumorigenesis. Understanding the role of AS in BC may be helpful to reveal new therapeutic targets for clinical interventions.

METHODS

RNA-seq, clinical and AS data of TCGA-BRCA were downloaded from TCGA and TCGA SpliceSeq databases. AS events associated with prognosis were filtered by univariate Cox regression. The AS risk model of BC was built by Lasso regression, random forest and multivariate Cox regression. The accuracy of the AS risk model and clinicopathological factors were evaluated by time-dependent receiver operating characteristic (ROC) curves. The significant factors were used to construct the nomogram model. Tumor microenvironment analysis, immune infiltration and immune checkpoint analysis were performed to show the differences between the high and low AS risk groups. The expression differences of genes of AS events constituting the risk model in tumor tissues and normal tissues were analyzed, the genes with significant differences were screened, and their relationship with prognosis, tumor microenvironment, immune infiltration and immune checkpoint were analyzed. Finally, Pearson correlation analysis was used to calculate the correlation coefficient between splicing factors (SF) and prognostic AS events in TCGA-BRCA. The results were imported into Cytoscape, and the associated network was constructed.

RESULTS

A total of 21,232 genes had 45,421 AS events occurring in TCGA-BRCA, while 1604 AS events were found to be significantly correlated with survival. The BRCA risk model consisted of 5 AS events, (TTC39C|44853|AT*- 2.67) + (HSPBP1|52052|AP*- 4.28) + (MAZ|35942|ES*2.34) + (ANK3|11845|AP*1.18) + (ZC3HAV1|81940|AT*1.59), which were confirmed to be valuable for predicting BRCA prognosis to a certain degree, including ROC curve, survival analysis, tumor microenvironment analysis, immune infiltration and immune checkpoint analysis. Based on this, we constructed a nomogram prediction model composed of clinicopathological features and the AS risk signature. Furthermore, we found that MAZ was a core gene indicating the connection of tumor prognosis and AS events. Ultimately, a network of SF-AS regulation was established to reveal the relationship between them.

CONCLUSIONS

We constructed a nomogram model combined with clinicopathological features and AS risk score to predict the prognosis of BC. The detailed analysis of tumor microenvironment and immune infiltration in the AS risk model may further reveal the potential mechanisms of BC recurrence and development.

ZAP isoforms regulate unfolded protein response and epithelial- mesenchymal transition

Human ZAP inhibits many viruses, including HIV and coronaviruses, by binding to viral RNAs to promote their degradation and/or translation suppression. However, the regulatory role of ZAP in host mRNAs is largely unknown. Two major alternatively spliced ZAP isoforms, the constitutively expressed ZAPL and the infection-inducible ZAPS, play overlapping yet different antiviral and other roles that need further characterization. We found that the splicing factors hnRNPA1/A2, PTBP1/2, and U1-snRNP inhibit ZAPS production and demonstrated the feasibility to modulate the ZAPL/S balance by splice-switching antisense oligonucleotides in human cells. Transcriptomic analysis of ZAP-isoform-specific knockout cells revealed uncharacterized host mRNAs targeted by ZAPL/S with broad cellular functions such as unfolded protein response (UPR), epithelial-mesenchymal transition (EMT), and innate immunity. We established that endogenous ZAPL and ZAPS localize to membrane compartments and cytosol, respectively, and that the differential localization correlates with their target-RNA specificity. We showed that the ZAP isoforms regulated different UPR branches under resting and stress conditions and affected cell viability during ER stress. We also provided evidence for a different function of the ZAP isoforms in EMT-related cell migration, with effects that are cell-type dependent. Overall, this study demonstrates that the competition between splicing and IPA is a potential target for the modulation of the ZAPL/S balance, and reports new cellular transcripts and processes regulated by the ZAP isoforms.

Multitranscriptome analyses of keloid fibroblasts reveal the role of the HIF-1α/HOXC6/ERK axis in keloid development

Background

A keloid is a disease of excessive fibrosis that is characterized by the aberrant proliferation of fibroblasts. However, the molecular mechanisms of fibroblasts during the development of keloids remain unclear. This study aims to identify new molecular targets that promote the proliferation and migration of keloid fibroblasts, providing new ideas for the prevention and treatment of keloids.

Methods

We utilized bioinformatics tools to analyze data from keloid fibroblasts (KFs) available in the Gene Expression Omnibus (GEO) database to identify the key genes involved in keloid development. Homeobox C6 (HOXC6) emerged as a hub gene in KFs from the GEO database was verified in keloid tissue samples and KFs using reverse transcription-quantitative polymerase chain reaction, western blot (WB) and immunohistochemistry. Subsequently, the effects of downregulated HOXC6 expression on the cellular behaviors of KFs were examined by performing Cell Counting Kit-8, flow cytometry, transwell migration and WB assays. Meanwhile, we performed transcriptome sequencing and gene set enrichment analysis (GSEA) to further explore HOXC6-related mechanisms and validated the signaling pathways by performing a series of experiments.

Results

HOXC6 was the top-ranking hub gene of KFs in microarray datasets from GEO and was upregulated in keloid tissue samples and KFs. Downregulation of HOXC6 inhibited proliferation, migration and extracellular matrix (ECM) accumulation and promoted KF apoptosis. GSEA predicted that the hypoxia signaling pathway was associated with HOXC6 in KFs. Transcriptome sequencing suggested that the extracellular regulated protein kinase (ERK) pathway was one of the downstream pathways of HOXC6 in KFs. Our experiments confirmed that hypoxia-inducible factor-1α (HIF-1α) upregulates HOXC6, contributing to KFs proliferation, migration, apoptosis inhibition and collagen accumulation through the ERK signaling pathway.

Conclusions

Our findings first revealed that HOXC6 acts as an oncogenic driver in the molecular mechanisms of fibroblasts in keloids. The HIF-1α/HOXC6/ERK axis promotes proliferation, migration and ECM production by KFs, contributing to the progression of keloids. Taken together, HOXC6 may serve as a promising novel therapeutic target and new focus for research designed to understand the pathogenesis of keloids.

Integrated analyses of an RNA binding protein-based signature related to tumor immune microenvironment and candidate drugs in osteosarcoma

OBJECTIVE

Osteosarcoma is the most frequent primary bone malignancy, associated with frequent recurrence and lung metastasis. RNA-binding proteins (RBPs) are pivotal in regulating several aspects of cancer biology. Nonetheless, interaction between RBPs and the osteosarcoma immune microenvironment is poorly understood. We investigated whether RBPs can predict prognosis and immunotherapy response in osteosarcoma patients.

METHODS

We constructed an RBP-related prognostic signature (RRPS) by univariate coupled with multivariate analyses and verified the independent prognostic efficacy of the signature. Single-sample Gene Set Enrichment Analysis (ssGSEA) along with ESTIMATE analysis were carried out to investigate the variations in immune characteristics between subgroups with various RRPS-scores. Furthermore, we investigatedpossible small molecule drugs using the connectivity map database and validated the expression of hub RBPs by qRT-PCR.

RESULTS

The RRPS, consisting of seven hub RBPs, was an independent prognostic factor compared to traditional clinical features. The RRPS could distinguish immune functions, immune score, stromal score, tumor purity and tumor infiltration by immune cells in different osteosarcoma subjects. Additionally, patients with high RRPS-scores had lower expression of immune checkpoint genes than patients with low RRPS-scores. We finally identified six small molecule drugs that may improve prognosis in osteosarcoma patients and substantiated notable differences in the contents of these RBPs.

CONCLUSION

We evaluated the prognostic value and clinical application of an RBPs-based prognostic signature and identified promising biomarkers to predict immune cell infiltration and immunotherapy response in osteosarcoma.