On Thursday 12 October 2006 23:34, Lisa Luo wrote:
> Dear List,
>
> I have a set of BAC array CGH data to analysis. I am new to this
kind of
> analysis. Could anyone please recommend me a good package?
>
> I read some papers on aCGH. For the size of BAC clones, is
averaging a few
> probes a good idea? I tried DNAcopy and I am not happy with the
results
> either. Another question: Can we use the ratio to call it
amplification
> or deletion? In my normal samples, I should expect ratio to be
about 1.
> But ratios in some probes/chromosomes tend to be high and ratios in
another
> probes/regions tend to be lower. So how to use the ratio?
Lisa,
Segmentation methods (GLAD, DNAcopy, aCGH, etc.) are the way to go
here, as
others have suggested. Averaging is not a good idea given the
availability
of segmentation methods, which accomplish the same goal of reducing
noise as
does averaging, but do so with minimal loss of resolution.
The ratio is the number of interest, yes. Your normal samples should
have a
ratio of approximately 1, yes. However, because of technical
variation, some
probes and regions may drift a bit away from 1. If you have these
normal/normal hybs, you can use these to help you to determine where
your
data (in terms of ratio) becomes believable. Another useful measure
is to
look at the ratios that you see on the X-chromosome when comparing two
samples with different genders in comparison to the ratio that you see
for
the autosomes.
Finally, were there DNA quality issues that could explain some
variation in
your results? Were the samples amplified or directly labeled?
Sean
On Friday 13 October 2006 10:04, Lisa Luo wrote:
> Thank you, Sean and all others who replied to me.
>
> First for segmentation method, I am not satisfied with DNAcopy. A
> extreme in this case is that when I have a break point near the
middle, it
> could not detect it if I used the SD trim method. Or it gave too
many
> breaking points.
You may have to change parameters multiple times to get the results
that you
like. Also, you might consider using a merging function like that
implemented in snapCGH rather than using trimming, but others may not
agree
here.
> Also, it can not detect full chromosome change that you
> have to rely on the ratio.
The goal of the segmentation methods is not to call gains and losses,
but to
determine regions of constant copy number. There is another step in
the
process to determine whether a copy number estimate actually
represents a
gain or loss.
> Is BAC array CGH single strand PCR arrays? If so, should repeat
affact
> the result? In my tumor samples with opposite sex for control, I
found
> that the ratio for X varied from sample to sample. In one sample,
it is
> even 1. The explaination I got is that chrX has more repeat
regions. Are
> you suggesting that DNA quality may play a role in this?
Of course quality plays a role in the array result--both quality of
the arrays
and quality of the DNA. If you have arrays on which you have a male
and a
female hybed and the mean of the X-chromosome is the same of the mean
of the
autosomes, there are a few explanations:
1) Sample mixup--the samples are actually the same sex
2) One of the samples has either gained/lost an X chromosome (cell
line
artifact, if the samples are from cell lines)
3) The processing/hybing of the array failed for some reason
4) Probes are somehow mismapped, either by you or by the person
supplying the
annotation.
5) There actually IS a difference between the X-chromosome and the
autosomes,
but the difference is small or difficult to notice given the noise;
for
practical purposes, the X-chromosome/autosome ratio represents signal
and
noise can be quantified in numerous ways (SD, MAD, etc). Arrays with
low
signal-noise-ratio for whatever reason will look like there is very
little
difference between the autosomes/X-chromosome.
> I am planning to look at each BAC at one time, comparing tumor vs
normal.
> Is this a good idea?
This sounds fine, but I would use the data calculated from the
segmentation
results as the input to your test; this will reduce the noise
associated with
individual probes. Also, note that the data are not likely to be
normally
distributed at a given probe across samples, since each sample is
potentially
drawn from a different population (mean copy number).
Sean
There are tho main packages: GLAD (plus MANOR for the qualit? test)
and
aCGH
Best regards
Dr Alberto Goldoni
Medical Genetics Unit
S. Orsola-Malpighi Hospital
Via Massarenti n.9, Pad 11
40100 Bologna, Italy
Mobile Phone: +39-338-4145970
Fax: +39-051-636-4004
alberto.goldoni at eurogene.org
www.eurogene.org
www.lagem.it
-----Messaggio originale-----
Da: bioconductor-bounces at stat.math.ethz.ch
[mailto:bioconductor-bounces at stat.math.ethz.ch] Per conto di Lisa
Luo
Inviato: venerd? 13 ottobre 2006 5.34
A: bioconductor
Oggetto: [BioC] array CGH
Dear List,
I have a set of BAC array CGH data to analysis. I am new to this kind
of analysis. Could anyone please recommend me a good package?
I read some papers on aCGH. For the size of BAC clones, is averaging
a
few probes a good idea? I tried DNAcopy and I am not happy with the
results either. Another question: Can we use the ratio to call it
amplification or deletion? In my normal samples, I should expect
ratio
to be about 1. But ratios in some probes/chromosomes tend to be high
and ratios in another probes/regions tend to be lower. So how to use
the ratio?
Thank you so much for your help!
Lisa
---------------------------------
[[alternative HTML version deleted]]
_______________________________________________
Bioconductor mailing list
Bioconductor at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
Hi Lisa,
The snapCGH package might be a good place to start because it allows
you to
try out methods from the most common CGH packages (GLAD, aCGH, DNAcopy
etc..) and displays the results in a common format, allowing for easy
comparison.
Regards,
Mark
On Oct 13 2006, Lisa Luo wrote:
>Dear List,
>
> I have a set of BAC array CGH data to analysis. I am new to this
kind of
> analysis. Could anyone please recommend me a good package?
>
> I read some papers on aCGH. For the size of BAC clones, is averaging
a
> few probes a good idea? I tried DNAcopy and I am not happy with the
> results either. Another question: Can we use the ratio to call it
> amplification or deletion? In my normal samples, I should expect
ratio to
> be about 1. But ratios in some probes/chromosomes tend to be high
and
> ratios in another probes/regions tend to be lower. So how to use the
> ratio?
>
>Thank you so much for your help!
>
>Lisa
>
>
>---------------------------------
>
> [[alternative HTML version deleted]]
>
> _______________________________________________ Bioconductor mailing
list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor Search the
archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
--
Mark Dunning
PhD Student
Computational Biology Group
Hutchison / MRC Research Centre
Department of Oncology
University of Cambridge
Hills Rd, Cambridge CB2 2XZ
Email : md392 at cam.ac.uk
Phone : +44 (0) 1223 763380
Hi Lisa,
if after applying segmentation methods to your data the normal samples
still indicate gains or losses you might want to have a look at
http://projects.tcag.ca/variation/ for a list of known large scale
copy
number polymorphims.
Hilmar
Lisa Luo schrieb:
> Dear List,
>
> I have a set of BAC array CGH data to analysis. I am new to this
kind of analysis. Could anyone please recommend me a good package?
>
> I read some papers on aCGH. For the size of BAC clones, is
averaging a few probes a good idea? I tried DNAcopy and I am not
happy with the results either. Another question: Can we use the
ratio to call it amplification or deletion? In my normal samples, I
should expect ratio to be about 1. But ratios in some
probes/chromosomes tend to be high and ratios in another
probes/regions tend to be lower. So how to use the ratio?
>
> Thank you so much for your help!
>
> Lisa
>
>
> ---------------------------------
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
>
--
Hilmar Berger
Studienkoordinator
Institut f?r medizinische Informatik, Statistik und Epidemiologie
Universit?t Leipzig
H?rtelstr. 16-18
D-04107 Leipzig
Tel. +49 341 97 16 101
Fax. +49 341 97 16 109
email: hilmar.berger at imise.uni-leipzig.de