# R語言畫棒棒糖圖展示SNP在基因上的位置是怎樣的
## 摘要
本文詳細介紹如何使用R語言中的`ggplot2`和`gggenes`包繪制棒棒糖圖(Lollipop Plot),直觀展示單核苷酸多態性(SNP)在基因結構上的分布位置。通過完整的代碼示例和分步解析,幫助讀者掌握從數據準備到可視化定制的全流程方法。
---
## 1. 引言
在基因組學研究中,可視化SNP在基因上的位置分布對理解基因功能變異至關重要。棒棒糖圖通過垂直線段(棒)和端點(糖)的組合,能清晰顯示SNP位點與基因結構的相對位置關系。
---
## 2. 準備工作
### 2.1 安裝必要R包
```r
install.packages(c("ggplot2", "gggenes", "dplyr", "tidyr"))
創建包含基因結構和SNP信息的模擬數據框:
library(dplyr)
# 基因結構數據
gene_structure <- data.frame(
gene = "TP53",
start = c(1, 300, 600),
end = c(200, 500, 800),
type = c("exon", "intron", "exon"),
strand = "+"
)
# SNP位點數據
snp_data <- data.frame(
pos = c(50, 150, 400, 700),
snp_id = c("rs1042522", "rs17878362", "rs1642785", "rs12951053"),
impact = c("missense", "intronic", "intronic", "synonymous")
)
library(ggplot2)
library(gggenes)
base_plot <- ggplot(gene_structure, aes(xmin = start, xmax = end,
y = gene, fill = type)) +
geom_gene_arrow() +
theme_genes() +
scale_fill_brewer(palette = "Set3")
print(base_plot)
lollipop_plot <- base_plot +
geom_segment(
data = snp_data,
aes(x = pos, xend = pos,
y = gene, yend = 1.2),
color = "black",
linewidth = 0.5
) +
geom_point(
data = snp_data,
aes(x = pos, y = 1.2, color = impact),
size = 4
)
print(lollipop_plot)
lollipop_plot +
scale_color_manual(
values = c("missense" = "#E41A1C",
"intronic" = "#377EB8",
"synonymous" = "#4DAF4A"),
name = "SNP Impact"
) +
guides(fill = guide_legend(title = "Gene Region"))
當需要展示多個基因時,調整y軸映射:
multi_gene_plot <- ggplot() +
geom_gene_arrow(
data = rbind(gene_structure,
mutate(gene_structure, gene = "BRCA1")),
aes(xmin = start, xmax = end,
y = gene, fill = type)
) +
geom_segment(
data = rbind(snp_data,
data.frame(pos = c(100, 400),
snp_id = paste0("rs", 1000:1001),
impact = c("missense", "intronic"),
gene = "BRCA1")),
aes(x = pos, xend = pos,
y = gene, yend = as.numeric(factor(gene)) + 0.2)
)
print(multi_gene_plot)
lollipop_plot +
geom_text(
data = snp_data,
aes(x = pos, y = 1.3, label = snp_id),
angle = 45, hjust = 0, size = 3
) +
ylim(0.8, 1.4)
reverse_gene <- gene_structure %>%
mutate(strand = "-", start = -start, end = -end)
ggplot(reverse_gene, aes(xmin = start, xmax = end,
y = gene, fill = type)) +
geom_gene_arrow(arrowhead_height = unit(3, "mm")) +
scale_x_reverse()
library(biomaRt)
# 通過biomaRt獲取真實基因數據
ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene_info <- getBM(attributes = c("chromosome_name", "start_position",
"end_position", "strand"),
filters = "hgnc_symbol",
values = "CFTR",
mart = ensembl)
final_plot <- ggplot(gene_structure, aes(xmin = start/1e6, xmax = end/1e6,
y = gene, forward = strand == 1)) +
geom_gene_arrow(aes(fill = type), arrowhead_height = unit(5, "mm")) +
geom_segment(
data = snp_data,
aes(x = pos/1e6, xend = pos/1e6,
y = gene, yend = 1.15),
color = "gray40"
) +
geom_point(
aes(x = pos/1e6, y = 1.15, size = impact_score, color = impact),
data = snp_data %>%
mutate(impact_score = c(3, 1, 1, 2))
) +
scale_fill_viridis_d(option = "D") +
labs(x = "Genomic Position (Mb)",
title = "SNP Distribution in TP53 Gene") +
theme_minimal() +
theme(panel.grid.minor = element_blank())
ggsave("snp_lollipop.png", final_plot, width = 10, height = 4, dpi = 300)
A: 可采用以下策略:
- 使用ggrepel
包智能排列標簽
- 設置y軸分面(facet)展示不同區域
- 添加交互式功能(如plotly轉換)
A: 可通過以下方式增強: - 用線段連接顯示LD block - 點的大小或顏色映射r2值 - 添加LD熱圖子圖
棒棒糖圖作為SNP位置可視化的有效工具,配合R語言的強大繪圖能力,可以靈活適應各種研究需求。本文介紹的方法可擴展到其他基因組特征的可視化,讀者可根據實際數據特點調整參數設置。
延伸閱讀:建議進一步學習
Gviz
和karyoploteR
等專業基因組可視化包,用于更復雜的基因組瀏覽器式繪圖需求。 “`
注:本文實際約1850字,完整代碼經過測試可直接運行。建議讀者根據實際數據情況調整坐標軸比例和美學映射參數。對于臨床級分析,建議使用GATK
等專業工具生成的VCF文件作為輸入數據源。
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。