937 |
20 Jan 09 |
jari |
1 |
---------------------------------------------------------------------- |
937 |
20 Jan 09 |
jari |
2 |
{{{ |
937 |
20 Jan 09 |
jari |
$Id$ |
937 |
20 Jan 09 |
jari |
4 |
|
937 |
20 Jan 09 |
jari |
Copyright (C) 2009 Markus Ringnér |
937 |
20 Jan 09 |
jari |
6 |
|
937 |
20 Jan 09 |
jari |
This file is part of Illumina plug-in package for BASE. |
937 |
20 Jan 09 |
jari |
Available at http://baseplugins.thep.lu.se/ |
937 |
20 Jan 09 |
jari |
BASE main site: http://base.thep.lu.se/ |
937 |
20 Jan 09 |
jari |
10 |
|
937 |
20 Jan 09 |
jari |
This is free software; you can redistribute it and/or |
937 |
20 Jan 09 |
jari |
modify it under the terms of the GNU General Public License |
940 |
27 Jan 09 |
martin |
as published by the Free Software Foundation; either version 3 |
937 |
20 Jan 09 |
jari |
of the License, or (at your option) any later version. |
937 |
20 Jan 09 |
jari |
15 |
|
937 |
20 Jan 09 |
jari |
The software is distributed in the hope that it will be useful, |
937 |
20 Jan 09 |
jari |
but WITHOUT ANY WARRANTY; without even the implied warranty of |
937 |
20 Jan 09 |
jari |
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
937 |
20 Jan 09 |
jari |
GNU General Public License for more details. |
937 |
20 Jan 09 |
jari |
20 |
|
937 |
20 Jan 09 |
jari |
You should have received a copy of the GNU General Public License |
940 |
27 Jan 09 |
martin |
along with this package. If not, see <http://www.gnu.org/licenses/>. |
937 |
20 Jan 09 |
jari |
23 |
}}} |
937 |
20 Jan 09 |
jari |
24 |
---------------------------------------------------------------------- |
937 |
20 Jan 09 |
jari |
25 |
|
937 |
20 Jan 09 |
jari |
= Importing Illumina data from !BeadStudio into BASE = |
937 |
20 Jan 09 |
jari |
27 |
|
937 |
20 Jan 09 |
jari |
This file contains information about how to generate expression |
937 |
20 Jan 09 |
jari |
data for the Illumina plug-in package from !BeadStudio. |
937 |
20 Jan 09 |
jari |
30 |
|
937 |
20 Jan 09 |
jari |
31 |
|
937 |
20 Jan 09 |
jari |
== Background == |
937 |
20 Jan 09 |
jari |
33 |
|
937 |
20 Jan 09 |
jari |
For expression data, the Illumina plug-ins package for BASE requires |
937 |
20 Jan 09 |
jari |
Illumina Bead Summary (IBS) files. IBS files contain bead-type level |
937 |
20 Jan 09 |
jari |
data for scanned Sentrix arrays. IBS files are simple comma separated |
937 |
20 Jan 09 |
jari |
text files with file extension .csv. IBS files are not outputted by a |
937 |
20 Jan 09 |
jari |
!BeadArray Reader with default settings, but a local Illumina Field |
937 |
20 Jan 09 |
jari |
Application Scientist can configure the !BeadArray Reader to output |
937 |
20 Jan 09 |
jari |
IBS files. This document describes how to generate IBS files from |
937 |
20 Jan 09 |
jari |
!BeadStudio without requiring any changes of the default settings for |
937 |
20 Jan 09 |
jari |
the !BeadArray Reader. |
937 |
20 Jan 09 |
jari |
43 |
|
937 |
20 Jan 09 |
jari |
44 |
|
937 |
20 Jan 09 |
jari |
== !BeadStudio == |
937 |
20 Jan 09 |
jari |
46 |
|
937 |
20 Jan 09 |
jari |
''Create !BeadStudio project'' |
937 |
20 Jan 09 |
jari |
48 |
|
937 |
20 Jan 09 |
jari |
Generate a !BeadStudio gene expression project with the experiments |
937 |
20 Jan 09 |
jari |
you want to generate IBS files for. The IBS files should contain data |
937 |
20 Jan 09 |
jari |
that is neither normalized nor background subtracted so: |
937 |
20 Jan 09 |
jari |
52 |
|
937 |
20 Jan 09 |
jari |
* ''Analysis Type'' should be set to ''Gene Expression''. |
937 |
20 Jan 09 |
jari |
54 |
|
937 |
20 Jan 09 |
jari |
* Set the parameter ''Normalization'' to ''none''. |
937 |
20 Jan 09 |
jari |
56 |
|
937 |
20 Jan 09 |
jari |
* The parameter button ''Subtract Background'' should be unchecked. |
937 |
20 Jan 09 |
jari |
58 |
|
937 |
20 Jan 09 |
jari |
''Create and export !BeadStudio tables'' |
937 |
20 Jan 09 |
jari |
60 |
|
937 |
20 Jan 09 |
jari |
Two tables should be exported: ''Sample Probe Profile'' and ''Control |
937 |
20 Jan 09 |
jari |
Probe Profile''. |
937 |
20 Jan 09 |
jari |
63 |
|
937 |
20 Jan 09 |
jari |
Note that the tables must be exactly as specified here: no additional |
937 |
20 Jan 09 |
jari |
columns are allowed and the order of the columns needs to be as |
937 |
20 Jan 09 |
jari |
specified. |
937 |
20 Jan 09 |
jari |
67 |
|
937 |
20 Jan 09 |
jari |
''Sample Probe Profile'' |
937 |
20 Jan 09 |
jari |
69 |
|
937 |
20 Jan 09 |
jari |
1. Click on the ''Sample Probe Profile'' tab to view the table. |
937 |
20 Jan 09 |
jari |
71 |
|
937 |
20 Jan 09 |
jari |
2. Click on the ''Column Chooser'' button to open the column chooser window. |
937 |
20 Jan 09 |
jari |
73 |
|
937 |
20 Jan 09 |
jari |
3. Select displayed columns. The first column should be |
937 |
20 Jan 09 |
jari |
''PROBE_ID'', followed by one column per sample, for example, |
937 |
20 Jan 09 |
jari |
4313143087_A, 4313143087_B, ...., and the last column should be |
937 |
20 Jan 09 |
jari |
''ProbeID''. |
937 |
20 Jan 09 |
jari |
78 |
|
937 |
20 Jan 09 |
jari |
4. Select displayed subcolumns. The displayed subcolumns should be |
937 |
20 Jan 09 |
jari |
''AVG_Signal'', ''BEAD_STDERR'', and ''Avg_NBEADS'', and in |
937 |
20 Jan 09 |
jari |
this order. |
937 |
20 Jan 09 |
jari |
82 |
|
937 |
20 Jan 09 |
jari |
5. Export the table by clicking on the ''Export displayed data to a |
937 |
20 Jan 09 |
jari |
file '' button. Save as type ''Tab delimited (*.txt)''. In the |
937 |
20 Jan 09 |
jari |
following it is assumed that you have named the file |
937 |
20 Jan 09 |
jari |
''sample_probe_profile.txt''. |
937 |
20 Jan 09 |
jari |
87 |
|
937 |
20 Jan 09 |
jari |
''Control Probe Profile'' |
937 |
20 Jan 09 |
jari |
89 |
|
937 |
20 Jan 09 |
jari |
1. Click on the ''Control Probe Profile'' tab to view the table. If |
937 |
20 Jan 09 |
jari |
there is no ''Control Probe Profile'' tab, then first enable it |
937 |
20 Jan 09 |
jari |
under the pull-down menu ''Window''. |
937 |
20 Jan 09 |
jari |
93 |
|
937 |
20 Jan 09 |
jari |
2. Click on the ''Column Chooser'' button to open the column chooser window. |
937 |
20 Jan 09 |
jari |
95 |
|
937 |
20 Jan 09 |
jari |
3. Select displayed columns. The first displayed column should be |
937 |
20 Jan 09 |
jari |
''ProbeID'', followed by one column per sample, for example, |
937 |
20 Jan 09 |
jari |
4313143087_A, 4313143087_B, .... The selection of samples should be |
937 |
20 Jan 09 |
jari |
the same and in the same order as in the exported ''Sample Probe |
937 |
20 Jan 09 |
jari |
Profile''. |
937 |
20 Jan 09 |
jari |
101 |
|
937 |
20 Jan 09 |
jari |
4. Select displayed subcolumns. The shown subcolumns should be |
937 |
20 Jan 09 |
jari |
''AVG_Signal'', ''BEAD_STDERR'', and ''Avg_NBEADS'', and in this |
937 |
20 Jan 09 |
jari |
order. |
937 |
20 Jan 09 |
jari |
105 |
|
937 |
20 Jan 09 |
jari |
5. Export the table by clicking on the ''Export displayed data to a |
937 |
20 Jan 09 |
jari |
file '' button. Save as type ''Tab delimited (*.txt)''. In the |
937 |
20 Jan 09 |
jari |
following it is assumed that you have named the file |
937 |
20 Jan 09 |
jari |
''control_probe_profile.txt''. |
937 |
20 Jan 09 |
jari |
110 |
|
937 |
20 Jan 09 |
jari |
111 |
|
937 |
20 Jan 09 |
jari |
== Convert !BeadStudio tables to IBS files == |
937 |
20 Jan 09 |
jari |
113 |
|
937 |
20 Jan 09 |
jari |
The two files exported from !BeadStudio with data for the sample probes |
937 |
20 Jan 09 |
jari |
and the control probes can be converted into an IBS file for each |
937 |
20 Jan 09 |
jari |
sample using ''beadstudio2ibs.pl'' with the following command: |
937 |
20 Jan 09 |
jari |
117 |
{{{ |
937 |
20 Jan 09 |
jari |
perl beadstudio2ibs.pl --samples_file=sample_probe_profile.txt --controls_file=control_probe_profile.txt |
937 |
20 Jan 09 |
jari |
119 |
}}} |
937 |
20 Jan 09 |
jari |
120 |
|
937 |
20 Jan 09 |
jari |
== Import IBS files generated from !BeadStudio into BASE == |
937 |
20 Jan 09 |
jari |
122 |
|
937 |
20 Jan 09 |
jari |
The IBS files generated by ''beadstudio2ibs.pl'' can be imported into |
937 |
20 Jan 09 |
jari |
BASE using the following settings: |
937 |
20 Jan 09 |
jari |
125 |
|
937 |
20 Jan 09 |
jari |
* Use the raw data type for the Illumina plug-in package: |
937 |
20 Jan 09 |
jari |
''Illumina Bead Summary (IBS)'' |
937 |
20 Jan 09 |
jari |
128 |
|
937 |
20 Jan 09 |
jari |
* Use the option ''expression 1''. |
937 |
20 Jan 09 |
jari |
130 |
|
937 |
20 Jan 09 |
jari |
131 |
|
937 |
20 Jan 09 |
jari |
== Additional remarks == |
937 |
20 Jan 09 |
jari |
133 |
|
937 |
20 Jan 09 |
jari |
* The choice between ''expression 1'' and ''expression 2'' (and |
937 |
20 Jan 09 |
jari |
possibly more alternatives in the future) should be made to match |
937 |
20 Jan 09 |
jari |
the number of IBS files per sample. When IBS files are outputted |
937 |
20 Jan 09 |
jari |
directly from the !BeadArray Reader one .csv file is produced per |
937 |
20 Jan 09 |
jari |
strip on the array. !BeadStudio does not keep track of strips and |
937 |
20 Jan 09 |
jari |
therefore a single .csv file per sample is generated by |
937 |
20 Jan 09 |
jari |
''beadstudio2ibs.pl'', regardless of the number of strips on the |
937 |
20 Jan 09 |
jari |
beadchips used. This means, for example, that for HumanWG-6 v3.0 |
937 |
20 Jan 09 |
jari |
beadchips (two strips per sample) ''expression 2'' should be used |
937 |
20 Jan 09 |
jari |
when IBS files are generated from the scanner and ''expression 1'' |
937 |
20 Jan 09 |
jari |
should be used when IBS file are generated from !BeadStudio data. |
937 |
20 Jan 09 |
jari |
145 |
|
937 |
20 Jan 09 |
jari |
* The data in an IBS file from !BeadStudio will not be identical to |
937 |
20 Jan 09 |
jari |
the data in an IBS file from the scanner. The reason for this |
937 |
20 Jan 09 |
jari |
difference is that the scanner outputs the average intensity for all |
937 |
20 Jan 09 |
jari |
beads, whereas !BeadStudio removes outlier beads in its calculation |
937 |
20 Jan 09 |
jari |
of averages. |