ยีนที่สนใจ

Bioconductor เบื้องต้นใน R

Paula Andrea Martinez, PhD.

Data Scientist

ตัวอย่างของ genomic intervals

Reads ที่ align กับ reference
ยีนที่สนใจ
บริเวณ Exon
Single nucleotide polymorphisms (SNPs)
บริเวณ transcription หรือ binding sites, RNA-seq หรือ ChIP-seq

Genomic Ranges

library(GenomicRanges)
(myGR <- GRanges("chr1:200-300"))

 GRanges object with 1 range and 0 metadata columns:
  seqnames     ranges strand
     <Rle>  <IRanges>  <Rle>
[1]     chr1 [200, 300]    *
 -----
seqinfo: 1 sequence from an unspecified genome; no seqlengths

คลาส GRanges ใช้เก็บ genomic intervals ตามโครโมโซม
อาร์กิวเมนต์ขั้นต่ำคือ chr1:200-300
GRanges มี seqnames และ seqinfo

# df a data.frame like structure
   seqnames start end strand score  GC
1     chrX    50 120      +     1 0.25 
2     chrX   130 140      +     2 0.25 
3     chrX   153 154      +     3 0.25 
4     chrY    30  40      *     4 0.25 
5     chrY    50  55      -     5 0.25

(myGR <- as(df, "GRanges")) # transform df into GRanges

GRanges object with 5 ranges and 2 metadata columns:
      seqnames     ranges strand |     score        GC     
         <Rle>  <IRanges>  <Rle> | <integer> <numeric>
  [1]     chrX [ 50, 120]      + |         1      0.25   
  [2]     chrX [130, 140]      + |         2      0.25   
  [3]     chrX [153, 154]      + |         3      0.25   
  [4]     chrY [ 30,  40]      * |         4      0.25   
  [5]     chrY [ 50,  55]      - |         5      0.25   
 -----
 seqinfo: 2 sequences from an unspecified genome; no seqlengths

Accessors ของ Genomic Ranges

methods(class = "GRanges") # to check available accessors

# used for chromosome names 
seqnames(gr)

# returns an IRanges object for ranges
ranges(gr)

# stores metadata columns
mcols(gr)

# generic function to store sequence information
seqinfo(gr)

# stores the genome name
genome(gr)

Accessors ทำหน้าที่ได้ทั้ง setter และ getter
Accessors รับช่วงต่อได้ผ่านนิยาม S4

ยีนที่สนใจ: ABCD1

ABCD1 อยู่บริเวณปลายแขนยาวของโครโมโซม X
เข้ารหัสโปรตีนที่มีความสำคัญต่อการทำงานของเซลล์สมองและปอดในสัตว์เลี้ยงลูกด้วยนม
chrX มีความยาว ~156 ล้าน bp
ตำแหน่งบน chrX ~153.70 ล้าน bp

https://www.ncbi.nlm.nih.gov/gene/215

ChrX-ABCD1

GRanges ของโครโมโซม X

library(TxDb.Hsapiens.UCSC.hg38.knownGene)
hg <- TxDb.Hsapiens.UCSC.hg38.knownGene

เลือกยีนจากโครโมโซม X

hg_chrXg <- genes(hg, filter = list(tx_chrom = c("chrX")))

GRanges object with 1192 ranges and 1 metadata column:
            seqnames              ranges strand |     gene_id
               <Rle>           <IRanges>  <Rle> | <character>
  100008586     chrX   49551278-49568218      + |   100008586
      10009     chrX 120250752-120258398      + |       10009
  100093698     chrX   13310652-13319933      + |   100093698
        ...      ...                 ...    ... .         ...
  -------
  seqinfo: 640 sequences (1 circular) from hg38 genome

มาฝึกค้นหายีนที่สนใจในจีโนมมนุษย์กัน!

Bioconductor เบื้องต้นใน R