Identification of candidate genes related to pancreatic cancer based on analysis of gene co-expression and protein-protein interaction network

A two-way additive model with unknown group-specific interactions applied to gene expression data

We propose a two-way additive model with group-specific interactions, where the group information is unknown. We treat the group membership as latent information and propose an EM algorithm for estimation. With a single observation matrix and under the situation of diverging row and column numbers, we rigorously establish the estimation consistency and asymptotic normality of our estimator. Extensive simulation studies are conducted to demonstrate the finite sample performance. We apply the model to the triple negative breast cancer (TNBC) gene expression data and provide a new way to classify patients into different subtypes. Our analysis detects the potential genes that may be associated with TNBC.

A Novel IGLC2 Gene Linked With Prognosis of Triple-Negative Breast Cancer

Background

Immunoglobulin-related genes are associated with the favorable prognosis of triple-negative breast cancer (TNBC) patients. We aimed to analyze the function and prognostic value of immunoglobulin lambda constant 2 (IGLC2) in TNBC patients.

Methods

We knocked down the gene expression of IGLC2 (IGLC2-KD) in MDA-MB-231 cells to evaluate the proliferation, migration, and invasion of tumors via 3-(4,5-Dimethythiazol-2-yl)-2,5-diphenyl tetrazolium bromide assay, wound healing, and transwell cell migration assay respectively. Relapse-free survival (RFS) and distant metastasis-free survival (DMFS) analyses were conducted using the KM plotter online tool. The GSE76275 data set was used to analyze the association of IGLC2 and clinical characteristics. A pathway enrichment analysis was conducted using the next-generation sequencing data of wild-type and IGLC2-KD MDA-MB-231 cells.

Results

The low gene expression of IGLC2 was related to unfavorable RFS, DMFS. The high expression of IGLC2 was exhibited in the basal-like immune-activated (BLIA) TNBC molecular subtype, which was immune-activated and showed excellent response to immune therapy. IGLC2 was positively correlated with programmed death-ligand 1 (PD-L1) as shown by Spearman correlation (r = 0.25, p < 0.0001). IGLC2 had a strong prognostic effect on lymph node-negative TNBC (RFS range: 0.31, q value= 8.2e-05; DMFS = 0.16, q value = 8.2e-05) but had no significance on lymph node-positive ones. The shRNA-mediated silencing of IGLC2 increased the proliferation, migration, and invasion of MDA-MB-231 cells. The results of pathway enrichment analysis showed that IGLC2 is related to the PI3K-Akt signaling pathway, MAPK signaling pathway, and extracellular matrix–receptor interaction. We confirmed that MDA-MB-231 tumor cells expressed IGLC2, subverting the traditional finding of generation by immune cells.

Conclusions

IGLC2 linked with the proliferation, migration, and invasion of MDA-MB-231 cells. A high expression of IGLC2 was related to favorable prognosis for TNBC patients. IGLC2 may serve as a biomarker for the identification of TNBC patients who can benefit the most from immune checkpoint blockade treatment.

Locally Adjust Networks Based on Connectivity and Semantic Similarities for Disease Module Detection

For studying the pathogenesis of complex diseases, it is important to identify the disease modules in the system level. Since the protein-protein interaction (PPI) networks contain a number of incomplete and incorrect interactome, most existing methods often lead to many disease proteins isolating from disease modules. In this paper, we propose an effective disease module identification method IDMCSS, where the used human PPI networks are obtained by adding some potential missing interactions from existing PPI networks, as well as removing some potential incorrect interactions. In IDMCSS, a network adjustment strategy is developed to add or remove links around disease proteins based on both topological and semantic information. Next, neighboring proteins of disease proteins are prioritized according to a suggested similarity between each of them and disease proteins, and the protein with the largest similarity with disease proteins is added into a candidate disease protein set one by one. The stopping criterion is set to the boundary of the disease proteins. Finally, the connected subnetwork having the largest number of disease proteins is selected as a disease module. Experimental results on asthma demonstrate the effectiveness of the method in comparison to existing algorithms for disease module identification. It is also shown that the proposed IDMCSS can obtain the disease modules having crucial biological processes of asthma and 12 targets for drug intervention can be predicted.

S100A14 promotes progression and gemcitabine resistance in pancreatic cancer

S100 calcium binding protein A14 (S100A14) plays an important role in the progression of several types of cancer. However, its roles in pancreatic ductal adenocarcinoma (PDAC) are largely unexplored. Here, we characterized the functional roles of S100A14 in the progression and chemoresistance of PDAC. Gene expression microarray identified that S100A14 was significantly highly expressed in four pairs of human PDAC tumor compared with corresponding non-tumor tissues genes. Quantitative reverse transcription PCR (qRT-PCR), western blotting and immunohistochemical staining (IHC) showed that S100A14 was frequently overexpressed in PDAC cell lines and tissues. Moreover, expression level of S100A14 was positively correlated to advanced cancer stages. Further, Kaplan-Meier survival analysis suggested that PDAC patients with low S100A14 expression had longer overall survival in TCGA PDAC datasets. Transient overexpressing of S100A14 promoted cell proliferation, anchorage-independent colony formation, cell migration and invasion in cell lines with low endogenous S100A14 levels, while transient silencing of S100A14 inhibited cell proliferation, anchorage-independent colony formation, cell migration and invasion in cell lines with high endogenous S100A14 levels. Persistent knockdown of S100A14 by transducing shRNAs carrying lentivirus inhibited subcutaneous tumor formation in nude mice, and sensitized the PDAC cells to gemcitabine treatment. Taken together, S100A14 exhibited oncogenic properties by promoting cell proliferation, transformation, migration and invasion, and enhanced in vivo tumor growth. More importantly, inhibition of S100A14 could effectively abrogate the cancerous properties of the PDAC cells. Our study indicated that S100A14 was a valuable target for the development of therapeutic strategy, as well as a diagnostic and prognosis biomarker for PDAC patients.

Transcriptome analysis of the procession from chronic pancreatitis to pancreatic cancer and metastatic pancreatic cancer

Exploring the underlying mechanisms of cancer development is useful for cancer treatment. In this paper, we analyzed the transcriptome profiles from the human normal pancreas, pancreatitis, pancreatic cancer and metastatic pancreatic cancer to study the intricate associations among pancreatic cancer progression. We clustered the transcriptome data, and analyzed the differential expressed genes. WGCNA was applied to construct co-expression networks and detect important modules. Importantly we selected the module in a different way. As the pancreatic disease deteriorates, the number of differentially expressed genes increases. The gene networks of T cells and interferon are upregulated in stages. In conclusion, the network-based study provides gradually activated gene networks in the disease progression of pancreatitis, pancreatic cancer, and metastatic pancreatic cancer. It may contribute to the rational design of anti-cancer drugs.

The Effects of Genetic and Epigenetic Alterations of BARD1 on the Development of Non-Breast and Non-Gynecological Cancers

Breast Cancer 1 (BRCA1) gene is a well-characterized tumor suppressor gene, mutations of which are primarily found in women with breast and ovarian cancers. BRCA1-associated RING domain 1 (BARD1) gene has also been identified as an important tumor suppressor gene in breast, ovarian, and uterine cancers. Underscoring the functional significance of the BRCA1 and BARD1 interactions, prevalent mutations in the BRCA1 gene are found in its RING domain, through which it binds the RING domain of BARD1. BARD1-BRCA1 heterodimer plays a crucial role in a variety of DNA damage response (DDR) pathways, including DNA damage checkpoint and homologous recombination (HR). However, many mutations in both BARD1 and BRCA1 also exist in other domains that significantly affect their biological functions. Intriguingly, recent genome-wide studies have identified various single nucleotide polymorphisms (SNPs), genetic alterations, and epigenetic modifications in or near the BARD1 gene that manifested profound effects on tumorigenesis in a variety of non-breast and non-gynecological cancers. In this review, we will briefly discuss the molecular functions of BARD1, including its BRCA1-dependent as well as BRCA1-independent functions. We will then focus on evaluating the common BARD1 related SNPs as well as genetic and epigenetic changes that occur in the non-BRCA1-dominant cancers, including neuroblastoma, lung, and gastrointestinal cancers. Furthermore, the pro- and anti-tumorigenic functions of different SNPs and BARD1 variants will also be discussed.

Effects of the MDM-2 inhibitor Nutlin-3a on PDAC cells containing and lacking WT-TP53 on sensitivity to chemotherapy, signal transduction inhibitors and nutraceuticals

Mutations at the TP53 gene are readily detected (approximately 50-75%) in pancreatic ductal adenocarcinoma (PDAC) patients. TP53 was previously thought to be a difficult target as it is often mutated, deleted or inactivated on both chromosomes in certain cancers. In the following study, the effects of restoration of wild-type (WT) TP53 activity on the sensitivities of MIA-PaCa-2 pancreatic cancer cells to the MDM2 inhibitor nutlin-3a in combination with chemotherapy, targeted therapy, as well as, nutraceuticals were examined. Upon introduction of the WT-TP53 gene into MIA-PaCa-2 cells, which contain a TP53 gain of function (GOF) mutation, the sensitivity to the MDM2 inhibitor increased. However, effects of nutlin-3a were also observed in MIA-PaCa-2 cells lacking WT-TP53, as upon co-treatment with nutlin-3a, the sensitivity to certain inhibitors, chemotherapeutic drugs and nutraceuticals increased. Interestingly, co-treatment with nutlin-3a and certain chemotherapeutic drug such as irinotecan and oxaliplatin resulted in antagonistic effects in cells both lacking and containing WT-TP53 activity. These studies indicate the sensitizing abilities that WT-TP53 activity can have in PDAC cells which normally lack WT-TP53, as well as, the effects that the MDM2 inhibitor nutlin-3a can have in both cells containing and lacking WT-TP53 to various therapeutic agents.

Introduction of WT-TP53 into pancreatic cancer cells alters sensitivity to chemotherapeutic drugs, targeted therapeutics and nutraceuticals

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive, highly metastatic malignancy and accounts for 85% of pancreatic cancers. PDAC patients have poor prognosis with a five-year survival of only 5-10%. Mutations at the TP53 gene are readily detected in pancreatic tumors isolated from PDAC patients. We have investigated the effects of restoration of wild-type (WT) TP53 activity on the sensitivity of pancreatic cancer cells to: chemotherapy, targeted therapy, as well as, nutraceuticals. Upon introduction of the WT-TP53 gene into the MIA-PaCa-2 pancreatic cancer cell line, the sensitivity to drugs used to treat pancreatic cancer cells such as: gemcitabine, fluorouracil (5FU), cisplatin, irinotecan, oxaliplatin, and paclitaxel increased significantly. Likewise, the sensitivity to drugs used to treat other cancers such as: doxorubicin, mitoxantrone, and 4 hydroxy tamoxifen (4HT) also increased upon introduction of WT-TP53 into MIA-PaCa-2 cells. Furthermore, the sensitivity to certain inhibitors which target: PI3K/mTORC1, PDK1, SRC, GSK-3, and biochemical processes such as proteasomal degradation and the nutraceutical berberine as increased upon introduction of WT-TP53. Furthermore, in some cases, cells with WT-TP53 were more sensitive to the combination of drugs and suboptimal doses of the MDM2 inhibitor nutlin-3a. However, TP53-independent effects of nutlin-3a were observed upon treatment with either a proteasomal or a PI3K/mTOR inhibitor. These studies indicate the sensitizing effects that WT-TP53 can have in PDAC cells which normally lack WT-TP53 to various therapeutic agents and suggest approaches to improve PDAC therapy.