Shoda a filtrování

Úvod do Bioconductoru v R

Paula Andrea Martinez, PhD.

Data Scientist

Duplicitní sekvence

Biologické duplikáty sekvencí se vyskytují přirozeně
Amplifikace při přípravě knihovny (PCR)
Vícenásobné sekvenování vzorku

Odstranit duplikáty nebo je označit

Sekvenování celého genomu nebo exomu

Označit duplikáty pomocí prahové hodnoty

RNA-seq a ChIP-seq

srduplicated

library(ShortRead)

# Counting duplicates TRUE is the number of duplicates
table(srduplicated(dfqsample))

FALSE  TRUE 
500   500

# Cleaning reads from duplicates x[fun(x)]
cleanReads <- mydReads[srduplicated(mydReads) == FALSE]

# Counting duplicates
table(srduplicated(cleanReads))

FALSE
500

Vytváření vlastních filtrů

srFilter pro filtrování dle podmínky x[fun(x)]

Příklad filtru

library(ShortRead)

# Use a custom filter to remove reads from fqsample
# This filter to remove reads shorter than a min number of bases
readWidthCutOff <- srFilter(function(x) {width(x) >= minWidth}, 
                            name = "MinWidth")

minWidth <- 51

fqsample[readWidthCutOff(fqsample)]

nFilter

library(ShortRead)

# save your filter, .name is optional
myFilter <- nFilter(threshold = 10, .name = "cleanNFilter")

# use the filter at reading point    
filtered <- readFastq(dirPath = "data", 
                      pattern = ".fastq",
                       filter = myFilter)

# you will retrieve only those reads that have a maximum of 10 N's
filtered

idFilter a polynFilter

library(ShortRead)

#id filter example
myFilterID <- idFilter(regex = ":3:1")
# will return only those ids that contain the regular expression

# optional parameters are .name, fixed and exclude
# use the filter at reading point    
filtered <- readFastq(dirPath = "data", pattern = ".fastq",
                       filter = myFilterID)

# filter to remove poly-A regions
myFilterPolyA <- polynFilter(threshold = 10, nuc = c("A"))
# will return the sequences that have a maximun number of 10 consecutive A's

# use the filter for subsetting    
filtered[myFilterPolyA(filtered)]

Pojďme procvičit filtry!

Úvod do Bioconductoru v R