Methods for identifying DNA copy number changes

• US 20060035258A1
• Filed: 08/08/2005
• Published: 02/16/2006
• Est. Priority Date: 08/06/2004
• Status: Abandoned Application

• Current Assignee: Affymetrix, Inc. (Thermo Fisher Scientific Incorporated)
• Original Assignee: Affymetrix, Inc. (Thermo Fisher Scientific Incorporated)
• Inventors: Wu, Kai, Tadakamalla, Rajeswari P.
• Application Number:
  US11/199,742
• Publication Number:
  US 20060035258A1
• Time in Patent Office: Days
• US Class Current: 435/6
• CPC Class Codes: C12Q 1/6809 Methods for determination or identification of nucleic acids involving differential detection
• CPC Class Codes: C12Q 2521/301 Endonuclease
• CPC Class Codes: C12Q 2531/119 Strand displacement amplification [SDA]
• CPC Class Codes: C12Q 2565/501 being an array of oligonucleotides
• CPC Class Codes: G16B 25/00 ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
• CPC Class Codes: G16B 25/10 Gene or protein expression profiling; Expression-ratio estimation or normalisation
• CPC Class Codes: G16B 25/20 Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
• CPC Class Codes: Y02A 90/10 Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

Methods and computer software products for identifying changes in genomic DNA copy number are disclosed. Methods for identifying homozygous deletions and genetic amplifications are disclosed. Genomic DNA is amplified generically and amplified sample is hybridized to an expression array. The expression array comprises probes to regions of genes that are expressed. The probes are complementary to genomic sequences found in mRNAs. Signal intensity is correlated to copy number. The methods may be used to detect copy number changes in cancerous tissue compared to normal tissue. The methods may be used to diagnose cancer and other diseases associated with chromosomal anomalies.

Claims

1. A method of estimating the copy number of a gene in a genomic sample comprising:

amplifying the genomic sample in an amplification reaction comprising random primers and a strand displacing DNA polymerase to obtain an amplified sample;

fragmenting and labeling the amplified sample;

hybridizing the labeled fragments to an array of probes comprising a plurality of probes, wherein said plurality of probes comprises at least 100,000 different sequence probes wherein each probe is present on the array at a feature of known or determinable location and wherein each probe in the plurality is complementary to an expressed region of a gene;

analyzing a resulting hybridization pattern to obtain a hybridization intensity measurement for each of a plurality of features of the array;

comparing the hybridization intensity measurement for each feature to an expected hybridization intensity measurement for that feature; and

estimating the copy number of one or more genes based on the hybridization intensity measurement.

2. The method of claim 1 wherein the random primers comprise random hexamers.

3. The method of claim 1 wherein the strand displacing DNA polymerase is a phi29 DNA polymerase.

4. The method of claim 1 wherein the strand displacing DNA polymerase is a Bst DNA polymerase.

5. The method of claim 1 wherein the strand displacing DNA polymerase is REPLI-g DNA polymerase.

6. The method of claim 1 wherein amplification of the genomic sample is by fragmentation with a restriction enzyme, ligation of an adaptor to the fragments and polymerase chain reaction amplification using a single primer that is complementary to the adaptor.

7. The method of claim 1 wherein the probes are between 15 and 60 nucleotides.

8. The method of claim 7 wherein the probes are 25 nucleotides.

9. The method of claim 7 wherein the array further comprises a plurality of control probes.

10. The method of claim 1 wherein the array is an expression array comprising more than 1,000,000 different features wherein each of said features has a different oligonucleotide probe sequence.

11. The method of claim 1 wherein the array comprises a plurality of probes to detect each of more than 30,000 different human mRNA transcripts.

12. A method for estimating genomic copy number at a plurality of genomic regions in a genomic sample comprising:

amplifying the genomic sample to obtain an amplified genomic sample;

fragmenting the amplified genomic sample to obtain fragments;

labeling the fragments;

hybridizing the fragments to an expression array to generate a hybridization pattern, wherein the expression array comprises at least 10,000 probe sets;

analyzing the hybridization pattern to obtain a plurality of probe set signals, wherein a probe set signal is a normalized measurement of the hybridization signal for a probe set;

calculating a Z-score for each probe set using a mean and standard deviation calculated from a training data set;

mapping the chromosomal location of each probe set to obtain a plurality of mapped probe sets that map to a single chromosomal location;

calculating a Stouffer Z-score for each mapped probe set; and

identifying chromosomal regions of amplification or deletion based on Stouffer Z-score.

13. The method of claim 12 wherein the genomic sample is amplified in a reaction comprising random primers and a strand displacing polymerase.

14. The method of claim 13 wherein the strand displacing polymerase is selected from the group consisting of phi29 DNA polymerase and Bst DNA polymerase.

15. The method of claim 12 wherein the fragments are end labeled with biotin in a reaction comprising terminal transferase.

16. The method of claim 12 wherein the training data set is obtained by analyzing the probe set signal from at least 30 control genomic samples.

17. The method of claim 12 wherein the genomic samples included in the training data set each have a mean normalized probe set signal of about 250 and a standard deviation of between 280 and 450.

18. The method of claim 16 wherein the control genomic samples are normal samples.

19. The method of claim 12 wherein probe sets with a Stouffer Z-score above a selected threshold are identified as being complementary to an amplified genomic region.

20. The method of claim 12 wherein probe sets with a Stouffer Z-score below a selected lower threshold are identified as being complementary to a deleted genomic region.

22. The computer software product of claim 18 wherein Stouffer Z-scores are obtained for a plurality of the probes sets and wherein the Stouffer Z-scores are plotted against chromosomal location and output as a graphical display.

21. A computer software product for analyzing hybridization data for a genomic sample hybridized to an expression array, comprising a computer readable medium having computer-executable instructions for performing logic steps comprising:

inputting probe intensities from probes designed to interrogate for the presence of mRNA transcripts;

obtaining a normalized signal for a plurality of probe sets;

partitioning the data into a training set and a test set;

generating a signal mean and a standard deviation from the training set;

generating a Z-score for a plurality of probe sets;

identifying probe sets that map to chromosomal locations and comparing Z-scores for probe sets to estimate copy number for selected genomic regions.