Using BioPAX-Parser (BiP) to enrich lists of genes or proteins with pathway data

BACKGROUND

Pathway enrichment analysis (PEA) is a well-established methodology for interpreting a list of genes and proteins of interest related to a condition under investigation. This paper aims to extend our previous work in which we introduced a preliminary comparative analysis of pathway enrichment analysis tools. We extended the earlier work by providing more case studies, comparing BiP enrichment performance with other well-known PEA software tools.

METHODS

PEA uses pathway information to discover connections between a list of genes and proteins as well as biological mechanisms, helping researchers to overcome the problem of explaining biological entity lists of interest disconnected from the biological context.

RESULTS

We compared the results of BiP with some existing pathway enrichment analysis tools comprising Centrality-based Pathway Enrichment, pathDIP, and Signaling Pathway Impact Analysis, considering three cancer types (colorectal, endometrial, and thyroid), for a total of six datasets (that is, two datasets per cancer type) obtained from the The Cancer Genome Atlas and Gene Expression Omnibus databases. We measured the similarities between the overlap of the enrichment results obtained using each couple of cancer datasets related to the same cancer.

CONCLUSION

As a result, BiP identified some well-known pathways related to the investigated cancer type, validated by the available literature. We also used the Jaccard and meet-min indices to evaluate the stability and the similarity between the enrichment results obtained from each couple of cancer datasets. The obtained results show that BiP provides more stable enrichment results than other tools.

Expression of genes and pathways associated with the B7-CD28 superfamily in response to irradiation of blood cells using 137Cs

PURPOSE

DNA damage is one of the main consequences of exposure to ionizing irradiation (IR). Recent studies indicate that IR can modulate the expression of immune system-related genes. However, the effects of IR on the expression of genes and pathways of the B7-CD28 superfamily remain poorly defined. The aim of this study was to evaluate the modulation of genes and pathways related to the B7-CD28 superfamily in response to IR.

MATERIALS AND METHODS

In this study, we used transcriptome data available from the Gene Expression Omnibus (GEO) database to investigate the modulation of the response of genes and pathways of samples of human peripheral blood irradiated with doses of 150, 300, and 600 cGy. The data were obtained at 6 and 24 h after irradiation. The relationship between genes and pathways was established through the Reactome database. The behavior of these pathways was analyzed using mathematical methods based on relative activity and diversity. Analysis of variance (ANOVA) followed by multiple comparisons tests (Bonferroni and Tamhanes) was used to identify differentially expressed genes. Data on transcriptomes were analyzed through ViaComplex V.1.0 and IBM SPSS Statistics 22.

RESULTS

For the pathways investigated in this study, we observed that the effects produced by these doses significantly modified the behavior of five pathways associated with the immune system. Also, the dose of 300 cGy might trigger signaling for the activation of T cells through the negative regulation (p < .05) of the co-inhibitory PDCD1LG2 gene. Positive regulation caused by 300 cGy (p < .05) of the CD80 receptor was observed by us, which might be related to a stimulatory signal. According to our findings, this dose induced the production of cytokines and genes that are associated with the activation and differentiation of T cells.

CONCLUSIONS

Our findings indicate that the irradiation modulated the organization of the biological system, suggesting that 300 cGy is more efficient in activating the immune system.

Expression deregulation of mir31 and CXCL12 in two types of oral precancers and cancer: importance in progression of precancer and cancer

Oral cancer generally progresses from precancerous lesions such as leukoplakia (LK), lichen planus (LP) and oral submucous fibrosis (OSMF). Since few of these precancers progress to cancers; it is worth to identify biological molecules that may play important roles in progression. Here, expression deregulation of 7 miRNAs (mir204, mir31, mir31*, mir133a, mir7, mir206 and mir1293) and their possible target genes in 23 cancers, 18 LK, 12 LP, 23 OSMF tissues compared to 20 healthy tissues was determined by qPCR method. Expression of mir7, mir31, mir31* and mir1293 was upregulated and that of mir133a, mir204 and mir206 was downregulated in cancer. Expression of most of these miRNAs was also upregulated in LK and LP tissues but not in OSMF. Expression deregulation of some of the target genes was also determined in cancer, LK and LP tissues. Significant upregulation of mir31 and downregulation of its target gene, CXCL12, in cancer, LK and LP tissues suggest their importance in progression of precancer to cancer. Expression upregulation of mir31 was also validated using GEO data sets. Although sample size is low, novelty of this work lies in studying expression deregulation of miRNAs and target genes in oral cancer and three types of precancerous lesions.

Expression of the γ-phosphorylated histone H2AX in gastric carcinoma and gastric precancerous lesions

The histone γH2AX is a marker of activated DNA damage that is overexpressed in various cancers and corresponding precursor lesions, indicating that γH2AX is a component in oncogenic transformation. The present study aimed to determine whether the immunohistochemical expression of γH2AX is involved in the progression between superficial gastritis (n=20), atrophic gastritis (n=24) and gastric carcinoma (n=79). There was no increase in γH2AX expression between superficial and atrophic gastritis, but there was a significant increase in γH2AX expression between these two conditions and gastric carcinoma (χ²=68.712; P<0.001). The expression of γH2AX in moderately-differentiated gastric adenocarcinoma (n=49) was evidently higher compared with poorly-differentiated gastric adenocarcinoma (n=26; χ²=14.241; P<0.01). Staining for γH2AX did not reveal a significant association between the expression of the histone and patient age, depth of invasion, lymph node metastasis or the tumor-node-metastasis stage of the gastric carcinoma. Overall, the present study demonstrated that enhanced γH2AX expression may be closely associated with gastric carcinoma, but is less likely to be involved in the genesis of gastric carcinoma.

Bufalin induces cell death in human lung cancer cells through disruption of DNA damage response pathways

Bufalin is a key component of a Chinese medicine (Chan Su) and has been proved effective in killing various cancer cells. Its role in inducing DNA damage and the inhibition of the DNA damage response (DDR) has been reported, but none have studied such action in lung cancer in detail. In this study, we demonstrated bufalin-induced DNA damage and condensation in NCI-H460 cells through a comet assay and DAPI staining, respectively. Western blotting indicated that bufalin suppressed the protein levels associated with DNA damage and repair, such as a DNA dependent serine/threonine protein kinase (DNA-PK), DNA repair proteins breast cancer 1, early onset (BRCA1), 14-3-3 σ (an important checkpoint keeper of DDR), mediator of DNA damage checkpoint 1 (MDC1), O6-methylguanine-DNA methyltransferase (MGMT) and p53 (tumor suppressor protein). Bufalin could activate phosphorylated p53 in NCI-H460 cells. DNA damage in NCI-H460 cells after treatment with bufalin up-regulated its ATM and ATR genes, which encode proteins functioning as sensors in DDR, and also up-regulated the gene expression (mRNA) of BRCA1 and DNA-PK. But bufalin suppressed the gene expression (mRNA) of p53 and 14-3-3 σ, however, bufalin did not significantly affect the mRNA of MGMT. In conclusion, bufalin induced DNA damage in NCI-H460 cells and also inhibited its DNA repair and checkpoint function.

Modelling the onset of senescence at the G1/S cell cycle checkpoint

BACKGROUND

DNA damage (single or double-strand breaks) triggers adapted cellular responses. These responses are elicited through signalling pathways, which activate cell cycle checkpoints and basically lead to three cellular fates: cycle arrest promoting DNA repair, senescence (permanent arrest) or cell death. Cellular senescence is known for having a tumour-suppressive function and its regulation arouses a growing scientific interest. Here, we advance a qualitative model covering DNA damage response pathways, focusing on G1/S checkpoint enforcement, supposedly more sensitive to arrest than G2/M checkpoint.

RESULTS

We define a discrete, logical model encompassing ATM (ataxia telangiectasia mutated) and ATR (ATM and Rad3-related) pathways activation upon DNA damage, as well as G1/S checkpoint main components. It also includes the stress responsive protein p38MAPK (mitogen-activated protein kinase 14) known to be involved in the regulation of senescence. The model has four outcomes that convey alternative cell fates: proliferation, (transient) cell cycle arrest, apoptosis and senescence. Different levels of DNA damage are considered, defined by distinct combinations of single and double-strand breaks. Each leads to a single stable state denoting the cell fate adopted upon this specific damage. A range of model perturbations corresponding to gene loss-of-function or gain-of-function is compared to experimental mutations.

CONCLUSIONS

As a step towards an integrative model of DNA-damage response pathways to better cover the onset of senescence, our model focuses on G1/S checkpoint enforcement. This model qualitatively agrees with most experimental observations, including experiments involving mutations. Furthermore, it provides some predictions.