PADI4‑mediated epithelial‑mesenchymal transition in lung cancer cells

Identifying Key Genes Involved in Axillary Lymph Node Metastasis in Breast Cancer Using Advanced RNA-Seq Analysis: A Methodological Approach with GLMQL and MAS

Our study aims to address the methodological challenges frequently encountered in RNA-Seq data analysis within cancer studies. Specifically, it enhances the identification of key genes involved in axillary lymph node metastasis (ALNM) in breast cancer. We employ Generalized Linear Models with Quasi-Likelihood (GLMQLs) to manage the inherently discrete and overdispersed nature of RNA-Seq data, marking a significant improvement over conventional methods such as the t-test, which assumes a normal distribution and equal variances across samples. We utilize the Trimmed Mean of M-values (TMMs) method for normalization to address library-specific compositional differences effectively. Our study focuses on a distinct cohort of 104 untreated patients from the TCGA Breast Invasive Carcinoma (BRCA) dataset to maintain an untainted genetic profile, thereby providing more accurate insights into the genetic underpinnings of lymph node metastasis. This strategic selection paves the way for developing early intervention strategies and targeted therapies. Our analysis is exclusively dedicated to protein-coding genes, enriched by the Magnitude Altitude Scoring (MAS) system, which rigorously identifies key genes that could serve as predictors in developing an ALNM predictive model. Our novel approach has pinpointed several genes significantly linked to ALNM in breast cancer, offering vital insights into the molecular dynamics of cancer development and metastasis. These genes, including ERBB2, CCNA1, FOXC2, LEFTY2, VTN, ACKR3, and PTGS2, are involved in key processes like apoptosis, epithelial-mesenchymal transition, angiogenesis, response to hypoxia, and KRAS signaling pathways, which are crucial for tumor virulence and the spread of metastases. Moreover, the approach has also emphasized the importance of the small proline-rich protein family (SPRR), including SPRR2B, SPRR2E, and SPRR2D, recognized for their significant involvement in cancer-related pathways and their potential as therapeutic targets. Important transcripts such as H3C10, H1-2, PADI4, and others have been highlighted as critical in modulating the chromatin structure and gene expression, fundamental for the progression and spread of cancer.

Mitosis targeting in non-small lung cancer cells by inhibition of PAD4

PAD4 expression and activity were significantly up-regulated in lung cancer tissues suggesting that PAD4 could be a possible target for lung cancer treatment. In this study we had demonstrated that PAD4 expression was higher in lung cancer patients whom with lymphnode metastasis and pleural invasion. Inhibiting PAD4 with a small molecular inhibitor could induce apoptosis and suppress growth in lung cancer cells. We used RNA-sequencing to further investigate transcriptional changes that induced by PAD4 inhibition, and results suggested its affected mostly on the cell cycle, mitotic cell cycle process, p53 signaling pathway. By using image flow cytometry analysis, we found that PAD4 inhibited by YW3-56 could accumulate cells in the G1/G0 phases and reducing the fraction of G2/M and S phase cells. Quantification of different phase of mitosis in cells treated with YW3-56 revealed an increasing trend of telophase and prophase cells. Taken together, our data indicated that PAD4 inhibitor could affect cell cycle and mitosis of lung cancer cells, and targeting PAD4 could be a promising strategy for discovery novel anti-NSCLC treatments.

Gene signatures associated with prognosis and chemotherapy resistance in glioblastoma treated with temozolomide

Background: Glioblastoma (GBM) prognosis remains extremely poor despite standard treatment that includes temozolomide (TMZ) chemotherapy. To discover new GBM drug targets and biomarkers, genes signatures associated with survival and TMZ resistance in GBM patients treated with TMZ were identified. Methods: GBM cases in The Cancer Genome Atlas who received TMZ (n = 221) were stratified into subgroups that differed by median overall survival (mOS) using network-based stratification to cluster patients whose somatic mutations affected genes in similar modules of a gene interaction network. Gene signatures formed from differentially mutated genes in the subgroup with the longest mOS were used to confirm their association with survival and TMZ resistance in independent datasets. Somatic mutations in these genes also were assessed for an association with OS in an independent group of 37 GBM cases. Results: Among the four subgroups identified, subgroup four (n = 71 subjects) exhibited the longest mOS at 18.3 months (95% confidence interval: 16.2, 34.1; p = 0.0324). Subsets of the 86 genes that were differentially mutated in this subgroup formed 20-gene and 8-gene signatures that predicted OS in two independent datasets (Spearman's rho of 0.64 and 0.58 between actual and predicted OS; p < 0.001). Patients with mutations in five of the 86 genes had longer OS in a small, independent sample of 37 GBM cases, but this association did not reach statistical significance (p = 0.07). Thirty-one of the 86 genes formed signatures that distinguished TMZ-resistant GBM samples from controls in three independent datasets (area under the curve ≥ 0.75). The prognostic and TMZ-resistance signatures had eight genes in common (ANG, BACH1, CDKN2C, HMGA1, IFI16, PADI4, SDF4, and TP53INP1). The latter three genes have not been associated with GBM previously. Conclusion: PADI4, SDF4, and TP53INP1 are novel therapy and biomarker candidates for GBM. Further investigation of their oncologic functions may provide new insight into GBM treatment resistance mechanisms.

Identification and Interaction Analysis of Molecular Markers in Pancreatic Ductal Adenocarcinoma by Bioinformatics and Next-Generation Sequencing Data Analysis

Background

Pancreatic ductal adenocarcinoma (PDAC) is one of the most common cancers worldwide. Intense efforts have been made to elucidate the molecular pathogenesis, but the molecular mechanisms of PDAC are still not well understood. The purpose of this study is to further explore the molecular mechanism of PDAC through integrated bioinformatics analysis.

Methods

To identify the candidate genes in the carcinogenesis and progression of PDAC, next-generation sequencing (NGS) data set GSE133684 was downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified, and Gene Ontology (GO) and pathway enrichment analyses were performed. The protein-protein interaction network (PPI) was constructed and the module analysis was performed using Integrated Interactions Database (IID) interactome database and Cytoscape. Subsequently, miRNA-DEG regulatory network and TF-DEG regulatory network were constructed using miRNet database, NetworkAnalyst database, and Cytoscape software. The expression levels of hub genes were validated based on Kaplan-Meier analysis, expression analysis, stage analysis, mutation analysis, protein expression analysis, immune infiltration analysis, and receiver operating characteristic (ROC) curve analysis.

Results

A total of 463 DEGs were identified, consisting of 232 upregulated genes and 233 downregulated genes. The enriched GO terms and pathways of the DEGs include vesicle organization, secretory vesicle, protein dimerization activity, lymphocyte activation, cell surface, transferase activity, transferring phosphorus-containing groups, hemostasis, and adaptive immune system. Four hub genes (namely, cathepsin B [CCNB1], four-and-a-half LIM domains 2 (FHL2), major histocompatibility complex, class II, DP alpha 1 (HLA-DPA1) and tubulin beta 1 class VI (TUBB1)) were obtained via taking interaction of different analysis results.

Conclusions

On the whole, the findings of this investigation enhance our understanding of the potential molecular mechanisms of PDAC and provide potential targets for further investigation.

PAD4-dependent citrullination of nuclear translocation of GSK3β promotes colorectal cancer progression via the degradation of nuclear CDKN1A

Peptidylarginine deiminase 4 (PAD4), a Ca²⁺-dependent enzyme, catalyzes the conversion of arginine to citrulline and has been strongly associated with many malignant tumors. However, the molecular mechanisms of PAD4 in the development and progression of colorectal cancer (CRC) remain unclearly defined. In our study, PAD4 expression was increased in CRC tissues and cells, and was closely related to tumor size, lymph node metastasis. Moreover, the transcription factor KLF9 directly bound to PADI4 gene promoter, leading to overexpression of PAD4 in CRC cells, which augmented cell growth and migration. We revealed that PAD4 interacted with and citrullinated glycogen synthase kinase-3β (GSK3β) in CRC cells, and GSK3β Arg-344 was the dominating PAD4-citrullination site. Furthermore, IgL2 and catalytic domains of PAD4 directly bound to the kinase domain of GSK3β in CRC cells. Mechanistically, PAD4 promoted the transport of GSK3β from the cytoplasm to the nucleus, thereby increasing the ubiquitin-dependent proteasome degradation of nuclear cyclin-dependent kinase inhibitor 1 (CDKN1A). Our study is the first to reveal the details of a critical PAD4/GSK3β/CDKN1A signaling axis for CRC progression, and provides evidence that PAD4 is a potential diagnosis biomarker and therapeutic target in CRC.

A Screening of Epigenetic Therapeutic Targets for Non-Small Cell Lung Cancer Reveals PADI4 and KDM6B as Promising Candidates

Despite advances in diagnostic and therapeutic approaches for lung cancer, new therapies targeting metastasis by the specific regulation of cancer genes are needed. In this study, we screened a small library of epigenetic inhibitors in non-small-cell lung cancer (NSCLC) cell lines and evaluated 38 epigenetic targets for their potential role in metastatic NSCLC. The potential candidates were ranked by a streamlined approach using in silico and in vitro experiments based on publicly available databases and evaluated by real-time qPCR target gene expression, cell viability and invasion assays, and transcriptomic analysis. The survival rate of patients with lung adenocarcinoma is inversely correlated with the gene expression of eight epigenetic targets, and a systematic review of the literature confirmed that four of them have already been identified as targets for the treatment of NSCLC. Using nontoxic doses of the remaining inhibitors, KDM6B and PADI4 were identified as potential targets affecting the invasion and migration of metastatic lung cancer cell lines. Transcriptomic analysis of KDM6B and PADI4 treated cells showed altered expression of important genes related to the metastatic process. In conclusion, we showed that KDM6B and PADI4 are promising targets for inhibiting the metastasis of lung adenocarcinoma cancer cells.

Dynamic proteomic change of tumor and immune organs in an immune-competent hepatocellular carcinoma mouse model

Subcutaneous implantation of a human cancer cell line in immune-deficient mice (CDX) is a commonly used tool in preclinical studies for the assessment of potential anti-cancer drugs. As immunotherapy is transforming cancer treatment, tumor models in immunocompetent mice are necessary for us to understand the immune aspects of tumor biology. However, the systemic immune response to the implantation of cancer cells at proteome level is unclear. In this study, we characterized the dynamic proteomic changes of subcutaneous tumors and 5 immune organs (draining lymph node, mesenteric lymph node, spleen, thymus and marrow) at six time points after implantation using a Hepa1-6 derived allograft mouse model. Our data suggest that interaction of the implanted tumor cells with mouse immune system followed the trajectory of "tumor rejection" to "immune evasion" in that the tumor gained the ability to evade the immune system for growth. Furthermore, anti-PDL2 antibody was validated here as an optional immunotherapy strategy to inhibit the growth of Hepa1-6 subcutaneous tumors. These findings from our study provided valuable information for the understanding of tumor and immune interaction and shed light on the rational design for clinical cancer treatment and other preclinical experiments.

Suppression of epithelial abnormalities by nintedanib in induced-rheumatoid arthritis-associated interstitial lung disease mouse model

Rheumatoid arthritis-associated interstitial lung disease (RA-ILD) is relevant for the prognosis in patients with RA. Nintedanib, which inhibits both receptor and non-receptor type tyrosine kinases, is an antifibrotic drug for the treatment of progressive fibrosing ILDs, such as idiopathic pulmonary fibrosis and systemic sclerosis-associated interstitial lung disease. Little is known about the effects of nintedanib on RA-ILD. We examined the characteristics of a novel induced RA-ILD (iRA-ILD) mouse model and the effects of nintedanib on the model. D1CC×D1BC mice are highly susceptible to arthritogenic antigens, such as bovine type II collagen, resulting in severe inflammatory arthritis. ILD develops after joint inflammation is alleviated. Serum surfactant protein D levels were monitored as an ILD marker. Nintedanib was orally administered to iRA-ILD mice for 2 months. The iRA-ILD model showed similar symptoms to those in patients with RA-ILD. The histopathological features of pulmonary disorder resembled nonspecific interstitial pneumonia, but with metaplastic epithelium. Histopathological analysis revealed that in addition to reducing fibrosis, nintedanib suppressed M2 macrophage polarisation and hyperplasia of Type 2 alveolar epithelial cells. The metaplastic epithelium acquired invasiveness because of the expression of E-cadherin, MMP7, Tgf-β, Col1a1, Padi2 and Padi4. Moreover, citrullinated peptides were detected in these invasive epithelial cells as well as in the bronchiolar epithelium. Administration of nintedanib reduced the expression of Pad4 and citrullinated peptides and eliminated invasive epithelial cells. The broad inhibitory effects of nintedanib on tyrosine kinases may contribute to the overall improvement in RA-ILD, including epithelial abnormalities associated with progressive lung fibrosis.

Gimap5 Inhibits Lung Cancer Growth by Interacting With M6PR

Objective

Several studies have demonstrated the impacts of GTPases of immunity-associated proteins (GIMAPs) on malignant cells. However, the mechanisms through which Gimap5 regulates lung cancer cells are yet to be thoroughly investigated in the literature. Our study aimed to investigate the function of Gimap5 in the development of lung cancer.

Methods

The expression levels of the GIMAP family were analyzed in lung cancer patients of various cancer databases and lung cancer cell lines. After the survival rates of the cells were analyzed, we constructed Gimap5 over-expressed lung cancer cell lines and assessed the effects of Gimap5 on cell migration, cell invasion, cell proliferation and the epithelial-mesenchymal transition (EMT). We later screened the interacting proteins of Gimap5 using Co-IP combined with mass spectrometry and then analyzed the expression and distribution of M6PR, including its impacts on protein-arginine deiminase type-4 (PADI4).

Results

Findings indicated that GIMAP family expression decreased significantly in lung cancer cell lines. We also noticed that the downregulation of the GIMAP family was related to the poor prognosis of lung cancer patients. Our experimental results showed that Gimap5 could inhibit the migration, invasion, proliferation and EMT of lung cancer cell lines. Moreover, we found that Gimap5 promoted the transport of M6PR from the cytoplasm to the cell membrane, thereby inhibiting the enhancement of EMT-related PADI4.

Conclusion

Our research suggested that Gimap5 could inhibit the growth of lung cancer by interacting with M6PR and that it could be a potential biomarker for the diagnosis and prognosis of lung cancer.

PADI1 contributes to EMT in PAAD by activating the ERK1/2-p38 signaling pathway

Background

Peptidylarginine deiminase 1 (PADI1) has been reported to promote tumorigenesis in breast cancer. However, the functional role of PADI1 in pancreatic ductal adenocarcinoma (PAAD) has remained elusive until now.

Methods

The expression pattern of PADI1 in PAAD tissues and normal tissues was analyzed using The Cancer Genome Atlas (TCGA) dataset. A Kaplan-Meier curve analysis was performed to evaluate the prognostic value of PADI1 in PAAD patients. PADI1 was knocked down in CFPAN-1 and HPAC cells, and overexpressed in PANC-1 and Bxpc-3 cells by RNA interference. A wound-healing assay was performed to analyze relative cell migration distance. Cell migration and invasion were assessed by a Transwell assay. Related protein expression levels were measured by western blot and immunofluorescence.

Results

The bioinformatics analysis showed that PADI1 was overexpressed in PAAD tissues and associated with a poor survival prognosis. The knockdown of PADI1 suppressed cell migration and invasion, and activated the ERK1/2-p38 signaling pathway in CFPAN-1 and HPAC cells. The overexpression of PADI1 produced the opposite results in PANC-1 and Bxpc-3 cells. Additionally, treatment with an MEK1/2 inhibitor significantly attenuated the effects of PADI1 knockdown on cell migration, invasion, the epithelial-mesenchymal transition (EMT) process, and p-ERK1/2 and p38 expression in CFPAN-1 and HPAC cells.

Conclusions

Our data suggested that PADI1 may function as an oncogene in regulating metastasis in vitro in PAAD.

Inhibition of PAD4 enhances radiosensitivity and inhibits aggressive phenotypes of nasopharyngeal carcinoma cells

BACKGROUND

Nasopharyngeal carcinoma (NPC) is a tumor deriving from nasopharyngeal epithelium. Peptidyl-arginine deiminase 4 (PAD4) is a vital mediator of histone citrullination and plays an essential role in regulating disease process. Radiotherapy is an essential method to treat NPC. In this research, we explored the effect of PAD4 on NPC radiosensitivity.

METHODS

We enrolled 50 NPC patients, established mice xenograft model, and purchased cell lines for this study. Statistical analysis and a series of experiments including RT-qPCR, clonogenic survival, EdU, Transwell, and wound healing assays were done.

RESULTS

Our data manifested that PAD4 (mRNA and protein) presented a high expression in NPC tissues and cells. GSK484, an inhibitor of PAD4, could inhibit activity of PAD4 in NPC cell lines. PAD4 overexpression promoted the radioresistance, survival, migration, and invasion of NPC cells, whereas treatment of GSK484 exerted inhibitory effects on radioresistance and aggressive phenotype of NPC cells. Additionally, GSK484 could attenuate the effect of PAD4 of NPC cell progression. More importantly, we found that GSK484 significantly inhibited tumor size, tumor weight and tumor volume in mice following irradiation.

CONCLUSIONS

PAD4 inhibitor GSK484 attenuated the radioresistance and cellular progression in NPC.

Biochemical and Computational Approaches for the Large-Scale Analysis of Protein Arginine Methylation by Mass Spectrometry

The absence of efficient mass spectrometry-based approaches for the large-scale analysis of protein arginine methylation has hindered the understanding of its biological role, beyond the transcriptional regulation occurring through histone modification. In the last decade, however, several technological advances of both the biochemical methods for methylated polypeptide enrichment and the computational pipelines for MS data analysis have considerably boosted this research field, generating novel insights about the extent and role of this post-translational modification. Here, we offer an overview of state-of-the-art approaches for the high-confidence identification and accurate quantification of protein arginine methylation by high-resolution mass spectrometry methods, which comprise the development of both biochemical and bioinformatics methods. The further optimization and systematic application of these analytical solutions will lead to ground-breaking discoveries on the role of protein methylation in biological processes.

Putative Roles for Peptidylarginine Deiminases in COVID-19

Peptidylarginine deiminases (PADs) are a family of calcium-regulated enzymes that are phylogenetically conserved and cause post-translational deimination/citrullination, contributing to protein moonlighting in health and disease. PADs are implicated in a range of inflammatory and autoimmune conditions, in the regulation of extracellular vesicle (EV) release, and their roles in infection and immunomodulation are known to some extent, including in viral infections. In the current study we describe putative roles for PADs in COVID-19, based on in silico analysis of BioProject transcriptome data (PRJNA615032 BioProject), including lung biopsies from healthy volunteers and SARS-CoV-2-infected patients, as well as SARS-CoV-2-infected, and mock human bronchial epithelial NHBE and adenocarcinoma alveolar basal epithelial A549 cell lines. In addition, BioProject Data PRJNA631753, analysing patients tissue biopsy data (n = 5), was utilised. We report a high individual variation observed for all PADI isozymes in the patients' tissue biopsies, including lung, in response to SARS-CoV-2 infection, while PADI2 and PADI4 mRNA showed most variability in lung tissue specifically. The other tissues assessed were heart, kidney, marrow, bowel, jejunum, skin and fat, which all varied with respect to mRNA levels for the different PADI isozymes. In vitro lung epithelial and adenocarcinoma alveolar cell models revealed that PADI1, PADI2 and PADI4 mRNA levels were elevated, but PADI3 and PADI6 mRNA levels were reduced in SARS-CoV-2-infected NHBE cells. In A549 cells, PADI2 mRNA was elevated, PADI3 and PADI6 mRNA was downregulated, and no effect was observed on the PADI4 or PADI6 mRNA levels in infected cells, compared with control mock cells. Our findings indicate a link between PADI expression changes, including modulation of PADI2 and PADI4, particularly in lung tissue, in response to SARS-CoV-2 infection. PADI isozyme 1-6 expression in other organ biopsies also reveals putative links to COVID-19 symptoms, including vascular, cardiac and cutaneous responses, kidney injury and stroke. KEGG and GO pathway analysis furthermore identified links between PADs and inflammatory pathways, in particular between PAD4 and viral infections, as well as identifying links for PADs with a range of comorbidities. The analysis presented here highlights roles for PADs in-host responses to SARS-CoV-2, and their potential as therapeutic targets in COVID-19.