US20060194224A1 - Computer-aided nucleic acid sequencing - Google Patents

Computer-aided nucleic acid sequencing Download PDF

Info

Publication number
US20060194224A1
US20060194224A1 US11/300,196 US30019605A US2006194224A1 US 20060194224 A1 US20060194224 A1 US 20060194224A1 US 30019605 A US30019605 A US 30019605A US 2006194224 A1 US2006194224 A1 US 2006194224A1
Authority
US
United States
Prior art keywords
nucleic acid
probes
probe
bases
acid probes
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US11/300,196
Inventor
Lubert Stryer
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Affymetrix Inc
Original Assignee
Affymetrix Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Affymetrix Inc filed Critical Affymetrix Inc
Priority to US11/300,196 priority Critical patent/US20060194224A1/en
Publication of US20060194224A1 publication Critical patent/US20060194224A1/en
Assigned to GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT reassignment GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT SECURITY AGREEMENT Assignors: AFFYMETRIX, INC.
Assigned to AFFYMETRIX, INC. reassignment AFFYMETRIX, INC. RELEASE BY SECURED PARTY (SEE DOCUMENT FOR DETAILS). Assignors: GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT
Abandoned legal-status Critical Current

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • C12Q1/6874Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/10Sequence alignment; Homology search

Definitions

  • the present invention relates to the field of computer systems. More specifically, the present invention relates to computer systems for sequencing biological molecules including nucleic acids.
  • an array of nucleic acid probes is fabricated at known locations on a chip or substrate.
  • a labeled nucleic acid is then brought into contact with the chip and a scanner generates an image file (also called a cell file) indicating the locations where the labeled nucleic acids are bound to the chip.
  • image file also called a cell file
  • Such systems have been used to form, for example, arrays of DNA that may be used to study and detect mutations relevant to genetic diseases, cancers, infectious diseases, HIV, and other genetic characteristics.
  • the VLSIPSTM technology provides methods of making very large arrays of oligonucleotide probes on very small chips. See U.S. Pat. No. 5,143,854 and PCT patent publication Nos. WO 90/15070 and 92/10092, each of which is incorporated by reference for all purposes.
  • the oligonucleotide probes on the DNA probe array are used to detect complementary nucleic acid sequences in a sample nucleic acid of interest (the “target” nucleic acid).
  • the chip may be tiled for a specific target nucleic acid sequence.
  • the chip may contain probes that are perfectly complementary to the target sequence and probes that differ from the target sequence by a single base mismatch. These probes are tiled on a chip in rows and columns of cells, where each cell includes multiple copies of a particular probe. Additionally, “blank” cells may be present on the chip which do not include any probes. As the blank cells contains no probes, labeled targets should not bind specifically to the chip in this area. Thus, a blank cell provides a measure of the background intensity.
  • the chip may include all the possible probes of a specific length. These probes are synthesized on the chip at known locations, typically with multiple copies of a particular probe in a cell. Blank cells may also be utilized to provide a measure of the background intensity.
  • the present invention provides an improved computer-aided system for sequencing sample nucleic acid sequences from nucleic acid hybridization information.
  • the accuracy of nucleic acid sequencing is increased by analyzing the hybridization strength of related probes, where the related probes are identified according to mismatch information among the probes.
  • the related probes may include single base mismatches or otherwise have identical subsequences.
  • the methods of the present invention allow sequencing under conditions that do not allow identification of all of the probes that are perfectly complementary to part of the target nucleic acid sequence.
  • a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes, the nucleic acid probes hybridizing with the nucleic acid sequence under conditions that do not allow identification of all of nucleic acid probes that are perfectly complementary to part of the nucleic acid sequence; and sequencing the nucleic acid sequence according to selected nucleic acid probes.
  • a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes; selecting nucleic acid probes with highest numbers of single base mismatch neighbors among the probes, a single base mismatch neighbor being another probe that has the same sequence except for a single base that is different; and sequencing the nucleic acid sequence according to the selected nucleic acid probes.
  • a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes; selecting nucleic acid probes that have fewer than a predetermined number of base mismatches with another probe; and sequencing the nucleic acid sequence according to the selected nucleic acid probes.
  • a nucleic acid is sequenced by a method including the steps of: contacting a set of oligonucleotide probes of predetermined sequence and length with the nucleic acid under hybridization conditions that do not allow differentiation between (i) those probes of the set which are perfectly complementary to part of the nucleic acid and (ii) those probes that are not perfectly complementary to part of the nucleic acid; selecting a subset of oligonucleotide probes that includes probes that are perfectly complementary to part of the nucleic acid and probes that are not perfectly complementary to part of the nucleic acid; and determining the sequence of the nucleic acid by compiling overlapping sequences of the subset of probes.
  • FIG. 1 illustrates an example of a computer system used to execute the software of the present invention
  • FIG. 2 shows a system block diagram of a typical computer system used to execute the software of the present invention
  • FIG. 3 illustrates an overall system for forming and analyzing arrays of biological materials such as DNA or RNA
  • FIG. 4 is an illustration: of the software for the overall system
  • FIG. 5 illustrates conceptually the binding of probes on chips
  • FIG. 6 shows a high level flow of sequencing utilizing mismatch information
  • FIG. 7 shows a high level flow of another embodiment of sequencing utilizing mismatch information
  • FIG. 8 shows a straight mismatch matrix for use with the process of FIG. 7 ;
  • FIG. 9 shows an offset mismatch matrix for use with the process of FIG. 7 .
  • FIG. 1 illustrates an example of a computer system used to execute the software of the present invention.
  • FIG. 1 shows a computer system 1 which includes a monitor 3 , screen 5 , cabinet 7 , keyboard 9 , and mouse 11 .
  • Mouse 11 may have one or more buttons such as mouse buttons 13 .
  • Cabinet 7 houses a floppy disk drive 14 and a hard drive (not shown) that may be utilized to store and retrieve software programs including computer or software code incorporating the present invention. Although a floppy disk 15 is shown as the computer readable storage medium, other computer readable storage media including CD-ROM, DRAM, hard drives, flash memory, tape, and the like may be utilized.
  • Cabinet 7 also houses familiar computer components (not shown) such as a processor, memory, and the like.
  • FIG. 2 shows a system block diagram of computer system 1 used to execute the software of the present invention.
  • computer system 1 includes monitor 3 and keyboard 9 .
  • Computer system 1 further includes subsystems such as a central processor 52 , system memory 54 , I/O controller 56 , display adapter 58 , serial port 62 , disk 64 , network interface 66 , and speaker 68 .
  • Disk 64 is representative of an internal hard drive, floppy drive, CD-ROM, flash memory, tape, or any other storage medium.
  • Other computer systems suitable for use with the present invention may include additional or fewer subsystems.
  • another computer system could include more than one processor 52 (i.e., a multi-processor system) or memory cache.
  • Arrows such as 70 represent the system bus architecture of computer system 1 . However, these arrows are illustrative of any interconnection scheme serving to link the subsystems. For example, speaker 68 could be connected to the other subsystems through a port or have an internal direct connection to central processor 52 .
  • Computer system 1 shown in FIG. 2 is but an example of a computer system suitable for use with the present invention. Other configurations of subsystems suitable for use with the present invention will be readily apparent to one of ordinary skill in the art.
  • the present invention provides methods of analyzing hybridization intensity files for a chip containing hybridized nucleic acid probes.
  • the files represent fluorescence data from a biological array, but the files may also represent other data such as radioactive intensity, light scattering, refractive index, conductivity, electroluminescence, or large molecule detection data. Therefore, the present invention is not limited to analyzing fluorescence measurements of hybridizations but may be readily utilized to analyze other measurements of hybridization.
  • the present invention is described as being part of a computer system that designs a chip mask, synthesizes the probes on the chip, labels the nucleic acids, and scans the hybridized nucleic acid probes.
  • a computer system that designs a chip mask, synthesizes the probes on the chip, labels the nucleic acids, and scans the hybridized nucleic acid probes.
  • Such a system is fully described in U.S. patent application Ser. No. 08/249,188 which has been incorporated by reference for all purposes.
  • the present invention may be used separately from the overall system for analyzing data generated by such systems.
  • FIG. 3 illustrates a computerized system for forming and analyzing arrays of biological materials such as RNA or DNA.
  • a computer 100 is used to design arrays of biological polymers such as RNA or DNA.
  • the computer 100 may be, for example, an appropriately programmed Sun Workstation or personal computer or workstation, such as an IBM PC equivalent, including appropriate memory and a CPU as shown in FIGS. 1 and 2 .
  • the computer system 100 obtains inputs from a user regarding characteristics of a gene of interest, and other inputs regarding the desired features of the array.
  • the computer system may obtain information regarding a specific genetic sequence of interest from an external or internal database 102 such as GenBank.
  • the output of the computer system 100 is a set of chip design computer files 104 in the form of, for example, a switch matrix, as described in PCT application WO 92/10092, and other associated computer files.
  • the chip design files are provided to a system 106 that designs the lithographic masks used in the fabrication of arrays of molecules such as DNA.
  • the system or process 106 may include the hardware necessary to manufacture masks 110 and also the necessary computer hardware and software 108 necessary to lay the mask patterns out on the mask in an efficient manner. As with the other features in FIG. 3 , such equipment may or may not be located at the same physical site, but is shown together for ease of illustration in FIG. 3 .
  • the system 106 generates masks 110 or other synthesis patterns such as chrome-on-glass masks for use in the fabrication of polymer arrays.
  • Synthesis system 112 includes the necessary hardware and software used to fabricate arrays of polymers on a substrate or chip 114 .
  • synthesizer 112 includes a light source 116 and a chemical flow cell 118 on which the substrate or chip 114 is placed.
  • Mask 110 is placed between the light source and the substrate/chip, and the two are translated relative to each other at appropriate times for deprotection of selected regions of the chip.
  • Selected chemical reagents are directed through flow cell 118 for coupling to deprotected regions, as well as for washing and other operations. All operations are preferably directed by an appropriately programmed computer 119 , which may or may not be the same computer as the computer(s) used in mask design and mask making;
  • the substrates fabricated by synthesis system 112 are optionally diced into smaller chips and exposed to marked targets.
  • the targets may or may not be complementary to one or more of the molecules on the substrate.
  • the targets are marked with a label such as a fluorescein label (indicated by an asterisk in FIG. 3 ) and placed in scanning system 120 .
  • Scanning system 120 again operates under the direction of an appropriately programmed digital computer 122 , which also may or may not be the same computer as the computers used in synthesis, mask making, and mask design.
  • the scanner 120 includes a detection device 124 such as a confocal microscope or CCD (charge-coupled device) that is used to detect the location where labeled target (*) has bound to the substrate.
  • a detection device 124 such as a confocal microscope or CCD (charge-coupled device) that is used to detect the location where labeled target (*) has bound to the substrate.
  • the output of scanner 120 is an image file(s) 124 indicating, in the case of fluorescein labeled target, the fluorescence intensity (photon counts or other related measurements, such as voltage) as a function of position on the substrate. Since higher photon counts will be observed where the labeled target has bound more strongly to the array of polymers (e.g., DNA probes on the substrate), and since the monomer sequence of the polymers on the substrate is known as a function of position, it becomes possible to determine the sequence(s) of polymer(s) on the substrate that are complementary to the target.
  • image file(s) 124 indicating, in the case of fluorescein labeled target, the fluorescence intensity (photon counts or other related measurements, such as voltage) as a function of position on the substrate. Since higher photon counts will be observed where the labeled target has bound more strongly to the array of polymers (e.g., DNA probes on the substrate), and since the monomer sequence of the polymers on the substrate is known as a function of position
  • the image file 124 is provided as input to an analysis system 126 that incorporates the visualization and analysis methods of the present invention.
  • the analysis system may be any one of a wide variety of computer system(s), but in a preferred embodiment the analysis system is based on a Sun Workstation or equivalent.
  • the present invention provides various methods of analyzing the chip design files and the image files, providing appropriate output 128 .
  • the present invention may further be used to identify specific mutations in a target such as DNA or RNA.
  • FIG. 4 provides a simplified illustration of the overall software system used in the operation of one embodiment of the invention.
  • the system first identifies the genetic sequence(s) or targets that would be of interest in a particular analysis at step 202 .
  • the sequences of interest may identify a virus, microorganism or individual. Additionally, the sequence of interest may provide information about genetic diseases, cancers or infectious diseases. Sequence selection may be provided via manual input of text files or may be from external sources such as GenBank. In a preferred embodiment that performs de novo sequencing of target nucleic acids, this steps is not necessary as the chip includes all the possible n-mer probes (where n represents the length of the nucleic acid probe).
  • a chip may be synthesized to include cells containing all the possible probes of a specific length. For example, a chip may be synthesized that includes all the possible 8-mer DNA probes. Such a chip would have 65,536 cells (4*4*4*4*4*4*4*4*4), with each cell corresponding to a particular probe. A chip may also include other probes including all the probes of other lengths.
  • the system determines which probes would be desirable on the chip, and provides an appropriate “layout” on the chip for the probes.
  • the layout implements desired characteristics such as an arrangement on the chip that permits “reading” of genetic sequence and/or minimization of edge effects, ease of synthesis, and the like.
  • the masks for the synthesis are designed.
  • the software utilizes the mask design and layout information to make the DNA or other polymer chips. This software 208 will control, among other things, relative translation of a substrate and the mask, the flow of desired reagents through a flow cell, the synthesis temperature of the flow cell, and other parameters.
  • another piece of software is used in scanning a chip thus synthesized and exposed to a labeled target. The software controls the scanning of the chip, and stores the data thus obtained in a file that may later be utilized to extract sequence information.
  • a computer system utilizes the layout information and the fluorescence information to evaluate the hybridized nucleic acid probes on the chip.
  • important pieces of information obtained from DNA probe arrays are the identification of mutant targets and determination of the genetic sequence of a particular target.
  • FIG. 5 illustrates the binding of a particular target DNA to an array of DNA probes 114 .
  • the following probes are formed in the array: 3′-AGAACGT AGACCGT AGAGCGT AGATCGT . . .
  • fluorescein-labeled (or otherwise marked) target 5′-TCTTGCA is exposed to the array, it is complementary only to the probe 3′-AGAACGT, and fluorescein will be primarily found on the surface of the chip where 3′-AGAACGT is located.
  • the chip contains cells that include multiple copies of a particular probe.
  • the image file will contain fluorescence intensities, one for each probe (or cell). By analyzing the fluorescence intensities associated with a specific probe, it becomes possible to extract sequence information from such arrays using the methods of the invention disclosed herein.
  • mismatch information relates to base mismatches between or among nucleic acid probes. Mismatch information may include the number of base mismatches, the location of the base mismatches and the base differences. The mismatch information maybe combined with information from the hybridization intensity to sequence the nucleic acid sequence with a high degree of accuracy. In a preferred embodiment, the present invention is utilized for de novo sequencing of nucleic acids.
  • the target nucleic acid is an 8-mer (meaning that the target nucleic acid consists of eight bases or nucleotides) and that the target is exposed to a chip including the complete set of 8-mer probes.
  • the 1 probe that is perfectly complementary to the target and the 24 probes that contain a single base mismatch have the highest hybridization intensities because they hybridize most strongly to the target.
  • the target is ACTGGTCT-3′
  • the following would be the probes having the highest measured intensities in this example: Perfect complement TGACCAGA-5′ One base mismatches GGACCAGA-5′ AGACCAGA-5′ CGACCAGA-5′ TAACCAGA-5′ TCACCAGA-5′ TTACCAGA-5′ . . . and so forth for the other six positions.
  • the set of these 25 probes may be analyzed to sequence the target nucleic acid. Although typically, the target nucleic acid is longer than the probes, the example provides a good illustration of aspects of the present invention.
  • probes that are perfectly or exactly complementary to the target may not have the highest hybridization intensities. Therefore, a probe that is perfectly complementary to the target often cannot be identified from the rank order of hybridization intensities.
  • the present invention utilizes mismatch information among the nucleic acid probes to sequence the target where all of the nucleic acid probes that are perfectly complementary to part of the target may not be readily identified.
  • the present invention identifies neighbor-rich probes which are then utilized to sequence the target nucleic acid.
  • a “neighbor-rich probe” is a probe that is related to many other probes in the probe space by a single base mismatch.
  • a probe that has a single base mismatch with another probe will be referred to as a “single base mismatch neighbor.”
  • Neighbor-rich probes may be identified according to mismatch information as follows.
  • each probe in the set is compared to the other probes to determine how the probe's sequence compares to the other probes.
  • this one probe is related to or has 24 single base mismatch neighbors.
  • the perfectly complementary probe may be identified as a neighbor-rich probe from the mismatch information because it has many single mismatch neighbor probes in the probe space.
  • the perfectly complementary probe had 8 times as many single mismatch neighbor probes as nearly-complementary probes.
  • mismatch information may be utilized for de novo sequencing of a target nucleic acid where oligonucleotide probes are contacted with the target under conditions that do not allow differentiation between those probes that are perfectly complementary to part of the target and those probes that are not.
  • the sequence of the target was known. However, in many applications including de novo sequencing the sequence of the target is unknown. Nevertheless, the example is useful in demonstrating how neighbor-rich probes may be identified.
  • FIG. 6 shows a high level flow of sequencing utilizing mismatch information.
  • hybridization intensities from probes or other data indicative of binding affinity are input into the system.
  • the system may receive the hybridization intensities many different ways.
  • the system may operate the scanning device directly, the system may receive the hybridization intensities from another computer system that measured the intensities, or an operator may manually enter the data.
  • There may be thousands or tens of thousands of hybridization intensities that correspond to nucleic acid probes on a chip.
  • the chip typically includes all possible probes of a specific length in order to sequence the target.
  • the system selects a set of probes associated with the highest hybridization intensities (i.e., that show the strongest binding affinity). Selecting the probes with the highest hybridization intensities may be done in any number of ways. For example, the system may use an intensity threshold value and select the probes whose hybridization intensities are higher than the intensity threshold (e.g., 100 photon counts). The system may select a specific number or percentage of probes (e.g., 50 probes or the top 10%) that have the highest hybridization intensities. Additionally, the system may select the probes that have a hybridization intensity greater than a specific percentage (e.g., 40%) of the highest hybridization intensity.
  • a specific percentage e.g. 40%
  • the system calculates the number of single base mismatch neighbors for each probe in the set at step 504 . For example, in one embodiment, a probe is first selected in order to compare the selected probe to the other probes. The system then determines how many of the other probes in the set are identical to the selected probe except for a single base mismatch at one base position. The number of single base mismatch neighbors is calculated for each of the probes having the highest hybridization intensities. Additionally, the system may calculate and utilize the number of double base mismatches in an extension of the concepts herein.
  • the system selects the probes in the set with the highest number of single base mismatch neighbors. Selecting the probes with the highest number of neighbors may be done in any number of ways including utilizing a threshold, a specific number of the probes, or greater than a specific percentage of the highest number of neighbors. In one embodiment, the system selects the probes in the set with the highest number of neighbors and the highest hybridization intensities. In other words, a second intensity threshold is utilized to further reduce the set of probes.
  • the selected probes with the highest number of single base mismatch neighbors are the neighbor-rich probes.
  • the neighbor-rich probes are then aligned at step 508 .
  • the neighbor-rich probes are aligned or compiled so that they have the most bases in common.
  • neighbor-rich probes that have a single base mismatch are aligned to form an aligned set of probes.
  • Aligned sets of probes are then aligned in an offset fashion in the way that reduces the number of base mismatches between sets of probes.
  • the aligned probes are utilized to sequence the target nucleic acid sequence.
  • the target may be sequenced in many different ways including the formation of a consensus sequence may be produced as described in the following example.
  • a target of 5′-AGTTGTAGTGGATGG was exposed to a chip containing 8-mer probes.
  • the highest hybridization intensity was 331 photon counts.
  • An intensity threshold of 90 photon counts was utilized and there were 133 probes that had a hybridization intensity greater than the intensity threshold of 90 .
  • each probe was compared to the other 132 probes to determine the number of single and double base mismatches the probe had with the other probes.
  • the highest number of single base mismatch neighbor probes was 12 and the probes are presented in decreasing order according to the number of single base mismatch neighbor probes.
  • the neighbor-rich probes were utilized to sequence the target nucleic acid sequence.
  • the system utilized the frequency of bases at each position to produce a consensus sequence, where a “consensus sequence” is a sequence generated by neighbor-rich probes to sequence the target.
  • the system aligned the neighbor-rich probes so that each probe had the highest number of bases in common with other probes.
  • the following are the aligned neighbor-rich probes with the complement of the target sequence shown for reference.
  • the target sequence is known in this example but the target sequence may be an unknown sequence or only partially known.
  • the system counts the number or frequency of each base (A, C, G and T) at each position. After the frequency of bases at each position is calculated, the system produces a consensus sequence. In one embodiment, the base that occurred the most at a position is utilized to produce the consensus sequence if the base occurred more than 2 times and the frequency that the base occurred is greater than 50% at that position.
  • the following is a matrix of base vs.
  • the consensus sequence is the complement of the target; thus, the target is sequenced according to the complement of the consensus sequence.
  • the target sequence was 5′-AGTTGTAGTGGATGG and it was correctly sequenced as 5′-GTTGTAGTGGAT (the terminal bases being ambiguous).
  • the parameters for producing the consensus sequence may be varied according to the experimental data. For example, if the consensus sequence was formed solely by the bases that occur most often at each position, the consensus sequence would be perfectly complementary to the target nucleic acid for this data. However, this will not always be the case.
  • Software Appendix A provides a BASIC source code illustration of this embodiment of the invention.
  • the source code is written in Quick BASIC for an IBM compatible personal computer.
  • FIG. 7 shows a high level flow of another embodiment of sequencing utilizing mismatch information.
  • hybridization intensities are input into the system.
  • the system may receive the hybridization intensities many different ways.
  • the system may operate the scanning device directly, the system may receive the hybridization intensities from another computer system that measured the intensities, or an operator may manually enter the data.
  • the system selects a set of probes associated with the highest hybridization intensities. Selecting the probes with the highest hybridization intensities may be done in any number of ways. For example, the system may use an intensity threshold value and select the probes whose hybridization intensities are higher than the intensity threshold. The system may select a specific number or percentage of probes that have the highest hybridization intensities. Alternatively, the system may select the probes that have a hybridization greater than a specific percentage of the highest hybridization intensity.
  • the system calculates the number of straight mismatches for each probe in the set at step 554 .
  • “Straight mismatches” refers to base mismatches between probes where the bases at corresponding positions are compared (e.g., as was done in the previous embodiment).
  • a probe 3′-AACAT is compared to a probe 3′-AACTT by aligning the probes as follows: 3′-AACAT 3′-AACTT
  • a matrix is formed to show the straight mismatches between probes.
  • FIG. 8 shows a straight mismatch matrix for 50 probes.
  • each probe is assigned an identification number from 1 to 50.
  • the numbers on the rows and columns of the matrix correspond to the identification number of the probe.
  • the value in the matrix represents the number of straight mismatches between the probes designated by the row and column. If there are more than 2 mismatches, the matrix shows a “.” at the appropriate matrix location. Since a diagonal 600 of the matrix shows the number of straight mismatches between the same probe, the diagonal contains 0's because the probe is being compared to itself (i.e., there are no straight mismatches). Also, the matrix is identical on each side of diagonal 600 .
  • Offset mismatches refers to base mismatches between probes where the bases at offset positions are compared.
  • a probe 3′-CGAATCAA is compared to a probe 3′-GCATCAAC by aligning the probes as follows: 3′-CGAATCAA 3′-GCATCCAC
  • there are two mismatches double mismatch
  • a matrix is formed to show the offset mismatches between probes.
  • FIG. 9 shows an offset mismatch matrix for 50 probes. As in FIG. 8 , the rows and columns of the matrix correspond to the identification number of the probe, which is 1 to 50. The value in the matrix represents the number of offset mismatches between the probes designated by the row and column. If there are more than 2 mismatches, the matrix shows a “.” at the appropriate matrix location. As shown, a diagonal will not contain O's and the matrix is not identical on each side of the diagonal. Although the probes were offset a single base position, the probes may be offset more positions when they are compared in an extension of the principles herein.
  • the system selects the probes with less than some small number straight mismatches and less than some small number of offset mismatches.
  • the system identifies the probes in the matrices that have less than 2 straight mismatches and less than 3 offset mismatches. The parameters for selecting these probes with few mismatches may be varied according to the experimental data.
  • the selected probes are then aligned at step 560 .
  • the probes are aligned so that they have the most bases in common.
  • the mismatch information concerning the straight mismatches and offset mismatches is utilized to align the probes so that the number of mismatches between the probes is reduced.
  • the aligned probes are utilized to sequence the target nucleic acid sequence.
  • the target may be sequenced in many different ways. For example, a consensus sequence may be produced as described in the following example.
  • FIGS. 11 and 12 show the straight and offset mismatch matrices for the fifty probes (steps 554 and 556 ). Forty-seven probes were selected that have less than 2 straight mismatches with at least one other probe and less than 3 offset mismatches with at least one other probe (step 558 ).
  • the straight and offset mismatch information was utilized to align the 47 probes (step 560 ).
  • FIG. 9 shows that the probe identified as 2 on the row had 0 offset mismatches with the probe identified as 1 on the column. Therefore, probes 2 and 1 align well if they are offset a single base position.
  • AACATCACCT CAACATCACC ACATCACCTA ACAACATCAC CAACATCACA ATCACCTACC AACATCACCA CACATCACCT AACATCACCG TCAACATCAC CATCACCTAC CCAACATCAC ACACCTACCA CAACATCACG GCAACATCAC ACATCACCTT AACATCACCC AGCACCTACC AAACATCACC ACCATCACCT ACATCACCAT CACCTACCAA ACATCACCTC ATCACCTACA ACATCACCTG ACACATCACC GAACATCACC TACATCACCT TAACATCACC CACCTACCAG GACATCACCT CACATCACCA ATCATCACCT ACCTACCATC ACAACATCAA CAACATCACT ACATCACCAA CACATCACCG ACATCACCCT CACCTACCAC CATCACCTAA TCACCTACCA CACCTACCAT CCAT CCATCACCTA ACATCACCCA ACATCACCGA ATCAACATCA
  • the system counts the number of occurrences of each base (A, C, G and T) at each position. After the frequency of bases at each position is calculated, the system produces a consensus sequence which should be complementary to the target sequence. If the system utilized bases that occurred more than 2 times and the frequency that the base occurred is greater than 50% at that position, the consensus sequence 3′-CAACATCACCTACCA is produced.
  • the consensus sequence is ideally the complement of the target; thus, the target is sequenced according to the complement of the consensus sequence (consensus' where the prime denotes the complement).
  • the target and consensus' sequence were as follows: Target AGTTGTAGTGGATGGT Consensus′ GTTGTAGTGGATGGT (one terminal base of the consensus sequence being ambiguous).
  • the target was sequenced with a high degree of accuracy utilizing mismatch information.
  • the parameters for producing the consensus sequence may be varied according to the experimental data.
  • Software Appendix B provides a BASIC source code illustration of this embodiment of the invention.
  • the source code is written in Quick BASIC for an IBM compatible personal computer.

Abstract

A computer system for sequencing nucleic acids is provided. The computer system may be used for de novo sequencing of a nucleic acid sequence by analyzing the fluorescence intensities of hybridized nucleic acid probes on biological chips. The probes with the highest intensities are utilized to sequence the nucleic acid and related probes are analyzed to increase the accuracy of nucleic acid sequencing. The sequence of the nucleic acid sequence may be determined from hybridization intensities that do not allow identification of perfectly complementary probes.

Description

    COPYRIGHT NOTICE
  • A portion of the disclosure of this patent document contains material which is subject to copyright protection. The copyright owner has no objection to the xerographic reproduction by anyone of the patent document or the patent disclosure in exactly the form it appears in the Patent and Trademark Office patent file or records, but otherwise reserves all copyright rights whatsoever.
  • SOFTWARE APPENDICES
  • Software Appendices A and B comprising six (6) sheets are included herewith.
  • BACKGROUND OF THE INVENTION
  • The present invention relates to the field of computer systems. More specifically, the present invention relates to computer systems for sequencing biological molecules including nucleic acids.
  • Devices and computer systems for forming and using arrays of materials on a substrate are known. For example, PCT applications WO92/10588 and 95/11995, incorporated herein by reference for all purposes, describe techniques for sequencing or sequence checking nucleic acids and other materials. Arrays for performing these operations may be formed in arrays according to the methods of, for example, the pioneering techniques disclosed in U.S. Pat. Nos. 5,445,934 and 5,384,261, and U.S. patent application Ser. No. 08/249,188, each incorporated herein by reference for all purposes.
  • According to one aspect of the techniques described therein, an array of nucleic acid probes is fabricated at known locations on a chip or substrate. A labeled nucleic acid is then brought into contact with the chip and a scanner generates an image file (also called a cell file) indicating the locations where the labeled nucleic acids are bound to the chip. Based upon the image file and identities of the probes at specific locations, it becomes possible to extract information such as the nucleotide or monomer sequence of DNA or RNA. Such systems have been used to form, for example, arrays of DNA that may be used to study and detect mutations relevant to genetic diseases, cancers, infectious diseases, HIV, and other genetic characteristics.
  • The VLSIPS™ technology provides methods of making very large arrays of oligonucleotide probes on very small chips. See U.S. Pat. No. 5,143,854 and PCT patent publication Nos. WO 90/15070 and 92/10092, each of which is incorporated by reference for all purposes. The oligonucleotide probes on the DNA probe array are used to detect complementary nucleic acid sequences in a sample nucleic acid of interest (the “target” nucleic acid).
  • For sequence checking applications, the chip may be tiled for a specific target nucleic acid sequence. For example, the chip may contain probes that are perfectly complementary to the target sequence and probes that differ from the target sequence by a single base mismatch. These probes are tiled on a chip in rows and columns of cells, where each cell includes multiple copies of a particular probe. Additionally, “blank” cells may be present on the chip which do not include any probes. As the blank cells contains no probes, labeled targets should not bind specifically to the chip in this area. Thus, a blank cell provides a measure of the background intensity.
  • For de novo sequencing applications, the chip may include all the possible probes of a specific length. These probes are synthesized on the chip at known locations, typically with multiple copies of a particular probe in a cell. Blank cells may also be utilized to provide a measure of the background intensity.
  • SUMMARY OF THE INVENTION
  • The present invention provides an improved computer-aided system for sequencing sample nucleic acid sequences from nucleic acid hybridization information. The accuracy of nucleic acid sequencing is increased by analyzing the hybridization strength of related probes, where the related probes are identified according to mismatch information among the probes. The related probes may include single base mismatches or otherwise have identical subsequences. The methods of the present invention allow sequencing under conditions that do not allow identification of all of the probes that are perfectly complementary to part of the target nucleic acid sequence.
  • According to one aspect of the present invention, a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes, the nucleic acid probes hybridizing with the nucleic acid sequence under conditions that do not allow identification of all of nucleic acid probes that are perfectly complementary to part of the nucleic acid sequence; and sequencing the nucleic acid sequence according to selected nucleic acid probes.
  • According to another aspect of the present invention, a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes; selecting nucleic acid probes with highest numbers of single base mismatch neighbors among the probes, a single base mismatch neighbor being another probe that has the same sequence except for a single base that is different; and sequencing the nucleic acid sequence according to the selected nucleic acid probes.
  • According to another aspect of the present invention, a computer system is used to sequence a nucleic acid by a method including the steps of: inputting hybridization intensities for a plurality of nucleic acid probes; selecting nucleic acid probes that have fewer than a predetermined number of base mismatches with another probe; and sequencing the nucleic acid sequence according to the selected nucleic acid probes.
  • According to another aspect of the present invention, a nucleic acid is sequenced by a method including the steps of: contacting a set of oligonucleotide probes of predetermined sequence and length with the nucleic acid under hybridization conditions that do not allow differentiation between (i) those probes of the set which are perfectly complementary to part of the nucleic acid and (ii) those probes that are not perfectly complementary to part of the nucleic acid; selecting a subset of oligonucleotide probes that includes probes that are perfectly complementary to part of the nucleic acid and probes that are not perfectly complementary to part of the nucleic acid; and determining the sequence of the nucleic acid by compiling overlapping sequences of the subset of probes.
  • A further understanding of the nature and advantages of the inventions herein may be realized by reference to the remaining portions of the specification and the attached drawings.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 illustrates an example of a computer system used to execute the software of the present invention;
  • FIG. 2 shows a system block diagram of a typical computer system used to execute the software of the present invention;
  • FIG. 3 illustrates an overall system for forming and analyzing arrays of biological materials such as DNA or RNA;
  • FIG. 4 is an illustration: of the software for the overall system;
  • FIG. 5 illustrates conceptually the binding of probes on chips;
  • FIG. 6 shows a high level flow of sequencing utilizing mismatch information;
  • FIG. 7 shows a high level flow of another embodiment of sequencing utilizing mismatch information;
  • FIG. 8 shows a straight mismatch matrix for use with the process of FIG. 7; and
  • FIG. 9 shows an offset mismatch matrix for use with the process of FIG. 7.
  • DESCRIPTION OF THE PREFERRED EMBODIMENTS
  • Introduction
  • In the description that follows, the present invention will be described in reference to a Sun Workstation in a UNIX environment. The present invention, however, is not limited to any particular hardware or operating system environment. Instead, those skilled in the art will find that the systems and methods of the present invention may be advantageously applied to a variety of systems, including IBM personal computers running MS-DOS or Microsoft Windows. Therefore, the following description of specific systems are for purposes of illustration and not limitation.
  • FIG. 1 illustrates an example of a computer system used to execute the software of the present invention. FIG. 1 shows a computer system 1 which includes a monitor 3, screen 5, cabinet 7, keyboard 9, and mouse 11. Mouse 11 may have one or more buttons such as mouse buttons 13. Cabinet 7 houses a floppy disk drive 14 and a hard drive (not shown) that may be utilized to store and retrieve software programs including computer or software code incorporating the present invention. Although a floppy disk 15 is shown as the computer readable storage medium, other computer readable storage media including CD-ROM, DRAM, hard drives, flash memory, tape, and the like may be utilized. Cabinet 7 also houses familiar computer components (not shown) such as a processor, memory, and the like.
  • FIG. 2 shows a system block diagram of computer system 1 used to execute the software of the present invention. As in FIG. 1, computer system 1 includes monitor 3 and keyboard 9. Computer system 1 further includes subsystems such as a central processor 52, system memory 54, I/O controller 56, display adapter 58, serial port 62, disk 64, network interface 66, and speaker 68. Disk 64 is representative of an internal hard drive, floppy drive, CD-ROM, flash memory, tape, or any other storage medium. Other computer systems suitable for use with the present invention may include additional or fewer subsystems. For example, another computer system could include more than one processor 52 (i.e., a multi-processor system) or memory cache.
  • Arrows such as 70 represent the system bus architecture of computer system 1. However, these arrows are illustrative of any interconnection scheme serving to link the subsystems. For example, speaker 68 could be connected to the other subsystems through a port or have an internal direct connection to central processor 52. Computer system 1 shown in FIG. 2 is but an example of a computer system suitable for use with the present invention. Other configurations of subsystems suitable for use with the present invention will be readily apparent to one of ordinary skill in the art.
  • The present invention provides methods of analyzing hybridization intensity files for a chip containing hybridized nucleic acid probes. In a representative embodiment, the files represent fluorescence data from a biological array, but the files may also represent other data such as radioactive intensity, light scattering, refractive index, conductivity, electroluminescence, or large molecule detection data. Therefore, the present invention is not limited to analyzing fluorescence measurements of hybridizations but may be readily utilized to analyze other measurements of hybridization.
  • For purposes of illustration, the present invention is described as being part of a computer system that designs a chip mask, synthesizes the probes on the chip, labels the nucleic acids, and scans the hybridized nucleic acid probes. Such a system is fully described in U.S. patent application Ser. No. 08/249,188 which has been incorporated by reference for all purposes. However, the present invention may be used separately from the overall system for analyzing data generated by such systems.
  • FIG. 3 illustrates a computerized system for forming and analyzing arrays of biological materials such as RNA or DNA. A computer 100 is used to design arrays of biological polymers such as RNA or DNA. The computer 100 may be, for example, an appropriately programmed Sun Workstation or personal computer or workstation, such as an IBM PC equivalent, including appropriate memory and a CPU as shown in FIGS. 1 and 2. The computer system 100 obtains inputs from a user regarding characteristics of a gene of interest, and other inputs regarding the desired features of the array. Optionally, the computer system may obtain information regarding a specific genetic sequence of interest from an external or internal database 102 such as GenBank. The output of the computer system 100 is a set of chip design computer files 104 in the form of, for example, a switch matrix, as described in PCT application WO 92/10092, and other associated computer files.
  • The chip design files are provided to a system 106 that designs the lithographic masks used in the fabrication of arrays of molecules such as DNA. The system or process 106 may include the hardware necessary to manufacture masks 110 and also the necessary computer hardware and software 108 necessary to lay the mask patterns out on the mask in an efficient manner. As with the other features in FIG. 3, such equipment may or may not be located at the same physical site, but is shown together for ease of illustration in FIG. 3. The system 106 generates masks 110 or other synthesis patterns such as chrome-on-glass masks for use in the fabrication of polymer arrays.
  • The masks 110, as well as selected information relating to the design of the chips from system 100, are used in a synthesis system 112. Synthesis system 112 includes the necessary hardware and software used to fabricate arrays of polymers on a substrate or chip 114. For example, synthesizer 112 includes a light source 116 and a chemical flow cell 118 on which the substrate or chip 114 is placed. Mask 110 is placed between the light source and the substrate/chip, and the two are translated relative to each other at appropriate times for deprotection of selected regions of the chip. Selected chemical reagents are directed through flow cell 118 for coupling to deprotected regions, as well as for washing and other operations. All operations are preferably directed by an appropriately programmed computer 119, which may or may not be the same computer as the computer(s) used in mask design and mask making;
  • The substrates fabricated by synthesis system 112 are optionally diced into smaller chips and exposed to marked targets. The targets may or may not be complementary to one or more of the molecules on the substrate. The targets are marked with a label such as a fluorescein label (indicated by an asterisk in FIG. 3) and placed in scanning system 120. Scanning system 120 again operates under the direction of an appropriately programmed digital computer 122, which also may or may not be the same computer as the computers used in synthesis, mask making, and mask design. The scanner 120 includes a detection device 124 such as a confocal microscope or CCD (charge-coupled device) that is used to detect the location where labeled target (*) has bound to the substrate. The output of scanner 120 is an image file(s) 124 indicating, in the case of fluorescein labeled target, the fluorescence intensity (photon counts or other related measurements, such as voltage) as a function of position on the substrate. Since higher photon counts will be observed where the labeled target has bound more strongly to the array of polymers (e.g., DNA probes on the substrate), and since the monomer sequence of the polymers on the substrate is known as a function of position, it becomes possible to determine the sequence(s) of polymer(s) on the substrate that are complementary to the target.
  • The image file 124 is provided as input to an analysis system 126 that incorporates the visualization and analysis methods of the present invention. Again, the analysis system may be any one of a wide variety of computer system(s), but in a preferred embodiment the analysis system is based on a Sun Workstation or equivalent. The present invention provides various methods of analyzing the chip design files and the image files, providing appropriate output 128. The present invention may further be used to identify specific mutations in a target such as DNA or RNA.
  • FIG. 4 provides a simplified illustration of the overall software system used in the operation of one embodiment of the invention. As shown in FIG. 4, in some cases (such as sequence checking systems) the system first identifies the genetic sequence(s) or targets that would be of interest in a particular analysis at step 202. The sequences of interest may identify a virus, microorganism or individual. Additionally, the sequence of interest may provide information about genetic diseases, cancers or infectious diseases. Sequence selection may be provided via manual input of text files or may be from external sources such as GenBank. In a preferred embodiment that performs de novo sequencing of target nucleic acids, this steps is not necessary as the chip includes all the possible n-mer probes (where n represents the length of the nucleic acid probe).
  • For de novo sequencing, a chip may be synthesized to include cells containing all the possible probes of a specific length. For example, a chip may be synthesized that includes all the possible 8-mer DNA probes. Such a chip would have 65,536 cells (4*4*4*4*4*4*4*4), with each cell corresponding to a particular probe. A chip may also include other probes including all the probes of other lengths.
  • At step 204 the system determines which probes would be desirable on the chip, and provides an appropriate “layout” on the chip for the probes. The layout implements desired characteristics such as an arrangement on the chip that permits “reading” of genetic sequence and/or minimization of edge effects, ease of synthesis, and the like.
  • Again referring to FIG. 4, at step 206 the masks for the synthesis are designed. At step 208 the software utilizes the mask design and layout information to make the DNA or other polymer chips. This software 208 will control, among other things, relative translation of a substrate and the mask, the flow of desired reagents through a flow cell, the synthesis temperature of the flow cell, and other parameters. At step 210, another piece of software is used in scanning a chip thus synthesized and exposed to a labeled target. The software controls the scanning of the chip, and stores the data thus obtained in a file that may later be utilized to extract sequence information.
  • At step 212 a computer system according to the present invention utilizes the layout information and the fluorescence information to evaluate the hybridized nucleic acid probes on the chip. Among the important pieces of information obtained from DNA probe arrays are the identification of mutant targets and determination of the genetic sequence of a particular target.
  • FIG. 5 illustrates the binding of a particular target DNA to an array of DNA probes 114. As shown in this simple example, the following probes are formed in the array:
    3′-AGAACGT
       AGACCGT
       AGAGCGT
       AGATCGT
    .
    .
    .

    As shown, when the fluorescein-labeled (or otherwise marked) target 5′-TCTTGCA is exposed to the array, it is complementary only to the probe 3′-AGAACGT, and fluorescein will be primarily found on the surface of the chip where 3′-AGAACGT is located. The chip contains cells that include multiple copies of a particular probe. Thus, the image file will contain fluorescence intensities, one for each probe (or cell). By analyzing the fluorescence intensities associated with a specific probe, it becomes possible to extract sequence information from such arrays using the methods of the invention disclosed herein.
  • For ease of reference, one may call bases by assigning the bases the following codes:
    Code Group Meaning
    A A Adenine
    C C Cytosine
    G G Guanine
    T T(U) Thymine (Uracil)
    M A or C aMino
    R A or G puRine
    W A or T(U) Weak interaction
    (2 H bonds)
    Y C or T(U) pYrimidine
    S C or G Strong interaction
    (3 H bonds)
    K G or T(U) Keto
    V A, C or G not T(U)
    H A, C or T(U) not G
    D A, G or T(U) not C
    B C, G or T(U) not A
    N A, C, G, or T(U) Insufficient intensity
    to call
    X A, C, G, or T(U) Insufficient
    discrimination to
    call

    Most of the codes conform to the IUPAC standard. However, code N has been redefined and code X has been added.
    Sequencing Utilizing Mismatch Information
  • The present invention provides methods of sequencing nucleic acid sequences utilizing mismatch information. When used herein, “mismatch information” relates to base mismatches between or among nucleic acid probes. Mismatch information may include the number of base mismatches, the location of the base mismatches and the base differences. The mismatch information maybe combined with information from the hybridization intensity to sequence the nucleic acid sequence with a high degree of accuracy. In a preferred embodiment, the present invention is utilized for de novo sequencing of nucleic acids.
  • In order to illustrate what mismatch information or the pattern of mismatches may include, a hypothetical example will be described. Suppose the target nucleic acid is an 8-mer (meaning that the target nucleic acid consists of eight bases or nucleotides) and that the target is exposed to a chip including the complete set of 8-mer probes. In order to simplify this example, further assume that the 1 probe that is perfectly complementary to the target and the 24 probes that contain a single base mismatch (i.e., perfectly complementary except for a single base mismatch) have the highest hybridization intensities because they hybridize most strongly to the target.
  • Thus, if the target is ACTGGTCT-3′, the following would be the probes having the highest measured intensities in this example:
    Perfect complement TGACCAGA-5′
    One base mismatches GGACCAGA-5′
    AGACCAGA-5′
    CGACCAGA-5′
    TAACCAGA-5′
    TCACCAGA-5′
    TTACCAGA-5′
    .
    .
    .

    and so forth for the other six positions. The set of these 25 probes may be analyzed to sequence the target nucleic acid. Although typically, the target nucleic acid is longer than the probes, the example provides a good illustration of aspects of the present invention.
  • For many reasons, probes that are perfectly or exactly complementary to the target may not have the highest hybridization intensities. Therefore, a probe that is perfectly complementary to the target often cannot be identified from the rank order of hybridization intensities. The present invention utilizes mismatch information among the nucleic acid probes to sequence the target where all of the nucleic acid probes that are perfectly complementary to part of the target may not be readily identified.
  • The present invention identifies neighbor-rich probes which are then utilized to sequence the target nucleic acid. A “neighbor-rich probe” is a probe that is related to many other probes in the probe space by a single base mismatch. A probe that has a single base mismatch with another probe will be referred to as a “single base mismatch neighbor.” Neighbor-rich probes may be identified according to mismatch information as follows.
  • After a set of probes is identified, each probe in the set is compared to the other probes to determine how the probe's sequence compares to the other probes. In the example above, one probe differs from the other 24 probes by a single base mismatch (i.e., m=1, where m is the number of mismatches). Thus, this one probe is related to or has 24 single base mismatch neighbors.
  • By contrast, twenty-four probes differ from 3 other probes in the set by a single base mismatch and from 21 other probes in the set by a double base mismatch (i.e., m=2). In this simple example, the perfectly complementary probe may be identified as a neighbor-rich probe from the mismatch information because it has many single mismatch neighbor probes in the probe space. The perfectly complementary probe had 8 times as many single mismatch neighbor probes as nearly-complementary probes. Although the hybridization conditions did not allow identification of the perfectly complementary probes, an analysis of mismatch information may be utilized to identify the perfectly complementary probe. In practice, mismatch information may be utilized for de novo sequencing of a target nucleic acid where oligonucleotide probes are contacted with the target under conditions that do not allow differentiation between those probes that are perfectly complementary to part of the target and those probes that are not.
  • In this example, the sequence of the target was known. However, in many applications including de novo sequencing the sequence of the target is unknown. Nevertheless, the example is useful in demonstrating how neighbor-rich probes may be identified.
  • A. One Embodiment
  • FIG. 6 shows a high level flow of sequencing utilizing mismatch information. At step 500, hybridization intensities from probes or other data indicative of binding affinity are input into the system. The system may receive the hybridization intensities many different ways. The system may operate the scanning device directly, the system may receive the hybridization intensities from another computer system that measured the intensities, or an operator may manually enter the data. There may be thousands or tens of thousands of hybridization intensities that correspond to nucleic acid probes on a chip. Typically, the chip includes all possible probes of a specific length in order to sequence the target.
  • At step 502, the system selects a set of probes associated with the highest hybridization intensities (i.e., that show the strongest binding affinity). Selecting the probes with the highest hybridization intensities may be done in any number of ways. For example, the system may use an intensity threshold value and select the probes whose hybridization intensities are higher than the intensity threshold (e.g., 100 photon counts). The system may select a specific number or percentage of probes (e.g., 50 probes or the top 10%) that have the highest hybridization intensities. Additionally, the system may select the probes that have a hybridization intensity greater than a specific percentage (e.g., 40%) of the highest hybridization intensity.
  • After the set of probes with the strongest binding affinity is selected, the system calculates the number of single base mismatch neighbors for each probe in the set at step 504. For example, in one embodiment, a probe is first selected in order to compare the selected probe to the other probes. The system then determines how many of the other probes in the set are identical to the selected probe except for a single base mismatch at one base position. The number of single base mismatch neighbors is calculated for each of the probes having the highest hybridization intensities. Additionally, the system may calculate and utilize the number of double base mismatches in an extension of the concepts herein.
  • At step 506, the system selects the probes in the set with the highest number of single base mismatch neighbors. Selecting the probes with the highest number of neighbors may be done in any number of ways including utilizing a threshold, a specific number of the probes, or greater than a specific percentage of the highest number of neighbors. In one embodiment, the system selects the probes in the set with the highest number of neighbors and the highest hybridization intensities. In other words, a second intensity threshold is utilized to further reduce the set of probes.
  • The selected probes with the highest number of single base mismatch neighbors are the neighbor-rich probes. The neighbor-rich probes are then aligned at step 508. The neighbor-rich probes are aligned or compiled so that they have the most bases in common. Thus, neighbor-rich probes that have a single base mismatch are aligned to form an aligned set of probes. Aligned sets of probes are then aligned in an offset fashion in the way that reduces the number of base mismatches between sets of probes. At step 510, the aligned probes are utilized to sequence the target nucleic acid sequence. The target may be sequenced in many different ways including the formation of a consensus sequence may be produced as described in the following example.
  • B. Example
  • A target of 5′-AGTTGTAGTGGATGG was exposed to a chip containing 8-mer probes. The highest hybridization intensity was 331 photon counts. An intensity threshold of 90 photon counts was utilized and there were 133 probes that had a hybridization intensity greater than the intensity threshold of 90. These 133 probes provided the set of probes with the highest hybridization intensities and are as follows:
    Probe intensity m = 1 m = 2
    3′-ACATCACC 331 12 10
    3′-CATCACCT 286 11 12
    3′-ATCACCTA 323 10 8
    3′-CATCACCA 253 8 12
    3′-ACAACATC 331 7 10
    3′-AATCACCT 131 7 14
    3′-ACATCACA 330 7 10
    3′-ACCTACCA 280 7 6
    3′-CACCTACC 204 7 2
    3′-ACTCACCT 188 7 11
    3′-CCATCACC 270 6 13
    3′-TTCACCTA 134 6 10
    3′-ACACCACC 98 6 14
    3′-TCAACATC 331 6 10
    3′-TCATCACC 238 5 17
    3′-CTCACCTA 203 5 10
    3′-ACACCAAC 122 5 10
    3′-ACACCTAC 272 5 8
    3′-ACCTACCC 108 5 7
    3′-TCATCACA 147 5 7
    3′-CAACATCA 275 5 5
    3′-CAACACCT 183 5 15
    3′-CACCACCT 113 5 15
    3′-ATCACCAC 157 5 6
    3′-TCCACCTA 112 5 7
    3′-TCACCTAC 248 5 6
    3′-TGCACCTA 105 5 6
    3′-TATCACCT 143 5 13
    3′-CCACCTAC 208 5 5
    3′-ACAACACC 147 5 18
    3′-CCAACATC 325 4 10
    3′-GCATCACC 262 4 15
    3′-GCACCTAC 199 4 6
    3′-AACATCAC 148 4 1
    3′-AGTCACCT 90 4 9
    3′-CAACATCT 101 4 5
    3′-CCATCACA 181 4 9
    3′-CGCACCTA 127 4 6
    3′-TCAACACC 129 4 13
    3′-ATCACCTT 129 4 9
    3′-ACACCTAA 155 4 6
    3′-CAACACCA 100 4 14
    3′-ACACACCT 305 4 11
    3′-TCACCTAA 176 4 4
    3′-ACACACCA 174 4 8
    3′-CAGCACCT 111 4 13
    3′-ATCACCAA 115 4 11
    3′-ATCACCTC 137 4 13
    3′-GCATCACA 156 4 7
    3′-TACACCTA 96 4 7
    3′-CCTCACCT 102 4 14
    3′-TCAACCTC 132 4 6
    3′-CACCACCA 91 4 13
    3′-CATCACCC 131 4 17
    3′-GCAACATC 319 4 9
    3′-CATCAACC 105 4 6
    3′-CACCTACA 187 4 5
    3′-ACACCATC 128 4 13
    3′-ACCTACCT 112 4 11
    3′-ATTCACCT 91 4 10
    3′-CCACCTAA 111 4 5
    3′-GCACCTAA 127 4 3
    3′-AGCACCTA 148 4 11
    3′-GTCACCTA 141 3 12
    3′-ACATCACT 164 3 13
    3′-CATCACCG 163 3 15
    3′-CCCTACCA 133 3 6
    3′-ACCTACCG 119 3 8
    3′-ACAGCACC 101 3 14
    3′-ATCACCCA 106 3 12
    3′-CGTCACCT 114 3 12
    3′-CAACATCC 148 3 7
    3′-ACAACCTC 114 3 10
    3′-ATCAACCT 120 3 7
    3′-ACCAACCA 104 3 12
    3′-GCCTACCA 111 3 4
    3′-CACCAACA 119 3 8
    3′-ACTCACCA 143 3 12
    3′-ACCACCTA 141 3 17
    3′-CTATCACC 100 3 4
    3′-CAACATCG 137 3 4
    3′-ACGCACCT 110 3 11
    3′-TCACCATC 102 3 10
    3′-CACACCTA 102 3 7
    3′-CACCTACT 106 3 6
    3′-CACCAACC 103 3 11
    3′-ATCACCTG 106 3 9
    3′-ACATCACG 149 3 13
    3′-GCAACCTC 93 3 4
    3′-AAGCACCT 92 3 11
    3′-ATCATCAC 90 3 6
    3′-TCCTACCA 91 3 4
    3′-ATCAACCA 103 3 5
    3′-ACCAACCT 97 3 10
    3′-GATCACCT 102 3 14
    3′-TCACCAAC 102 3 10
    3′-ACCTACTC 102 3 3
    3′-CACCTACG 99 3 5
    3′-CCTCACCA 91 3 12
    3′-ATCACCAT 125 3 8
    3′-TCAACCTA 104 2 8
    3′-ACCATCAC 113 2 5
    3′-CATCTACC 94 2 8
    3′-CAATCACC 94 2 6
    3′-ACATCAAC 154 2 15
    3′-ACCTACAC 113 2 5
    3′-ACACATCA 128 2 6
    3′-CCACATCA 90 2 7
    3′-TATCACCA 97 2 12
    3′-CACATCAC 154 2 3
    3′-TCAACACA 97 2 8
    3′-TACCACCT 91 2 9
    3′-ATCCACCT 105 2 11
    3′-ACACACCG 122 2 5
    3′-ACACCACA 90 2 14
    3′-ATATCACC 96 2 12
    3′-TACATCAC 128 2 2
    3′-CAACCTAC 116 1 6
    3′-CATCACAA 107 1 7
    3′-ACCTCACC 102 1 13
    3′-ACCAACTC 96 1 9
    3′-TATCAACC 94 1 6
    3′-TACCTACC 99 1 8
    3′-ACCACATC 128 1 12
    3′-ATCACAAC 153 1 5
    3′-CCTACATC 93 1 4
    3′-CACCTAAC 95 1 7
    3′-CCTACCAA 128 0 0
    3′-TACACACC 91 0 2
    3′-CAACCATC 93 0 5
    3′-GTTAAGAG 329 0 0
    3′-AGCAACAT 94 0 3
    3′-TCTATGCG 33 0 0

    where the columns denoted m=1 and m=2 indicate the number of single and double base mismatch neighbors, respectively. Thus, each probe was compared to the other 132 probes to determine the number of single and double base mismatches the probe had with the other probes. The highest number of single base mismatch neighbor probes was 12 and the probes are presented in decreasing order according to the number of single base mismatch neighbor probes.
  • A set of neighbor-rich probes was identified by selecting the probes that had a hybridization intensity greater than 40% of 331 (0.40*331=132.4), and the number of single base mismatch neighbors greater than 40% of 12 (0.40*12=4.8). The following is the list of neighbor-rich probes selected in this manner:
    Probe intensity m = 1 m = 2
    3′-ACACCTAC 272 5 8
    3′-ACATCACC 331 12 10
    3′-ACATCACA 330 7 10
    3′-ACAACACC 147 5 18
    3′-ACAACATC 331 7 10
    3′-ACCTACCA 280 7 6
    3′-ATCACCAC 157 5 6
    3′-ATCACCTA 323 10 8
    3′-ACTCACCT 188 7 11
    3′-CAACATCA 275 5 5
    3′-CAACACCT 183 5 15
    3′-CCACCTAC 208 5 5
    3′-CCATCACC 270 6 13
    3′-CTCACCTA 203 5 10
    3′-CATCACCT 286 11 12
    3′-CATCACCA 253 8 12
    3′-TCACCTAC 248 5 6
    3′-TCATCACC 238 5 17
    3′-TCATCACA 147 5 7
    3′-TCAACATC 331 6 10
    3′-TTCACCTA 134 6 10
    3′-TATCACCT 143 5 13

    where again m=1 and m=2 indicates the number of single and double base mismatch neighbors, respectively, with other probes in the set of probes with a hybridization intensity greater than 90.
  • Once the neighbor-rich probes having a high hybridization intensity and high number of single base mismatch neighbors have been selected, the neighbor-rich probes were utilized to sequence the target nucleic acid sequence. The system utilized the frequency of bases at each position to produce a consensus sequence, where a “consensus sequence” is a sequence generated by neighbor-rich probes to sequence the target.
  • In order to produce a consensus sequence, the system aligned the neighbor-rich probes so that each probe had the highest number of bases in common with other probes. The following are the aligned neighbor-rich probes with the complement of the target sequence shown for reference. The target sequence is known in this example but the target sequence may be an unknown sequence or only partially known.
          A C A C C T A C
          A C A T C A C C
          A C A T C A C A
    A C A A C A C C
    A C A A C A T C
            A C C T A C C A
              A T C A C C A C
              A T C A C C T A
            A C T C A C C T
      C A A C A T C A
            C A A C A C C T
            C C A C C T A C
          C C A T C A C C
          C T C A C C T A
              C A T C A C C T
            C A T C A C C A
                T C A C C T A C
          T C A T C A C C
          T C A T C C A
    T C A A C A T C
              T T C A C C T A
            T A T C A C C T
    Target T C A A C A T C A C C T A C C
    complement-
  • After the neighbor-rich probes are aligned, the system counts the number or frequency of each base (A, C, G and T) at each position. After the frequency of bases at each position is calculated, the system produces a consensus sequence. In one embodiment, the base that occurred the most at a position is utilized to produce the consensus sequence if the base occurred more than 2 times and the frequency that the base occurred is greater than 50% at that position. The following is a matrix of base vs. frequency that was used to produce the consensus sequence in this manner:
    Frequency
    Base A
    2 0 5 8 2 17 2 0 17 2 2 3 5 0 0
    C 0 5 0 1 15 4 3 22 2 16 16 1 1 2 1
    G 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    T 1 0 0 2 1 1 18 1 1 1 0 9 0 0 0
    Consensus C A A C A T C A C C T A

    Therefore, for the first position in the consensus sequence (left-most in the matrix), base A occurred 66% (2 divided by 3) of the time which is greater than 50%, however, the base did not occur more than 2 times so the base is called as ambiguous (i.e., “−”).
  • The consensus sequence is the complement of the target; thus, the target is sequenced according to the complement of the consensus sequence. In this example, the target sequence was 5′-AGTTGTAGTGGATGG and it was correctly sequenced as 5′-GTTGTAGTGGAT (the terminal bases being ambiguous). The parameters for producing the consensus sequence may be varied according to the experimental data. For example, if the consensus sequence was formed solely by the bases that occur most often at each position, the consensus sequence would be perfectly complementary to the target nucleic acid for this data. However, this will not always be the case.
  • Software Appendix A provides a BASIC source code illustration of this embodiment of the invention. The source code is written in Quick BASIC for an IBM compatible personal computer.
  • C. Alternate Embodiment
  • FIG. 7 shows a high level flow of another embodiment of sequencing utilizing mismatch information. At step 550, hybridization intensities are input into the system. The system may receive the hybridization intensities many different ways. The system may operate the scanning device directly, the system may receive the hybridization intensities from another computer system that measured the intensities, or an operator may manually enter the data.
  • At step 552, the system selects a set of probes associated with the highest hybridization intensities. Selecting the probes with the highest hybridization intensities may be done in any number of ways. For example, the system may use an intensity threshold value and select the probes whose hybridization intensities are higher than the intensity threshold. The system may select a specific number or percentage of probes that have the highest hybridization intensities. Alternatively, the system may select the probes that have a hybridization greater than a specific percentage of the highest hybridization intensity.
  • After the set of probes with the highest hybridization intensities is selected, the system calculates the number of straight mismatches for each probe in the set at step 554. “Straight mismatches” refers to base mismatches between probes where the bases at corresponding positions are compared (e.g., as was done in the previous embodiment). For example, a probe 3′-AACAT is compared to a probe 3′-AACTT by aligning the probes as follows:
    3′-AACAT
    3′-AACTT

    Each probe has n bases, where n=5. If the base at the 3′ end is at position 1, there is a single mismatch at the fourth position where the A and T do not match. Thus, straight mismatches are determined by comparing bases at the same position in each probe (i.e., z=0, where z indicates the number of bases one of the probes was offset from the other probe before comparing).
  • In one embodiment, a matrix is formed to show the straight mismatches between probes. FIG. 8 shows a straight mismatch matrix for 50 probes. For simplicity, each probe is assigned an identification number from 1 to 50. The numbers on the rows and columns of the matrix correspond to the identification number of the probe. The value in the matrix represents the number of straight mismatches between the probes designated by the row and column. If there are more than 2 mismatches, the matrix shows a “.” at the appropriate matrix location. Since a diagonal 600 of the matrix shows the number of straight mismatches between the same probe, the diagonal contains 0's because the probe is being compared to itself (i.e., there are no straight mismatches). Also, the matrix is identical on each side of diagonal 600.
  • The system calculates the number of offset mismatches for each probe in the set at step 556. “Offset mismatches” refers to base mismatches between probes where the bases at offset positions are compared. For example, a probe 3′-CGAATCAA is compared to a probe 3′-GCATCAAC by aligning the probes as follows:
    3′-CGAATCAA
    3′-GCATCCAC

    Each probe has n bases, where n=8. If the base at the 3′ end is at position 1, bases at position 1 through 7 (or n−1) of the first probe is compared to bases at position 2 through 8 (or n). As shown, there are two mismatches (double mismatch) when the probes are offset a single base position (i.e., z=1, where z indicates the number of bases one of the probes is offset from the other probe before comparing).
  • In one embodiment, a matrix is formed to show the offset mismatches between probes. FIG. 9 shows an offset mismatch matrix for 50 probes. As in FIG. 8, the rows and columns of the matrix correspond to the identification number of the probe, which is 1 to 50. The value in the matrix represents the number of offset mismatches between the probes designated by the row and column. If there are more than 2 mismatches, the matrix shows a “.” at the appropriate matrix location. As shown, a diagonal will not contain O's and the matrix is not identical on each side of the diagonal. Although the probes were offset a single base position, the probes may be offset more positions when they are compared in an extension of the principles herein.
  • At step 558, the system selects the probes with less than some small number straight mismatches and less than some small number of offset mismatches. In one embodiment, the system identifies the probes in the matrices that have less than 2 straight mismatches and less than 3 offset mismatches. The parameters for selecting these probes with few mismatches may be varied according to the experimental data.
  • The selected probes are then aligned at step 560. The probes are aligned so that they have the most bases in common. The mismatch information concerning the straight mismatches and offset mismatches is utilized to align the probes so that the number of mismatches between the probes is reduced. At step 562, the aligned probes are utilized to sequence the target nucleic acid sequence. The target may be sequenced in many different ways. For example, a consensus sequence may be produced as described in the following example.
  • D. Example
  • A target of 5′-AGTTGTAGTGGATGGT was exposed to a chip containing 10-mer probes. Fifty probes were selected that have the highest hybridization intensities (step 552). FIGS. 11 and 12 show the straight and offset mismatch matrices for the fifty probes (steps 554 and 556). Forty-seven probes were selected that have less than 2 straight mismatches with at least one other probe and less than 3 offset mismatches with at least one other probe (step 558).
  • The straight and offset mismatch information was utilized to align the 47 probes (step 560). For example, FIG. 9 shows that the probe identified as 2 on the row had 0 offset mismatches with the probe identified as 1 on the column. Therefore, probes 2 and 1 align well if they are offset a single base position. The following are the aligned 47 probes:
      AACATCACCT
     CAACATCACC
       ACATCACCTA
    ACAACATCAC
     CAACATCACA
         ATCACCTACC
      AACATCACCA
      CACATCACCT
      AACATCACCG
    TCAACATCAC
        CATCACCTAC
    CCAACATCAC
          ACACCTACCA
     CAACATCACG
    GCAACATCAC
       ACATCACCTT
      AACATCACCC
         AGCACCTACC
     AAACATCACC
      ACCATCACCT
       ACATCACCAT
           CACCTACCAA
       ACATCACCTC
         ATCACCTACA
       ACATCACCTG
     ACACATCACC
     GAACATCACC
      TACATCACCT
     TAACATCACC
           CACCTACCAG
       GACATCACCT
       CACATCACCA
       ATCATCACCT
             ACCTACCATC
     ACAACATCAA
      CAACATCACT
        ACATCACCAA
       CACATCACCG
        ACATCACCCT
            CACCTACCAC
         CATCACCTAA
           TCACCTACCA
            CACCTACCAT
        CCATCACCTA
        ACATCACCCA
        ACATCACCGA
    ATCAACATCA
  • After the selected probes were aligned, the system counts the number of occurrences of each base (A, C, G and T) at each position. After the frequency of bases at each position is calculated, the system produces a consensus sequence which should be complementary to the target sequence. If the system utilized bases that occurred more than 2 times and the frequency that the base occurred is greater than 50% at that position, the consensus sequence 3′-CAACATCACCTACCA is produced.
  • The consensus sequence is ideally the complement of the target; thus, the target is sequenced according to the complement of the consensus sequence (consensus' where the prime denotes the complement). In this example, the target and consensus' sequence were as follows:
    Target AGTTGTAGTGGATGGT
    Consensus′ GTTGTAGTGGATGGT

    (one terminal base of the consensus sequence being ambiguous). Thus, the target was sequenced with a high degree of accuracy utilizing mismatch information. The parameters for producing the consensus sequence may be varied according to the experimental data.
  • Software Appendix B provides a BASIC source code illustration of this embodiment of the invention. The source code is written in Quick BASIC for an IBM compatible personal computer.
  • CONCLUSION
  • The above description is illustrative and not restrictive. Many variations of the invention will become apparent to those of skill in the art upon review of this disclosure. Merely by way of example, while the invention is illustrated with particular reference to the evaluation of DNA (natural or unnatural), the methods can be used in the analysis from chips with other materials synthesized thereon, such as RNA. The scope of the invention should, therefore, be determined not with reference to the above description, but instead should be determined with reference to the appended claims along with their full scope of equivalents.
  • APPENDIX A
  • ’SEARCH.BAS 11/21/94
    ’Finds pixels with intensities greater than fmin
    ’Calculates the match score of bright pixels
    ’The input file eightmer.dat is based on fs8mer.exe
    ’The output file score.dat lists the bright pixels
    ’  and gives the number of 1 and 2 mismatch-related
    ’  pixels
    DIM a$(1000), f%(1000), m1%(1000), m2%(1000)
    tstart = TIMER
    inf$ = “eightmer.dat”
    outf$ = “score.dat”
    OPEN inf$ FOR INPUT AS #1
    OPEN outf$ FOR OUTPUT AS #2
    CLS
    ’Read the input file and store the bright pixels
    fmin = 90 ’threshold value for inclusion
    n = 0: fmax = 0
    WHILE NOT EOF(1)
    LINE INPUT #1, g$
    seq$ = MID$(g$, 1, 8)
    intens = VAL(MID$(g$, 9, 6))
    IF intens >= 90 THEN
       n = n + 1
       a$(n) = seq$
       f%(n) = intens
       IF intens > fmax THEN fmax = intens
       PRINT n;
    END IF
    WEND
    PRINT
    PRINT USING “#### intensity values above ####”; n; fmin
    PRINT USING “Highest intensity is ####”; fmax
    ’Calculate m1 and m2 for each bright pixel
    ’  m1 and m2 are the number of other pixels that are related
    ’  by 1 and 2 mismatches, respectively
    m1max = 0 ’Keep track of highest m1 score
    FOR j = 1 TO n
    PRINT j;
    FOR i = 1 TO n
    m = 0
    FOR k = 1 TO 8
    IF MID$(a$(j), k, l) <> MID$(a$(i), k, l) THEN m = m + 1
    NEXT k
    IF m = 1 THEN m1%(j) = m1%(j) + 1
    IF m = 2 THEN m2%(j) = m2%(j) + 1
    NEXT i
    IF m1%(j) > m1max THEN m1max = m1%(j)
    NEXT j
    PRINT #2, USING “SEARCH.BAS  &  &”; DATE$; TIME$
    PRINT #2, USING “Input file: &   Output file: &”; inf$; outf$
    PRINT #2, USING “#### intensity values above ####”; n; fmin
    PRINT #2, USING “Highest intensity is ####”; fmax
    PRINT #2, USING “Greatest number of 1-mismatch relations
    is ##”; m1max
    PRINT #2,
    PRINT #2, “List of probes with highest intensity and best matching”
    PRINT #2, “ f m1 m2 sequence”
    Figure US20060194224A1-20060831-P00801
    k = 1 TO n
    IF f%(k) > .4 * fmax AND m1%(k) > .4 * m1max THEN
       PRINT #2, USING “#### ### ### &”;
       f%(k); m1%(k); m2%(k); a$(k)
    END IF
    NEXT k
    PRINT #2, CHR$(12)
    ’Sort according to f
    s% = n \ 2
    DO WHILE s% > 0
       FOR i% = s% TO n − 1
       j% = i% − s% + 1
       FOR j% = (i% − s% + 1) TO 1 STEP −s%
          IF f%(j%) >= f%(j% + s%) THEN EXIT FOR
          SWAP f%(j%), f%(j% + s%)
          SWAP m1%(j%), m1%(j% + s%)
          SWAP m2%(j%), m2%(j% + s%)
          SWAP a$(j%), a$(j% + s%)
       NEXT j%
       NEXT i%
       s% = s% \ 2
    LOOP
    PRINT #2,
    PRINT #2, “ f m1 m2 sequence”
    FOR k = 1 TO n
    PRINT #2, USING “#### ### ### &”; f%(k); m1%(k); m2%(k); a$(k)
    NEXT k
    PRINT CHR$(12)
    ’Sort according to m1
    s% = n \ 2
    DO WHILE s% > 0
       FOR i% = s% TO n − 1
       j% = i% − s% + 1
       FOR j% = (i% − s% + 1) TO 1 STEP −s%
          IF m1%(j%) >= m1%(j% + s%) THEN EXIT FOR
          SWAP f%(j%), f%(j% + s%)
          SWAP m1%(j%), m1%(j% + s%)
          SWAP m2%(j%), m2%(j% + s%)
          SWAP a$(j%), a$(j% + s%)
       NEXT j%
       NEXT i%
       s% = s% \ 2
    LOOP
    PRINT #2,
    PRINT #2, “ f m1 m2 sequence”
    FOR k = 1 TO n
    PRINT #2, USING “#### ### ### &”; f%(k); m1%(k); m2%(k); a$(k)
    NEXT
    PRINT USING “Time= ####.# seconds”; TIMER − tstart
  • APPENDIX B
  • ‘CONSENS.BAS 1/8/95
    ‘Derive a consensus sequence from the highest scoring probes
    DIM a$(70), m%(1, 70, 70), f(70), s(−20 TO 20, 4)
    CLS
    INPUT “Input file: ”, inf$
    INPUT “Output file: ”, outf$
    OPEN inf$ FOR INPUT AS #1
    OPEN outf$ FOR OUTPUT AS #2
    LINE INPUT #1, descr$ ‘File description
    INPUT #1, pl ‘Probe length
    INPUT #1, n ‘Number of sequences
    FOR j = 1 TO n
    LINE INPUT #1, a$(j)
    NEXT j
    CLOSE #1
    ‘Initialize the mismatch matrix
    FOR z = 0 TO 1: FOR i = 1 TO n: FOR j = 1 TO n
    m%(z, i, j) = 100
    NEXT j: NEXT i: NEXT z
    PRINT #2,
    PRINT #2, “CONSENS.BAS ”; DATE$; “ ”; TIME$
    PRINT #2, : PRINT #2,
    PRINT #2, “Input file: ”; inf$; “  Output file: ”; outf$
    PRINT #2, descr$
    PRINT #2, USING “The ##  ##-mer sequences with the
    highest scores are:”;
    n; pl
    PRINT #2,
    FOR j = 1 TO n
    PRINT #2, USING “##  &”; j; a$(j)
    NEXT j
    PRINT #2, : PRINT #2,
    z = 0
    PRINT #2, USING “z=##”; z
    PRINT #2, “ ”;
    FOR k = 1 TO n: PRINT #2, USING “##”; k; : NEXT k
    FOR i = 1 TO n
    PRINT #2,
    PRINT #2, USING “##  ”; i;
    FOR j = 1 TO n
       m = 0
       FOR k = 1 TO pl
       IF MID$(a$(j), k, l) <> MIDS(a$(i), k, l) THEN m = m + 1
       NEXT k
       m%(0, i, j) = m
       IF m <= 2 THEN PRINT #2, USING “ #”; m; ELSE
       PRINT #2, “ .”;
    NEXT j
    NEXT i
    PRINT #2, : PRINT #2,
    z = 1
    PRINT #2, USING “z=##”;
    PRINT #2, “ ”;
    FOR k = 1 TO n: PRINT #2, USING “##”; k; : NEXT k
    FOR i = 1 TO n
    PRINT #2,
    PRINT #2, USING “##  ”; i;
    FOR j = 1 TO n
       m = 0
       FOR k = 1 TO pl − 1
       IF MID$(a$(j), k, 1) <> MID$(a$(i), k + 1, 1) THEN m = m + 1
       NEXT k
       m%(1, i, j) = m
       IF m <= 2 THEN PRINT #2, USING “ #”; m; ELSE
       PRINT #2, “ .”;
    NEXT j
    NEXT i
    PRINT #2, : PRINT #2,
    ‘Mark all sequences with a 100 tag
    FOR i = 1 TO n: f(i) = 100: NEXT i
    ‘Designate the first sequence as the origin
    f(1) = 0
    ‘Find the frames of sequences that can be aligned
    FOR i = 1 TO n
    FOR j = 1 TO n
    IF m%(1, i, j) <= 2 AND f(i) <> 100 THEN
       f(j) = f(i) + 1
    END IF
    NEXT j
    NEXT i
    FOR i = 1 TO n
    FOR j = 1 TO n
    IF m%(l, j, i) <= 2 AND f(i) <> 100 THEN
       f(j) = f(i) − 1
    END IF
    NEXT j
    NEXT i
    FOR i = 1 TO n
    FOR j = i + 1 TO n
    IF m%(0, i, j) <= 1 AND f(i) <> 100 THEN
       f(j) = f(i)
    END IF
    NEXT j
    NEXT i
    PRINT #2, : PRINT #2,
    PRINT #2, “Alignment criteria: <= 1 mismatch allowed for z=0”
    PRINT #2, “  <=2 mismatches for z=1”
    PRINT #2,
    PRINT #2, “The aligned sequences are:”
    ‘Print the aligned sequences
    FOR i = 1 TO n
    IF f(i) <> 100 THEN
       PRINT #2, SPACE$(15 + f(i)); a$(i)
    END IF
    NEXT i
    PRINT #2, : PRINT #2,
    ‘Accumulate the sequence scores
    offset = 0
    FOR i = 1 TO n
    IF f(i) <> 100 THEN
       FOR k = 1 TO pl
       g = INSTR(“ACGT”, MID$(a$(i), k, 1))
       s(offset + k + f(i), g) = s(offset + k + f(i), g) + 1
       ‘PRINT offset + k + f(i); g; “   ”;
       NEXT k
    END IF
    NEXT i
    PRINT #2, CHR$(12)
    PRINT #2, “CONSENS.BAS ”; DATE$; “ ”; TIME$
    PRINT #2, USING “Input file: &  Output file: &”; inf$; outf$
    PRINT #2, USING “### ##mer sequences”; n; pl
    PRINT #2, descr$
    PRINT #2,
    PRINT #2, “The frequencies of bases in the aligned sequences are:”
    PRINT #2,
    ‘Print the scores
    FOR g = 1 TO 4
    FOR j = −10 TO 18
    PRINT #2, USING “## ”; s(j, g);
    ‘PRINT USING “## ”; s(j, g);
    NEXT j
    PRINT #2,
    NEXT g
    ‘Find and print the, consensus
    c$(0) = “−”: c$(1) = “A”: c$(2) = “C”: c$(3) = “G”: c$(4) = “T”
    FOR j = −10 TO 18
    most 0: mg = 0: sum = 0: b$ = “−”
    FOR g = 1 TO 4
    IF s(j, g) > most THEN most = s(j, g): mg = g
    sum = sum + s(j, g)
    NEXT g
    ‘A base is defined if present in at least 2 sequences
    ‘  and 55% of those aligned at that position
    IF most >= 3 THEN
       IF most / sum > .5 THEN b$ = c$(mg)
    END IF
    PRINT #2, USING “ & ”; b$;
    cons$ = cons$ + b$
    NEXT j
    PRINT #2, : PRINT #2, : PRINT #2, “The consensus sequence
    is: ”; cons$:
    PRINT #2,
    PRINT cons$
    PRINT #2, : PRINT #2,
    PRINT #2, “The correct sequence is TCAACATCACCTACCA”
    PRINT #2,
    PRINT #2, “The stray sequences are:”
    FOR i = 1 TO n
    IF f(i) = 100 THEN PRINT #2, SPACE$(5); a$(i)
    NEXT i

Claims (29)

1. In a computer system, a method of sequencing a nucleic acid sequence, the method comprising the steps of:
inputting hybridization intensities for a plurality of nucleic acid probes, the nucleic acid probes hybridizing with the nucleic acid sequence under conditions that do not allow identification of all of nucleic acid probes that are perfectly complementary to part of the nucleic acid sequence; and
sequencing the nucleic acid sequence according to selected nucleic acid probes.
2. The method of claim 1, further comprising the step of selecting nucleic acid probes utilizing mismatch information among the nucleic acid probes.
3. The method of claim 1, further comprising the step of selecting nucleic acid probes with highest numbers of single base mismatch neighbors among the probes, a single base mismatch neighbor being another probe that has the same sequence except for a single base that is different.
4. The method of claim 1, further comprising the step of selecting nucleic acid probes that have fewer than a predetermined number of base mismatches with another probe.
5. The method of claim 1, further comprising the step of selecting nucleic acid probes with highest hybridization intensities.
6. The method of claim 5, wherein the nucleic acid probes with the highest hybridization intensities are above an intensity threshold.
7. The method of claim 1, wherein the sequencing step includes the steps of:
aligning the selected nucleic acid probes so that bases that are common overlap; and
sequencing the nucleic acid sequence according to bases that occur most often at each position in the aligned nucleic acid probes.
8. The method of claim 7, wherein the nucleic acid sequence is sequenced as complementary to the bases that occur most often.
9. A computer program product that sequences a nucleic acid sequence, comprising:
code that receives as input hybridization intensities for a plurality of nucleic acid probes, the nucleic acid probes hybridizing with the nucleic acid sequence under conditions that do not allow identification of all nucleic acid probes that are perfectly complementary to part of the nucleic acid sequence;
code that sequences the nucleic acid sequence according to selected nucleic acid probes; and
a computer readable storage medium that stores the codes.
10. In a computer system, a method of sequencing a nucleic acid sequence, the method comprising the steps of:
inputting hybridization intensities for a plurality of nucleic acid probes;
selecting nucleic acid probes with highest numbers of single base mismatch neighbors among the probes, a single base mismatch neighbor being another probe that has the same sequence except for a single base that is different; and
sequencing the nucleic acid sequence according to the selected nucleic acid probes.
11. The method of claim 10, further comprising the step of selecting nucleic acid probes with highest hybridization intensities.
12. The method of claim 11, wherein the nucleic acid probes with the highest hybridization intensities are above an intensity threshold.
13. The method of claim 10, wherein the step of selecting includes the step of calculating a number of single base mismatch neighbors for each nucleic acid probe.
14. The method of claim 10, wherein the step of selecting includes the step of utilizing double base mismatch neighbors, a double mismatch neighbor being another probe that has the same sequence except for two bases that are different.
15. The method of claim 10, wherein the sequencing step includes the steps of:
aligning the selected nucleic acid probes so that bases that are common overlap; and
sequencing the nucleic acid sequence according to bases that occur most often at each position in the aligned nucleic acid probes.
16. The method of claim 15, wherein the nucleic acid sequence is sequenced as complementary to the bases that occur most often.
17. A computer program product that sequences a nucleic acid sequence, comprising:
code that receives as input hybridization intensities for a plurality of nucleic acid probes;
code that selects nucleic acid probes with highest numbers of single base mismatch neighbors among the probes, a single base mismatch neighbor being another probe that has the same sequence except for a single base that is different;
code that sequences the nucleic acid sequence according to the selected nucleic acid probes; and
a computer readable storage medium that stores the codes.
18. In a computer system, a method of sequencing a nucleic acid sequence, the method comprising the steps of:
inputting hybridization intensities for a plurality of nucleic acid probes;
selecting nucleic acid probes that have fewer than a predetermined number of base mismatches with another probe; and
sequencing the nucleic acid sequence according to the selected nucleic acid probes.
19. The method of claim 18, further comprising the step of selecting nucleic acid probes with highest hybridization intensities.
20. The method of claim 19, wherein the nucleic acid probes with the highest hybridization intensities are above an intensity threshold.
21. The method of claim 18, wherein each nucleic acid probe includes n bases.
22. The method of claim 21, further comprising the steps of:
comparing bases at positions 1 through n of a first nucleic acid probe to bases at positions 1 through n of a second nucleic acid probe; and
counting base mismatches between the first and second nucleic acid probe.
23. The method of claim 21, further comprising the steps of:
comparing bases at positions 1 through n−1 of a first nucleic acid probe to bases at positions 2 through n of a second nucleic acid probe; and
counting base mismatches between the first and second nucleic acid probe.
24. The method of claim 18, wherein the sequencing step includes the steps of:
aligning the selected nucleic acid probes so that bases that are common overlap; and
sequencing the nucleic acid sequence according to bases that occur most often at each position in the aligned nucleic acid probes.
25. The method of claim 24, wherein the nucleic acid sequence is sequenced as complementary to the bases that occur most often.
26. A computer program product that sequences a nucleic acid sequence, comprising:
code that receives as input hybridization intensities for a plurality of nucleic acid probes;
code that selects nucleic acid probes that have fewer than a predetermined number of base mismatches with another probe;
code that sequences the nucleic acid sequence according to the selected nucleic acid probes; and
a computer readable storage medium that stores the codes.
27. A method of sequencing a nucleic acid, comprising the steps of:
contacting a set of oligonucleotide probes of predetermined sequence and length with the nucleic acid under hybridization conditions that do not allow differentiation between (i) those probes of the set which are perfectly complementary to part of the nucleic acid and (ii) those probes that are not perfectly complementary to part of the nucleic acid;
selecting a subset of oligonucleotide probes that includes probes that are perfectly complementary to part of the nucleic acid and probes that are not perfectly complementary to part of the nucleic acid; and
determining the sequence of the nucleic acid by compiling overlapping sequences of the subset of probes.
28. The method of claim 27, further comprising the step of fragmenting the nucleic acid before contacting the nucleic acid with the set of oligonucleotide probes.
29. The method of claim 28, further comprising the step of labeling the nucleic acid before contacting the nucleic acid with the set of oligonucleotide probes.
US11/300,196 1996-03-19 2005-12-13 Computer-aided nucleic acid sequencing Abandoned US20060194224A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US11/300,196 US20060194224A1 (en) 1996-03-19 2005-12-13 Computer-aided nucleic acid sequencing

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US08/618,834 US6361937B1 (en) 1996-03-19 1996-03-19 Computer-aided nucleic acid sequencing
US1495301A 2001-12-11 2001-12-11
US11/300,196 US20060194224A1 (en) 1996-03-19 2005-12-13 Computer-aided nucleic acid sequencing

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US1495301A Continuation 1996-03-19 2001-12-11

Publications (1)

Publication Number Publication Date
US20060194224A1 true US20060194224A1 (en) 2006-08-31

Family

ID=24479320

Family Applications (2)

Application Number Title Priority Date Filing Date
US08/618,834 Expired - Lifetime US6361937B1 (en) 1996-03-19 1996-03-19 Computer-aided nucleic acid sequencing
US11/300,196 Abandoned US20060194224A1 (en) 1996-03-19 2005-12-13 Computer-aided nucleic acid sequencing

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US08/618,834 Expired - Lifetime US6361937B1 (en) 1996-03-19 1996-03-19 Computer-aided nucleic acid sequencing

Country Status (1)

Country Link
US (2) US6361937B1 (en)

Families Citing this family (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6361937B1 (en) * 1996-03-19 2002-03-26 Affymetrix, Incorporated Computer-aided nucleic acid sequencing
US6780591B2 (en) * 1998-05-01 2004-08-24 Arizona Board Of Regents Method of determining the nucleotide sequence of oligonucleotides and DNA molecules
US7875440B2 (en) 1998-05-01 2011-01-25 Arizona Board Of Regents Method of determining the nucleotide sequence of oligonucleotides and DNA molecules
US6142681A (en) 1999-02-22 2000-11-07 Vialogy Corporation Method and apparatus for interpreting hybridized bioelectronic DNA microarray patterns using self-scaling convergent reverberant dynamics
US20040111219A1 (en) * 1999-02-22 2004-06-10 Sandeep Gulati Active interferometric signal analysis in software
US6245511B1 (en) * 1999-02-22 2001-06-12 Vialogy Corp Method and apparatus for exponentially convergent therapy effectiveness monitoring using DNA microarray based viral load measurements
US6136541A (en) 1999-02-22 2000-10-24 Vialogy Corporation Method and apparatus for analyzing hybridized biochip patterns using resonance interactions employing quantum expressor functions
US7501245B2 (en) * 1999-06-28 2009-03-10 Helicos Biosciences Corp. Methods and apparatuses for analyzing polynucleotide sequences
US6818395B1 (en) * 1999-06-28 2004-11-16 California Institute Of Technology Methods and apparatus for analyzing polynucleotide sequences
US7013221B1 (en) * 1999-07-16 2006-03-14 Rosetta Inpharmatics Llc Iterative probe design and detailed expression profiling with flexible in-situ synthesis arrays
JP2004523243A (en) * 2001-03-12 2004-08-05 カリフォルニア インスティチュート オブ テクノロジー Method and apparatus for analyzing polynucleotide sequences by asynchronous base extension
WO2002083952A1 (en) * 2001-04-12 2002-10-24 Caliper Technologies Corp. Systems and methods for high throughput genetic analysis
EP1502098B1 (en) * 2002-05-03 2008-09-17 ViaLogy, LLC Method for characterizing microarray output data
US20050170367A1 (en) * 2003-06-10 2005-08-04 Quake Stephen R. Fluorescently labeled nucleoside triphosphates and analogs thereof for sequencing nucleic acids
US7169560B2 (en) 2003-11-12 2007-01-30 Helicos Biosciences Corporation Short cycle methods for sequencing polynucleotides
US20060172408A1 (en) * 2003-12-01 2006-08-03 Quake Steven R Device for immobilizing chemical and biochemical species and methods of using same
US7981604B2 (en) 2004-02-19 2011-07-19 California Institute Of Technology Methods and kits for analyzing polynucleotide sequences
US20060046258A1 (en) * 2004-02-27 2006-03-02 Lapidus Stanley N Applications of single molecule sequencing
US20050239085A1 (en) * 2004-04-23 2005-10-27 Buzby Philip R Methods for nucleic acid sequence determination
US20050260609A1 (en) * 2004-05-24 2005-11-24 Lapidus Stanley N Methods and devices for sequencing nucleic acids
US20070117103A1 (en) * 2005-11-22 2007-05-24 Buzby Philip R Nucleotide analogs
US7476734B2 (en) * 2005-12-06 2009-01-13 Helicos Biosciences Corporation Nucleotide analogs
US20070117104A1 (en) * 2005-11-22 2007-05-24 Buzby Philip R Nucleotide analogs
EP1766090B1 (en) * 2004-05-25 2011-04-27 Helicos Biosciences Corporation Methods for nucleic acid immobilisation
US20060024678A1 (en) * 2004-07-28 2006-02-02 Helicos Biosciences Corporation Use of single-stranded nucleic acid binding proteins in sequencing
US20060099612A1 (en) * 2004-09-02 2006-05-11 Suntory Limited Method for analyzing genes of industrial yeasts
EP1792263A2 (en) * 2004-09-02 2007-06-06 Vialogy Corporation Detecting events of interest using quantum resonance interferometry
US20060046253A1 (en) * 2004-09-02 2006-03-02 Suntory Limited Method for analyzing genes of industrial yeasts
US7682782B2 (en) 2004-10-29 2010-03-23 Affymetrix, Inc. System, method, and product for multiple wavelength detection using single source excitation
US20060118754A1 (en) * 2004-12-08 2006-06-08 Lapen Daniel C Stabilizing a polyelectrolyte multilayer
US7220549B2 (en) * 2004-12-30 2007-05-22 Helicos Biosciences Corporation Stabilizing a nucleic acid for nucleic acid sequencing
US20060172328A1 (en) * 2005-01-05 2006-08-03 Buzby Philip R Methods and compositions for correcting misincorporation in a nucleic acid synthesis reaction
US7482120B2 (en) * 2005-01-28 2009-01-27 Helicos Biosciences Corporation Methods and compositions for improving fidelity in a nucleic acid synthesis reaction
US20060263790A1 (en) * 2005-05-20 2006-11-23 Timothy Harris Methods for improving fidelity in a nucleic acid synthesis reaction
US7666593B2 (en) 2005-08-26 2010-02-23 Helicos Biosciences Corporation Single molecule sequencing of captured nucleic acids
US20070117102A1 (en) * 2005-11-22 2007-05-24 Buzby Philip R Nucleotide analogs
US20070128610A1 (en) * 2005-12-02 2007-06-07 Buzby Philip R Sample preparation method and apparatus for nucleic acid sequencing
US20090062138A1 (en) * 2007-08-31 2009-03-05 Curry Bo U Array-based method for performing SNP analysis
US20100034444A1 (en) * 2008-08-07 2010-02-11 Helicos Biosciences Corporation Image analysis

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4741043A (en) * 1985-11-04 1988-04-26 Cell Analysis Systems, Inc. Method of and an apparatus for image analyses of biological specimens
US4965725A (en) * 1988-04-08 1990-10-23 Nueromedical Systems, Inc. Neural network based automated cytological specimen classification system and method
US5002867A (en) * 1988-04-25 1991-03-26 Macevicz Stephen C Nucleic acid sequence determination by multiple mixed oligonucleotide probes
US5143854A (en) * 1989-06-07 1992-09-01 Affymax Technologies N.V. Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof
US5200313A (en) * 1983-08-05 1993-04-06 Miles Inc. Nucleic acid hybridization assay employing detectable anti-hybrid antibodies
US5202231A (en) * 1987-04-01 1993-04-13 Drmanac Radoje T Method of sequencing of genomes by hybridization of oligonucleotide probes
US5235626A (en) * 1991-10-22 1993-08-10 International Business Machines Corporation Segmented mask and exposure system for x-ray lithography
US5288514A (en) * 1992-09-14 1994-02-22 The Regents Of The University Of California Solid phase and combinatorial synthesis of benzodiazepine compounds on a solid support
US5384261A (en) * 1991-11-22 1995-01-24 Affymax Technologies N.V. Very large scale immobilized polymer synthesis using mechanically directed flow paths
US5470710A (en) * 1993-10-22 1995-11-28 University Of Utah Automated hybridization/imaging device for fluorescent multiplex DNA sequencing
US5503980A (en) * 1992-11-06 1996-04-02 Trustees Of Boston University Positional sequencing by hybridization
US5525464A (en) * 1987-04-01 1996-06-11 Hyseq, Inc. Method of sequencing by hybridization of oligonucleotide probes
US5527681A (en) * 1989-06-07 1996-06-18 Affymax Technologies N.V. Immobilized molecular synthesis of systematically substituted compounds
US5733729A (en) * 1995-09-14 1998-03-31 Affymetrix, Inc. Computer-aided probability base calling for arrays of nucleic acid probes on chips
US6361937B1 (en) * 1996-03-19 2002-03-26 Affymetrix, Incorporated Computer-aided nucleic acid sequencing

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8810400D0 (en) 1988-05-03 1988-06-08 Southern E Analysing polynucleotide sequences
DE69132905T2 (en) 1990-12-06 2002-08-01 Affymetrix Inc N D Ges D Staat Detection of nucleic acid sequences
CA2097708A1 (en) 1990-12-06 1992-06-07 Stephen P. A. Fodor Very large scale immobilized polymer synthesis
US5273632A (en) 1992-11-19 1993-12-28 University Of Utah Research Foundation Methods and apparatus for analysis of chromatographic migration patterns
AU8126694A (en) 1993-10-26 1995-05-22 Affymax Technologies N.V. Arrays of nucleic acid probes on biological chips
US5807522A (en) 1994-06-17 1998-09-15 The Board Of Trustees Of The Leland Stanford Junior University Methods for fabricating microarrays of biological samples

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5200313A (en) * 1983-08-05 1993-04-06 Miles Inc. Nucleic acid hybridization assay employing detectable anti-hybrid antibodies
US4741043B1 (en) * 1985-11-04 1994-08-09 Cell Analysis Systems Inc Method of and apparatus for image analyses of biological specimens
US4741043A (en) * 1985-11-04 1988-04-26 Cell Analysis Systems, Inc. Method of and an apparatus for image analyses of biological specimens
US5525464A (en) * 1987-04-01 1996-06-11 Hyseq, Inc. Method of sequencing by hybridization of oligonucleotide probes
US5202231A (en) * 1987-04-01 1993-04-13 Drmanac Radoje T Method of sequencing of genomes by hybridization of oligonucleotide probes
US5492806A (en) * 1987-04-01 1996-02-20 Hyseq, Inc. Method of determining an ordered sequence of subfragments of a nucleic acid fragment by hybridization of oligonucleotide probes
US4965725B1 (en) * 1988-04-08 1996-05-07 Neuromedical Systems Inc Neural network based automated cytological specimen classification system and method
US4965725A (en) * 1988-04-08 1990-10-23 Nueromedical Systems, Inc. Neural network based automated cytological specimen classification system and method
US5002867A (en) * 1988-04-25 1991-03-26 Macevicz Stephen C Nucleic acid sequence determination by multiple mixed oligonucleotide probes
US5143854A (en) * 1989-06-07 1992-09-01 Affymax Technologies N.V. Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof
US5445934A (en) * 1989-06-07 1995-08-29 Affymax Technologies N.V. Array of oligonucleotides on a solid substrate
US5527681A (en) * 1989-06-07 1996-06-18 Affymax Technologies N.V. Immobilized molecular synthesis of systematically substituted compounds
US5235626A (en) * 1991-10-22 1993-08-10 International Business Machines Corporation Segmented mask and exposure system for x-ray lithography
US5384261A (en) * 1991-11-22 1995-01-24 Affymax Technologies N.V. Very large scale immobilized polymer synthesis using mechanically directed flow paths
US5288514A (en) * 1992-09-14 1994-02-22 The Regents Of The University Of California Solid phase and combinatorial synthesis of benzodiazepine compounds on a solid support
US5503980A (en) * 1992-11-06 1996-04-02 Trustees Of Boston University Positional sequencing by hybridization
US5470710A (en) * 1993-10-22 1995-11-28 University Of Utah Automated hybridization/imaging device for fluorescent multiplex DNA sequencing
US5733729A (en) * 1995-09-14 1998-03-31 Affymetrix, Inc. Computer-aided probability base calling for arrays of nucleic acid probes on chips
US6361937B1 (en) * 1996-03-19 2002-03-26 Affymetrix, Incorporated Computer-aided nucleic acid sequencing

Also Published As

Publication number Publication date
US6361937B1 (en) 2002-03-26

Similar Documents

Publication Publication Date Title
US6361937B1 (en) Computer-aided nucleic acid sequencing
US6957149B2 (en) Computer-aided probability base calling for arrays of nucleic acid probes on chips
EP0923050B1 (en) Scanned image alignment systems and methods
US6242180B1 (en) Computer-aided visualization and analysis system for sequence evaluation
US6600996B2 (en) Computer-aided techniques for analyzing biological sequences
US6223127B1 (en) Polymorphism detection utilizing clustering analysis
US8204695B2 (en) Techniques for identifying, confirming, mapping, and categorizing polymers
US20070067111A1 (en) Computer-aided visualization of expression comparison
US6699659B2 (en) Products and methods for analyzing nucleic acids including identification of substitutions, insertions and deletions
US20040175718A1 (en) Computer-aided visualization and analysis system for sequence evaluation
KR100442839B1 (en) Method for scoring and selection for optimum probes in probes design
US20030220748A1 (en) Computer-aided techniques for analyzing biological sequences
CA2302193A1 (en) Techniques for identifying, confirming, mapping and categorizing polymers
EP1632579A2 (en) Computer-aided techniques for analyzing biological sequences

Legal Events

Date Code Title Description
AS Assignment

Owner name: GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT, MA

Free format text: SECURITY AGREEMENT;ASSIGNOR:AFFYMETRIX, INC.;REEL/FRAME:028465/0541

Effective date: 20120625

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION

AS Assignment

Owner name: AFFYMETRIX, INC., CALIFORNIA

Free format text: RELEASE BY SECURED PARTY;ASSIGNOR:GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT;REEL/FRAME:037109/0132

Effective date: 20151028