首先查看可供选择的物种信息
> msigdbr_species()
# A tibble: 20 x 2
species_name species_common_name
<chr> <chr>
1 Anolis carolinensis Carolina anole, green an~
2 Bos taurus bovine, cattle, cow, dai~
3 Caenorhabditis elegans NA
4 Canis lupus familiaris dog, dogs
5 Danio rerio leopard danio, zebra dan~
6 Drosophila melanogaster fruit fly
7 Equus caballus domestic horse, equine, ~
8 Felis catus cat, cats, domestic cat
9 Gallus gallus bantam, chicken, chicken~
10 Homo sapiens human
11 Macaca mulatta rhesus macaque, rhesus m~
12 Monodelphis domestica gray short-tailed opossum
13 Mus musculus house mouse, mouse
14 Ornithorhynchus anatinus duck-billed platypus, du~
15 Pan troglodytes chimpanzee
16 Rattus norvegicus brown rat, Norway rat, r~
17 Saccharomyces cerevisiae baker's yeast, brewer's ~
18 Schizosaccharomyces pombe 972h- NA
19 Sus scrofa pig, pigs, swine, wild b~
20 Xenopus tropicalis tropical clawed frog, we~
使用msigdbr()函数可以得到想要的基因集，其中category参数可以指定基因集的类别，比如“H”或者“C1”
> m_df = msigdbr(
+ species = "Mus musculus",
+ category = "H")
之后可以利用split函数把得到的基因集分组建成一个list
> m_list = m_df %>% split(x = .$gene_symbol, f = .$gs_name)
> str(m_list)
List of 50
$ HALLMARK_ADIPOGENESIS : chr [1:200] "Abca1" "Abcb8" "Acaa2" "Acadl" ...
$ HALLMARK_ALLOGRAFT_REJECTION : chr [1:204] "Aars" "Abce1" "Abi1" "Ache" ...
$ HALLMARK_ANDROGEN_RESPONSE : chr [1:126] "Abcc4" "Abhd2" "Acsl3" "Actn1" ...
$ HALLMARK_ANGIOGENESIS : chr [1:36] "Apoh" "App" "Ccnd2" "Col3a1" ...
$ HALLMARK_APICAL_JUNCTION : chr [1:199] "Acta1" "Actb" "Actc1" "Actg1" ...
$ HALLMARK_APICAL_SURFACE : chr [1:44] "Adam10" "Adipor2" "Afap1l2" "Akap7" ...
$ HALLMARK_APOPTOSIS : chr [1:161] "Add1" "Aifm3" "Ank" "Anxa1" ...
$ HALLMARK_BILE_ACID_METABOLISM : chr [1:112] "Abca1" "Abca2" "Abca3" "Abca4" ...
$ HALLMARK_CHOLESTEROL_HOMEOSTASIS : chr [1:75] "Abca2" "Acat3" "Acss2" "Actg1" ...
$ HALLMARK_COAGULATION : chr [1:139] "A2m" "Acox2" "Adam9" "Ang" ...
$ HALLMARK_COMPLEMENT : chr [1:198] "Actn2" "Adam9" "Adra2b" "Akap10" ...
$ HALLMARK_DNA_REPAIR : chr [1:150] "Aaas" "Ada" "Adcy6" "Adrm1" ...
$ HALLMARK_E2F_TARGETS : chr [1:200] "Ak2" "Anp32e" "Asf1a" "Asf1b" ...
$ HALLMARK_EPITHELIAL_MESENCHYMAL_TRANSITION: chr [1:206] "Abi3bp" "Acta2" "Adam12" "Anpep" ...
$ HALLMARK_ESTROGEN_RESPONSE_EARLY : chr [1:200] "Abat" "Abca3" "Abhd2" "Ablim1" ...
$ HALLMARK_ESTROGEN_RESPONSE_LATE : chr [1:206] "Abca3" "Abhd2" "Acox2" "Add3" ...
$ HALLMARK_FATTY_ACID_METABOLISM : chr [1:158] "Aadat" "Acaa1a" "Acaa2" "Acadl" ...
$ HALLMARK_G2M_CHECKPOINT : chr [1:199] "Abl1" "Amd1" "Arid4a" "Atf5" ...
$ HALLMARK_GLYCOLYSIS : chr [1:200] "Abcb6" "Adora2b" "Agl" "Agrn" ...
$ HALLMARK_HEDGEHOG_SIGNALING : chr [1:36] "Ache" "Adgrg1" "Amot" "Cdk5r1" ...
$ HALLMARK_HEME_METABOLISM : chr [1:199] "Abcb6" "Abcg2" "Ackr1" "Acp5" ...
$ HALLMARK_HYPOXIA : chr [1:200] "Ackr3" "Adm" "Adora2b" "Ak4" ...
$ HALLMARK_IL2_STAT5_SIGNALING : chr [1:199] "Abcb1a" "Adam19" "Ager" "Ahcy" ...
$ HALLMARK_IL6_JAK_STAT3_SIGNALING : chr [1:89] "A2m" "Acvr1b" "Acvrl1" "Bak1" ...
$ HALLMARK_INFLAMMATORY_RESPONSE : chr [1:202] "Abca1" "Abi1" "Acvr1b" "Acvr2a" ...
$ HALLMARK_INTERFERON_ALPHA_RESPONSE : chr [1:100] "Adar" "B2m" "Batf2" "Bst2" ...
$ HALLMARK_INTERFERON_GAMMA_RESPONSE : chr [1:211] "Adar" "Apol6" "Arid5b" "Arl4a" ...
$ HALLMARK_KRAS_SIGNALING_DN : chr [1:202] "Abcb11" "Abcg4" "Actc1" "Adra2c" ...
$ HALLMARK_KRAS_SIGNALING_UP : chr [1:203] "Abcb1a" "Ace" "Adam17" "Adam8" ...
$ HALLMARK_MITOTIC_SPINDLE : chr [1:199] "Abi1" "Abl1" "Abr" "Actn4" ...
$ HALLMARK_MTORC1_SIGNALING : chr [1:201] "Abcf2" "Acaca" "Acly" "Acsl3" ...
$ HALLMARK_MYC_TARGETS_V1 : chr [1:201] "Abce1" "Acp1" "Aimp2" "Ap3s1" ...
$ HALLMARK_MYC_TARGETS_V2 : chr [1:58] "Aimp2" "Bysl" "Cbx3" "Cdk4" ...
$ HALLMARK_MYOGENESIS : chr [1:201] "Ablim1" "Ache" "Acsl1" "Acta1" ...
$ HALLMARK_NOTCH_SIGNALING : chr [1:32] "Aph1a" "Arrb1" "Ccnd1" "Cul1" ...
$ HALLMARK_OXIDATIVE_PHOSPHORYLATION : chr [1:199] "Abcb7" "Acaa1a" "Acaa2" "Acadm" ...
$ HALLMARK_P53_PATHWAY : chr [1:199] "Abat" "Abcc5" "Abhd4" "Acvr1b" ...
$ HALLMARK_PANCREAS_BETA_CELLS : chr [1:40] "Abcc8" "Akt3" "Chga" "Dcx" ...
$ HALLMARK_PEROXISOME : chr [1:103] "Abcb1a" "Abcb4" "Abcb9" "Abcc5" ...
$ HALLMARK_PI3K_AKT_MTOR_SIGNALING : chr [1:105] "Acaca" "Actr2" "Actr3" "Adcy2" ...
$ HALLMARK_PROTEIN_SECRETION : chr [1:96] "Abca1" "Adam10" "Anp32e" "Ap1g1" ...
$ HALLMARK_REACTIVE_OXYGEN_SPECIES_PATHWAY : chr [1:49] "Abcc1" "Atox1" "Cat" "Cdkn2d" ...
$ HALLMARK_SPERMATOGENESIS : chr [1:135] "Ace" "Acrbp" "Acrv1" "Actl7b" ...
$ HALLMARK_TGF_BETA_SIGNALING : chr [1:54] "Acvr1" "Apc" "Arid4b" "Bcar3" ...
$ HALLMARK_TNFA_SIGNALING_VIA_NFKB : chr [1:203] "Abca1" "Ackr3" "Areg" "Atf3" ...
$ HALLMARK_UNFOLDED_PROTEIN_RESPONSE : chr [1:114] "Aldh18a1" "Arfgap1" "Asns" "Atf3" ...
$ HALLMARK_UV_RESPONSE_DN : chr [1:144] "Abcc1" "Acvr2a" "Add3" "Adgrl2" ...
$ HALLMARK_UV_RESPONSE_UP : chr [1:159] "Abcb1a" "Acaa1a" "Ago2" "Alas1" ...
$ HALLMARK_WNT_BETA_CATENIN_SIGNALING : chr [1:42] "Adam17" "Axin1" "Axin2" "Ccnd2" ...
$ HALLMARK_XENOBIOTIC_METABOLISM : chr [1:205] "Abcc2" "Abcc3" "Abcd2" "Abhd6" ...
之后就可以利用这个list为所欲为了！
网友评论