文章预览
大豆的数据来源论文 https://www.sciencedirect.com/science/article/pii/S0092867420306188 Pan-Genome of Wild and Cultivated Soybeans 大豆基因组数据下载链接 https://ngdc.cncb.ac.cn/soyomics/download 下载基因组fasta和对应的蛋白注释文件，用gffread提取cds序列和蛋白序列 作图代码 dat.family.group %>% pull(group) %>% table() %>% as.data.frame() %>% magrittr::set_colnames(c("group","Freq")) %>% mutate(group=factor(group,levels=c("Core","SoftCore", "Dispensable","Private"))) %>% arrange(group) %>% mutate(x=cumsum(Freq)+0.5) -> segment.df library(ggrastr) pdf(file = "cell.soybean.PanGenome/fig2c01.pdf", width = 12,height = 6) read_tsv("cell.soybean.PanGenome/Orthogroups.GeneCount.tsv") %>% dplyr::select(-Total) %>% column_to_rownames("Orthogroup") %>% mutate(across(everything(),~ifelse(.>0,1,0))) %>% rownames_to_column("familyID") %>% left_join(dat.family.group) %>%
………………………………
