Probe design methods and microarrays for comparative genomic hybridization and location analysis

• US 8,036,835 B2
• Filed: 06/09/2010
• Issued: 10/11/2011
• Est. Priority Date: 11/23/2004
• Status: Active Grant

• Current Assignee: Agilent Technologies, Inc.
• Original Assignee: Agilent Technologies, Inc.
• Inventors: Sampas, Nicholas M., Lipson, Doron, Yakhini, Zohar H., Curry, Bo, Tsang, Peter
• Primary Examiner(s): Zeman; Mary
• Application Number:
  US12/797,521
• Publication Number:
  US 20100279883A1
• Time in Patent Office: 489 Days
• Field of Search: None
• US Class Current: 702/19
• CPC Class Codes: B01J 2219/00689 using computers
• CPC Class Codes: B01J 2219/00695 Synthesis control routines, e.g. using computer programs
• CPC Class Codes: B01J 2219/00722 Nucleotides
• CPC Class Codes: C12Q 1/6883 for diseases caused by alterations of genetic material
• CPC Class Codes: G16B 25/00 ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
• CPC Class Codes: G16B 25/20 Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation

Abstract

Methods and systems for identifying and selecting nucleic acid probes for detecting a target with a nucleic acid probe array or comparative genome hybridization microarray, comprising selecting a plurality of potential target sequences, generating a plurality of candidate probes from the target sequences, filtering the plurality of candidate probes by analyzing candidate probes for selected probe properties in silico. Microarrays comprising probes selected by the methods of the invention are particularly useful for comparative genome hybridization and location analysis.

Claims

1. A computer-implemented method for generating a set of probe nucleic acid sequences, the method comprising:

(a) sorting a plurality of candidate probe nucleic acid sequences, for a genomic region of interest, from smallest genomic distance to largest genomic distance between neighboring candidate probe nucleic acid sequences to produce a sorted plurality of candidate probe nucleic acid sequences;

(b) evaluating a probe parameter for a neighboring pair of candidate probe nucleic acid sequences from said sorted plurality to identify a first member of said neighboring pair with a more desirable probe parameter than a second pair member of said neighboring pair;

(c) removing said second pair member from said plurality;

(d) reiterating said sorting, evaluating and removing steps at least once to generate a set of probe nucleic acid sequences; and

(e) outputting said set of probe nucleic acid sequences

wherein said method is performed by a computer that is programmed to perform said method.

2. The method according to claim 1, wherein said neighboring pair evaluated in step (b) is a pair that is closest to each other in terms of genomic distance in said sorted plurality.

3. The method according to claim 2, wherein said probe parameter is an in silico probe parameter.

4. The method of claim 3, wherein said in silico probe parameter is selected from a group consisting of duplex melting temperature, hairpin stability, GC content, probe is within an exon, probe is within a gene, probe is within an intron, probe is within a intergenic region, and homology score.

5. The method according to claim 1, comprising synthesizing at least one probe nucleic acid having a sequence of a member probe nucleic acid sequence of said set of probe nucleic acid sequences.

6. The method according to claim 5, wherein said method further comprises assaying said probe nucleic acid in a hybridization assay.

7. The method according to claim 1, comprising generating a nucleic acid array comprising at least one probe nucleic acid having a sequence of a member probe nucleic acid sequence of said set of probe nucleic acid sequences.

8. The method according to claim 7, wherein said method further comprises contacting said nucleic acid array with a genomic sample.

9. The method according to claim 4, wherein said evaluating comprises evaluating duplex melting temperature and wherein said first member has a lower duplex melting temperature than said second pair member.

10. The method according to claim 4, wherein said evaluating comprises evaluating homology score and wherein said first member has a higher homology score than said second pair member.

11. The method according to claim 1, wherein said plurality of candidate probe nucleic acid sequences are generated by:

(i) selecting target sequences from said genomic region of interest;

(ii) repeat-masking said target sequences to form non-repeat masked regions;

(iii) tiling sequences across said non-repeat masked regions to generate said candidate nucleic acid sequences; and

(iv) screening said candidate probes nucleic acid sequences according to at least one in silico parameter.

12. The method according to claim 11, comprising identifying restriction enzyme sites in the genomic region of interest, and selecting target sequences that exclude said restriction enzyme sites.

13. A non-transitory computer readable medium comprising instructions for performing the following method:

(a) sorting a plurality of candidate probe nucleic acid sequences, for a genomic region of interest, from smallest genomic distance to largest genomic distance between neighboring candidate probe nucleic acid sequences to produce a sorted plurality of candidate probe nucleic acid sequences;

(b) evaluating a probe parameter for a neighboring pair of candidate probe nucleic acid sequences from said sorted plurality to identify a first member of said neighboring pair with a more desirable probe parameter than a second pair member of said neighboring pair;

(c) removing said second pair member from said plurality;

(d) reiterating said sorting, evaluating and removing steps at least once to generate a set of probe nucleic acid sequences; and

(e) outputting said set of probe nucleic acid sequences.

14. The non-transitory computer readable medium of claim 13, wherein said neighboring pair evaluated in step (b) is a pair that is closest to each other in terms of genomic distance in said sorted plurality.

15. The non-transitory computer readable medium of claim 14, wherein said probe parameter is an in silico probe parameter.

16. The non-transitory computer readable medium of claim 15, wherein said in silico probe parameter is selected from a group consisting of duplex melting temperature, hairpin stability, GC content, probe is within an exon, probe is within a gene, probe is within an intron, probe is within a intergenic region, and homology score.

17. The non-transitory computer readable medium of claim 16, wherein said evaluating comprises evaluating duplex melting temperature and wherein said first member has a lower duplex melting temperature than said second pair member.

18. The non-transitory computer readable medium of claim 16, wherein said evaluating comprises evaluating homology score and wherein said first member has a higher homology score than said second pair member.

19. The non-transitory computer readable medium of claim 13, where said plurality of candidate probe nucleic acid sequences are generated by:

(i) selecting target sequences from said genomic region of interest;

(ii) repeat-masking said target sequences to form non-repeat masked regions;

(iii) tiling sequences across said non-repeat masked regions to generate said candidate nucleic acid sequences; and

(iv) screening said candidate probes nucleic acid sequences according to at least one in silico parameter.

20. The non-transitory computer readable medium of claim 19, wherein said method comprises identifying restriction enzyme sites in the genomic region of interest, and selecting target sequences that exclude said restriction enzyme sites.