농축 검정

R에서 limma로 하는 차등 발현 분석

John Blischak

Instructor

결과 해석

results <- decideTests(fit2)
summary(results)

topTable(fit2, number = 3)

            symbol entrez  chrom    logFC  AveExpr        t
205225_at     ESR1   2099 6q25.1 3.762901 11.37774 22.68392
209603_at    GATA3   2625  10p15 3.052348  9.94199 18.98154
209604_s_at  GATA3   2625  10p15 2.431309 13.18533 17.59968
                 P.Value    adj.P.Val        B
205225_at   2.001001e-70 4.458832e-66 149.1987
209603_at   1.486522e-55 1.656209e-51 115.4641
209604_s_at 5.839050e-50 4.337052e-46 102.7571

생물학 데이터베이스

KEGG: 교토 유전자·게놈 백과사전
- https://www.genome.jp/kegg/
- 예: 광합성, 단백질 수송
유전자 온톨로지 컨소시엄 (GO)
- http://geneontology.org/
- 예: 스트레스 반응, 발생 과정

농축 검정

	유전자 집합 포함	유전자 집합 미포함
DE	10	90
all	100	900

fisher.test(matrix(c(10, 100, 90, 900), nrow = 2))

Fisher's Exact Test for Count Data

data:  matrix(c(10, 100, 90, 900), nrow = 2)
p-value = 1
alternative hypothesis: true odds ratio is 
                        not equal to 1
95 percent confidence interval:
 0.4490765 2.0076377
sample estimates:
odds ratio
         1

농축 검정

	유전자 집합 포함	유전자 집합 미포함
DE	30	70
all	100	900

fisher.test(matrix(c(30, 100, 70, 900), nrow = 2))

Fisher's Exact Test for Count Data

data:  matrix(c(30, 100, 70, 900), nrow = 2)
p-value = 1.88e-07
alternative hypothesis: true odds ratio is 
                        not equal to 1
95 percent confidence interval:
 2.306911 6.320992
sample estimates:
odds ratio
  3.850476

head(fit2$genes, 3)

          symbol entrez   chrom
1007_s_at   DDR1    780  6p21.3
1053_at     RFC2   5982 7q11.23
117_at     HSPA6   3310    1q23

entrez <- fit2$genes[, "entrez"]

enrich_kegg <- kegga(fit2, geneid = entrez, species = "Hs")

topKEGG(enrich_kegg, number = 3)

                         Pathway    N  Up Down         P.Up       P.Down
path:hsa04110         Cell cycle  115  30   82 6.192773e-01 5.081518e-12
path:hsa05166   HTLV-I infection  233  55  135 8.959082e-01 9.285167e-09
path:hsa01100 Metabolic pathways 1033 350  373 3.175782e-08 9.969693e-01

GO 농축 검정

enrich_go <- goana(fit2, geneid = entrez, species = "Hs")

topGO(enrich_go, ontology = "BP", number = 3)

                             Term Ont    N  Up Down P.Up       P.Down
GO:0002376  immune system process  BP 1935 426  914    1 7.925179e-32
GO:0006955        immune response  BP 1236 230  619    1 3.625368e-29
GO:0045087 innate immune response  BP  645 113  346    1 1.635833e-22

주의 사항

결과를 과도하게 해석하지 마십시오
상향/하향 조절 결과에 대해 신중하게 접근하십시오
배경 유전자 집합은 검정된 유전자만 포함해야 합니다
limma 함수 camera 및 roast 등 더 고급 방법도 있습니다

연습해 봅시다!

R에서 limma로 하는 차등 발현 분석