Question

Failed to annotated my wells

0

Entering edit mode

wstoo1 • 0

@wstoo1-8468

Last seen 9.3 years ago

European Union

Hi,

I would like to ask everyone some questions.

I am doing a genome wide screening of siRNA.

At the end of data analysis, I manage to generate the correct data but fail to annotate the wells.

Even though the state is TRUE for annotation but they are not in html report.

State:
configured=TRUE, annotated=TRUE

x <- readPlateList("platelist.txt", name = expName, path=dataPath)
x<-configure(x,confFile = confFile, logFile = logFile,descripFile = desFile, path = dataPath)
x<-annotate(x,geneIDFile="annotationST1.txt",path=dataPath)

xn<-summarizeChannels(x,fun = function (r1,r2,r3,r4,r5,r6,r7) r5/r3)

   xnp <- normalizePlates(xn, scale = "multiplicative", log = TRUE, method = "negatives", varianceAdjust = "none")
   xsc <- scoreReplicates(xnp, sign = "+", method = "zscore")
   xf <- summarizeReplicates(xsc, summary = "mean")

I tried to annotate the x (raw) and xf (summarized) one, but they generated the same results.

Plate   Well   GeneID
1   A01   DMPK
1   A02   empty
1   A03   TTK
1   A04   PIK3CA
1   A05   SYK
1   A06   KIAA1804
1   A07   EPHA8
1   A08   NAGK
1   A09   empty
1   A10   IRAK1
total 384.

Please someone help me with this.

Thank you.

cellhts2 • 2.0k views

ADD COMMENT • link updated 9.3 years ago by Joseph Barry ▴ 160 • written 9.3 years ago by wstoo1 • 0

score 0 · Answer 1 · 2015-07-24

Hi, Can you print the feature data before and after the annotate step to clarify whether or not the annotate function did its job? e.g. the equivalent to > x <- readPlateList("Platelist.txt", path=datadir, name="KcViabSmall") KcViabSmall: found data in 16 x 24 (384 well) format. Read 6 plates. > head(fData(x)) plate well controlStatus 1 1 A01 unknown 2 1 A02 unknown 3 1 A03 unknown 4 1 A04 unknown 5 1 A05 unknown 6 1 A06 unknown > x <- annotate(x, "GeneIDs_Dm_HFAsubset_1.1.txt", path=datadir) > head(fData(x)) plate well controlStatus HFAID GeneID 1 1 A01 unknown <na> <na> 2 1 A02 unknown <na> <na> 3 1 A03 unknown HFA00274 CG11371 4 1 A04 unknown HFA00646 CG31671 5 1 A05 unknown HFA00307 CG11376 6 1 A06 unknown HFA00324 CG11723 If the annotate step worked you should see an extra column. What you expect to see in the report is potentially a separate issue. Joe On Fri, Jul 24, 2015 at 11:00 AM, wstoo1 [bioc] <noreply@bioconductor.org> wrote: > Activity on a post you are following on support.bioconductor.org > > User wstoo1 <https: support.bioconductor.org="" u="" 8468=""/> wrote Question: > Failed to annotated my wells <https: support.bioconductor.org="" p="" 70344=""/>: > > Hi, > > I would like to ask everyone some questions. > > I am doing a genome wide screening of siRNA. > > At the end of data analysis, I manage to generate the correct data but > fail to annotate the wells. > > Even though the state is TRUE for annotation but they are not in html > report. > > > State: > configured=TRUE, annotated=TRUE > > x <- readPlateList("platelist.txt", name = expName, path=dataPath) > x<-configure(x,confFile = confFile, logFile = logFile,descripFile = > desFile, path = dataPath) > x<-annotate(x,geneIDFile="annotationST1.txt",path=dataPath) > > xn<-summarizeChannels(x,fun = function (r1,r2,r3,r4,r5,r6,r7) r5/r3) > > > xnp <- normalizePlates(xn, scale = "multiplicative", log = TRUE, > method = "negatives", varianceAdjust = "none") > xsc <- scoreReplicates(xnp, sign = "+", method = "zscore") > xf <- summarizeReplicates(xsc, summary = "mean") > > I tried to annotate the x (raw) and xf (summarized) one, but they > generated the same results. > > Plate Well GeneID > 1 A01 DMPK > 1 A02 empty > 1 A03 TTK > 1 A04 PIK3CA > 1 A05 SYK > 1 A06 KIAA1804 > 1 A07 EPHA8 > 1 A08 NAGK > 1 A09 empty > 1 A10 IRAK1 > total 384. > > Please someone help me with this. > > Thank you. > > > > > > ------------------------------ > > Post tags: cellhts2 > > You may reply via email or visit Failed to annotated my wells >

score 0 · Answer 2 · 2015-07-25

Hi,

I tried it.

> head(fData(xsc))
plate well controlStatus
1     1 A01        sample
2     1 A02           neg
3     1 A03        sample
4     1 A04        sample
5     1 A05        sample
6     1 A06        sample
>
> xsc<-annotate(xsc,geneIDFile="annotationST1.txt",path=dataPath)
> head(fData(xsc))
plate well controlStatus   GeneID X
1     1 A01        sample     DMPK NA
2     1 A02           neg    empty NA
3     1 A03        sample      TTK NA
4     1 A04        sample   PIK3CA NA
5     1 A05        sample      SYK NA
6     1 A06        sample KIAA1804 NA

It seems like annotation function is working. But there is another column exist called "x".

At first, x was assigned as my raw data.

However, If i assigned xfw to scored, it wont work since xfw doesnt have the scored data.

With the same function, my html result doesnt show annotations to my wells,

score 0 · Answer 3 · 2015-07-26

0

Entering edit mode

Joseph Barry ▴ 160

@joseph-barry-5000

Last seen 7.9 years ago

Dana-Farber Cancer Institute, Boston, U…

Great, so the annotate function is working. Can you try moving the annotation step to earlier in your pipeline e.g. after configure. The annotation will get carried along. Try running

state(x)
state(xnp)
state(xf)

and so on to follow what objects have been properly annotated. Since xf comes later in the pipeline, running annotate on that will not cause the previously saved objects to also inherit the annotation. Please let me know how this affects your report.

ADD COMMENT • link 9.3 years ago Joseph Barry ▴ 160

0

Entering edit mode

Yes. I tried to annotate my raw data x but nothing change in the html report.

The annotation status at the end was true. The html report didnt show the annotation correctly.

It showed:

WellAnno finalWellAnno

sample sample

It seems like the annotation function is working normally, but failed to link the write report function.

ADD REPLY • link 9.3 years ago wstoo1 • 0

0

Entering edit mode

Okay, at this point it is hard for me to pinpoint the problem without seeing the data. Would you mind privately sending me example data? joseph.barry@... please get the rest from my user profile

ADD REPLY • link 9.3 years ago Joseph Barry ▴ 160

score 0 · Answer 4 · 2015-08-06

An update to the thread for completeness. User wstoo1 sent me an example dataset privately. We figured out that he was assuming that the 'wellAnno' column should contain gene identifiers, which is not a valid expectation. The 'wellAnno' column contains information about the control status of different wells, while gene identifiers are stored in a separate column, 'GeneID'. The package is annotating the screen properly and outputting screen annotation to the report correctly.