Analysis Of RNA-seq Data Using R

Nadia Atallah Lanman
November 19, 2019

Central Dogma of Molecular Biology

  • By measuring how many mRNA molecules there are for each gene, we can quantify which genes are “on” or “off”

Why do gene expression analyses?

  • Now feasible and simple to measure the expression levels of thousands of genes simultaneously
  • Gene Annotation
  • Functional studies
  • Some molecular features are only seen at RNA level

Several types of gene expression analyses

  • qRT-PCR
  • microarray
  • RNA-seq
    • bulk
    • single-cell

Perform a GO enrichment analysis

mf.go <- enrichGO(gene          = deg_ann.tb$entrezgene_id,
                  OrgDb         = org.Hs.eg.db,
                  ont           = "MF",
                  pAdjustMethod = "BH",
                  qvalueCutoff  = 0.05,
                  readable      = TRUE)
head(mf.go)
                   ID                                 Description GeneRatio
GO:0050839 GO:0050839              cell adhesion molecule binding 493/15137
GO:0045296 GO:0045296                            cadherin binding 328/15137
GO:0004674 GO:0004674    protein serine/threonine kinase activity 427/15137
GO:0003779 GO:0003779                               actin binding 416/15137
GO:0031267 GO:0031267                        small GTPase binding 425/15137
GO:0019787 GO:0019787 ubiquitin-like protein transferase activity 392/15137
             BgRatio       pvalue     p.adjust       qvalue
GO:0050839 499/17697 2.594115e-26 3.107749e-23 2.348357e-23
GO:0045296 331/17697 6.634980e-19 3.974353e-16 3.003202e-16
GO:0004674 439/17697 5.250043e-17 2.096517e-14 1.584224e-14
GO:0003779 431/17697 2.206193e-14 6.607548e-12 4.992963e-12
GO:0031267 443/17697 4.257279e-13 8.551020e-11 6.461539e-11
GO:0019787 407/17697 4.282648e-13 8.551020e-11 6.461539e-11
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  geneID
GO:0050839 LASP1/ITGAL/ITGA3/CX3CL1/BAIAP2L1/PROM1/ADAM22/TENM1/CD9/NISCH/LYPLA2/ANLN/GPRC5A/IGF1/MRE11/NRXN3/SNX1/IBSP/VCL/TIMP2/CASR/CDH1/CDH10/CAPG/CTNNA1/HSPA5/DSG2/LIMA1/LAMA3/PKP2/CDH3/CCN5/CNN2/SLC9A3R2/PKN2/ZC3H15/SLK/CTNNA2/PFKP/PKM/DHX29/CDH19/ADGRL1/ACTN1/EPN2/MARK2/PVR/ADAM11/PICALM/NOTCH3/CDHR2/ENO1/FSCN1/DLG1/ICAM3/ACTN2/CAPZB/FAP/CDH17/PTPRH/PSEN1/CDH7/PKP1/BZW1/COBLL1/ITGB5/ERC1/CHMP2B/COL16A1/RAB10/CTTN/EPS15/CHMP5/ADD1/DOCK9/SNX5/RPL6/GCN1/PXN/FXYD5/LTBP4/NUDC/ICAM1/GOLGA3/LAMB1/ITGA6/EZR/TBC1D2/HSP90AB1/DSP/SH3GLB1/CDHR5/MADCAM1/ARVCF/RANBP1/CRKL/TBC1D10A/MICALL1/PACSIN2/MYH9/RANGAP1/AHSA1/EIF5/NOP56/MAPRE1/CHMP4B/CDH20/VAPA/EMD/STK24/OLFM4/UBFD1/EHD4/TJP1/EIF3E/CCN4/NDRG1/EEF1D/ADAM2/PPP1R13L/FCER2/PLIN3/ICAM5/COMP/ITGB8/PTN/ZC3HAV1/PTPRZ1/YKT6/EIF4H/ECM2/PDLIM1/CXCL12/PFN1/ICAM2/LRRC59/YWHAE/VTN/CPE/RPL34/HSPA8/EHD1/NRXN2/EIF4G2/DDX6/CBL/NECTIN1/CD81/VWF/PPFIBP1/KRT18/PTPN6/STK38/CCN6/CLINT1/THBS4/CDH6/FGF1/RARS1/DBN1/RPL24/ADAM23/ACVR1/ITGB6/ITGA4/SPTBN1/RTN4/FN1/STAT1/TXNDC9/HDLBP/PARK7/CAPZA1/LGALS8/IGSF21/CD46/PRDX1/CNN3/PRDX6/ESYT2/CCN2/SPP1/PPL/TJP2/ITGB1BP1/CTNNAL1/EPCAM/PCMT1/TNN/TGFBI/TMPO/TRIM25/SEPTIN7/CALD1/GIPC1/VAPB/CDH26/LRRFIP1/AHNAK/IL1B/VASP/MMP24/KTN1/TSPAN8/PTPRB/EPS15L1/MACF1/PAICS/KDR/CDC42EP1/CDHR3/MYO1B/CDH15/NECTIN2/SNX9/AFDN/ACTN4/COL5A1/PAK4/LAMA5/EIF2S3/ARHGEF16/EPS8L1/LILRB2/DNAJB1/CC2D1A/ARFIP2/RAN/ITGB4/ANGPTL3/SLC14A2/MPRIP/TRPC4/POSTN/SWAP70/NUMB/LDHA/ARGLU1/ADAMTS8/ANXA1/CEMIP2/P2RX4/TES/FLNB/LCP1/GPNMB/BZW2/DBNL/CIB2/ABI1/NIBAN2/DAB2IP/CCN3/UBAP2/TLN1/TRIM29/RDX/THBS1/PAK6/ADAM10/FNBP1L/RAB1A/EMILIN1/ANXA7/ATIC/IDH1/OLA1/ITGAV/USP8/TMOD3/SEMA7A/FGF2/USO1/ITGB7/DIAPH3/CDH24/FBLN5/MFGE8/UNC45A/IQGAP1/CDH11/CDH13/RPS2/CSNK1D/SH3GL1/SCYL1/PSMB6/EPHA2/EFHD2/SERBP1/CCN1/PTPRF/NECTIN4/HMCN1/CGN/ADAM15/PKP4/PHLDB2/EIF2A/SFRP2/CDH18/TENM2/ARHGAP18/EGFR/NLGN4X/MSN/ITGB1BP2/RPL7A/SH3GLB2/LRRC4C/TNKS1BP1/PTPRJ/TENM4/CDH22/FERMT3/PLCB3/ALDOA/MPP7/ITGB1/CDH8/CD226/ANK3/TWF1/PTPRO/ADAM8/ADAM17/DST/BAG3/UTRN/CAST/ASAP1/PTPRD/THY1/CDH12/PRKCA/CXADR/ADAMTS5/LARP1/CCT8/NPTN/NRG1/MMP14/CCNB2/FMNL2/NCK1/WASF2/CD1D/TAGLN2/F11R/LAD1/ITGB2/ADAMTS13/JAML/ABCF3/ITGA5/STX5/RPL29/NTNG1/VCAM1/PDLIM5/S100A11/CIP2A/S100P/ITGA2/EDIL3/ESM1/CAMLG/YWHAZ/SYK/HNRNPK/GAPVD1/NLGN4Y/JAM3/FBN1/ILK/YWHAB/GOLGA2/IGF2/EEF2/GLOD4/EVPL/SLC3A2/CTNNB1/VASN/SEPTIN2/TNXB/ATXN2L/COL3A1/ADAM9/NPNT/STXBP6/COL4A3/PCBP1/FASN/NLGN1/CTNND2/NLGN2/CDH2/KIF5B/RSL1D1/FGG/FGA/FGB/FRMD5/LAMB2/BSG/PPP1CA/HCFC1/CORO1B/PTPRM/JUP/SPTBN2/CCS/RPL15/KLC2/CKAP5/ARHGAP1/PTPN2/RUVBL1/SFN/BAIAP2/CDK5R1/EPS8L2/CD151/NECTIN3/PLEC/CALR/CDH4/PTPN11/CDH5/PDXDC1/NRXN1/PUF60/EXOC3/PAK2/MB21D2/SCRIB/IST1/ANXA2/CTNNA3/TACSTD2/PKP3/SEPTIN9/H1-10/RAB11B/P4HB/MRTFB/SHTN1/ISG15/RPL14/HMGB1/PTPRT/NLGN3/PTPN1/FLNA/SRC/H3C12/SND1/H3C4/DCHS2/SPTAN1/PARVA/RPS26/CD2AP/RPL23A/ITGBL1/CTNND1/EGFL6/GIGYF2/BMPR2/HSPA1A/RACK1/CD177/SNX2/S1PR3/CLIC1/EMP2/TSPAN4/PPME1/DDX3X/TENM3/PI4KA/TWF2/LYN/EEF1G/ITGB3/S1PR2/H3C8/H3C6/H3C11/H3C1/H3C7/H3C10
GO:0045296                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               LASP1/BAIAP2L1/PROM1/LYPLA2/ANLN/GPRC5A/MRE11/SNX1/VCL/CDH1/CDH10/CAPG/CTNNA1/HSPA5/LIMA1/PKP2/CDH3/CNN2/SLC9A3R2/PKN2/ZC3H15/SLK/CTNNA2/PFKP/PKM/DHX29/CDH19/EPN2/MARK2/PICALM/NOTCH3/ENO1/FSCN1/DLG1/CAPZB/CDH17/PTPRH/PSEN1/CDH7/PKP1/BZW1/COBLL1/ERC1/CHMP2B/RAB10/CTTN/EPS15/CHMP5/ADD1/DOCK9/SNX5/RPL6/GCN1/FXYD5/NUDC/GOLGA3/ITGA6/EZR/TBC1D2/HSP90AB1/SH3GLB1/ARVCF/RANBP1/CRKL/TBC1D10A/MICALL1/PACSIN2/MYH9/RANGAP1/AHSA1/EIF5/NOP56/MAPRE1/CHMP4B/CDH20/VAPA/EMD/STK24/OLFM4/UBFD1/EHD4/TJP1/EIF3E/NDRG1/EEF1D/PPP1R13L/PLIN3/ZC3HAV1/YKT6/EIF4H/PDLIM1/PFN1/LRRC59/YWHAE/RPL34/HSPA8/EHD1/EIF4G2/DDX6/CBL/PPFIBP1/KRT18/STK38/CLINT1/CDH6/RARS1/DBN1/RPL24/ACVR1/SPTBN1/RTN4/STAT1/TXNDC9/HDLBP/PARK7/CAPZA1/CD46/PRDX1/CNN3/PRDX6/ESYT2/PPL/TJP2/CTNNAL1/EPCAM/PCMT1/TMPO/TRIM25/SEPTIN7/CALD1/GIPC1/VAPB/CDH26/LRRFIP1/AHNAK/VASP/MMP24/KTN1/PTPRB/EPS15L1/MACF1/PAICS/KDR/CDC42EP1/CDHR3/MYO1B/CDH15/SNX9/AFDN/PAK4/EIF2S3/ARHGEF16/EPS8L1/DNAJB1/CC2D1A/ARFIP2/RAN/MPRIP/TRPC4/SWAP70/NUMB/LDHA/ARGLU1/ANXA1/CEMIP2/P2RX4/TES/FLNB/BZW2/DBNL/ABI1/NIBAN2/DAB2IP/UBAP2/TLN1/TRIM29/RDX/PAK6/FNBP1L/RAB1A/ATIC/IDH1/OLA1/USP8/TMOD3/USO1/DIAPH3/CDH24/UNC45A/IQGAP1/CDH11/CDH13/RPS2/CSNK1D/SH3GL1/SCYL1/PSMB6/EPHA2/EFHD2/SERBP1/CGN/PKP4/PHLDB2/EIF2A/CDH18/ARHGAP18/EGFR/RPL7A/SH3GLB2/TNKS1BP1/PTPRJ/CDH22/PLCB3/ALDOA/MPP7/ITGB1/CDH8/ANK3/TWF1/PTPRO/BAG3/CAST/ASAP1/CDH12/LARP1/CCT8/CCNB2/FMNL2/NCK1/WASF2/TAGLN2/F11R/LAD1/ABCF3/STX5/RPL29/PDLIM5/S100A11/CIP2A/S100P/YWHAZ/HNRNPK/GAPVD1/YWHAB/GOLGA2/EEF2/GLOD4/EVPL/SLC3A2/CTNNB1/VASN/SEPTIN2/ATXN2L/STXBP6/PCBP1/FASN/CTNND2/CDH2/KIF5B/RSL1D1/BSG/PPP1CA/HCFC1/CORO1B/PTPRM/JUP/SPTBN2/CCS/RPL15/KLC2/CKAP5/ARHGAP1/RUVBL1/SFN/BAIAP2/CDK5R1/EPS8L2/PLEC/CDH4/CDH5/PDXDC1/PUF60/EXOC3/PAK2/MB21D2/SCRIB/IST1/ANXA2/CTNNA3/TACSTD2/PKP3/SEPTIN9/H1-10/RAB11B/MRTFB/SHTN1/RPL14/PTPRT/PTPN1/FLNA/SRC/H3C12/SND1/H3C4/DCHS2/SPTAN1/PARVA/RPS26/CD2AP/RPL23A/CTNND1/GIGYF2/BMPR2/HSPA1A/RACK1/SNX2/CLIC1/PPME1/DDX3X/PI4KA/TWF2/EEF1G/H3C8/H3C6/H3C11/H3C1/H3C7/H3C10
GO:0004674                                                                                                                                                                                                                                                                                                                                                                 CAMKK1/MAP3K14/MAP3K9/CDKL3/MARK4/CDKL5/CAMK1G/CDK11A/DYRK4/MAP4K3/MAP4K5/CLK1/MNAT1/PRKCH/VRK2/MAP2K3/DAPK2/TRIO/PHKA2/LTBP1/MAPK9/EIF2AK2/CDK14/CAMK2B/RIOK2/CDK17/WNK1/TNK2/ADCK1/HIPK2/PKN2/MYLK/MAP2K4/SLK/PRKCQ/CDK13/PHKA1/PRKCZ/ROCK1/PDK3/MAST4/TGFBR3/MAPK6/TESK2/CSNK2A2/CAMK2A/MAP4K4/RPS6KA2/MYO3B/PRKACA/RPS6KA6/SPEG/MARK2/STK10/ALPK1/MAP3K13/NUAK1/MARK3/MAP2K7/PAK3/ARAF/MKNK1/MOK/STK17B/STRADB/GSK3B/ULK2/MAP3K4/CPNE3/MAST2/EIF2AK1/NLK/AURKA/MAPKAPK5/LTBP4/CCNK/IRAK3/MAP3K20/TRPM7/MAP3K1/MYO3A/SRPK1/CDC7/MAST3/MKNK2/MAPK1/GRK3/CDKL1/VRK1/RPS6KA5/SGK2/STK4/CSNK2A1/MYLK2/PAK5/RIOK3/PIM2/CDK16/CAB39L/STK24/MAPK3/PHKB/EEF2K/BCKDK/SGK3/C8orf44-SGK3/RIPK2/IKBKB/STK3/MAP4K1/ERCC2/DMPK/VRK3/AURKC/DYRK1B/AKT2/PRKD2/MAST1/GSK3A/CDK6/PIK3CG/PRKAG2/LIMK1/TGFBR1/TESK1/MAPK8/BMPR1A/MAP3K8/RPS6KB1/MAP2K6/MAPK10/HIPK3/GTF2H1/CAMKK2/GTF2H3/PRKAB1/MAK/MAPK14/STK38/ICK/CCND3/PRPF4B/TTK/CLK4/CSNK1A1/NEK11/MAPKAPK3/ACVR2B/NEK4/ACVR1/PRKAG3/STK16/PASK/STK25/PRKD3/AAK1/MARK1/AKT3/CDK18/NEK2/RPS6KA1/STK11/SGK1/TRPM6/NEK6/NEK9/MASTL/PTK2B/PIK3CA/ACVR2A/CIT/PKN1/CDK2/NRK/ACVR1C/RIOK1/GRK4/STK35/WNK4/PRKCG/MAP2K2/DYRK2/PKMYT1/EIF2AK4/TTBK2/RIPK3/STK33/PAK4/MAP3K10/PNCK/LATS1/PRKAB2/RAF1/PRKAA1/CDK8/DSTYK/DCLK1/CSNK1G2/ADCK2/CDK7/IRAK2/CAMK1/ROCK2/ERN2/CCNH/STK26/TAOK3/SRPK2/MAP3K7/AMHR2/CDK4/ACVR1B/CAB39/NEK3/ALPK3/RPS6KC1/CDK9/PIM1/RIPK1/NEK1/MAP2K5/ITPKA/PAK6/CDK15/PRKG2/BMPR1B/BMP2K/CDKL2/ACVRL1/MAP3K12/TSSK4/ULK3/MYLK3/PDPK1/MINK1/CSNK1D/MAPK4/BRD4/HUNK/SIK1/AKT1/LMTK3/PLK4/MAP3K6/PRKACB/DYRK3/MAP3K21/CDC42BPA/DCAF1/CAMK2D/PLK2/GTF2H2/MYLK4/TTBK1/EGFR/TLK2/CASK/TAF1/CAMK2G/PAK1/ATM/CHEK1/TAOK2/LATS2/CSNK1G3/NEK7/UHMK1/CAMK4/STK32B/PRKCA/LRRK1/TNIK/OBSCN/CDK19/TTN/CDK20/MAPK13/PHKG2/BUB1B/SMG1/DYRK1A/BRAF/PINK1/PSKH1/KALRN/HIPK4/PKN3/BRSK1/TAOK1/SIK3/NEK8/SQSTM1/RPS6KA4/PRKAA2/TSSK3/MAPKAPK2/ERCC3/HIPK1/STK36/NEK10/TGFBR2/NUAK2/PRKCI/DCLK3/SNRK/PRKCD/CAMKV/STK17A/LMTK2/PHKG1/CDK5/FASTK/SYK/WNK2/MELK/STK32C/ILK/MAPK7/PRKCB/PLK1/ATP23/CDK12/RSKR/DAPK3/ULK4/MAP4K2/PBK/MLKL/HJV/MAP2K1/CSNK1G1/STK32A/BUB1/MAP3K2/SIK2/ANKK1/CDK1/DCLK2/AKAP13/PPM1D/PRKCE/CDC42BPG/SOSTDC1/KSR2/EIF2AK3/TP53RK/OXSR1/GRK2/ADCK5/MAP3K11/PLK3/BRSK2/ATR/PDIK1L/RPS6KB2/CLK2/MAP3K19/ULK1/RPS6KA3/NIM1K/HASPIN/FAM20C/TSSK6/ERN1/GAK/AURKB/CLK3/PAK2/MAP3K15/MAPK15/AATK/PRKAG1/LIMK2/CAMK1D/RIPK4/TBK1/CHEK2/PRKX/IRAK1/PRKD1/SRPK3/NRBP2/CDK10/MAPK11/PRKG1/BCR/SBK2/MAPK12/SBK1/LRRK2/STK40/STK31/PIK3R4/WNK3/DAPK1/NEK5/MAP3K5/IRAK4/GRK6/RPS6KL1/PIM3/TLK1/STK39/CDC42BPB/MTOR/ALPK2/STKLD1/GRK5/TOP1/MAP3K3/BMPR2/STK19/CSNK2B/CDKL4/STK38L/CHUK/GTF2H4/CSNK1E/SBK3/NME2/CDK11B/CDK3/PRKDC/IKBKE/STRADA/SIK1B/CCL3
GO:0003779                                                                                                                                                                                                                                                                                                                                                                                                                                                          LASP1/BAIAP2L1/USH1C/SCIN/MYH13/GAS7/SYN1/ANLN/WAS/TMSB10/VCL/MYOM2/MYO16/CAPG/CTNNA1/LIMA1/SYNE2/CAMK2B/CROCC/LIMCH1/CNN2/MYLK/CTNNA2/MYO9A/PFN2/SPTB/WDR1/VASH1/MYO3B/ACTN1/TRPC5/FERMT2/ADD2/FSCN1/KLHL20/CAMSAP3/ACTN2/CAPZB/NEBL/MYH7B/TNS1/CEACAM1/OPHN1/EPB41L2/EPB41L3/COBLL1/SSH1/CTTN/MYBPC2/ADD1/PHACTR3/EPB41L1/FXYD5/MAEA/MYO15A/MYH7/TRPM7/EZR/HDAC6/SORBS1/MYO3A/ABL1/ABLIM1/MYO9B/MISP/TRIOBP/PICK1/MYH9/DAAM1/FERMT1/SNTA1/TNNC2/MYOM1/PLS3/EMD/KLHL4/CORO1A/COTL1/MEFV/SLC6A2/CORO2B/HOMER2/VPS18/BLOC1S6/NCALD/MYH14/COBL/WASL/ACTR3C/BCL7B/CORO2A/PDLIM1/PFN1/SLC6A4/PPP1R9B/MYH1/MYH3/KLHL2/FRG1/KLHL5/PANX1/CORO1C/TRPV4/ARPC3/MYL2/TULP1/PHACTR1/CAP2/WASF1/PHACTR2/PACRG/DBN1/TNNC1/ACTR3/SPTBN1/CCDC88A/MLPH/WIPF1/CAPZA1/GBP1/CNN3/KPTN/TNNT2/MED28/FKBP15/CTNNAL1/ABITRAM/MYOT/PLS1/PDLIM2/KIF18A/HPCA/CXCR4/LDB3/WIPF3/CALD1/GIPC1/OBSL1/POF1B/MYH2/VASP/DSTN/AIF1L/MACF1/VIL1/HIP1/YWHAH/FLNC/MYO1B/MYO5C/TMOD2/INO80/PALLD/CSRP3/AJUBA/TNNI3/CNN1/AFDN/ACTN4/ARPC1B/MAP1S/LSP1/TNNT3/TNNI2/GMFG/HIP1R/SYNE1/EPS8L1/CAP1/PDLIM4/DIAPH1/MAP1B/TNS4/MTSS2/WASF3/MYH8/MYH10/MPRIP/MYH11/MYO18B/ACTR3B/MICAL2/SPIRE1/MYBPC3/HOOK1/FHOD3/AVIL/MICAL1/FHOD1/VILL/FLNB/LCP1/DBNL/MYO1G/ERMN/BIN1/TMOD1/ARPC5L/TLN1/FHDC1/MYO7A/TRPC6/RDX/SPTBN5/ACTR2/MYPN/ITPRID2/TMOD3/SHROOM3/PARVG/FGD4/GAS2L3/DIAPH3/PSTPIP1/TLNRD1/TPM1/SSH2/MYO1F/MYOM3/ABL2/CGN/VASH2/TPM3/ESPNL/MYH15/GC/SNCA/MYO10/IQGAP2/KLHL3/PPP1R18/DAAM2/EGFR/SHROOM2/MSN/DIAPH2/FBXO25/ZNF185/GSN/MSRB2/ADD3/TAGLN/ALDOA/ITGB1/KLHL1/DIXDC1/TWF1/EPS8/DST/PSTPIP2/JMY/PLEKHH2/UTRN/HNRNPU/IMPACT/ENAH/PDLIM3/TMSB4Y/TTN/FMN2/ADCY8/KCNMA1/WHAMM/MYO1E/FMNL2/TMSB15A/WASF2/SHROOM4/PPP1R9A/DMTN/EPB41/TNNI1/PIP/SPTBN4/RUSC1/CCR5/MYL3/ALKBH4/FMNL3/NEXN/ARPC5/XIRP2/PDLIM5/LMOD3/LMOD1/ARPC2/ABLIM2/SHROOM1/MYOZ3/NOS3/MICALL2/CFL2/INPPL1/PKNOX2/CLMN/CACNB2/SYNPO2L/MYO1A/MAP1A/SAMD14/NOD2/MYO5B/TPM4/CORO6/EEF2/PXK/FAM107A/XIRP1/SNTB2/PTK2/ANTXR1/MYO7B/MYRIP/LMOD2/MTSS1/PRKCE/CDK5R2/WIPF2/ENC1/TLN2/SYNPO/SNTB1/MYOZ2/SYNPO2/SNTG2/CORO1B/CFL1/SSH3/ABLIM3/DAG1/SPTBN2/ABRA/GCSAM/MYO1H/MARCKSL1/SYNE3/MYO1D/PFN4/CDK5R1/EPS8L2/PAWR/FLII/MYOZ1/PPP1R42/PLEC/IGSF22/TIGD5/HCLS1/WASHC1/NEB/CTNNA3/FHL3/SETD3/SMTN/SPATA32/FMNL1/ADSS1/PRKN/NF2/FSCN2/MAPT/ESPN/SHTN1/KLHL17/PARVB/TMEM201/LRRK2/MYBPC1/S100A4/EVL/AFAP1/MYO18A/MYO6/MRTFA/PDLIM7/FLNA/ANXA6/GMFB/SVIL/IPP/MIB2/MYO5A/MYH6/SPTAN1/PARVA/MYO1C/NRAP/MYL4/TPM2/MSRB1/EPS8L3/RCSD1/CAPZA2/DMD/INF2/PHACTR4/MACO1/AIF1/SPIRE2/TMSB4X/STK38L/DNASE1/ANG/VPS16/ARPC4/ARPC1A/MICAL3/TWF2/FMN1/SHANK3/CORO7/ANXA8/CYFIP1/MARCKS/PRICKLE4/MYO19
GO:0031267                                                                                              ALS2/DVL2/MYCBP2/FARP2/ARHGAP44/NOX1/CDKL5/PLEKHG6/RANBP9/ANLN/RABGAP1/WAS/NPC1L1/RALBP1/RANBP3/TRIO/CLEC16A/ADRB1/PREX2/ARHGEF5/SIKE1/MCF2L2/TRAPPC3/TBC1D22A/GDI2/RASGRF1/YIPF1/YBX3/NCKAP1/SPHK2/NGFR/IPO5/PKN2/TBC1D22B/TBC1D1/NGEF/GOLGA5/EVI5/ROCK1/TBC1D25/EXOC5/ATP6AP1/LLGL2/PICALM/NSF/ARHGEF10L/RAP1GAP/ARHGEF1/USP33/PAK3/ADCYAP1R1/RAPGEF3/RIMS1/ERC1/XPO1/TNPO1/HACE1/IPO11/DNM1L/DOCK3/RPH3A/KIF16B/BIRC5/ARHGAP4/RAB11FIP3/PLEKHG2/EXOC1/RAPGEF4/UNC13D/NUP50/TBC1D2/MYO9B/RANBP1/LZTR1/TBC1D10A/GGA1/SH3BP1/HPS4/TRIOBP/MICALL1/SBF1/SGSM3/RANGAP1/SOS2/DAAM1/RIN3/RABGGTA/RIMS4/PAK5/KIF3B/MCF2/FGD1/SYTL4/ARHGEF7/NUTF2/PARD6A/TSC2/GGA2/DMXL2/NDRG1/ARHGEF10/ARHGEF18/TBC1D17/PPP6R1/ANKRD27/DENND3/TNPO2/CAV1/TBC1D13/RIC1/RGP1/RAPGEF1/DVL1/DNMBP/RAB11FIP2/TBC1D12/KPNB1/PFN1/SLC6A4/RANGRF/TNFAIP1/IFT20/RAB34/TBC1D9/RAPGEF2/EHD1/ARHGEF17/MADD/EXPH5/CORO1C/KCTD10/ARHGDIB/TBC1D30/WASF1/EXOC2/RASGRF2/ACAP2/ECT2/PEX5L/ARHGEF26/RTKN/MLPH/RAB3GAP1/SOS1/ERRFI1/RAP1A/ARHGEF2/DOCK7/NCF2/EXOC8/RIMS3/RAB29/IPO13/RAB3GAP2/MAPKAP1/DENND1A/PLEKHG1/TBC1D15/OCRL/PKN1/OPTN/NCKAP1L/PREX1/PARD6B/CSE1L/XPO5/GGA3/TBC1D20/MCF2L/PLEKHG3/FGD3/RAB3IP/CDC42EP1/DOCK4/SERGEF/USP6/ARHGEF6/SH3BP4/XPO7/AFDN/PAK4/YIPF2/ARHGEF16/EPS8L1/ARHGEF9/RAB11FIP4/SH3BP5/TBC1D5/RBSN/KIF3A/DIAPH1/EXOC4/WDR44/LLGL1/ARFIP2/TBC1D14/MTSS2/NXT1/RIN2/ARHGEF11/TBC1D8B/IPO8/SBF2/VAV3/ROCK2/DOCK2/BICDL1/MICAL1/RAB11FIP5/TMEM127/ARHGEF4/TBC1D4/RCBTB2/RAPGEF5/RAC1/RALGPS1/ARHGEF39/DENND4C/SYTL2/SORL1/ITPKA/PAK6/RABGGTB/PLCE1/ABI2/RASGEF1B/EGF/FGD4/PEX5/DIAPH3/GRTP1/IQGAP1/ARHGAP17/GAS8/SGSM2/ARHGDIA/VAV1/SOD1/EVI5L/ARHGEF19/SYTL1/STRIP1/RGL1/ADPRH/STXBP5L/TBCK/IQGAP2/TBC1D7/TBC1D7-LOC100130357/DAAM2/FGD2/TIAM2/DENND2A/SYTL5/DIAPH2/DOCK11/C9orf72/USP6NL/PAK1/CDC42EP2/DOCK1/PIH1D2/RILPL2/EPS8/BICD1/RABGAP1L/RASGRP3/FARP1/RGPD3/RANBP2/PLEKHG4B/OBSCN/RABGEF1/FGD5/WHAMM/TIAM1/RAB11FIP1/PDE6D/DGKI/FMNL2/CDC42SE2/RAPGEF6/ABR/KALRN/VAV2/PKN3/DVL3/FMNL3/GRASP/BICDL2/DENND2D/CDC42EP3/ARHGEF3/RANBP3L/WDR41/RHOBTB3/STXBP5/SYTL3/MICALL2/ABCA1/ATP7A/ARHGEF40/HPS6/DENND2B/AP1G1/SGSM1/TBC1D2B/TBC1D16/MYO5B/RAB8A/CDC42EP5/DAPK3/RILP/RAB3IL1/RHOH/STXBP6/XPO6/TBC1D10B/RGPD8/MYRIP/TRAPPC1/DENND5B/AKAP13/PLEKHG5/RASGRP4/KNDC1/SRGAP2C/RASGRP1/MAP3K11/NET1/PIFO/DENND4A/RIN1/DENND6A/KCTD13/SH3BP5L/ARHGAP1/TBC1D10C/DENND2C/RIMS2/RNF152/SMCR8/EPS8L2/ULK1/PARD6G/RCC2/CDC42EP4/RCC1/FGD6/PAK2/TRAPPC5/RPH3AL/RNF41/ANXA2/SPATA13/ARHGEF37/IQGAP3/DENND5A/XPOT/FMNL1/AP3M1/CIB1/RGPD2/ANKFY1/BICD2/BCR/PLEKHG7/SBK2/RILPL1/C15orf62/CHM/NOXA1/LRRK2/PLEKHG4/SRGAP3/SRGAP2B/IPO4/TRAPPC4/STRN3/ARHGEF12/FLNA/SRGAP1/TBC1D9B/MYO5A/CDC42SE1/MYO1C/ITSN2/IPO9/UNC13B/CDC42BPB/EPS8L3/DENND4B/ARHGEF15/RUSC2/RUNDC1/RASGEF1A/INF2/CHML/ECT2L/GDI1/STK19/LSM2/TBC1D8/RANBP17/IPO7/RGL3/DENND6B/ITSN1/DENND1C/DENND1B/ARHGEF33/ARHGEF28/VPS52/ARHGEF38/RGL2/ARHGEF25/ARHGDIG/MICAL3/BRK1/SRGAP2/RASSF5/CYFIP1/TBC1D3D/TBC1D3/TBC1D3E
GO:0019787                                                                                                                                                                                                                                                                                                                                                                                                                                                                                ANKIB1/HECW1/KLHL13/MYCBP2/FBXL3/ASB4/MYLIP/UBE3C/RNF216/BRCA1/UBR7/RNF14/UFL1/RNF10/BIRC3/UBR2/PIAS1/RNF19A/CUL3/TSPAN17/ERCC8/NEDD4L/ANAPC4/RC3H2/RNF4/ASB1/PDZD4/HDAC4/NEDD4/RNF126/CDC42/RAD18/SEL1L/HLTF/LNX1/UBE2D1/CHFR/FBXW11/KLHL20/UBE2T/UBE2A/PIAS2/UBE2K/ITCH/UBE2D4/CNOT4/RNF13/HACE1/ATG16L1/NFX1/HUWE1/KLHL42/ANAPC5/BRAP/MAEA/MUL1/RNF31/G2E3/HECTD1/RFFL/CDC23/MARCHF2/CDC34/RNF215/PPIL2/FBXO7/RBX1/TRIM9/ASB2/BIRC7/RNF24/SEL1L2/RNF125/MIB1/XIAP/UBL4A/MGRN1/STUB1/UBE2I/BFAR/RNF40/HERC1/UBE2W/UBR5/PIAS4/DMAC2/CBLL1/FBXO24/UBE2R2/FBXL15/NEURL1/UBE2S/ZMIZ1/FBXL20/RNF43/TRIM37/RNF167/MED31/PEX12/SMURF2/TNFAIP1/UBE2D3/TRIM2/DTX4/PRPF19/TRIM3/RNF141/BIRC2/UBE4A/CBL/FBXO3/KCTD10/RNF8/FBXO9/ZNF451/FBXL4/RNF130/SKP1/TRIM23/UBE3A/RNF7/CBLB/FANCL/BIRC6/RNF19B/TRIM62/FBXO2/FBXO6/UBE3D/FBXO30/TNFAIP3/RNF146/FBXL5/UBE2B/TRIM32/FBXW2/IRF2BPL/AREL1/RNF170/TRIM25/TRIM6/PDZRN3/RNF2/RBBP6/ZMIZ2/TRIM24/EGR2/RNF11/WWP1/NEURL2/MED20/LRRC29/RNF113A/MED1/RBCK1/HECTD3/TRAF2/FBXL12/UBR4/FBXL16/RNF6/RNF112/HERC2/UBE2M/TRIM28/UBE4B/RLIM/TRAF3/UBE2D2/TRAF7/PIAS3/TRIM21/TRIM5/UBE2G1/TRIM47/FBXO44/RNF128/MED10/MKRN1/RNF122/VHL/DDB2/RNF138/FBXO21/MDM2/FBXL8/RC3H1/LMO7/RNF38/RNF144B/RNF121/FBXO11/BARD1/HECW2/HERC3/HERC6/HERC5/RNF185/PELI2/UBE2Q2/GID4/ANAPC11/RNF157/CBX4/RNF165/FEM1A/CBLC/NOSIP/COP1/UFC1/DTL/UBR3/MARCHF4/UBA3/ATG3/MARCHF1/RNF175/MARCHF6/SKP2/FBXL17/RNF145/RMND5B/RNF217/SHPRH/TAF1/FBXO25/UHRF2/FBXO10/LRSAM1/HERC4/UBE3B/RNF144A/AMN1/FBXO4/ATG10/TRIM36/MED21/ANAPC1/RANBP2/FBXL2/RMND5A/TRIP12/TRIM11/SH3RF1/ELOC/RABGEF1/FBXL18/RNF20/MED7/SH3RF2/UBE2L6/NSMCE2/RNF111/TRIM63/UBE2Z/UBR1/AMFR/UBE2J2/MED27/UBE2Q1/FBXL13/FBXO27/MED11/SYVN1/KLHL21/TRIM58/TRIM17/NEURL3/RNF149/DCST1/RNF25/RCHY1/FBXO40/DTX3L/RNF168/RNF123/RNF180/MED30/RNF183/HECTD2/MARCHF8/PDZRN4/BTRC/ARIH1/CUL5/RAG1/FBXO22/RNF214/TRIM68/ZNF598/KIAA1586/RNF187/RFWD3/RNF181/TMEM129/FEM1B/UBE2V2/NSMCE1/TRIM56/UBE2E3/UBE2E1/RNF150/IRF2BP1/RNF34/RNF139/FBXL14/MALT1/ZNF738/SH3RF3/HECTD4/RNF26/RNF213/MARCHF3/PELI3/MSL2/KCTD13/FBXW8/UBE2C/TRAF6/UBE2O/KCMF1/RNF152/TRIM72/ARIH2/UBE2N/RNF212/RNF186/SHARPIN/ZNRF2/RNF182/PJA1/RNF135/SIAH2/RNF41/UBA7/UBE2E2/FBXL6/TTC3/BCOR/ZNRF3/FBXL7/TRAIP/UBE2F/MED12/UBE2G2/UBOX5/PRKN/UBE2L3/TRIM69/ZNRF1/UBE2H/ZFP91/RNF220/NHLRC1/NCCRP1/NHLRC3/RNFT1/WDSUB1/SIAH1/TRIM33/PELI1/FBXL22/MIB2/TOPORS/MARCHF5/WWP2/KLHL9/SMURF1/UBE2J1/LTN1/ASB12/DZIP3/PJA2/RING1/RNF5/PPP1R11/TRIM27/TRIM13/RNF208/TRIM59/NEURL1B/ZBED1/NEURL4/UBE2QL1/SIAH3/RNF148/RNF103/UBE2V1/MAGEL2/RNF115/FBXO17/UHRF1
           Count
GO:0050839   493
GO:0045296   328
GO:0004674   427
GO:0003779   416
GO:0031267   425
GO:0019787   392
#write.csv(bp.go,file="mf.go.csv")

Perform a KEGG pathway enrichment analysis

kk <- enrichKEGG(gene         = deg_ann.tb$entrezgene_id,
                 organism     = 'hsa',
                 pAdjustMethod = "BH",
                 qvalueCutoff  = 0.05)
head(kk)
               ID                    Description GeneRatio  BgRatio
hsa04510 hsa04510                 Focal adhesion  198/6809 199/7946
hsa04360 hsa04360                  Axon guidance  180/6809 181/7946
hsa04144 hsa04144                    Endocytosis  243/6809 249/7946
hsa04010 hsa04010         MAPK signaling pathway  284/6809 295/7946
hsa04722 hsa04722 Neurotrophin signaling pathway  119/6809 119/7946
hsa04150 hsa04150         mTOR signaling pathway  151/6809 153/7946
               pvalue     p.adjust       qvalue
hsa04510 1.041735e-12 3.354385e-10 1.546153e-10
hsa04360 1.643594e-11 2.646186e-09 1.219720e-09
hsa04144 9.927879e-11 1.065592e-08 4.911687e-09
hsa04010 1.086962e-09 8.750041e-08 4.033200e-08
hsa04722 9.004526e-09 5.798915e-07 2.672922e-07
hsa04150 1.566294e-08 8.405780e-07 3.874518e-07
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    geneID
hsa04510                                                                                                                                                                                                                                                                                                                                                                                                                                                     572/3675/3674/2534/3479/3082/330/3381/7414/2324/3371/1298/10319/5601/9564/5291/3909/3918/4659/5923/4638/387/6093/998/87/60/4660/5063/8516/1286/2932/3693/5829/22798/3912/3655/1299/1399/5594/5155/6655/394/85366/10398/57144/10627/284217/331/2321/5595/208/5296/1311/3696/858/857/4233/93408/2889/5599/1277/7448/5602/595/329/7450/4633/1297/896/7422/3910/7060/5159/3694/3676/2335/6654/7143/5906/10000/8503/103910/6696/894/5228/63923/5290/54776/7408/5582/2002/5908/3791/5880/2318/25759/81/10298/3911/1729/5894/3691/10451/9475/5156/1284/3679/3915/2317/5879/1101/7094/7057/22801/56924/106821730/3685/1950/64098/1280/3695/3480/91807/5170/2064/7409/1291/1292/207/8515/3680/56034/5295/340156/1956/53358/5058/3688/7424/1793/5578/673/7791/2909/7410/6464/3678/824/1293/3673/1278/2277/3611/5579/1398/1499/7148/1285/5604/7059/5747/5881/80310/5293/596/5728/83660/3913/5499/256076/7423/3265/3725/2885/29895/5062/859/71/399694/23396/5501/7058/1282/1287/29780/5649/3908/3914/2316/6714/5154/1288/55742/131873/5500/3672/4636/3690
hsa04360                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      6405/56920/2534/7869/2042/5291/816/2773/9423/25791/4756/387/5590/6093/815/998/2046/8440/7224/4775/59277/10371/10512/5362/5063/285220/54437/2932/54434/1947/57556/1808/57715/25/3983/8829/1943/5594/4776/4773/655/10398/57144/10627/5595/50855/7976/5296/4233/2051/3984/6387/219699/5532/1949/81029/9037/90249/7474/2771/5361/2043/8503/103910/5533/5290/7852/84612/5335/1948/6237/2770/5880/64101/6608/10298/55558/5894/7223/2048/3845/9475/2045/10505/10154/5879/7225/56924/106821730/80031/8482/658/7222/5530/5998/85464/1969/10500/1944/151449/91653/7220/9353/2044/817/5295/5921/2041/56288/818/57689/5058/3688/5163/223117/64221/5578/55740/2047/137970/56896/4690/54361/22854/6259/84448/5364/6469/1020/1073/3611/10501/54910/1942/5747/5881/6091/9723/5293/27289/1072/54961/29984/22885/3265/84552/5781/5062/8633/56963/2242/3985/2049/6586/1946/6092/10509/5727/6585/1630/10507/5535/64218/9901/84628/2050/23654/23365/57522/6714/5336/5365/3897/659/4893/4636/5534/91584/1945/23380
hsa04144                                                                                                           381/163/9265/51534/30845/30846/6642/9135/9765/1175/116984/55040/23327/5662/29924/2263/10565/387/5590/2261/55737/4734/8218/998/22905/7037/9744/7251/5337/7879/832/273/83737/1785/25978/64744/10890/2060/51510/55616/27131/51324/9727/51100/387680/27128/157/3560/55738/128866/8411/50855/8766/30844/9266/858/857/8976/7046/1759/8395/6456/22841/8729/28964/9267/5878/64750/3312/10938/2348/867/10094/51699/5869/11021/60682/8724/23527/868/8723/10254/51652/7456/829/9829/9525/64411/440073/51019/7852/644150/1211/9559/11059/84612/10564/57403/23096/10617/2868/23550/3306/58513/5338/58533/3949/10095/27243/5119/84440/64145/84313/116983/27183/23111/116988/377/3559/5156/116986/26056/4193/23325/51028/11311/274/8027/81873/408/9101/9815/79720/3480/5371/6457/1213/409/93343/6455/23624/8394/375/5868/9922/84249/1956/30011/29934/3561/91782/56904/56288/84364/112936/8853/1601/50807/118813/137492/3798/23362/80223/9372/26119/116987/11267/1234/2264/51160/1173/10092/3577/10109/7048/5584/92421/9897/2350/382/4088/138429/4218/55048/5867/3800/11031/10015/7314/3799/147179/729092/253725/254122/156/3310/3265/7189/4087/1212/79643/155382/80230/84552/3579/100287171/10193/9798/859/161/390243/9230/9146/23396/116985/119016/89853/160/3482/6714/26052/2870/57154/2869/830/414189/642517/3304/3303/3305/3107/3133/3135/3134/6643/3105/3106/10093/10552/26286/29082/57132/119385
hsa04010 9020/8913/9254/23542/5536/8491/3479/3082/355/5606/2324/4254/5601/23118/5923/4804/2065/6416/2263/7132/782/2261/10235/2255/27006/998/9448/6196/5566/27330/4775/9175/2249/27092/5609/4791/2260/369/26281/8569/4208/779/2122/4216/51701/8550/2323/51776/285/7042/4214/1943/4616/2872/1399/5594/5155/10454/8911/6655/5494/9252/6789/994/51378/778/2254/2321/8717/5595/3551/6788/11184/5971/208/7040/8605/4233/3315/7046/5599/5532/1326/1845/786/5608/8649/4790/374/5602/9693/3312/84867/80824/1432/6722/7422/5924/2246/5159/7867/3552/3554/6654/1386/5906/5321/1647/10000/356/3925/6195/5228/7043/1843/7010/1846/5533/9479/2322/3164/2069/3553/4149/6237/5582/2002/3306/5605/7186/5908/3791/468/5880/2318/51295/10912/59285/3727/1852/4772/5894/3845/22800/4803/5156/51347/6885/2317/5879/4609/408/5607/5495/1847/2247/1950/5530/23162/1848/7786/2252/3480/409/7157/2064/773/207/59283/1969/9064/5567/11221/1944/5778/56034/5921/2768/1956/4915/774/5058/9344/7424/93589/775/25780/5801/781/284/5578/5922/8817/5603/776/3815/55799/673/1844/22808/2005/8822/283748/57551/2264/8986/9965/9261/7039/7048/1849/836/2277/783/5598/5579/1398/3481/784/5871/255189/8681/5604/1942/5881/10746/2353/929/80310/3643/115727/55970/10125/4615/5970/3310/4296/7423/3265/7189/1649/627/6197/3725/2885/2066/5062/785/1436/3654/1946/1435/1850/5600/4908/4137/2248/123745/6300/5535/3556/8912/4763/2316/4217/5154/51135/777/4914/4215/1616/3304/3303/3305/4893/1147/5534/4909/100506012/7124/1945/100137049/3630/8517
hsa04722                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     572/4145/5601/5291/816/4804/387/815/998/6196/27330/5609/7161/5663/2932/581/11213/4214/25/1399/5594/6655/9252/4792/5595/8767/3551/4793/208/5296/2889/5599/7531/4790/5602/2549/11108/10019/397/1432/6654/5906/10000/8503/356/6195/2309/5290/5335/5605/5908/468/91860/25759/5894/3845/3656/6272/4803/57498/5879/5607/4916/5170/7157/396/207/5664/805/817/5295/4794/4915/53358/818/814/5603/673/808/6464/10603/9261/5580/10818/5598/27018/1398/5604/3667/163688/10782/5293/596/5970/3265/7189/627/6197/3725/2885/25970/810/51806/9500/5781/3654/5600/399694/4908/6300/4217/5336/51135/4914/801/4215/4893/4909/398
hsa04150                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            51384/1856/3479/58528/51606/5291/57600/1975/7132/387/4040/6196/27330/7479/55437/2932/9706/10325/81929/7481/5594/9681/6655/51382/9894/9675/81617/5595/79726/8131/7249/8140/7976/3551/208/5296/7472/2887/6009/1855/7473/6198/8649/81029/11211/7474/10641/523/7475/6654/525/28956/64121/10000/8503/6195/6794/6446/79109/5290/7471/5582/5605/84219/9296/83667/6249/529/5894/5562/3845/7482/10542/23175/80326/9470/51719/9550/6194/3480/5170/57521/207/7483/245973/6502/5295/51256/107080638/526/55004/1977/127124/5578/89780/7476/201163/528/8324/64798/10670/6396/8321/673/7484/1857/4041/5563/54361/7855/253260/54468/8323/7248/5579/64223/6520/54541/5604/3667/7480/3643/5293/5728/3265/8322/6199/220441/153129/8408/6197/8325/2885/2535/55615/1978/7477/389541/8326/2475/84335/4893/1147/534/96459/7124/652968/90423/3630/92335/729438
         Count
hsa04510   198
hsa04360   180
hsa04144   243
hsa04010   284
hsa04722   119
hsa04150   151
#write.csv(kk,file="enrichedKeggPathways.csv")

Perform a KEGG pathway enrichment analysis

dotplot(kk, showCategory=30) + ggtitle("dotplot for KEGG")

plot of chunk unnamed-chunk-37