0.1 Introduction

0.1.1 libraries required

0.1.2 CARDIA Project:: Proteiomics pathway analysis_Feb14-2024

Steps followed:

Identify the variables, model in use and separate the significant protein markers based on its p. value<0.05.

Pathway Analysis

Identify the UniProt ID column

UniProt IDs to Entrez Gene IDs

Perform pathway enrichment analysis

Print the pathway enrichment results

Plot pathway enrichment results

The format of CARDIA proteiomics dataset should look familiar to anyone who has worked with differential expression results. It contains columns of gene identifiers (UNIPORT ID ), Protein symbols, log2foldchange values, P-Values and adjusted P-Values.

Identify the Uniport ID and convert them to Entrez Gene ID as its required for KEGG pathay analysis.

Now let’s prepare significant protein lists using some conventional criteria (Pvalue<0.05).

We have Uniport IDs, but the package we are going to use wants Entrez IDs. Fortunately, the package provides its own converter called bitr. This is the first function in the clusterProfiler package that we’ll be using:

Pathway Enrichment Analysis:: we are ready to perform enrichment analysis. Let’s start with Gene Ontology.

Convert Enrichment result to table

knitr::opts_chunk$set(collapse = TRUE, include = TRUE)

enrich_result_df<- data.frame(enrich_result)

enrich_result_2 <- enrichKEGG(gene         = gene_ids_df$gene_ids,
                            organism     = 'hsa',
                            pvalueCutoff = 0.05)

kegg <- enrichKEGG(  gene_ids_df$gene_ids,  organism = "hsa",  keyType = "kegg",  pvalueCutoff = 0.05,
                     pAdjustMethod = "BH", qvalueCutoff = 0.05)

p<- dotplot(kegg, showCategory=20) + ggtitle("CARDIA_protiomics ")
p + theme(
  plot.title = element_text(color="black", size=14, face="bold", hjust=0.5))