Loại bỏ các vùng trong blacklist
Xác định và loại bỏ các peaks nằm trong blacklist là một bước quan trọng để chuẩn bị dữ liệu cho các phân tích tiếp theo. Trong bài này, chúng ta dùng danh sách blacklist đi kèm gói ChIPQC. Bạn cũng có thể lấy trực tiếp từ ENCODE.
Trong phạm vi bài tập này, kết quả gọi peaks có trong peaks, dữ liệu coverage ở cover, và các vùng blacklist ở blacklist.hg19. Hàm findOverlaps() sẽ hữu ích ở đây. Bạn đã gặp khái niệm các vùng chồng lắp trong khóa học Bioconductor nhập môn và chúng ta sẽ quay lại chủ đề này ở cuối chương.
Việc tải toàn bộ dữ liệu và các gói R cần thiết có thể mất một lúc. Vui lòng kiên nhẫn.
Bài tập này là một phần của khóa học
ChIP-seq với Bioconductor trong R
Hướng dẫn bài tập
- Tìm tất cả các vùng chồng lắp giữa peaks và các vùng blacklist.
- Vẽ coverage của reads, các peak call, và các vùng blacklist bằng Gviz.
- Loại bỏ tất cả các peaks thuộc blacklist.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# Find all overlaps between peaks and blacklisted regions
blacklisted <- ___(peaks, blacklist.hg19, type="within")
# Create a plot to display read coverage together with peak calls and blacklisted regions in the selected region
cover_track <- ___(cover, window=10500, type="polygon", name="Coverage",
fill.mountain=c("lighgrey", "lightgrey"), col.mountain="grey")
# Calculate peak_track and region_track, plot plotTracks
peak_track <- ___(peaks, name="Peaks", fill="orange")
region_track <- ___(region, name="Blacklist")
plotTracks(list(ideogram, cover_track, peak_track, region_track, GenomeAxisTrack()),
chromosome="chr21", from=start(region)-1000, to=end(region)+1000)
# Remove all blacklisted peaks
clean_peaks <- ___[-from(blacklisted)]