US20020106673A1 - Method of sequencing by hybridization of oligonucleotide probes - Google Patents
Method of sequencing by hybridization of oligonucleotide probes Download PDFInfo
- Publication number
- US20020106673A1 US20020106673A1 US09/954,426 US95442601A US2002106673A1 US 20020106673 A1 US20020106673 A1 US 20020106673A1 US 95442601 A US95442601 A US 95442601A US 2002106673 A1 US2002106673 A1 US 2002106673A1
- Authority
- US
- United States
- Prior art keywords
- sequence
- hybridization
- subfragments
- nucleic acid
- target
- 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
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6813—Hybridisation assays
- C12Q1/6827—Hybridisation assays for detection of mutation or polymorphism
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
- C12Q1/6874—Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/20—Sequence assembly
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
Definitions
- the present invention belongs to the field of molecular biology. It involves a novel method of sequencing of a target nucleic acid sequence by hybridization of short oligonucleotide probes to a nucleic acid target.
- the oligonucleotide probes can comprise all known combinations of the four nucleotides of a given length, i.e. oligonucleotides of base composition adenine (A), thymine (T), guanine (G), and cytosine (C) for DNA and A, G, C, and uridine (U) for RNA. Conditions are described which allow hybridization discrimination between oligonucleotides which are as short as six nucleotides long and have a single base end-mismatch with the target sequence.
- Hybridization depends on the pairing of complementary bases in nucleic acids and is a specific tool useful for the general recognition of informational polymers.
- Diverse research problems using hybridization of synthetic oligonucleotide probes of known sequence include, amongst others, the different techniques of identification of specific clones from cDNA and genomic libraries; detecting single base pair polymorphisms in DNA; generation of mutations by oligonucleotide mutagenesis; and the amplification of nucleic acids in vitro from a single sperm, an extinct organism, or a single virus infecting a single cell.
- Hybridization conditions that discriminate between (1) a perfectly complementary hybridizing pair of nucleic acid sequences where one partner of the pair is a short oligonucleotide, and (2) a pair wherein a mismatch of one nucleotide occurs on the 5′ or 3′ end of the oligonucleotide, provide a more stringent environment than is required for internal mismatches because hybrid stability is affected less by a mismatch at the end of a hybridizing pair of complementary nucleic acids than for an internal mismatch.
- Nucleic acid sequencing methods where the position of each base in a nucleic acid molecule in relation to its neighbors is determined to define its primary structure, were developed in the early 1960's for RNA molecules and in the late 1970's for DNA.
- the two major methods for DNA sequencing i.e. chemical degradation and dideoxy-chain termination, involve identification and characterization of 1-500 nucleotide long DNA fragments, specific for each one of at least four nucleotide bases, on polyacrylamide gels.
- the polyacrylamide gels must be able to distinguish single base pair differences in length between fragments. The fragments are generated either by chemical degradation [Maxam-Gilbert, Proc. Natl. Acad. Sci.
- the genome of higher eucaryotes has up to a million times greater physical complexity than is the complexity of individual genes it encodes, giving it a a corresponding huge informational complexity. From the present knowledge of genome organization and biochemical, biophysical and biological functions, the following approximate scale of the informational complexity for higher eucaryotes can be proposed: 10,000 gene families—100,000 genes—1,000,000 biological functions. The number of basic biochemical functions represented by a single gene family is probably not significantly increased compared to procaryotic and lower eucaryotic genomes.
- sequencing projects are extended to include the determination of the genomic sequences of characteristic individuals or species of organisms, especially those that have economic, social or medical importance.
- Such sequencing projects would advance not only our understanding of the evolution of organisms and the evolution of biochemical processes, but would also further the detection, treatment and understanding of disease, and would aid agriculture, the food industry and biotechnology in general.
- beneficial the results of such projects would be, their successful completion requires the development of a new, rapid, reproducible and reliable sequencing method such as those described in this invention.
- the present invention provides a new method of sequencing that is ideally suited to the sequencing of large complex genomes because it avoids the intensive manual labor involved in resolving gel fragments by size on polyacrylamide gels.
- the present invention provides methods for sequencing a target nucleic acid by hybridization of overlapping short oligonucleotide probes of known or predicted sequence to the nucleic acid target serially or simultaneously.
- the oligomer probes of a given size can contain all or most existing combinations of nucleotides for complete sequencing and a part of all possible variants for partial sequencing. Probes can also be composed of oligomers of different sizes as well as comprising all known combinations of nucleotides that are possible for that size oligonucleotide. As the size of the probes that are used decreases, hybridization conditions that are still able to distinguish between mismatched and perfectly matched short oligonucleotides must be used.
- multiple oligonucleotides that are 11 nucleotides long or longer are hybridized to the target sequence. hybridization occurs using conditions which are controlled and varied to ensure discrimination between perfectly matched oligonucleotides and oligonucleotides having a one base pair mismatch with the target sequence where the mismatch is located at either one of two ends of the oligonucleotide.
- the sequence of the target nucleic acid is generated by an algorithm using the principle of maximal nonidentical overlap of probe.
- oligonucleotides are prepared, target fragments are prepared appropriate for the length of oligonucleotide used for hybridization, and hybridization of the target with all the oligonucleotides occurs under defined conditions that allow discrimination in binding of perfectly matched complementary oligonucleotides and mismatched oligonucleotides.
- the relationship of probe size and target length is defined and allows complete sequencing of genomes. The novel theoretical basis of the relationship between oligonucleotide probe size and target length is described infra.
- the number of target fragments that compose the entire sequence is multiplied by the number of different oligonucleotides required to define the sequence of the target fragment. The shorter the size of the oligonucleotides that are hybridized, the more target fragments that must be analyzed. Similarly, as the oligonucleotide size increases, fewer target fragments must be examined.
- Hybridization reactions can be performed in separate reaction vessels or by binding one of the two components (oligomers and DNA fragments) to a solid surface, like nylon filters etc. Since the described method does not require macromolecular separation like gel-based sequencing methods, the surface, bound with either an oligomer or nucleic acid fragment can have microdimensions.
- a adenine bp base pair C cytosine G guanine IF an M13 clone containing a 921 bp EcoR1-Bg1II human ⁇ 1 interferon fragment kD kilo Dalton nG nanogram nM nenomolar pmol picomole sc subclone SF subfragment SOH short oligonucleotide hybridization T thymine CCD Charge Coupled Device DNA Deoxyribonucleic acid DP Discrete particle HA Hybridization area LAR Ligation-amplification reaction ON Oligonucleotide ONP Oligonucleotide probe ONS Oligonucleotide sequence PCR Polymerase chain reaction RE Restriction Enzyme RFLP Restriction fragment length polymorphism RNA Ribonucleic acid SBH Sequencing by hybridization
- FIG. 1 shows the generation of subfragments in sequence by hybridization and their ordering.
- A The sequence of a hypothetical clone. NNNNNNN-ends of vector sequence. AGTCCCT and TTGGCTG are the only oligonucleotides 7 bp or longer repeated within the depicted sequence.
- B Formation of subfragments. Assuming that the content of 8-mers for the depicted sequence is known, these 8-mers are ordered by maximal overlap, in this case 7 bp. Starting from the 5′ 8-mer (NNNNNNNc), ordering is unambiguous up to gAGTCCCT, which on its 3′ end contains a repeated 7-mer.
- AGTCCCTc and AGTCCCTg can be overlapped with gAGTCCCT preventing further ordering.
- Each of the two sequences serves as a starting point for new ordering (not shown). Therefore, each repeated sequence 7 bp or longer represents a branching point. Unambiguous sequences are obtained between two consecutive branching points only.
- C Listing of subfragments formed from 8-mers of depicted sequence. Subfragments are horizontally displaced to indicate overlap; the orientation is 5′ to 3′ and end subfragments are identifiable.
- D The subfragments cannot be unambiguously ordered into a starting sequence without additional information.
- FIG. 2 presents the average number of SFs (N sf ) as a function of the length of DNA fragment (L f ) for various values of the length of the overlapping sequence (N-1, in bp), or average distance of two consecutive identical N-1 sequences in DNA subjected to sequencing by hybridization (A o ), in kb.
- the curves are obtained using equation one as described below in section 5.2.
- FIG. 3 describes the kinetic stability of a fully matched hybrid obtained with a probe 8 nucleotides in length. Stability is expressed as a fraction of the hybrid dissociated in unit time (minutes) as a function of temperature.
- the equal amounts of hybrid were incubated at the indicated temperatures for a short time in a large volume of buffer and the remaining hybrid measured. Each point represents the average value for four dots.
- FIG. 4 indicates the properties of short oligonucleotide hybridization.
- FIG. 4 a non-optimized discrimination with probes 6, 7, and 8 nucleotides in length is illustrated.
- the probe GCTCAT was hybridized to the target sequence NC ATGAGC ANN which contains the perfectly matching sequence (underlined).
- the NNCATGAG T TN target sequence contains an end mismatch (double underlined). 1.4 pM of each target was applied to the filter.
- the probe GCTCATG, and the prose TGCTCATG were used against 50 ng of IF and M13 DNA. The probe concentration was 4 nM.
- FIG. 4 b limits of signal detection are examined.
- the indicated volumes of IF culture supernatants of average titer of 6 ⁇ 10 11 pfu/ml were mixed with an equal volume of 1 M NaOH, 3M NaCl and spotted on a filter as described in a above. Hybridization was at 2° C. with TGCTCATG as the probe.
- FIG. 4 c the time course of hybridization at 13° C. is shown.
- the IF-M13 system was used with 50 ng of phage DNA per dot, and the probe was TGCTCATG.
- the 3 hr IF dot contained 18020 cpm measured with 20% efficiency.
- FIG. 5 the effect of the washing step on discrimination is indicated.
- FIG. 5 a inversion of the signal in IF-M13 pair upon washing is shown. 10 ng of IF and 500 ng M13 DNA were applied, and the probe was TGCTCATG. The top row was not washed, the other rows were washed at 7, 13 and 25° C., respectively for the indicated times. A DNA control is included in the top row also. Hybridization with the M13 specific probe AGCTGCTC measures amounts of DNA in the two dots.
- FIG. 5 b the change of discrimination with time of washing at 0° C. (full circles) and 13° C. (open circles) is depicted.
- FIG. 6 demonstrates the effects of complexity of target sequences on discrimination. 50 ng each of IF and M13 were hybridized with the indicated probes at a concentration of 4 nM. No wash was performed. The number of matched and end base mismatched targets in IF and M13 is indicated for each probe.
- FIG. 7 examines an array of clones for the presence of an oligonucleotide sequence.
- 51 recombinant plasmid DNAs (10+5 ng) were spotted in rows B to II, columns 1 to 8 (except row H).
- Line A and column 9 contained control DNAs of known sequence.
- Unknown clones were taken from human brain cDNA library in Bluescript vector (BS) (Stratogene Cat. No. 935205).
- Controls of known sequence in lines A1 to A8 and A9 to G9 are: IF(M13), M13, Alu(M3), IF(BS), BS , 1M(pUC 9), pUC 9, 2M(pUC 9), respectively except that in the vertical row Alu(M13) was omitted.
- 1M and 2M are rat ⁇ -globin gene subclones.
- the probe concentration was 8 nM.
- BS specific probe CTCCCTTT was also contained in IF and 2M inserts but not in M13 and pUC vectors.
- the sequence of probe CCAGTTTT was contained in the IF insert but not in either vector.
- the sequence of probe GCCTTCTC was contained in the 1M insert only.
- FIG. 8 sequencing 100 bp of 921 bp ⁇ 1 -human interferon gene fragment. (IF) by hybridization.
- Hybridization results A. Hybridization with 93 probes (72) octamers and 21 nonamers with the full match in IF. IF and controls rat globin clones pHEA and pHI were PCR amplified while M13 mp18 and pUC18 were in linearized double stranded form. Base denatured DNA (20 ng of IF and equimolar amounts of control DNA) were spotted on Gene Screen membranes (N.E.N.). Hybridization was according to Drmanac et al., described in ⁇ 6 below.
- ⁇ - 32 end labeled probes (3.3 pm, 10 mCI, Amersham 3000 C/mM in concentration of 10 ng/ml were hybridized at 12° C. in 0.5 M Na 2 HPO 4 pH 7.2, 7% Na-lauryl Sarcosine for 3 hours. All probes were made by Genesys, Inc., Houston. Hybrids were washed in 6 ⁇ SSC at 0° for 40 minutes and autoradiographed for 4-48 hours. Test dot signal intensity, Hp, and discrimination as ratio of signals of test over control dot, D, were visually estimated. For probes 34 and 74, dot radioactivity was measured in a scintillation counter.
- Hp was 6,000 and 300 cpm
- D was 20 and 4, and a film was exposed for 4 and 48 hr respectively.
- B Hybridization with 12 probes (11 octamers and 1 nonamer) which have end mismatch in IF fragment. Control DNAs having single full match targets were pHEA for probes 97., 98., 102., pUC18 for 95., 100., 104., 105., and M13 for 94., 96., 99., 101., and 103. Probes 104 and 105 have 3 end-mismatched targets in IF. Hybridization procedures were as described in A.
- C DNA Calibration. 1. and 2. IF and pHEA, probe CTGATAT. 3.
- the ratios of DNA amounts in IF and control dot were 1:1 in panel 1., 3., 4., and 1:3 in panel 2., respectively.
- the remainder of probes were used on filters of the type shown on panel 1.
- FIG. 9 an example of the discriminative hybridization with 8-mers.
- the upper dot contains 15 ng of EcoRI-BgIII 920 bp genomic fragment with interferon beta-1 gene (IF) and the lower dot contains 25 ng of 1500 bp of rat beta-globin gene (GL) prepared by PCR from clones in bluescript vector using universal and reverse sequencing primers.
- Each probe was hybridized and washed under the same conditions.
- 32 p labeled probes (by T4-polynucleotide kinase) at a concentration of 4 nM in 0.5 M Na 2 HPO 4 pH 7.2, 7% Sodium lauroyl sarcosine were hybridized at 12° C. for 3 h.
- Hybrids were washed in 6 ⁇ SSC at 2° C. for 40 minutes and autoradiographed for 6h. Panel on top represents hybridization with probe complementary to amplified vector end and reveals the molar ratio of the two target DNAS.
- First three 8-mer probes are designed to have perfect target in interferon sequence and target with end mismatch in globin sequence, and probe 4 opposite. Over 5 fold discrimination is obtained with each probe, having 3 to 6 (C+G). In the case of the most stable probe (probe 1) similar discrimination is obtained but both dots are stronger than in case of other probes.
- FIG. 10 k-tuple overlap rule. S-start of 30 subfragment assembly. Filled bars—the identical parts of k-tuples. Empty bars—non-identical parts of k-tuples. 0, 1 non occurrence or occurrence, respectively in k-tuple set. (1-3)—situations violating the rule. (4) situation recognized by the overlapping rule.
- FIG. 11 Examples of generation of short pSFs. Bars are k-tuples. Filled parts of are bars overlapping sequences which length are indicated on the right site. Open bars are correct k-tuples. Bars with dashed line are false negative k-tuples. Bars with (+) are false positive k-tuples; positions of (+) indicate the wrong base. Horizontal lines are extensions of subfragments in the considered site.
- FIG. 12 Process of formation of misconnected subfragments.
- Column I Misconnection due to the presence of false negative k-tuples in the set.
- Column II Misconnection due to-t ⁇ t simultaneous presence of both false positive and false negative k-tuples.
- Oblique lines show the wrong position of a part of a sequence within the misconnected subfragment. Dashed line indicate the position in subfragment of sequence used for overlap, causing the misconnection.
- the method of the invention of sequencing by hybridization is based on an entirely different logical principle than previous methods of sequencing which have been described above in ⁇ 2.3.
- Known methods of sequencing define a nucleic acid sequence in a base by base fashion that is read from the position of DNA fragments in polyacrylamide gels where the fragments are produced by base specific chemical degradation or chain termination techniques.
- the method of sequencing of the present invention is a process comprising hybridizing oligonucleotides to a target sequence, detecting the hybridized oligonucleotides, and assembling the sequence of the target from the sequence of the hybridizing oligonucleotides.
- the method of sequencing according to the present invention is most suitable for determining the sequence of complex nucleic acid fragments, molecules, or genomes and especially those containing one million or more nucleotide base pairs.
- the number of oligonucleotide probes, the lengths and types of target nucleic acid sequences, the numbers of subfragments of the target nucleic acid sequences and the numbers of hybridizations necessary for sequencing have been calculated and show that sequencing by hybridization is more efficient than sequencing by traditional methods.
- a mismatch-free hybridization of oligonucleotides to an unknown target nucleic acid represents, in essence, a sequencing of complementary target. Probability calculations and, in part, computer simulations estimate the types and numbers of oligonucleotides that would have to be synthesized in order to sequence a megabase plus segment of DNA.
- the number of probes can be reduced depending on the number of fragments used and vice versa.
- the optimal ratio of the two numbers is defined by the technological requirements of a particular sequencing by hybridization process.
- the useful probes are those longer than 6 bases, since the shorter ones would require use of unsuitable and unpractical nucleic acid fragments shorter than a few hundred bases long.
- probe length There are two parameters which influence the choice of probe length. The first is the success in obtaining hybridization results that show the required degree of discrimination. The second is the technological feasibility of synthesis of the required number of probes. The requirement of obtaining sufficient hybridization discrimination with practical and useful amounts of target nucleic acid limits the probe length from both sides. It is difficult to obtain a sufficient amount of hybrid with short probes, and to discriminate end mismatch with long probes. There is no evidence for use of probes shorter than 11-mers in the literature, except for very stable ones [Estivill et al., Nucl. Acids Res. 15: 1415 (1987)]. On the other hand, probes longer than 15 bases discriminate end mismatch with difficulty [Wood et al., Proc. Nat'l. Acad. Sci. USA 82: 1585 (1985)].
- One solution for the problem of unstable probes and end mismatch discrimination is the use of a group of longer probes representing a single shorter probe in an informational sense. Groups of sixty-four 11-mers can be used instead of single 8-mers. Every member of the group has a common core 8-mer and one of three possible variations on outer positions with two variations at each end.
- the probe can be represented as 5′(A,T,C,G)(A,T,C,G)N8(A,T,C,G) 3′. With this type of probe one does not need to discriminate the non-informative end bases (two on 5′ end, and one on 3′ end) since only the internal 8-mer is read. This solution requires the use of higher mass amounts of probes and label in hybridization reactions.
- the number of hybridization reactions is dependent on the number of discrete labeled probes. Therefore in the cases of sequencing shorter nucleic acids using a smaller number of fragments-clones than the number of oligonucleotides, it is better to use oligomers as the target and nucleic acid fragments as probes. For example, one needs 20,000 clone. and up to five fold more oligonucleotides (50-100,000) to sequence 10 million bp of DNA with 8-mer oligonucleotides.
- Target nucleic acids which have undefined sequences can be produced as a mixture of representative libraries in a phage or plasmid vector having inserts of genomic fragments of different sizes or in samples prepared by PCR. Inevitable gaps and uncertainties in alignment of sequenced fragments arise from nonrandom or repetitive sequence organization of complex genomes and difficulties in cloning “poisonous” sequences in Escherichia coli. These problems are inherent in sequencing large complex molecules using any method. In this invention such problems are minimized by the choice of libraries and number of subclones used for hybridization. Alternatively such difficulties can be overcome through the use of amplified target sequences, e.g. by PCR amplification, ligation-amplified reactions etc. Because this invention is based on simple biochemical procedures, the method of the invention is inherently easier to automate than existing sequencing methods. The sequence can be derived from simple primary data using only extensive computing.
- Sequencing by hybridization (SBH) of a target nucleic acid can be visualized as consisting of two steps; 1) a process of dissolving the target nucleic acid into all its constituent oligonucleotide N-mers, and 2) the back assembly of N-mers detected by hybridization and assembled by overlap into an extended sequence.
- hybridization of all possible N-mer oligonucleotide probes to the target nucleic acid determines the N-mer oligonucleotide subset contained in the primary sequence of the target nucleic acid and is the first step in the process of the invention.
- a nucleic acid fragment can be dissolved into all constituent oligonucleotides.
- Positively hybridizing N-mer oligonucleotide probes are ordered and the sequence of the target DNA is determined using (N-1)mer overlapping frames between the oligonucleotide probes.
- sequence assembly will be interrupted where ever a given overlapping (N-1)mer is duplicated two or more times. Then either of the two N-mers differing in the last nucleotide can be used in extending the sequence. This branching point limits unambiguous assembly of sequence.
- sequence subfragment results if any part of the sequence of a target nucleic acid starts sand ends with an (N-1)mer that is repeated two or more times within the target sequence.
- subfragments are sequences generated between two points of branching in the process of assembly of the sequences in the method of the invention. The sum of all subfragments is longer than the actual target nucleic acid because of overlapping short ends.
- Nsf 1 +A 0 ⁇ K ⁇ P ( K, L f )
- P(K, L f ) represents the probability of an N-mer occurring K-times on a fragment L f base long.
- FIG. 2 the number of subfragments increases with the increase of lengths of fragments for a given length of probe. Obtained subfragments may not be uniquely ordered among themselves. Although not complete, this information is very useful for comparative sequence analysis and the recognition of functional squence characteristics. This type of information we call partial sequence. Another way of obtaining partial sequence is the use of only a subset of oligonucleotide probes of a given length.
- This number consists of 700 random 7 kb clones (basic library), 1250 pools of 20 clones of 500 bp (subfragments ordering library) and 150 clones from jumping (or similar) library.
- the developed algorithm can regenerate sequence using hybridization data of these described samples.
- oligonucleotides were either purchased from Genetic Designs, Inc. Houston, Texas or made on an Applied Biosystems 381A DNA synthesizer. Most of the probes used were not purified by HPLC or gel electrophoresis. Probes were designed to have both a single perfectly complementary target in IF, a M13 clone containing a 921 bp Eco RI-Bgl II human ⁇ 1 -interferon fragment [Ohno and Tangiuchi, Proc. Natl. Acad. Sci. 74: 4370-4374 (1981)], and at least one target with an end base mismatch in M13 vector itself.
- oligonucleotides End labeling of oligonucleotides was performed as described [Maniatis et al., in Molecular Cloning: a Laboratory Manual, Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y. (1982)] in 10 ⁇ l containing T 4 -polynucleotide-kinase (5 units Amersham), 7 - 32 p -ATP (3.3 pM, 10 ⁇ Ci Amersham 3000 Ci/mM) and oligonucleotide (4 pM, 10 ng). Specific activities of the probes were 2.5-5 ⁇ 10 9 cpm/nM.
- a drop of hybridization solution (10 ⁇ l, 0.5 M Na 2 HPO 4 pH 7.2, 7% sodium lauroyl sarcosine) with a 32 P end labeled oligomer probe at 4 nM concentration was placed over 1-6 dots, overlayed with a square piece of polyethylene (approximately 1 ⁇ 1 cm), and incubated in a moist chamber at the indicated temperatures for 3 hr.
- Hybridization was stopped by placing the filter in 6 ⁇ SSC washing solution for 3 ⁇ 5 minute at 0° C. to remove unhybridized probe.
- the filter was either dried, or further washed for the indicated times and temperatures, and autoradiographed.
- the dots were excised from the dried filters after autoradiography, placed in liquid scintillation cocktail and counted. The uncorrected ratio of cpms for IF and M13 dots is given as D.
- S t and S o are target sequence concentrations at time t and t o , respectively.
- (OP) is probe concentration.
- Kh increases only slightly in the 0° C. to 30° C. range [Porschke and Eigen, J. Mol. Biol. 62: 361 (1971); Craig et al., J. Mol. Biol. 62: 383 (1971)].
- Hybrid melting is a first order reaction with respect to hybrid concentration (here replaced by mass due to filter bound state) as shown in equation three:
- H t and H o are hybrid concentrations at times t and t o , respectively;
- k m is a rate constant for hybrid melting which is dependant on temperature and salt concentration [Ikuta et al., Nucl. Acids Res. 15: 797 (1987); Porschke and Eigen, J. Mol. Biol. 62: 361 (1971); Craig et al., J. Mol. Biol. 62: 383 (1971)].
- the back melting, or strand dissociation, reaction takes place as well.
- the amount of hybrid formed in time t is a result of forward and back reactions.
- the equilibrium can be moved towards hybrid formation by increasing probe concentration and/or decreasing temperature. However, during washing cycles in large volumes of buffer, the melting reaction is dominant and the back reaction of hybridization is insignificant, since the probe is absent.
- SOH short oligonucleotide hybridization
- H p (t w ) and H i (tW) are the amounts of hybrids remaining after a washing time, t w , for the identical amounts of perfectly and imperfectly complementary duplex, respectively.
- the background, B represents the lowest hybridization signal detectable in the system. Since any further decrease of H i cannot be examined, D would increase upon continued washing. Washing past t w just decreases H p relative to B, and is seen as a decrease in D.
- the optimal washing time, t w for imperfect hybrids, from equation three and equation five is Equation six:
- Equation seven Since H p is being washed for the same t w , combining Equation three for Hp, Equation four and Equation six, one obtains the optimal discrimination function in Equation seven:
- Equation nine for E ⁇ ,i less than E ⁇ ,p is that discrimination more strongly decreases with T, if H i (t o ) increases relative to H p (t o ).
- Equation nine shows that D at lower temperatures depends to a higher degree on the H p (t o )/B ratio than on the H p (t o )/H i (t o ) ratio. This result indicates that it is better to obtain a sufficient quantity of H p in the hybridization regardless of the discrimination that can be achieved in this step. Better discrimination can then be obtained by washing, since the higher amounts of perfect hybrid allow more time for differential melting to show an effect. Similarly, using larger amounts of target nucleic acid a necessary discrimination can be obtained even with small differences between K m,p and K m,i .
- the detergent in these buffers is essential for obtaining tolerable background with up to 40 nM concentrations of labeled probe.
- Preliminary characterization of the thermal stability of short oligonucleotide hybrids was determined on a prototype octamer with 50% GC content, i.e. probe of sequence TGCTCATG. The theoretical expectation is that this probe is among the less stable octamers. Its transition enthalpy is similar to those of more stable heptamers or, even to probes 6 nucleotides in length [Bresslauer et al., Proc. Nati. Acad. Sci. U.S.A. 83:.3746 (1986)].
- T d the temperature at which 50% of the hybrid is melted in unit time of a minute is 18° C. The result shows that T d is 15° C. lower for the 8 bp hybrid than for an 11 bp duplex [Wallace et al., Nucleic Acids Res. 6: 3543 (1979)].
- an M13 vector was chosen as a system for a practical demonstration of short oligonucleotide hybridization.
- the main aim was to show useful end-mismatch discrimination with a target similar to the ones which will be used in various applications of the method of the invention.
- Oligonucleotide probes for the M13 model were chosen in such a way that the M13 vector itself contains the end mismatched base.
- Vector IF an M13 recombinant containing a 921 bp human interferon gene insert, carries a single perfectly matched target. Thus, IF has either the identical or a higher number of mismatched targets in comparison to the M13 vector itself.
- the hybridization signal depends on the amount of target available on the filter for reaction with the probe. A necessary control is to show that the difference in signal intensity is not a reflection of varying amounts of nucleic acid in the two dots.
- Hybridization with a probe that has the same number and kind of targets in both IF and M13 shows that there is an equal amount of DNA in the dots (FIG. 4A). Since the efficiency of hybrid formation increases with hybrid length, the signal for a duplex having six nucleotides was best detected with a high mass of oligonucleotide target bound to the filter (FIG. 4A). Due to their lower molecular weight, a larger number of oligonucleotide target molecules can be bound to a given surface area when compared to large molecules of nucleic acid that serves as target.
- the equilibrium yield of hybrid depends on probe concentration and/or temperature of reaction. For instance, the signal level for the same amount of target with 4 nM octamer at 13° C. is 3 times lower than with a probe concentration of 40 nM, and is decreased 4.5 times by raising the hybridization temperature to 25° C.
- the M13 system has the advantage of showing the effects of target DNA complexity on the levels of discrimination.
- FIG. 6 shows the discrimination levels prior to washing for several probes with the IF-M13 system. The number of mismatched targets present in the insert and M13 vector are also shown. The general trend of a decrease in measured discrimination with the number of end-mismatched targets is obvious. Thus for two octamers having either none or five mismatched targets and differing in only one GC pair the observed discriminations were 18.3 and 1.7, respectively (FIG. 6).
- FIG. 7 shows the variation in dots from sample to sample, since the differences in hybridization with the vector specific probe should reflect the differences in DNA concentration from dot to dot.
- This probe sequence also has a complementary target in the interferon insert. Since the M13 dot (A2 FIG.
- FIG. 7 a is negative while the interferon insert in either M13 or Bluescript (Al and A4, respectively, FIG. 7 a ) is positive, the hybridization is sequence specific.
- probes that detect the target sequence in only one of 51 inserts (FIG. 7 b ), or in none of the examined inserts (FIG. 7 c ) are shown along with controls that confirm that hybridization would have ocurred if the appropriate targets were present in the clones.
- scoring positive hybridization comparison of signal intensities was accomplished not between dots on the same panel, but comparisons were also made for the same dot on panels b and c with the dot on panel a in order to correct for differences in DNA amounts in the dots and, in the case of controls, in complexity.
- Thermal stability curves for very short oligonucleotide hybrids that are 6-8 nucleotides in length are at least 15° C. lower than for hybrids 11-12 nucleotides in length [FIG. 1 and Wallace et al., Nucleic Acids Res. 6: 3543-3557 (1979)].
- performing the hybridization reaction at a low temperature and with a very practical 4-40 nM concentration of oligonucleotide probe allows the detection of complementary sequence in a known or unknown nucleic acid target.
- an entire set containing 65,535 8-mer probes can be used.
- nucleic acid for this purpose are present in convenient biological samples such as a few microliters of M13 culture, a plasmid prep from 10 ⁇ l of bacterial culture or a single colony of bacteria, or less than 1 ⁇ l of a standard PCR reaction.
- probes having E ⁇ less than 45 Kcal/mol e.g. for many heptamers and a majority of hexamers modified oligonucleotides give a more stable hybrid [Asseline, et al., Proc. Natl. Acad. Sci. 81: 3297 (1984)] than their unmodified counterparts.
- hybridization conditions described in this invention for short oligonucleotide hybridization using low temperatures give better discrimination for all sequences and duplex hybrid inputs.
- the only price paid in achieving uniformity in hybridization conditions for different sequences is an increase in washing time from minutes to up to 24 hours depending on the sequence. Moreover, the washing time can be further reduced by decreasing the salt concentration.
- mismatch oligonucleotides wherein internal or end mismatches would occur in the duplex formed by the oligonucleotide and the target were examined.
- the lowest practical temperature was used to maximize hybridization formation. Washes were accomplished at the same or lower temperatures to ensure maximal discrimination by utilizing the greater dissociation rate of mismatch versus matched oligonucleotide/target hybridization.
- the least distabilizing mismatch that can be postulated is a simple end mismatch, so that the test of sequencing by hybridization is the ability to discriminate perfectly matched oligonucleotide/target duplexes from end-mismatched oligonucleotide/target duplexes.
- the set of 93 probes provided consecutive overlapping frames of the target sequence displaced by one or two bases as described in FIG. 8.
- hybridization was examined for 12-additional probes that contained at least one end mismatch when hybridized to the 100 bp test target sequence. Also tested was the hybridization of twelve probes with target end-mismatch to four other control nucleic acid sequences chosen so that the 12 oligonucleotides could form perfectly matched duplex hybrids with the four control DNAs. Thus, the hybridization of internal mismatched, end-mismatched and perfectly matched duplex pairs of oligonucleotide and target were evaluated for each oligonucleotide used in the experiment.
- Hp defines the amount of hybrid duplex formed between a test target and an oligonucleotide probe.
- D Discrimination
- D was defined as the ratio of signal intensities between 1) the dot containing a perfect matched duplex formed between test oligonucleotide and target or control nucleic acid and 2) the dot containing a mismatch duplex formed between the same oligonucleotide and a different site within the target or control nucleic acid.
- Variations in the value of D result from either 1) perturbations in the hybridization efficiency which allows visualization of signal over background, or 2) the type of mismatch found between the test oligonucleotide and the target.
- the D values obtained in this experiment were between 2 and 40 for 102 of the 105 oligonucleotide probes examined. Calculations of D for the group of 102 oligonucleotides as a whole showed the average D was 10.6.
- Low hybrid formation values of lip could result from hybrid instability or from an inability to form hybrid duplex.
- An inability to form hybrid duplexes would result from either 1) self complementarity of the chosen probe or 2) target/target self hybridization.
- Oligonucleotide/oligonucleotide duplex formation could be favored over oligonucleotide/target hybrid duplex formation if the probe was self-complementary.
- target/target association could be favored if the target was self-complementary or could form internal palindromes. In evaluating these possibilities, it was apparent from probe analysis that the questionable probes were unable to form hybrids with themselves.
- the sequence which is generated by computer using the algorithm described below is of high fidelity.
- the algorithm can tolerate false positive signals from the hybridization dots as is indicated from the fact the sequence generated from the 105 hybridization values, which included four less reliable values, was correct.
- This fidelity in sequencing by hybridization is due to the “all or none” kinetics of short oligonucleotide hybridization and the difference in duplex stability that exists between perfectly matched duplexes and mismatched duplexes.
- the ratio of duplex stability of matched and end-mismatched duplexes increases with decreasing duplex length. Moreover binding energy decreases with decreasing duplex length resulting in a lower hybridization efficiency.
- results presented in other sections show that oligonucleotides that are 6, 7 or 8 nucleotides can be effectively used to generate reliable sequence on targets that are 0.5 kb (for hexamers) 2 kB (for septamers) and 6 kB long (for octamers).
- the sequence of long fragments can be overlapped to generate a complete genome sequence.
- This appendix describes an algorithm for generation of a long sequence written in a four letter alphabet from constituent k-tuple words in a minimal number of separate, randomly defined fragments of a starting nucleic acid sequence where K is the length of an oligonucleotide probe.
- the algorithm is primarily intended for use in the sequencing by hybridization (SBH) process.
- the algorithm is based on subfragments (SF), informative fragments (IF) and the possibility of using pools of physical nucleic sequences for defining informative fragments.
- subfragments can be caused by branch points in the assembly process resulting from the repetition of a K-1 oligomer sequence in a target nucleic acid.
- Subfragments are sequence fragments found between any two repetitive words of the length K-1 that occur in a sequence. Multiple occurrences of K-1 words are the cause of interruption of ordering the overlap of K-words in the process of sequence generation. Interruption leads to a sequence remaining in the form of subfragments. Thus, the unambiguous segments between branching points whose order can not be uniquely determined are called sequence subfragments.
- Appendix A informative fragments are defined as fragments of a sequence that are determined by the nearest ends of overlapped physical sequence fragments.
- a certain number of physical fragments can be pooled without losing the possibility of defining informative fragments.
- the total length of randomly pooled fragments depends on the length of k-tuples that are used in the sequencing process.
- the algorithm consists of two main units.
- the first part is used for generation of subfragments from the set of k-tuples contained in a sequence.
- Subfragments can be generated within the coding region of physical nucleic acid sequence of certain sizes, or within the informative fragments defined within long nucleic acid sequences. Both types of fragments are members of the basic library.
- This algorithm does not describe the determination of the content of the k-tuples of the informative fragments of the basic library, i.e. the step of preparation of informative fragments to be used in the sequence generation process.
- the second part of the algorithm determines the linear order of obtained subfragments with the purpose of regenerating the complete sequence of the nucleic acid fragments of the basic library.
- a second, ordering library is used, made of randomly pooled fragments of the starting sequence.
- the algorithm does not include the step of combining sequences of basic fragments to regenerate an entire, megabase plus sequence. This can be accomplished using the link-up of fragments of the basic library which is a prerequisite for informative fragment generation. Alternatively it can be accomplished after generation of sequences of fragments of the basic library by this algorithm, using search for their overlap, based. on the presence of common end-sequences.
- the algorithm requires neither knowledge of the number of appearances of a given K-tuple in a nucleic acid sequence of the basic and ordering libraries, nor does it require the information of which k-tuple words are present on the ends of a fragment.
- the algorithm operates with the mixed content of k-tuples of various length.
- the concept of the algorithm enables operations with the k-tuple sets that contain false positive and false negative k-tuples. only in specific cases does the content of the false k-tuples primarily influence the completeness and correctness of the generated sequence.
- the algorithm can be used for optimization of parameters in simulation experiments, as well as for sequence generation in the actual SB2I experiments e.g. generation of the genomic DNA sequence. In optimization of parameters, the choice of the oligonucleotide probes (k-tuples) for practical and convenient fragments and/or the choice of the optimal lengths and the number of fragments for the defined probes are especially important.
- the algorithm is implemented on an IBM PC/AT compatible computer running under MSDOS 3.30 System works on the clock frequency-of-6 or 10 MHz.
- the workstation is equipped with one 360 kb-1.2 Mb floppy disk drive, and a 40 Mb hard disk, formatted as two 6Mb and 32 Mb partitions.
- This part of the algorithm has a central role in the process of the generation of the sequence from the content of k-tuples. It is based on the unique ordering of k-tuples by means of maximal overlap. The main obstacles in sequence generation are specific repeated sequences and false positive and/or negative k-tuples. The aim of this part of the algorithm is to obtain the minimal number of the longest possible subfragments, with correct sequence. This part of the algorithm consists of one basic, and several control steps. A two-stage process is necessary since certain information can be used only after generation of all primary subfragments.
- the main problem of sequence generation is obtaining a repeated sequence from word contents that by definition do not carry information on the number of occurrences of the particular k-tuples.
- the concept of the entire algorithm depends on the basis on which this problem is solved.
- pSFs contain an excess of sequences and in the second case, they contain a deficit of sequences.
- the first approach requires elimination of the excess sequences generated, and the second requires permitting multiple use of some of the subfragments in the process of the final assembling of the sequence.
- k-tuple X is unambiguously-maximally overlapped with k-tuple Y if and only if, the rightmost K-1 end of k-tuple X is present only on the leftmost end of k-tuple Y and if the leftmost k-1 end of k-tuple Y is not present on the rightmost end of any other k-tuple.
- the algorithm based on the stricter rule is simpler, and is described herein.
- FIG. 10 show the minimal conditions for unique overlap, or conversely stopping the process of elongation of a given subfragment.
- the process is stopped when the right k-1 end of the last k-tuple included is not present on the left end of any k-tuple (FIG. 10. 1 ) or is present on two or more k-tuples (FIG. 10. 2 ). If it is present on only one k-tuple the second part of the rule is tested (FIG. 10. 3 and 10 . 4 ). If in addition there is a k-tuple which differs from the previously included one, the assembly of the given subfragment is terminated only on the first leftmost position (FIG. 10. 3 ). If this additional k-tuple does not Sexist, the conditions are met for unique k-1 overlap and a given subfragment is extended to the right by one element.
- a supplementary one is used to allow the usage of k-tuples of different lengths.
- the maximal overlap is the length of k-1 of the shorter k-tuple of the overlapping pair.
- Generation of the pSFs is performed starting from the first k-tuple from the file in which k-tuples are displayed randomly and independently from their order in a nucleic acid sequence.
- the first k-tuple in the file is not necessarily on the beginning of the sequence, nor on the start of the particular subfragment.
- the process of subfragment generation is performed by ordering the k-tuples by means of unique overlap, which is defined by the described rule. Each used k-tuple is erased from the file.
- the building of subfragment is terminated and the buildup of another pSF is started. Since generation of a majority of subfragments does not begin from their actual starts, the formed pSF are added to the k-tuple file and are considered as a longer k-tuple. Another possibility is to form subfragments going in both directions from the starting k-tuple. The process ends when further overlap, i.e. the extension of any of the subfragments, is not possible.
- the pSFs can be divided in three groups: 1) Subfragments of the maximal length and correct sequence in cases of exact k-tuple set; 2) short subfragments, formed due to the use of the maximal and unambiguous overlap rule on the incomplete set, and/or the set with some false Positive k-tuples; and 3) pSFs of an incorrect sequence.
- the incompleteness of the set in 2) is caused by false negative results of a hybridization experiment, as well as by using an incorrect set of k-tuples. These are formed due to the false positive and false negative k-tuples and can be: a) misconnected subfragments; b) subfragments with the wrong end; and c) false positive k-tuples which appears as false minimal subfragments.
- FIG. 11 Several more common examples of obtaining short subfragments with or without a wrong base at the end because of errors in the k-tuple set are shown in FIG. 11.
- pSFs are generated because of the impossibility of maximal overlapping (FIG. 10.1).
- pSFs are generated because of the impossibility of unambiguous overlapping (FIGS. 10, 11 and 12 ).
- pSFs are generated, and one of these pSFs contains the wrong k-tuple at the relevant end (Fig. 11C).
- the second possibility requires the presence of both false negative and false positive k-tuples in The file.
- the second possibility requires the presence of both false negative and false positive k-tuples in The file.
- the lack of only one k-tuple is sufficient when either end is repeated internally as shown in FIG. 12-A and B.
- the lack of two is needed for strictly internal repetition (FIG. 12-C and D).
- the reason is that the end of a sequence can be considered informatically as an endless linear array of false negative k-tuples. From the “smaller than K-1 case”, only the repeated sequence of the length of k-2, which requires two (FIG. 12-E) or three (FIG. 12-F) specific erroneous k-tuples, will be considered. It is very likely that these will be the only cases which will be detected in a real experiment, the others being much less frequent.
- a fragment beside at least two identical k-1 sequences, contains any k-2 sequence from k-1 or a fragment contains k-2 sequence repeated at least twice and at least one false negative k-tuple containing given k-2 sequence in the middle, etc.
- the aim of this part of the algorithm is to reduce the number of pSFs to a minimal number of longer subfragments with correct sequence.
- the generation of unique longer subfragments or a complete sequence is possible in two situations.
- the first situation concerns the specific order of repeated k-1 words.
- the merging of subfragments that can be uniquely connected is accomplished in the second step.
- the rule for connection is: two subfragments can be unambiguously connected if, and only if, the overlapping sequence at the relevant end or start of two subfragments is not present at the start and/or end of any other subfragment.
- the short subfragments that are obtained are of two kinds. In the common case, these subfragmentss can be unambiguously connected among themselves because of the distribution of repeated k-1 sequences. This can be done only after the process of generation of pSFs and is a good example of the necessity for two steps in the process of pSF generation.
- short pSFs are obtained on the sites of nonrepeated K-1 sequences. Considering false positive k-tuples, there is possibility of presence of the k-tuple containing more than one wrong base (or containing one wrong base somewhere in the middle), as well as k-tuple with wrong base on the end. Generation of short and erroneous (or misconnected) subfragments is caused by the latter k-tuples.
- the k-tuples of the former kind represent wrong pSFs with length equal to k-tuple length.
- the aim of merging pSF part of the algorithm is the reduction of the number of pSFs to the minimal number of longer subfragments with the correct sequence. All k-tuple subfragments that do not have an overlap on either end, of the length of longer than k-a on one, and longer than k-b on the other end, are eliminated to enable the maximal number of connections. In this way, the majority of false positive k-tuples are discarded.
- the rule for connection is: two subfragments can be unambiguously connected if, and only if the overlapping sequence of the relevant end or start of two subfragments is not present on the start and/or end of any other subfragment. The exception is a subfragment with the identical beginning and end.
- connection is permitted, provided that there is another subfragment with the same end present in the file.
- the main problem here is of precise definition of overlapping sequence.
- the presence of at least two specific false negative k-tuples on the points of repetition of k-1 or k-2 sequences, as well as combining of the false positive and false negative k-tuples can destroy or “mask” some overlapping sequences and can produce an unambiguous, but wrong connection of pSFs.
- completeness must be sacrificed on account of exactness: the connection is not permitted on the end-sequences shorter than k-2, and in the presence of an extra overlapping sequence longer than k-4.
- the overlapping sequences are defined from the end of the pSFs, or omitting one, or few last bases.
- the process of ordering of subfragments is similar to the process of their generation. If one considers subfragments as longer k-tuples, ordering is performed by their unambiguous connection via overlapping ends.
- the informational basis for unambiguous connection is the division of subfragments generated in fragments of the basic library into groups representing segments of those fragments.
- the method is analogous to the biochemical solution of this problem based on hybridization with longer oligonucleotides with relevant connecting sequence.
- the connecting sequences are generated as subfragments using the k-tuple sets of the appropriate segments of basic library fragments.
- Relevant segments are defined by the fragments of the ordering library that overlap with the respective fragments of the basic library.
- the shortest segments are informative fragments of the ordering library.
- the longer ones are several neighboring informative fragments or total overlapping portions of fragments corresponding of the ordering and basic libraries.
- fragments of the ordering library are randomly pooled, and the unique k-tuple content is determined.
- the process of ordering consists of three steps: (1) generation of the k-tuple contents of each segment; (2) generation of subfragments in each segment; and (3) connection of the subfragments of the fragment of the basic library using subfragments of the segments.
- Primary segments are defined as significant intersections and differences of k-tuple contents of a given fragment of the basic library with the k-tuple contents of the pools of the ordering library.
- Secondary (shorter) segments are defined as intersections and differences of the k-tuple contents of the primary segments.
- the set of false negatives from the original sequences are enlarged for false positives from intersections and the set of false positives for those k-tuples which are not included in the intersection by error, i.e. are false negatives in the intersection. If the starting sequences contains 10% false negative data, the primary and secondary intersections will contain 19% and 28% false negative k-tuples, respectively. On the other hand, a mathematical expectation of 77 false positives can be predicted if the basic fragment and the pools have lengths of 500 bp and 10,000 bp, respectively. However, there is a possibility of recovering most of the “lost” k-tuples and of eliminating most of the false positive k-tuples.
- intersection a set of common k-tuples—is defined for each pair (a basic fragment) ⁇ (a pool of ordering library). If the number of k-tuples in the set is significant it is enlarged with the false negatives according to the described rule.
- the primary difference set is obtained by subtracting from a given basic fragment the obtained intersection set.
- the false negative k-tuples are appended to the difference set by borrowing from the intersection set according to the described rule and, at the same time, removed from the intersection set as false positive k-tuples.
- this difference can represent the two separate segments which somewhat reduces its utility in further steps.
- the primary segments are all generated intersections and differences of pairs (a basic fragment) ⁇ (a pool of ordering library) containing the-significant number of k-tuples.
- k-tuple sets of secondary segments are obtained by comparison of k-tuple sets of all possible pairs of primary segments. The two differences are defined for each pair which produces the intersection with the significant number of k-tuples. The majority of available information from overlapped fragments is recovered in this step so that there is little to be gained from the third round of forming intersections and differences.
- connection of subfragments consists of sequentially determining the correctly linked pairs of subfragments among the subfragments from a given basic library fragment which have some overlapped ends.
- 4 relevant subfragments two of which contain the same beginning and two having the same end, there are 4 different pairs of subfragments that can be connected.
- 2 are correct and 2 are wrong.
- the length and the position of the connecting sequence are chosen to avoid interference with sequences which occur by chance. They are k+2 or longer, and include at least one element 2beside overlapping sequence in both subfragments of a given pair.
- the connection is permitted only if the two connecting sequences are found and the remaining two do not exist.
- the two linked subfragments replace former subfragments in the file and the process is cyclically repeated.
- Lsf is the length of subfragment
- Ck is the number of common k-tuples for a given subfragment/segment pair
- F is the parameter that includes relations between lengths of k-tuples, fragments of basic library, the size of the pool, and the error percentage.
- Subfragments attributed to a particular segment are treated as redundant short pSFs and are submitted to a process of unambiguous connection.
- the definition of unambiguous connection is slightly different in this case, since it is based on a probability that subfragments with overlapping end(s) belong to the segment considered.
- the accuracy of unambiguous connection is controlled by following the connection of these subfragments in other segments. After the connection in different segments, all of the obtained subfragments are merged together, shorter subfragments included within longer ones are eliminated, and the remaining ones are submitted to the ordinary connecting process. If the sequence is not regenerated completely, the process of partition and connection of subfragments is repeated with the same or less severe criterions of probability of belonging to the particular segment, followed by unambiguous connection.
- the described algorithm is tested on a randomly generated, 50 kb sequence, containing 40% GC to simulate the GC content of the human genome. In the middle part of this sequence were inserted various Alu, and some other repetitive sequences, of a total length of about 4 kb. To simulate an in vitro SBH experiment, the following operations were performed to prepare appropriate data.
- Positions of 1000 500 bp “clones” were randomly determined, to simulate making the ordering library. These fragments were extracted from the sequence. Random pools of fragments were made, and k-tuple sets of pools were determined and stored on the hard disk. These data are used in the subfragment ordering phase: For the same density of clones 4 million clones in basic library and 3 million clones in ordering library will be needed for the entire human genome. The total number of 7 million clones is several fold smaller than the number of clones a few kb long that are needed for random cloning of almost all of genomic DNA and its sequencing by a gel-based method.
- Elongation of the end of the sequence caused by false positive k-tuples is the special case of “insertions”, since the end of the sequence can be considered as the endless linear array of false negative k-tuples.
- An insertion,-or insertion in place of a deletion can arise as a result of specific combinations of false positive and false negative k-tuples. In the first case the number of consecutive false negatives has to be smaller than k. Both cases require several overlapping false positive k-tuples.
- the insertions and deletions are mostly theoretical possibilities without sizable practical repercussions since the requirements in the number and specificity of false k-tuples are simply too high.
Abstract
The conditions under which oligonucleotide probes hybridize preferentially with entirely complementary and homologous nucleic acid targets are described. Using these hybridization conditions, overlapping oligonucleotide probes associate with a target nucleic acid. Following washes, positive hybridization signals are used to assemble the sequence of a given nucleic acid fragment. Representative target nucleic acids are applied as dots. Up to to 100,000 probes of the type (A,T,C,G) (A,T,C,G)N8(A,T,C,G) are used to determine sequence information by simultaneous hybridization with nucleic acid molecules bound to a filter. Additional hybridization conditions are provided that allow stringent hybridization of 6-10 nucleotide long oligomers which extends the utility of the invention. A computer process determines the information sequence of the target nucleic acid which can include targets with the complexity of mammalian genomes. Sequence generation can be obtained for a large complex mammalian genome in a single process.
Description
- The present invention belongs to the field of molecular biology. It involves a novel method of sequencing of a target nucleic acid sequence by hybridization of short oligonucleotide probes to a nucleic acid target. The oligonucleotide probes can comprise all known combinations of the four nucleotides of a given length, i.e. oligonucleotides of base composition adenine (A), thymine (T), guanine (G), and cytosine (C) for DNA and A, G, C, and uridine (U) for RNA. Conditions are described which allow hybridization discrimination between oligonucleotides which are as short as six nucleotides long and have a single base end-mismatch with the target sequence.
- The invention is demonstrated by way of examples in which sequence information is generated using the method of the invention.
- Hybridization depends on the pairing of complementary bases in nucleic acids and is a specific tool useful for the general recognition of informational polymers. Diverse research problems using hybridization of synthetic oligonucleotide probes of known sequence include, amongst others, the different techniques of identification of specific clones from cDNA and genomic libraries; detecting single base pair polymorphisms in DNA; generation of mutations by oligonucleotide mutagenesis; and the amplification of nucleic acids in vitro from a single sperm, an extinct organism, or a single virus infecting a single cell.
- It is possible to discriminate perfect hybrids from those hybrids containing a single internal
mismatch using oligonucleotides 11 to 20 nucleotides in length [Wallace et al., Nucl. Acids Res. 6: 3543 (1979)]. Mismatched hybrids are distinguished on the basis of the difference in the amount of hybrid formed in the hybridization step and/or the amount remaining after the washing steps [Ikuta et al., Nucl. Acids Res. 15: 797 (1987); Thein and Wallace, in Human Genetic Diseases: A Practical Approach, ed. by J. Davies, IRL Press Ltd., Oxford, pp. 33-50 (1986)]. - The reproducible hybridization of different and diverse short oligonucleotides less than 11 nucleotides long has not been well characterized previously. Detailed hybridization data that allows a constant set of conditions for all predictable oligonucleotides is not available [Besmer et al., J. Mol. Biol. 72: 503 (1972); Smith, in Methods of DNA and RNA Sequencing, ed. S. Weissman, Praeger Publishers, New York, N.Y., pp. 23-68 (1983); Estivill et al., Nucl. Acids Res. 15: 1415 (1987).
- Information is also not available on the effects of a single noncomplementary base pair located at the 5′ or 3′ end of a hybridizing oligonucleotide that produces a mismatched hybrid when associated with a target nucleic acid. Hybridization conditions that discriminate between (1) a perfectly complementary hybridizing pair of nucleic acid sequences where one partner of the pair is a short oligonucleotide, and (2) a pair wherein a mismatch of one nucleotide occurs on the 5′ or 3′ end of the oligonucleotide, provide a more stringent environment than is required for internal mismatches because hybrid stability is affected less by a mismatch at the end of a hybridizing pair of complementary nucleic acids than for an internal mismatch.
- The length of nucleotides that can distinguish a unique sequence in a nucleic acid of defined size has been predicted (Smith in Methods of DNA and RNA Sequencing, ed. S. Weissman, Praeger Publishers, New York, N.Y., pp. 23-68 (1983)]. Thus random oligonucleotide sequences 16-17 long are expected to occur only once in random DNA of 3×109 bp, the size of the human genome. However, with decreasing probe length, e.g. for
oligonucleotides 5 to 10 nucleotides in length, there is an exponential increase in the frequency of occurance within a random DNA of a given size and complexity. Thus, the purposes for which oligonucleotide probe ore employed can impact on the length of the oligonucleotides that are used experimentally. - Wallace et al. [Nucl. Acids Res. 6: 3543 (1979)] describe conditions that differentiate the hybridization of 11 to 17 base long oligonucleotide probes that match perfectly and are completely homologous to the target nucleic acid as compared to similar oligonucleotide probes that contain a single internal base pair mismatch. Wood et al. [Proc. Natl. Acad. Sci. 82: 1585 (1985)] describe conditions for hybridization of 11 to 20 base long oligonucleotides using 3M tetramethyl ammonium chloride wherein the melting point of the hybrid depends only on the length of the oligonucleotide probe, regardless of its GC content. However, as disclosed in these references eleven mer oligonucleotides are the shortest ones that generally can be hybridized successfully, reliably and reproducibly using known hybridization conditions.
- Nucleic acid sequencing methods, where the position of each base in a nucleic acid molecule in relation to its neighbors is determined to define its primary structure, were developed in the early 1960's for RNA molecules and in the late 1970's for DNA. The two major methods for DNA sequencing, i.e. chemical degradation and dideoxy-chain termination, involve identification and characterization of 1-500 nucleotide long DNA fragments, specific for each one of at least four nucleotide bases, on polyacrylamide gels. The polyacrylamide gels must be able to distinguish single base pair differences in length between fragments. The fragments are generated either by chemical degradation [Maxam-Gilbert, Proc. Natl. Acad. Sci. 74: 560 (1977)] or by dideoxy-chain termination of DNA fragments synthesized by DNA polymerase [Sanger et al., Proc. Natl. Acad. Sci. 74: 5463 (1977)]. A sufficient quantity of isolated fragments is ensured by recombinant DNA technology methods which include cloning, restriction enzyme digestion, gel electrophoresis, and polymerase chain reaction amongst others. These methods allow the identification and amplification of the target DNA to provide material for sequencing.
- An intensive amount of manual labor is required in the preparation of appropriate polyacrylamide gels to resolve small differences in fragment size. The speed of sequencing in experienced laboratories throughout the world is approximately 100 bp per person daily. Although the use of electronic robots and computers allows acceleration of the number of base pairs actually determined, preparation of polyacrylamide gels, application of sample, electrophoresis and the subsequent manipulations necessary to obtain high quality autoradiograms that can be read by machines still involve significant intensive, skilled, manual labor for which no substitutes have been found.
- The genome of higher eucaryotes has up to a million times greater physical complexity than is the complexity of individual genes it encodes, giving it a a corresponding huge informational complexity. From the present knowledge of genome organization and biochemical, biophysical and biological functions, the following approximate scale of the informational complexity for higher eucaryotes can be proposed: 10,000 gene families—100,000 genes—1,000,000 biological functions. The number of basic biochemical functions represented by a single gene family is probably not significantly increased compared to procaryotic and lower eucaryotic genomes.
- Recently, there has been a surge of interest in mapping and sequencing the entire human genome [Lewin, Science 232: 1598 (1986); Wada, Nature 325: 771 (1987); Smith and Hood, Bio/Technology 5: 933-939 (1987)]. This stems from the fact that only 1 in about 75 human genes is either cloned or mapped (Human Gene Mapping 9, 1987). Unknown genes will have much to tell us about human biology. In the future, the progress of studies on molecular evolution may depend on the sequencing of genomes of Species besides humans.
- Because sequence information has already provided accelerated knowledge and potential resolution of diverse biological, medical and therapeutic research problems, it is not surprising that ideas of sequencing the whole human genome were discussed at various scientific meetings during the early and mid-1980's [Research News in Science 232: 1598 (1986)]. Such massive sequencing projects envision the final determination of approximately 3 billion base pairs of information encoded in the DNA of humans and are expected to take at least 10 years at a cost of at least $3 billion dollars using current technology. However, in practice, actual sequencing of at least three times that number of base pairs is required to obtain a reliable sequence for the human genome, thus requiring even more money and time.
- Such endeavors present a challenge to the technology of the twentieth century. Further challenges arise if sequencing projects are extended to include the determination of the genomic sequences of characteristic individuals or species of organisms, especially those that have economic, social or medical importance. Such sequencing projects would advance not only our understanding of the evolution of organisms and the evolution of biochemical processes, but would also further the detection, treatment and understanding of disease, and would aid agriculture, the food industry and biotechnology in general. However beneficial the results of such projects would be, their successful completion requires the development of a new, rapid, reproducible and reliable sequencing method such as those described in this invention.
- Although the ultimate goal of human genome characterization is the determination of sequence information, progress in characterizing portions of the human genome or the genome of other organisms have been achieved in several areas. A linkage map of the human genome based on cloned DNA probes detecting RFLPs has been obtained [Donis-Keller et al., Cell 51: 319-337 (1987)]. Once mapped, a gene can be approached from a neighboring DNA marker not only by walking [Cross et al., Trends Genet. 2: 174 (1986)] but also by the use of jumping [Collins+Weissman, Proc. Natl'l. Acad. Sci. USA 81: 6812 (1984); Poustka et al., Nature 325: 353 (1987)] and linking [Poustka et al., Trends Genet, 2: 174 (1986)] libraries. The task of going from a marker to a mapped gene is facilitated immensely if an ordered collection of overlapping cosmid or phage clones representing individual chromosomes is available. Attempts to provide a library of overlapping clones using similarities in their patterns of restriction digests have been tried [Coulson et al., Proc. Natl. Acad. Sci. USA 83: 7821 (1986); Olson et al., Proc. Natl. Acad. Sci. USA 83: 7826 (1986); Kohara et al., Cell 50: 495 (1987)]. Alternatively, the hybridization of a collection of 100 specific oligonucleotides to an array of 3-10×106 cosmid-containing colonies on filters has been proposed. The resulting patterns of hybridization identify specific regions along the genome to which a small collection-of cosmids from chromosome libraries can be fitted in the second step [Poustka et al., Cold Spring Harbor Symp. Quant. Biol. 51: 351 (1986); Craig et al., in Human Genetics, Proceedings of the 7th International Congress, Berlin, (1986); Michiels et al., CABIOS 3: 203 (1987)]. Such identification however does not provide desired and useful sequence information of the DNA in a particular identified fragment.
- In the area of human genetics, the emphasis is on an individual's DNA and the methods to detect patterns of its variation and inheritance which may influence the determination of a patient's chances for health or disease. The number of genetic regions to be scored in the DNA of an individual requires a large number of polymorphic probes and makes the use of traditional Southern blotting unpractical. However, a method that is capable of amplifying 1000-bp stretches of DNA starting from two flanking oligonucleotide primers and that requires DNA from only 150 cells of an individual has been described recently as well as oligonucleotide probes that can detect mutants in amplified DNA in dot blot hybridization [Saiki et al., Science 239: 487 (1986)]. Both the method of ordering cosmid libraries and the method of amplifying DNA use the work of Wallace for conditions of hybridization that only allowed oligonucleotides of almost perfect homology to their target DNA to hybridize at all [Wallace et al., Nucl. Acids Res.: 3543 (1979)]. In these conditions, almost perfect homology means that the perfect homology has to exist at least in the central part of the hybridizing oligonucleotide/target duplex.
- The present invention provides a new method of sequencing that is ideally suited to the sequencing of large complex genomes because it avoids the intensive manual labor involved in resolving gel fragments by size on polyacrylamide gels. The present invention provides methods for sequencing a target nucleic acid by hybridization of overlapping short oligonucleotide probes of known or predicted sequence to the nucleic acid target serially or simultaneously. The oligomer probes of a given size can contain all or most existing combinations of nucleotides for complete sequencing and a part of all possible variants for partial sequencing. Probes can also be composed of oligomers of different sizes as well as comprising all known combinations of nucleotides that are possible for that size oligonucleotide. As the size of the probes that are used decreases, hybridization conditions that are still able to distinguish between mismatched and perfectly matched short oligonucleotides must be used.
- In one embodiment of the invention, multiple oligonucleotides that are 11 nucleotides long or longer are hybridized to the target sequence. hybridization occurs using conditions which are controlled and varied to ensure discrimination between perfectly matched oligonucleotides and oligonucleotides having a one base pair mismatch with the target sequence where the mismatch is located at either one of two ends of the oligonucleotide.
- In another embodiment of the invention, as an alternative to previous numerous conditions each specific for different sizes and sequences of probes, a single, or few, sets of conditions is invented for all lengths and sequence of probes. These hybridization conditions allow discrimination between perfectly matched and mismatched oligonucleotides that are as short as six nucleotides long. The conditions allow discrimination between a perfectly matching oligonucleotide and one that has a single base mismatch as compared to the target sequence, the mismatch being located at one of the ends of the oligonucleotide.
- Following the detection of hybridization of perfectly matched oligonucleotides of known sequence, the sequence of the target nucleic acid is generated by an algorithm using the principle of maximal nonidentical overlap of probe.
- In determining sequence by hybridization, oligonucleotides are prepared, target fragments are prepared appropriate for the length of oligonucleotide used for hybridization, and hybridization of the target with all the oligonucleotides occurs under defined conditions that allow discrimination in binding of perfectly matched complementary oligonucleotides and mismatched oligonucleotides. The relationship of probe size and target length is defined and allows complete sequencing of genomes. The novel theoretical basis of the relationship between oligonucleotide probe size and target length is described infra.
- To determine the amount of hybridization data that is needed for sequence determination, the number of target fragments that compose the entire sequence is multiplied by the number of different oligonucleotides required to define the sequence of the target fragment. The shorter the size of the oligonucleotides that are hybridized, the more target fragments that must be analyzed. Similarly, as the oligonucleotide size increases, fewer target fragments must be examined.
- Hybridization reactions can be performed in separate reaction vessels or by binding one of the two components (oligomers and DNA fragments) to a solid surface, like nylon filters etc. Since the described method does not require macromolecular separation like gel-based sequencing methods, the surface, bound with either an oligomer or nucleic acid fragment can have microdimensions.
- Some of the advantages of the method of the present invention include the following:
- (1) rapidity, resulting in time effectiveness: (2) elimination of polyacrylamide gel electrophoresis and the intensive manual labor it requires; (3) reliability of the predicted base within the determined sequence due to the hybridization of multiple oligonucleotides to the same base within a target sequence; (4) the possibility of substantial miniaturization of the process; (5) ease of automation; (6) resulting cost effectiveness.
- The following terms and abbreviations will have the meanings indicated:
A adenine bp base pair C cytosine G guanine IF an M13 clone containing a 921 bp EcoR1-Bg1II human β1 interferon fragment kD kilo Dalton nG nanogram nM nenomolar pmol picomole sc subclone SF subfragment SOH short oligonucleotide hybridization T thymine CCD Charge Coupled Device DNA Deoxyribonucleic acid DP Discrete particle HA Hybridization area LAR Ligation-amplification reaction ON Oligonucleotide ONP Oligonucleotide probe ONS Oligonucleotide sequence PCR Polymerase chain reaction RE Restriction Enzyme RFLP Restriction fragment length polymorphism RNA Ribonucleic acid SBH Sequencing by hybridization - FIG. 1 shows the generation of subfragments in sequence by hybridization and their ordering. (A) The sequence of a hypothetical clone. NNNNNNN-ends of vector sequence. AGTCCCT and TTGGCTG are the
only oligonucleotides 7 bp or longer repeated within the depicted sequence. (B) Formation of subfragments. Assuming that the content of 8-mers for the depicted sequence is known, these 8-mers are ordered by maximal overlap, in thiscase 7 bp. Starting from the 5′ 8-mer (NNNNNNNc), ordering is unambiguous up to gAGTCCCT, which on its 3′ end contains a repeated 7-mer. The large capitals denote overlapping sequences shared by different oligonucleotides, while the small letters denote unshared bases. Both AGTCCCTc and AGTCCCTg can be overlapped with gAGTCCCT preventing further ordering. Each of the two sequences serves as a starting point for new ordering (not shown). Therefore, each repeatedsequence 7 bp or longer represents a branching point. Unambiguous sequences are obtained between two consecutive branching points only. (C) Listing of subfragments formed from 8-mers of depicted sequence. Subfragments are horizontally displaced to indicate overlap; the orientation is 5′ to 3′ and end subfragments are identifiable. (D) The subfragments cannot be unambiguously ordered into a starting sequence without additional information. Both arrangements shown are possible since subfragments AGTCCCTcggTTGGCTG and AGTCCCTgatTTGGCTG have the same 7-mers at their 5′ and 3′ ends, respectively. (E) Building the sequence from oligonucleotide blocks. The left box represents all 8-mer oligosequences which occur in 15 base long DNA molecule of unknown sequence (NNN . . . NNN), 8-mers can be ordered by 7 base overlap (right box). Each 7 oligomer extends the sequence of the starting 8-mer ACCGTAAA by one base. Thus, the sequence is generated by uniquely overlapped oligomer blocks. - FIG. 2 presents the average number of SFs (Nsf) as a function of the length of DNA fragment (Lf) for various values of the length of the overlapping sequence (N-1, in bp), or average distance of two consecutive identical N-1 sequences in DNA subjected to sequencing by hybridization (Ao), in kb. The curves are obtained using equation one as described below in section 5.2.
- FIG. 3 describes the kinetic stability of a fully matched hybrid obtained with a
probe 8 nucleotides in length. Stability is expressed as a fraction of the hybrid dissociated in unit time (minutes) as a function of temperature. 1.4 pmol of NCATGAGCANN as applied to each dot and hybridized with TGCTCATG as probe in a concentration of 4 nM. The equal amounts of hybrid were incubated at the indicated temperatures for a short time in a large volume of buffer and the remaining hybrid measured. Each point represents the average value for four dots. The curve is computer fitted with Eα=47.3 Kcal/mol obtained from the experimental points by the least squares method. - FIG. 4 indicates the properties of short oligonucleotide hybridization. In FIG. 4a, non-optimized discrimination with
probes - In FIG. 4b, limits of signal detection are examined. The indicated volumes of IF culture supernatants of average titer of 6×1011 pfu/ml were mixed with an equal volume of 1 M NaOH, 3M NaCl and spotted on a filter as described in a above. Hybridization was at 2° C. with TGCTCATG as the probe.
- In FIG. 4c, the time course of hybridization at 13° C. is shown. The IF-M13 system was used with 50 ng of phage DNA per dot, and the probe was TGCTCATG. The 3 hr IF dot contained 18020 cpm measured with 20% efficiency.
- In FIG. 5 the effect of the washing step on discrimination is indicated. In FIG. 5a, inversion of the signal in IF-M13 pair upon washing is shown. 10 ng of IF and 500 ng M13 DNA were applied, and the probe was TGCTCATG. The top row was not washed, the other rows were washed at 7, 13 and 25° C., respectively for the indicated times. A DNA control is included in the top row also. Hybridization with the M13 specific probe AGCTGCTC measures amounts of DNA in the two dots. In FIG. 5b, the change of discrimination with time of washing at 0° C. (full circles) and 13° C. (open circles) is depicted. 100 ng each of IF and M13 were applied to form dots. The dots were hybridized to probe TGCTCATG and probe AGCTGCTC was used in the control DNA hybridization (see top row, on the right, panel a). The dots were then washed at the indicated temperatures. At each time point the pairs of dots were removed and the ratio of radioactivity remaining in the each dot was measured. The D or discrimination was calculated as the mean value of the ratios for the duplicate pairs of dots.
- FIG. 6 demonstrates the effects of complexity of target sequences on discrimination. 50 ng each of IF and M13 were hybridized with the indicated probes at a concentration of 4 nM. No wash was performed. The number of matched and end base mismatched targets in IF and M13 is indicated for each probe.
- FIG. 7 examines an array of clones for the presence of an oligonucleotide sequence. 51 recombinant plasmid DNAs (10+5 ng) were spotted in rows B to II,
columns 1 to 8 (except row H). Line A andcolumn 9 contained control DNAs of known sequence. Unknown clones were taken from human brain cDNA library in Bluescript vector (BS) (Stratogene Cat. No. 935205). Controls of known sequence in lines A1 to A8 and A9 to G9 are: IF(M13), M13, Alu(M3), IF(BS), BS, 1M(pUC 9),pUC 9, 2M(pUC 9), respectively except that in the vertical row Alu(M13) was omitted. 1M and 2M are rat β-globin gene subclones. The probe concentration was 8 nM. In FIG. 7a, BS specific probe CTCCCTTT was also contained in IF and 2M inserts but not in M13 and pUC vectors. In FIG. 7b, the sequence of probe CCAGTTTT was contained in the IF insert but not in either vector. In FIG. 7c, the sequence of probe GCCTTCTC was contained in the 1M insert only. - FIG. 8
sequencing 100 bp of 921 bp β1-human interferon gene fragment. (IF) by hybridization. -
Part 1. Hybridization results. A. Hybridization with 93 probes (72) octamers and 21 nonamers with the full match in IF. IF and controls rat globin clones pHEA and pHI were PCR amplified while M13 mp18 and pUC18 were in linearized double stranded form. Base denatured DNA (20 ng of IF and equimolar amounts of control DNA) were spotted on Gene Screen membranes (N.E.N.). Hybridization was according to Drmanac et al., described in § 6 below. Briefly, γ-32 end labeled probes (3.3 pm, 10 mCI, Amersham 3000 C/mM in concentration of 10 ng/ml were hybridized at 12° C. in 0.5 M Na2HPO4 pH 7.2, 7% Na-lauryl Sarcosine for 3 hours. All probes were made by Genesys, Inc., Houston. Hybrids were washed in 6×SSC at 0° for 40 minutes and autoradiographed for 4-48 hours. Test dot signal intensity, Hp, and discrimination as ratio of signals of test over control dot, D, were visually estimated. For probes 34 and 74, dot radioactivity was measured in a scintillation counter. Hp was 6,000 and 300 cpm, D was 20 and 4, and a film was exposed for 4 and 48 hr respectively. B. Hybridization with 12 probes (11 octamers and 1 nonamer) which have end mismatch in IF fragment. Control DNAs having single full match targets were pHEA for probes 97., 98., 102., pUC18 for 95., 100., 104., 105., and M13 for 94., 96., 99., 101., and 103. Probes 104 and 105 have 3 end-mismatched targets in IF. Hybridization procedures were as described in A. C. DNA Calibration. 1. and 2. IF and pHEA, probe CTGATAT. 3. IF and pUC18 probe CAGATGGT. 4 IF and M13mp18, probe GACTGTCT. The ratios of DNA amounts in IF and control dot were 1:1 in panel 1., 3., 4., and 1:3 in panel 2., respectively. Filters with IF and pH had 1:2 ratio with probe CTGATGAT. Filters show in 2. were used with probes 1., 3., 4., 6. to 8., 10. to 13., containing; pUC18 with probes 31 and 85.; containing M13 with probes 53 and 74; and containing pH with probes 22, 54., 55., 69., 70., 83., and 84. The remainder of probes were used on filters of the type shown onpanel 1. -
Part 2. 10 bp sequence, position 625-726 in Eco RI β1-interferon fragment. The locations of the first 10 probes in the sequence are indicated. -
Part 3. List of probes used. The numbers after the probe sequences are ΔH, Hp. and D respectively. ΔH was calculated Hp and D were visually estimated on the basis of the radioactivity measurements of the selected examples ranging from 300-10,000 cpm per dot. Each symbol after the probe indicates the presence of one end mismatched target in its control DNA dot. - FIG. 9 an example of the discriminative hybridization with 8-mers. On each filter the upper dot contains 15 ng of EcoRI-BgIII 920 bp genomic fragment with interferon beta-1 gene (IF) and the lower dot contains 25 ng of 1500 bp of rat beta-globin gene (GL) prepared by PCR from clones in bluescript vector using universal and reverse sequencing primers. Each probe was hybridized and washed under the same conditions.32p labeled probes (by T4-polynucleotide kinase) at a concentration of 4 nM in 0.5 M Na2HPO4 pH 7.2, 7% Sodium lauroyl sarcosine were hybridized at 12° C. for 3 h. Hybrids were washed in 6×SSC at 2° C. for 40 minutes and autoradiographed for 6h. Panel on top represents hybridization with probe complementary to amplified vector end and reveals the molar ratio of the two target DNAS. First three 8-mer probes are designed to have perfect target in interferon sequence and target with end mismatch in globin sequence, and
probe 4 opposite. Over 5 fold discrimination is obtained with each probe, having 3 to 6 (C+G). In the case of the most stable probe (probe 1) similar discrimination is obtained but both dots are stronger than in case of other probes. - FIG. 10 k-tuple overlap rule. S-start of30subfragment assembly. Filled bars—the identical parts of k-tuples. Empty bars—non-identical parts of k-tuples. 0, 1 non occurrence or occurrence, respectively in k-tuple set. (1-3)—situations violating the rule. (4) situation recognized by the overlapping rule.
- FIG. 11 Examples of generation of short pSFs. Bars are k-tuples. Filled parts of are bars overlapping sequences which length are indicated on the right site. Open bars are correct k-tuples. Bars with dashed line are false negative k-tuples. Bars with (+) are false positive k-tuples; positions of (+) indicate the wrong base. Horizontal lines are extensions of subfragments in the considered site.
- FIG. 12 Process of formation of misconnected subfragments. Column I: Misconnection due to the presence of false negative k-tuples in the set. Column II: Misconnection due to-t~t simultaneous presence of both false positive and false negative k-tuples. Top lines-actual sequences. Immediately below-false k-tuples. Lower lines-generated subfragments representing the sequence. Filled bars-overlapping sequences present at the point of misconnection. Open bar-any other repeated sequences. Oblique lines show the wrong position of a part of a sequence within the misconnected subfragment. Dashed line indicate the position in subfragment of sequence used for overlap, causing the misconnection.
- The method of the invention of sequencing by hybridization is based on an entirely different logical principle than previous methods of sequencing which have been described above in § 2.3. Known methods of sequencing define a nucleic acid sequence in a base by base fashion that is read from the position of DNA fragments in polyacrylamide gels where the fragments are produced by base specific chemical degradation or chain termination techniques.
- In contrast, the method of sequencing of the present invention is a process comprising hybridizing oligonucleotides to a target sequence, detecting the hybridized oligonucleotides, and assembling the sequence of the target from the sequence of the hybridizing oligonucleotides. The method of sequencing according to the present invention is most suitable for determining the sequence of complex nucleic acid fragments, molecules, or genomes and especially those containing one million or more nucleotide base pairs.
- According to the method of the invention, the number of oligonucleotide probes, the lengths and types of target nucleic acid sequences, the numbers of subfragments of the target nucleic acid sequences and the numbers of hybridizations necessary for sequencing have been calculated and show that sequencing by hybridization is more efficient than sequencing by traditional methods.
- A mismatch-free hybridization of oligonucleotides to an unknown target nucleic acid represents, in essence, a sequencing of complementary target. Probability calculations and, in part, computer simulations estimate the types and numbers of oligonucleotides that would have to be synthesized in order to sequence a megabase plus segment of DNA.
- In order to obtain the required amount of experimental data defined as a matrix of (number of fragments-clones)×(number of probes), the number of probes can be reduced depending on the number of fragments used and vice versa. The optimal ratio of the two numbers is defined by the technological requirements of a particular sequencing by hybridization process. The useful probes are those longer than 6 bases, since the shorter ones would require use of unsuitable and unpractical nucleic acid fragments shorter than a few hundred bases long.
- There are two parameters which influence the choice of probe length. The first is the success in obtaining hybridization results that show the required degree of discrimination. The second is the technological feasibility of synthesis of the required number of probes. The requirement of obtaining sufficient hybridization discrimination with practical and useful amounts of target nucleic acid limits the probe length from both sides. It is difficult to obtain a sufficient amount of hybrid with short probes, and to discriminate end mismatch with long probes. There is no evidence for use of probes shorter than 11-mers in the literature, except for very stable ones [Estivill et al., Nucl. Acids Res. 15: 1415 (1987)]. On the other hand, probes longer than 15 bases discriminate end mismatch with difficulty [Wood et al., Proc. Nat'l. Acad. Sci. USA 82: 1585 (1985)].
- One solution for the problem of unstable probes and end mismatch discrimination is the use of a group of longer probes representing a single shorter probe in an informational sense. Groups of sixty-four 11-mers can be used instead of single 8-mers. Every member of the group has a common core 8-mer and one of three possible variations on outer positions with two variations at each end. The probe can be represented as 5′(A,T,C,G)(A,T,C,G)N8(A,T,C,G) 3′. With this type of probe one does not need to discriminate the non-informative end bases (two on 5′ end, and one on 3′ end) since only the internal 8-mer is read. This solution requires the use of higher mass amounts of probes and label in hybridization reactions.
- These disadvantages are eliminated by the use of a single condition of discriminative hybridization for oligomer probes as short as 6-mers (see Examples 6.3 and 6.4). The length of probes are restricted by the technological problem of having a practical number of probes that can be synthesized and utilized. At this time it may be difficult for some to synthesize more than about 100,000 8-mer and some additional 9-mer probes or groups of 11-mers of
types 5′(A,T,C,G)(A,T,C,G)N8(A,T,C,G)3′ and 5′(A,T,C,G)N9(A,T,C,G)3′ by the current methods for synthesis of individual oligonucleotides. - The number of hybridization reactions is dependent on the number of discrete labeled probes. Therefore in the cases of sequencing shorter nucleic acids using a smaller number of fragments-clones than the number of oligonucleotides, it is better to use oligomers as the target and nucleic acid fragments as probes. For example, one needs 20,000 clone. and up to five fold more oligonucleotides (50-100,000) to sequence 10 million bp of DNA with 8-mer oligonucleotides.
- Target nucleic acids which have undefined sequences can be produced as a mixture of representative libraries in a phage or plasmid vector having inserts of genomic fragments of different sizes or in samples prepared by PCR. Inevitable gaps and uncertainties in alignment of sequenced fragments arise from nonrandom or repetitive sequence organization of complex genomes and difficulties in cloning “poisonous” sequences inEscherichia coli. These problems are inherent in sequencing large complex molecules using any method. In this invention such problems are minimized by the choice of libraries and number of subclones used for hybridization. Alternatively such difficulties can be overcome through the use of amplified target sequences, e.g. by PCR amplification, ligation-amplified reactions etc. Because this invention is based on simple biochemical procedures, the method of the invention is inherently easier to automate than existing sequencing methods. The sequence can be derived from simple primary data using only extensive computing.
- Sequencing by hybridization (SBH) of a target nucleic acid can be visualized as consisting of two steps; 1) a process of dissolving the target nucleic acid into all its constituent oligonucleotide N-mers, and 2) the back assembly of N-mers detected by hybridization and assembled by overlap into an extended sequence. In the invention, hybridization of all possible N-mer oligonucleotide probes to the target nucleic acid determines the N-mer oligonucleotide subset contained in the primary sequence of the target nucleic acid and is the first step in the process of the invention. As indicated in FIG. 1a nucleic acid fragment can be dissolved into all constituent oligonucleotides. Positively hybridizing N-mer oligonucleotide probes are ordered and the sequence of the target DNA is determined using (N-1)mer overlapping frames between the oligonucleotide probes.
- Reassembling the sequence of known oligonucleotides that hybridize to the target nucleic acid to generate the sequence of the target nucleic acid cannot be accomplished in some cases. This is because some information is lost if the target nucleic acid is not in fragments of appropriate size (see FIG. 2) in relation to the size of oligonucleotide that is used for hybridizing. The quantity of information lost is proportional to the length of a target being sequenced. However, if sufficiently short targets are used, their sequence can be unambiguously determined.
- The sequence assembly will be interrupted where ever a given overlapping (N-1)mer is duplicated two or more times. Then either of the two N-mers differing in the last nucleotide can be used in extending the sequence. This branching point limits unambiguous assembly of sequence.
- The probable frequency of duplicated sequences that would interfere with sequence assembly which is distributed along a certain length of DNA can be calculated. This derivation requires the introduction of the definition of a parameter having to do with sequence organization: the sequence subfragment (SF). A sequence subfragment results if any part of the sequence of a target nucleic acid starts sand ends with an (N-1)mer that is repeated two or more times within the target sequence. Thus, subfragments are sequences generated between two points of branching in the process of assembly of the sequences in the method of the invention. The sum of all subfragments is longer than the actual target nucleic acid because of overlapping short ends. Generally, surfragments cannot be assembled in a linear order without additional information since they have shared (N-1)mers at their ends and starts. Different numbers of subfragments are obtained for each nucleic acid target depending on the number of its repeated (N-1)mers. The number depends on the value of N-1 and the length of the target.
- Probability calculations can estimate the interrelationship of the two factors which are given in graphical form in FIG. 2. If the ordering of positive N-mers is accomplished by using overlapping sequences of length N-1 or at an average distance of Ao, the N-1 of a fragment Lf bases long is given by equation one:
- Nsf=1+A 0 ×ΣK×P(K, L f)
- Where K greater than or =2, and P(K, Lf) represents the probability of an N-mer occurring K-times on a fragment Lf base long. Also, a computer program that is able to form subfragments from the content of N-mers for any given sequence is described below in Appendix A.
- As shown (FIG. 2) the number of subfragments increases with the increase of lengths of fragments for a given length of probe. Obtained subfragments may not be uniquely ordered among themselves. Although not complete, this information is very useful for comparative sequence analysis and the recognition of functional squence characteristics. This type of information we call partial sequence. Another way of obtaining partial sequence is the use of only a subset of oligonucleotide probes of a given length.
- There is a relatively good agreement between predicted sequence according to theory and a computer simulation for a random DNA sequence. For instance, for N-1=7, [using an 8-mer or groups of sixty-four 11-mers of
type 5′(A,T,C,G)(A,T,C,G)N8(A,T,C,G)3′] a target nucleic acid of 200 bases will have an average of three subfragments. However, because of the dispersion around the mean, a library of target nucleic acid should have inserts of 500 bp so that less than 1 in 2000 targets have more than three subfragments. Thus in an ideal case of sequence determination of a long nucleic acid of random sequence, it is only necessary to use a representative library with sufficiently short inserts of target nucleic acid. For such inserts it is possible to reconstruct the individual target by the method of the invention. The entire sequence of a large nucleic acid is then obtained by overlapping of the defined individual insert sequences. - We found a solution to reduce the need for very short fragments, e.g. 50 bases for 8-mer probes. The basis for the use of longer fragments lies in the use of information contained in the overlapped fragments present in every random DNA fragmentation process like cloning, or random PCR. It is also possible to use pools of short physical nucleic acid fragments. Using 8-mers or 11-mers like 5′(A,T,C,G) (A,T,C,G)N8(A,T,C,G)3′ for
sequencing 1 megabase, instead of needing 20,000 50 bp fragments only 2,100 samples are sufficient. This number consists of 700 random 7 kb clones (basic library), 1250 pools of 20 clones of 500 bp (subfragments ordering library) and 150 clones from jumping (or similar) library. The developed algorithm (see Appendix) can regenerate sequence using hybridization data of these described samples. - oligonucleotides were either purchased from Genetic Designs, Inc. Houston, Texas or made on an Applied Biosystems 381A DNA synthesizer. Most of the probes used were not purified by HPLC or gel electrophoresis. Probes were designed to have both a single perfectly complementary target in IF, a M13 clone containing a 921 bp Eco RI-Bgl II human β1-interferon fragment [Ohno and Tangiuchi, Proc. Natl. Acad. Sci. 74: 4370-4374 (1981)], and at least one target with an end base mismatch in M13 vector itself.
- End labeling of oligonucleotides was performed as described [Maniatis et al., in Molecular Cloning: a Laboratory Manual, Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y. (1982)] in 10 μl containing T4-polynucleotide-kinase (5 units Amersham), 7-32p -ATP (3.3 pM, 10 μCi Amersham 3000 Ci/mM) and oligonucleotide (4 pM, 10 ng). Specific activities of the probes were 2.5-5×109 cpm/nM.
- Single stranded DNA (2 to 4 μl in 0.5 NaOH, 1.5 M NaCl) was spotted on a Gene Screen membrane wetted with the same solution, the filters were neutralized in 0.05 M Na2HPO4 pH 6.5, baked in an oven at 80° C. for 60 min and UV irradiated for 1 min. Then, the filters were incubated in hybridization solution (0.5 M Na2HPO4 pH 7.2, 7% sodium lauroyl sarcosine for 5 min. at room temperature and placed on the surface of a plastic Petri dish. A drop of hybridization solution (10 μl, 0.5 M Na2 HPO4 pH 7.2, 7% sodium lauroyl sarcosine) with a 32 P end labeled oligomer probe at 4 nM concentration was placed over 1-6 dots, overlayed with a square piece of polyethylene (approximately 1×1 cm), and incubated in a moist chamber at the indicated temperatures for 3 hr. Hybridization was stopped by placing the filter in 6×SSC washing solution for 3×5 minute at 0° C. to remove unhybridized probe. The filter was either dried, or further washed for the indicated times and temperatures, and autoradiographed. For discrimination measurements, the dots were excised from the dried filters after autoradiography, placed in liquid scintillation cocktail and counted. The uncorrected ratio of cpms for IF and M13 dots is given as D.
- The conditions reported herein allow very short oligonucleotide hybridization but ensure discriminations between matched and mismatched oligonucleotides that-are complementary to and therefore bind to a target nuclei acid. Factors which influence the efficient detection of hybridization of specific short sequences based on the degree of discrimination (D) between a perfectly complementary target and an imperfectly complementary target with a single mismatch in the hybrid are defined. Experimental tests used dot blot hybridization of twenty-eight probes that were 6 to 8 nucleotides in length to two M13 clones or to model oligonucleotides bound to membrane filters. The principles guiding the experimental procedures are given below.
- Oligonucleotide hybridization to filter bound target nucleic acids only a few nucleotides longer than the probe in conditions of probe excess is a pseudo-first order reaction with respect to target concentration. This reaction is defined by equation two:
- S t /S o =e −Kh op T
- St and So are target sequence concentrations at time t and to, respectively. (OP) is probe concentration. The rate constant for hybrid formation, Kh increases only slightly in the 0° C. to 30° C. range [Porschke and Eigen, J. Mol. Biol. 62: 361 (1971); Craig et al., J. Mol. Biol. 62: 383 (1971)]. Hybrid melting is a first order reaction with respect to hybrid concentration (here replaced by mass due to filter bound state) as shown in equation three:
- H t /H o =e −k m t
- In equation three, Ht and Ho are hybrid concentrations at times t and to, respectively; km is a rate constant for hybrid melting which is dependant on temperature and salt concentration [Ikuta et al., Nucl. Acids Res. 15: 797 (1987); Porschke and Eigen, J. Mol. Biol. 62: 361 (1971); Craig et al., J. Mol. Biol. 62: 383 (1971)]. During hybridization, which is a strand association process, the back melting, or strand dissociation, reaction takes place as well. Thus, the amount of hybrid formed in time t is a result of forward and back reactions. The equilibrium can be moved towards hybrid formation by increasing probe concentration and/or decreasing temperature. However, during washing cycles in large volumes of buffer, the melting reaction is dominant and the back reaction of hybridization is insignificant, since the probe is absent. This analysis indicates workable short oligonucleotide hybridization (SOH) conditions can be varied for probe concentration or temperature.
- D or discrimination is defined in equation four:
- D=H p(t w)/H i(t w)
- Hp (tw) and Hi(tW) are the amounts of hybrids remaining after a washing time, tw, for the identical amounts of perfectly and imperfectly complementary duplex, respectively. For a given temperature, the discrimination D changes with the length of washing time and reaches the maximal value when Hi=B which is equation five.
- The background, B, represents the lowest hybridization signal detectable in the system. Since any further decrease of Hi cannot be examined, D would increase upon continued washing. Washing past tw just decreases Hp relative to B, and is seen as a decrease in D. The optimal washing time, tw, for imperfect hybrids, from equation three and equation five is Equation six:
- t w=−1n(B/H i(t o))/Km.i
- Since Hp is being washed for the same tw, combining Equation three for Hp, Equation four and Equation six, one obtains the optimal discrimination function in Equation seven:
- D=e 1n(B/H i (t o ))k m,p /k m,i ×H p(t o)/B
- The change of D as a function of T is important because of the choice of an optimal washing temperature. It is obtained by substituting the Arhenius equation which is Equation eight:
- K−=Ae −Eα/RT
-
- B less than Hi(to).
- Since the activation energy for perfect hybrids, Eαp, and the activation energy for imperfect hybrids, Eα, i, can be either equal, or Eα,i less than Eα,p, D is temperature independent, or decreases with increasing temperature, respectively. This result implies that the search for stringent temperature conditions for good discrimination in SOH is unjustified. By washing at lower temperatures, one obtains equal or better discrimination, but the time of washing exponentially increases with the decrease of temperature, as can be seen by substituting Equation eight in Equation six. A further property of Equation nine for Eα,i less than Eα,p is that discrimination more strongly decreases with T, if Hi(to) increases relative to Hp(to).
- Finally, the analysis of Equation nine shows that D at lower temperatures depends to a higher degree on the Hp(to)/B ratio than on the Hp(to)/Hi(to) ratio. This result indicates that it is better to obtain a sufficient quantity of Hp in the hybridization regardless of the discrimination that can be achieved in this step. Better discrimination can then be obtained by washing, since the higher amounts of perfect hybrid allow more time for differential melting to show an effect. Similarly, using larger amounts of target nucleic acid a necessary discrimination can be obtained even with small differences between Km,p and Km,i.
- Extrapolated to a more complex situation than covered in this simple model, the result is that washing at lower temperatures is even more important for obtaining discrimination in the case of hybridization of a probe having many end-mismatches within a given nucleic acid target.
- Using the described theoretical principles as a guide for experiments, reliable hybridizations have been obtained with probes six to eight nucleotides in length. All experiments were performed with a floating plastic sheet providing a film of hybridization solution above the filter. This procedure allows maximal reduction in the amount of probe, and thus reduced label costs in dot blot hybridizations. The high concentration of sodium lauroyl sarcosine instead of sodium lauroyl sulfate in the
phosphate 10 hybridization buffer allows dropping the reaction from room temperature down to 12° C. Similarly, the 4-6×SSC, 10% sodium lauroyl sarcosine buffer allows hybridization at temperatures as low as 2° C. The detergent in these buffers is essential for obtaining tolerable background with up to 40 nM concentrations of labeled probe. Preliminary characterization of the thermal stability of short oligonucleotide hybrids was determined on a prototype octamer with 50% GC content, i.e. probe of sequence TGCTCATG. The theoretical expectation is that this probe is among the less stable octamers. Its transition enthalpy is similar to those of more stable heptamers or, even toprobes 6 nucleotides in length [Bresslauer et al., Proc. Nati. Acad. Sci. U.S.A. 83:.3746 (1986)]. The stability of the 8 bp oligonucleotide duplex hybrid as a function of temperature is shown in FIG. 3. Parameter Td, the temperature at which 50% of the hybrid is melted in unit time of a minute is 18° C. The result shows that Td is 15° C. lower for the 8 bp hybrid than for an 11 bp duplex [Wallace et al., Nucleic Acids Res. 6: 3543 (1979)]. - In addition to experiments with model oligonucleotides, an M13 vector was chosen as a system for a practical demonstration of short oligonucleotide hybridization. The main aim was to show useful end-mismatch discrimination with a target similar to the ones which will be used in various applications of the method of the invention. Oligonucleotide probes for the M13 model were chosen in such a way that the M13 vector itself contains the end mismatched base. Vector IF, an M13 recombinant containing a 921 bp human interferon gene insert, carries a single perfectly matched target. Thus, IF has either the identical or a higher number of mismatched targets in comparison to the M13 vector itself.
- Using low temperature conditions and dot blots, sufficient differences in hybridization signals were obtained between the dot containing the perfect and the mismatched targets and the dot containing the mismatched targets only. This was true for the 6 nucleotides long oligonucleotides and was also true for the 7 and 8 nucleotide long oligonucleotides hybridized to the large IF-M13 pair of nucleic acids (FIG. 4A).
- The hybridization signal depends on the amount of target available on the filter for reaction with the probe. A necessary control is to show that the difference in signal intensity is not a reflection of varying amounts of nucleic acid in the two dots. Hybridization with a probe that has the same number and kind of targets in both IF and M13 shows that there is an equal amount of DNA in the dots (FIG. 4A). Since the efficiency of hybrid formation increases with hybrid length, the signal for a duplex having six nucleotides was best detected with a high mass of oligonucleotide target bound to the filter (FIG. 4A). Due to their lower molecular weight, a larger number of oligonucleotide target molecules can be bound to a given surface area when compared to large molecules of nucleic acid that serves as target.
- To measure the sensitivity of detection with unpurified DNA, various amounts of phage supernatants were spotted on the filter and hybridized with a32P-labeled octamer. As little as 50 million unpurified phage containing no more than 0.5 ng of DNA gave a detectable signal (FIG. 4B) indicating that sensitivity of the short oligonucleotide hybridization method is sufficient. Reaction time is short, adding to the practicality. For instance, the time of hybridization at 15° C. with 4 nM probe is 3 hours, after which equilibrium is reached (FIG. 4C).
- As mentioned in the theoretical section above, the equilibrium yield of hybrid depends on probe concentration and/or temperature of reaction. For instance, the signal level for the same amount of target with 4 nM octamer at 13° C. is 3 times lower than with a probe concentration of 40 nM, and is decreased 4.5 times by raising the hybridization temperature to 25° C.
- The utility of the low temperature wash for achieving maximal discrimination is demonstrated in FIG. 5. To make the phenomenon visually obvious, 50 times more DNA was put in the M13 dot than in the IF dot as shown using hybridization with a vector specific probe (FIG. 5 DNA control). In this way, the signal after the hybridization step with the actual probe was.made stronger in the mismatched than in the matched case. The HP/Hi ratio was 1:4. Inversion of signal intensities after prolonged-washing at 7° C. was achieved without a massive loss of perfect hybrid, resulting in a ratio of 2:1 (FIG. 5). In contrast, it is impossible to achieve any discrimination at 25° C., since the matched target signal is already brought down to the background level with 2 minute washing; at the same time, the signal from the mismatched hybrid is still detectable. The loss of discrimination at 13° C. compared to 7° C. is not so great but is clearly visible. If one considers the 90 minute point at 7° C. and the 15 minute point at 13° C. when the mismatched hybrid signal is near the background level, which represents optimal washing.times for the respective conditions, it is obvious that the amount of matched hybrid is several times greater at 7° C. than at 13° C. To illustrate this further, the time course of the change in discrimination with washing of the same amount of starting hybrid at the two temperatures shows the higher maximal D at the lower temperature (FIG. 5b). These results confirm the trend in the change of D with temperature and the ratio of amounts of the two types of hybrids at the start of the washing step predicted by Equation five.
- In order to show the general utility of the short oligonucleotide hybridization conditions, we have looked at hybridization of 4 heptamers, 10 octamers and an additional 14 probes up to 12 nucleotides in length in our simple M13 system. These include-the nonamer GTTTTTTAA and octamer GGCAGGCG representing the two extremes of GC content. Although GC content and sequence are expected to influence the stability of short hybrids [Bresslauer et al., Proc. Natl. Acad. Sci. U.S.A. 83: 3746 (1986)], the low temperature short oligonucleotide conditions were applicable to all tested probes in achieving sufficent discrimination. Since the best discrimination value obtained with
probes 8 nucleotides in length was 20, a several fold drop due to sequence variation is easily tolerated. - The M13 system has the advantage of showing the effects of target DNA complexity on the levels of discrimination. FIG. 6 shows the discrimination levels prior to washing for several probes with the IF-M13 system. The number of mismatched targets present in the insert and M13 vector are also shown. The general trend of a decrease in measured discrimination with the number of end-mismatched targets is obvious. Thus for two octamers having either none or five mismatched targets and differing in only one GC pair the observed discriminations were 18.3 and 1.7, respectively (FIG. 6).
- In order to show the utility of this method, three
probes 8 nucleotides in length were tested on a collection of 51 plasmid DNA dots made from a library in Bluescript vector (FIG. 7). One probe was present and specific for Bluescript vector but was absent in M13, while the other two probes had targets that were inserts of known sequence. This system allowed the use of hybridization negative or positive control DNAs with each probe. FIG. 7a shows the variation in dots from sample to sample, since the differences in hybridization with the vector specific probe should reflect the differences in DNA concentration from dot to dot. This probe sequence (CTCCCTTT) also has a complementary target in the interferon insert. Since the M13 dot (A2 FIG. 7a) is negative while the interferon insert in either M13 or Bluescript (Al and A4, respectively, FIG. 7a) is positive, the hybridization is sequence specific. Similarly, probes that detect the target sequence in only one of 51 inserts (FIG. 7b), or in none of the examined inserts (FIG. 7c) are shown along with controls that confirm that hybridization would have ocurred if the appropriate targets were present in the clones. In scoring positive hybridization, comparison of signal intensities was accomplished not between dots on the same panel, but comparisons were also made for the same dot on panels b and c with the dot on panel a in order to correct for differences in DNA amounts in the dots and, in the case of controls, in complexity. - Thermal stability curves for very short oligonucleotide hybrids that are 6-8 nucleotides in length are at least 15° C. lower than for hybrids 11-12 nucleotides in length [FIG. 1 and Wallace et al., Nucleic Acids Res. 6: 3543-3557 (1979)]. However, performing the hybridization reaction at a low temperature and with a very practical 4-40 nM concentration of oligonucleotide probe allows the detection of complementary sequence in a known or unknown nucleic acid target. To determine an unknown nucleic acid sequence completely, an entire set containing 65,535 8-mer probes can be used. Sufficient amounts of nucleic acid for this purpose are present in convenient biological samples such as a few microliters of M13 culture, a plasmid prep from 10 μl of bacterial culture or a single colony of bacteria, or less than 1 μl of a standard PCR reaction.
- Short oligonucleotides 6-10 nucleotides long give excellent discrimination. The relative decrease in hybrid stability with a single end mismatch is greater than for longer probes. Our results with the octamer TGCTCATG (FIG. 5) support this conclusion. In the experiments illustrated in FIG. 5, the target-with a G/T end mismatch is discriminated. The hybridization to the target of this type of mismatch is the most stable of all other types of mismatches which occur in more internal positions of an oligonucleotide. The discrimination achieved is the same as or greater than an internal G/T mismatch in a 19 base paired duplex [Ikuta et al., Nucl. Acids Res. 15: 797 (1987)]. Exploiting these discrimination properties using the described hybridization conditions for short oligonucleotide hybridization allows a very precise determination of oligonucleotide targets in an unknown sequence.
- In contrast to the ease of detecting discrimination between perfect and imperfect hybrids, a problem that may exist with using very short oligonucleotides is the preparation of sufficient amounts of hybrids. In practice, the need to discriminate Hp and Hi is aided by increasing the amount of DNA in the dot and/or the probe concentration, or by decreasing the hybridization temperature. However, higher probe concentrations usually increases background. Moreover, there are limits to the amounts of target nucleic acid that are practical to use. This problem was solved by the high concentration of the detergent Sarcosyl in our experiments which gave an effective background with 4 nM of probe. Further improvements can be effected either in the use of competitors for unspecific binding of probe to filter, or by changing the hybridization support material as will be recognized by the skilled artisan. Moreover, for probes having Eα less than 45 Kcal/mol (e.g. for many heptamers and a majority of hexamers modified oligonucleotides give a more stable hybrid [Asseline, et al., Proc. Natl. Acad. Sci. 81: 3297 (1984)] than their unmodified counterparts.
- The hybridization conditions described in this invention for short oligonucleotide hybridization using low temperatures give better discrimination for all sequences and duplex hybrid inputs. The only price paid in achieving uniformity in hybridization conditions for different sequences is an increase in washing time from minutes to up to 24 hours depending on the sequence. Moreover, the washing time can be further reduced by decreasing the salt concentration.
- We have shown excellent discrimination of one matched hybrid over a mismatched hybrids. However, in short oligonucleotide hybridization, signals from mismatched hybrids will always exist, with the majority of the mismatch hybrids resulting from end mismatch. This sets the limit on insert sizes that can be effectively examined by a probe of a certain length.
- The influence of sequence complexity on discrimination cannot be ignored. However, the complexity effects are more significant when defining sequence information by short oligonucleotide hybridization for specific, nonrandom sequences, and can be overcome by using an appropriate probe to target length ratio. The length ratio is chosen to make unlikely, on statistical grounds, the occurence of specific sequences which have a number of end-mismatches which would be able to eliminate or falsely invert discrimination. our results suggest the use of
oligonucleotides - In this example, hybridization conditions that were used are described supra in ↑ 6. Data resulting from the hybridization of octamer and nonomer oligonucleotides shows that sequencing by hybridization provides an extremely high degree of accuracy. In this experiment, a known sequence was used to predict a series of contiguous overlapping component octamer and nonamer oligonucleotides.
- In addition to the perfectly matching oligonucleotides, mismatch oligonucleotides wherein internal or end mismatches would occur in the duplex formed by the oligonucleotide and the target were examined. In these analyses, the lowest practical temperature was used to maximize hybridization formation. Washes were accomplished at the same or lower temperatures to ensure maximal discrimination by utilizing the greater dissociation rate of mismatch versus matched oligonucleotide/target hybridization. These conditions are shown to be applicable to all sequences although the absolute hybridization yield is shown to be sequence dependent.
- The least distabilizing mismatch that can be postulated is a simple end mismatch, so that the test of sequencing by hybridization is the ability to discriminate perfectly matched oligonucleotide/target duplexes from end-mismatched oligonucleotide/target duplexes.
- The discrimination values for 102 of 105 hybridizing oligonucleotides in a dot blot format were greater than 2 allowing a highly accurate generation of the sequence. This system also allowed an analysis of the effect of sequence on hybridization formation and hybridization instability.
- One hundred base pairs of a known portion of a human β-interferon genes prepared by PCR i.e. the 100 bp target sequence, was generated with data resulting from the hybridization of 105 oligonucleotides probes of known sequence to the target nucleic acid. The oligonucleotide probes used included 72 octamer and 21 nonamer oligonucleotides whose sequence was perfectly complementary to the target. The set of 93 probes provided consecutive overlapping frames of the target sequence displaced by one or two bases as described in FIG. 8.
- To evaluate the effect of mismatches, hybridization was examined for 12-additional probes that contained at least one end mismatch when hybridized to the 100 bp test target sequence. Also tested was the hybridization of twelve probes with target end-mismatch to four other control nucleic acid sequences chosen so that the 12 oligonucleotides could form perfectly matched duplex hybrids with the four control DNAs. Thus, the hybridization of internal mismatched, end-mismatched and perfectly matched duplex pairs of oligonucleotide and target were evaluated for each oligonucleotide used in the experiment.
- The effect of absolute DNA target concentration on the hybridization with the test octamer and nonamer oligonucleotides was determined by defining target DNA concentration by detecting hybridization of a different oligonucleotide probe to a single occurence non-target site within the co-amplified plasmid DNA.
- The results of this experiment showed that all oligonucleotides containing perfect matching complementary sequence to the target or control DNA hybridized more strongly than those oligonucleotides having mismatches. To come to this conclusion we examined Hp and D values for each probe. Hp defines the amount of hybrid duplex formed between a test target and an oligonucleotide probe. By assigning values of between 0 and 10 to the hybridization obtained for the 105 probes, it was apparent that 68.5% of the 105 probes had an HP greater than 2.
- Discrimination (D) values were obtained where D was defined as the ratio of signal intensities between 1) the dot containing a perfect matched duplex formed between test oligonucleotide and target or control nucleic acid and 2) the dot containing a mismatch duplex formed between the same oligonucleotide and a different site within the target or control nucleic acid. Variations in the value of D result from either 1) perturbations in the hybridization efficiency which allows visualization of signal over background, or 2) the type of mismatch found between the test oligonucleotide and the target. The D values obtained in this experiment were between 2 and 40 for 102 of the 105 oligonucleotide probes examined. Calculations of D for the group of 102 oligonucleotides as a whole showed the average D was 10.6.
- There were 20 cases where oligonucleotide/target duplexes exhibited an end-mismatch. In five of these, V was greater than 10. The large D value in these cases is most likely due to hybridization destabilization caused by other than the most stable (G/T and G/A) end mismatches. The other possibility is there was an error in the sequence of either the oligonucleotides or the target.
- Error in the target for probes with low Hp was excluded as a possibility because such an error would have affected the hybridization of each of the other eight overlapping oligonucleotides. There was no apparent instability due to sequence mismatch for the other overlapping oligonucleotides, indicating the target sequence was correct.
- Error in the oligonucleotide sequence was excluded as a possibility after the hybridization of seven newly synthesized oligonucleotides was re-examined. Only 1 of the seven oligonucleotides resulted in a better D value.
- Low hybrid formation values of lip could result from hybrid instability or from an inability to form hybrid duplex. An inability to form hybrid duplexes would result from either 1) self complementarity of the chosen probe or 2) target/target self hybridization. Oligonucleotide/oligonucleotide duplex formation could be favored over oligonucleotide/target hybrid duplex formation if the probe was self-complementary. Similarly, target/target association could be favored if the target was self-complementary or could form internal palindromes. In evaluating these possibilities, it was apparent from probe analysis that the questionable probes were unable to form hybrids with themselves. Moreover, in examining the contribution of target/target hybridization, it was determined that one of the questionable oligonucleotide probes hybridized inefficiently with two different DNAs containing the same target. The low probability that two different DNAs would have a self-complementary region for the same target sequence lead to the conclusion that target/target hybridization did not contribute to low hybridization formation. Thus, these results indicate that hybrid instability and not the inability to form hybrids was the cause of the low hybrid formation observed for specific oligonucleotides. The results also indicate that low hybrid formation is due to the specific sequences of certain oligonucleotides. Moreover, the results indicate that reliable results can be obtained to generate sequences if octamer and nonamer oligonucleotides are used.
- These results show that using the methods described long sequences of any specific target nucleic acid can be generated by maximal and unique overlap of constituent oligonucleotides. Such sequencing methods are dependent on the content of the individual component oligomers regardless of their frequency and their position.
- The sequence which is generated by computer using the algorithm described below is of high fidelity. The algorithm can tolerate false positive signals from the hybridization dots as is indicated from the fact the sequence generated from the 105 hybridization values, which included four less reliable values, was correct. This fidelity in sequencing by hybridization is due to the “all or none” kinetics of short oligonucleotide hybridization and the difference in duplex stability that exists between perfectly matched duplexes and mismatched duplexes. The ratio of duplex stability of matched and end-mismatched duplexes increases with decreasing duplex length. Moreover binding energy decreases with decreasing duplex length resulting in a lower hybridization efficiency. However, the results provided show that octamer hybridization allows the balancing of the factors affecting duplex stability and discrimination to produce a highly accurate method of sequencing by hybridization. Results presented in other sections show that oligonucleotides that are 6, 7 or 8 nucleotides can be effectively used to generate reliable sequence on targets that are 0.5 kb (for hexamers) 2 kB (for septamers) and 6 kB long (for octamers). The sequence of long fragments can be overlapped to generate a complete genome sequence.
- The computer algorithm to determine sequence by hybridization is described in Appendix A.
- This appendix describes an algorithm for generation of a long sequence written in a four letter alphabet from constituent k-tuple words in a minimal number of separate, randomly defined fragments of a starting nucleic acid sequence where K is the length of an oligonucleotide probe. The algorithm is primarily intended for use in the sequencing by hybridization (SBH) process. The algorithm is based on subfragments (SF), informative fragments (IF) and the possibility of using pools of physical nucleic sequences for defining informative fragments.
- As described in Appendix A, subfragments can be caused by branch points in the assembly process resulting from the repetition of a K-1 oligomer sequence in a target nucleic acid. Subfragments are sequence fragments found between any two repetitive words of the length K-1 that occur in a sequence. Multiple occurrences of K-1 words are the cause of interruption of ordering the overlap of K-words in the process of sequence generation. Interruption leads to a sequence remaining in the form of subfragments. Thus, the unambiguous segments between branching points whose order can not be uniquely determined are called sequence subfragments.
- In Appendix A informative fragments are defined as fragments of a sequence that are determined by the nearest ends of overlapped physical sequence fragments.
- A certain number of physical fragments can be pooled without losing the possibility of defining informative fragments. The total length of randomly pooled fragments depends on the length of k-tuples that are used in the sequencing process.
- The algorithm consists of two main units. The first part is used for generation of subfragments from the set of k-tuples contained in a sequence. Subfragments can be generated within the coding region of physical nucleic acid sequence of certain sizes, or within the informative fragments defined within long nucleic acid sequences. Both types of fragments are members of the basic library. This algorithm does not describe the determination of the content of the k-tuples of the informative fragments of the basic library, i.e. the step of preparation of informative fragments to be used in the sequence generation process.
- The second part of the algorithm determines the linear order of obtained subfragments with the purpose of regenerating the complete sequence of the nucleic acid fragments of the basic library. For this purpose a second, ordering library is used, made of randomly pooled fragments of the starting sequence. The algorithm does not include the step of combining sequences of basic fragments to regenerate an entire, megabase plus sequence. This can be accomplished using the link-up of fragments of the basic library which is a prerequisite for informative fragment generation. Alternatively it can be accomplished after generation of sequences of fragments of the basic library by this algorithm, using search for their overlap, based. on the presence of common end-sequences.
- The algorithm requires neither knowledge of the number of appearances of a given K-tuple in a nucleic acid sequence of the basic and ordering libraries, nor does it require the information of which k-tuple words are present on the ends of a fragment. The algorithm operates with the mixed content of k-tuples of various length. The concept of the algorithm enables operations with the k-tuple sets that contain false positive and false negative k-tuples. only in specific cases does the content of the false k-tuples primarily influence the completeness and correctness of the generated sequence. The algorithm can be used for optimization of parameters in simulation experiments, as well as for sequence generation in the actual SB2I experiments e.g. generation of the genomic DNA sequence. In optimization of parameters, the choice of the oligonucleotide probes (k-tuples) for practical and convenient fragments and/or the choice of the optimal lengths and the number of fragments for the defined probes are especially important.
- The algorithm is implemented on an IBM PC/AT compatible computer running under MSDOS 3.30 System works on the clock frequency-of-6 or 10 MHz. The workstation is equipped with one 360 kb-1.2 Mb floppy disk drive, and a 40 Mb hard disk, formatted as two 6Mb and 32 Mb partitions.
- The software is written using Microsoft's QuickBasic Compiler, version 2.0. It does not require any graphics facilities.
- 1. Generation of Subfragments
- This part of the algorithm has a central role in the process of the generation of the sequence from the content of k-tuples. It is based on the unique ordering of k-tuples by means of maximal overlap. The main obstacles in sequence generation are specific repeated sequences and false positive and/or negative k-tuples. The aim of this part of the algorithm is to obtain the minimal number of the longest possible subfragments, with correct sequence. This part of the algorithm consists of one basic, and several control steps. A two-stage process is necessary since certain information can be used only after generation of all primary subfragments.
- 1.1. Generation of Primary subfragments or pSFs
- The main problem of sequence generation is obtaining a repeated sequence from word contents that by definition do not carry information on the number of occurrences of the particular k-tuples. The concept of the entire algorithm depends on the basis on which this problem is solved. In principle, there are two opposite approaches: 1. repeated sequences can be obtained at the beginning, in the process of generation of pSFs, or 2. repeated sequences can be obtained later, in the process of the final ordering of the subfragments. In the first case, pSFs contain an excess of sequences and in the second case, they contain a deficit of sequences. The first approach requires elimination of the excess sequences generated, and the second requires permitting multiple use of some of the subfragments in the process of the final assembling of the sequence.
- The difference in the two approaches is in the degree of strictness of the rule of unique overlap of k-tuples. The less severe rule is: k-tuple X is unambiguously maximally overlapped with k-tuple Y if and only if, the rightmost R-1 end of k-tuple X is present only on the leftmost end of k-tuple Y. This rule allows the generation of repetitive sequences and the formation of surplus sequences.
- A stricter rule which is used in the second approach has an additional caveat: k-tuple X is unambiguously-maximally overlapped with k-tuple Y if and only if, the rightmost K-1 end of k-tuple X is present only on the leftmost end of k-tuple Y and if the leftmost k-1 end of k-tuple Y is not present on the rightmost end of any other k-tuple. The algorithm based on the stricter rule is simpler, and is described herein.
- FIG. 10 show the minimal conditions for unique overlap, or conversely stopping the process of elongation of a given subfragment. The process is stopped when the right k-1 end of the last k-tuple included is not present on the left end of any k-tuple (FIG. 10.1) or is present on two or more k-tuples (FIG. 10.2). If it is present on only one k-tuple the second part of the rule is tested (FIG. 10.3 and 10.4). If in addition there is a k-tuple which differs from the previously included one, the assembly of the given subfragment is terminated only on the first leftmost position (FIG. 10.3). If this additional k-tuple does not Sexist, the conditions are met for unique k-1 overlap and a given subfragment is extended to the right by one element.
- Beside the basic rule, a supplementary one is used to allow the usage of k-tuples of different lengths. The maximal overlap is the length of k-1 of the shorter k-tuple of the overlapping pair. Generation of the pSFs is performed starting from the first k-tuple from the file in which k-tuples are displayed randomly and independently from their order in a nucleic acid sequence. Thus, the first k-tuple in the file is not necessarily on the beginning of the sequence, nor on the start of the particular subfragment. The process of subfragment generation is performed by ordering the k-tuples by means of unique overlap, which is defined by the described rule. Each used k-tuple is erased from the file. At the point when there are no further k-tuples unambiguously overlapping with the last one included, the building of subfragment is terminated and the buildup of another pSF is started. Since generation of a majority of subfragments does not begin from their actual starts, the formed pSF are added to the k-tuple file and are considered as a longer k-tuple. Another possibility is to form subfragments going in both directions from the starting k-tuple. The process ends when further overlap, i.e. the extension of any of the subfragments, is not possible.
- The pSFs can be divided in three groups: 1) Subfragments of the maximal length and correct sequence in cases of exact k-tuple set; 2) short subfragments, formed due to the use of the maximal and unambiguous overlap rule on the incomplete set, and/or the set with some false Positive k-tuples; and 3) pSFs of an incorrect sequence. The incompleteness of the set in 2) is caused by false negative results of a hybridization experiment, as well as by using an incorrect set of k-tuples. These are formed due to the false positive and false negative k-tuples and can be: a) misconnected subfragments; b) subfragments with the wrong end; and c) false positive k-tuples which appears as false minimal subfragments.
- Considering false positive k-tuples, there is the possibility for the presence of a k-tuple containing more than one wrong base or containing one wrong base somewhere in the middle, as well as the possibility for a k-tuple with a wrong base on the-end. Generation of short, erroneous or misconnected subfragments is caused by the latter k-tuples. The k-tuples of the former two kinds represent wrong pSFs with length equal to k-tuple length.
- Several more common examples of obtaining short subfragments with or without a wrong base at the end because of errors in the k-tuple set are shown in FIG. 11. In the case of one false negative k-tuple (as shown in. FIG. 11A), pSFs are generated because of the impossibility of maximal overlapping (FIG. 10.1). In the case of the presence of one false positive k-tuple with the wrong base on its leftmost or rightmost end (either end in FIG. 11.B), pSFs are generated because of the impossibility of unambiguous overlapping (FIGS. 10, 11 and12). When both false positive and false negative k-tuples with a common K-1 sequence are present in the file, pSFs are generated, and one of these pSFs contains the wrong k-tuple at the relevant end (Fig. 11C).
- The process of correcting subfragments with errors in sequence and the linking of unambiguously connected pSF is performed after subfragment generation and in the process of subfragment ordering. The first step which consists of cutting the misconnected pSFs and obtaining the final subfragments by unambiguous connection of pSFs is described below.
- 1.2. Recognition and Treatment of Misconnected Subfragments
- We have defined two possibilities for the formation of misconnected subfragments. All variations of these two possibilities are shown in FIG. 12. In the first (demonstrated as column I in FIG. 12) a mistake occurs when an erroneous k-tuple appears on the points of assembly of the repeated sequences of lengths k-1. In the second, the repeated sequences are shorter than k-1 (FIG. 12, column II). These situations can occur in two variants each. In the first variant, one of the repeated sequences represents either end of a fragment (FIG. 12-A and B for column I and FIG. 12-E for column II). In the second variant, the repeated sequence occurs at any position within the fragment (FIG. 12-C and D for column I and FIG. 12-F for column II). For the first possibility, the absence of some k-tuples from the file (false negatives) is required to generate a misconnection. The second possibility requires the presence of both false negative and false positive k-tuples in The file. Considering the repetitions of K-1 sequence, the lack of only one k-tuple is sufficient when either end is repeated internally as shown in FIG. 12-A and B. The lack of two is needed for strictly internal repetition (FIG. 12-C and D). The reason is that the end of a sequence can be considered informatically as an endless linear array of false negative k-tuples. From the “smaller than K-1 case”, only the repeated sequence of the length of k-2, which requires two (FIG. 12-E) or three (FIG. 12-F) specific erroneous k-tuples, will be considered. It is very likely that these will be the only cases which will be detected in a real experiment, the others being much less frequent.
- Recognition of the misconnected subfragments is more strictly defined when a repeated sequence does not appear at the end of the fragment. In this situation, one can detect further two subfragments, one of which contains on its leftmost (FIG. 12-D), and the other on its rightmost end (FIG. 12-C and D and FIG. 12-F) k-2 sequences which are also present in the misconnected subfragment. When the repeated sequence is on the end of the fragment, there is only one subfragment which contains k-2 sequence causing the mistake in subfragment formation on its leftmost or rightmost end (as shown in FIG. 12-A and B and FIG. 12-E).
- The removal of-misconnected subfragments by their cutting is performed according to the common rule: If the leftmost or rightmost sequence of the length of k-2 of any subfragments is present in any other subfragment, the subfragment is to be cut into two subfragments, each of them containing k-2 sequence. This rule does not cover rarer situations of a repeated end when there are more than one false negative k-tuple on the point of repeated k-1 sequence. Misconnected subfragments of this kind can be recognized by using the information from the overlapped fragments, or informative fragments of both the basic and ordering libraries. In addition, the misconnected subfragment will remain when two or more false negative k-tuples occur on both positions which contain the identical k-1 sequence. This is a very rare situation since it requires at least 4 specific false k-tuples. An additional rule can be introduced to cut these subfragments on sequences of length k if the given sequence can be obtained by combination of sequences shorter than k-2 from the end of one subfragment and the start of another.
- By strict application of the described rule, some completeness is lost to ensure the accuracy of the output. Some of the subfragments will be cut although they are not misconnected, since they fit into the pattern of a misconnected subfragment. There are several situations of this kind. For example, a fragment, beside at least two identical k-1 sequences, contains any k-2 sequence from k-1 or a fragment contains k-2 sequence repeated at least twice and at least one false negative k-tuple containing given k-2 sequence in the middle, etc.
- 1.3. Merging of the pSFs into longer subfragments
- The aim of this part of the algorithm is to reduce the number of pSFs to a minimal number of longer subfragments with correct sequence. The generation of unique longer subfragments or a complete sequence is possible in two situations. The first situation concerns the specific order of repeated k-1 words. There are cases in which some or all maximally extended pSFs (the first group of pSFs) can be uniquely ordered.- For example in fragment S-R1-a-R2-b-R1-c-R2-E where S and E are the start and end of a fragment, a, b and c are different sequences specific to respective subfragments and R1 and R2 are two k-1 sequences that are tandemly repeated, five subfragments are generated (S-R1, R1-a-R2, R2-b-R1, R1-c-R2 and R2-E). They can be ordered in two ways; the original sequence above or S-R1-c-R2-b-R1-a-R2-E. In contrast, in a fragment with the same number and types of repeated sequences but ordered differently, i.e. S-R1-a-R1-b-R2-c-R2-E, there is no other sequence which includes all subfragments. Examples of this type can be recognized only after the process of generation of pSFs. They represent the necessity for two steps in the process of pSF generation. The second situation of generation of false short subfragments on positions of nonrepeated k-1 sequences when the files contain false negative and/or positive k-tuples is more important (second and third group of pSF, see FIG. 4.
- The solution for both pSF groups consists of two parts. First, the false positive k-tuples appearing as the nonexisting minimal subfragments are eliminated. All k-tuple subfragments of length k which do not have an overlap on either end, of the length of longer than k-a on one end and longer than k-b on the other end, are eliminated to enable formation of the maximal number of connections. In our experiments, the values for a and b of 2 and 3, respectively, appeared to be adequate to eliminate a sufficient number of false positive k-tuples.
- The merging of subfragments that can be uniquely connected is accomplished in the second step. The rule for connection is: two subfragments can be unambiguously connected if, and only if, the overlapping sequence at the relevant end or start of two subfragments is not present at the start and/or end of any other subfragment.
- The exception is if one subfragment from the considered pair has the identical beginning and end. In that case connection is permitted, even if there is another subfragment with the same end present in the file. The main problem here is the precise definition of overlapping sequence. The connection is not permitted if the overlapping sequence unique for only one pair of subfragments is shorter than k-2, of if it is k-2 or longer but an additional subfragment exists with the overlapping sequence of any length longer than k-4. Also, both the conanical ends of pSFs and the ends after omitting one (or few) last bases are considered as the overlapping sequences.
- After this step some false positive k-tuples (as minimal subfragments) and some subfragments with a wrong end can survive. In addition, in very rare occasions where a certain number of some specific false k-tuples are simultaneously present, an erroneous connection can take place. These cases will be detected and solved in the subfragment ordering process, and in the additional control steps along with the handling of uncut “misconnected” subfragments.
- The short subfragments that are obtained are of two kinds. In the common case, these subfragmentss can be unambiguously connected among themselves because of the distribution of repeated k-1 sequences. This can be done only after the process of generation of pSFs and is a good example of the necessity for two steps in the process of pSF generation. In the case of using the file containing false positive and/or false negative k-tuples, short pSFs are obtained on the sites of nonrepeated K-1 sequences. Considering false positive k-tuples, there is possibility of presence of the k-tuple containing more than one wrong base (or containing one wrong base somewhere in the middle), as well as k-tuple with wrong base on the end. Generation of short and erroneous (or misconnected) subfragments is caused by the latter k-tuples. The k-tuples of the former kind represent wrong pSFs with length equal to k-tuple length.
- The aim of merging pSF part of the algorithm is the reduction of the number of pSFs to the minimal number of longer subfragments with the correct sequence. All k-tuple subfragments that do not have an overlap on either end, of the length of longer than k-a on one, and longer than k-b on the other end, are eliminated to enable the maximal number of connections. In this way, the majority of false positive k-tuples are discarded. The rule for connection is: two subfragments can be unambiguously connected if, and only if the overlapping sequence of the relevant end or start of two subfragments is not present on the start and/or end of any other subfragment. The exception is a subfragment with the identical beginning and end. In that case connection is permitted, provided that there is another subfragment with the same end present in the file. The main problem here is of precise definition of overlapping sequence. The presence of at least two specific false negative k-tuples on the points of repetition of k-1 or k-2 sequences, as well as combining of the false positive and false negative k-tuples can destroy or “mask” some overlapping sequences and can produce an unambiguous, but wrong connection of pSFs. To prevent this, completeness must be sacrificed on account of exactness: the connection is not permitted on the end-sequences shorter than k-2, and in the presence of an extra overlapping sequence longer than k-4. The overlapping sequences are defined from the end of the pSFs, or omitting one, or few last bases.
- In the very rare situations, with the presence of a certain number of some specific false positive and false negative k-tuples, some subfragments with the wrong end can survive, some false positive k-tuples (as minimal subfragments) can remain, or the erroneous connection can take place. These cases are detected and solved in the subfragments ordering process, and in the additional control steps along with the handling of uncut, misconnected subfragments.
- 2. Generation of the Sequence by the Ordering of Subfragments
- The process of ordering of subfragments is similar to the process of their generation. If one considers subfragments as longer k-tuples, ordering is performed by their unambiguous connection via overlapping ends. The informational basis for unambiguous connection is the division of subfragments generated in fragments of the basic library into groups representing segments of those fragments. The method is analogous to the biochemical solution of this problem based on hybridization with longer oligonucleotides with relevant connecting sequence. The connecting sequences are generated as subfragments using the k-tuple sets of the appropriate segments of basic library fragments. Relevant segments are defined by the fragments of the ordering library that overlap with the respective fragments of the basic library. The shortest segments are informative fragments of the ordering library. The longer ones are several neighboring informative fragments or total overlapping portions of fragments corresponding of the ordering and basic libraries. In order to decrease the number of separate samples, fragments of the ordering library are randomly pooled, and the unique k-tuple content is determined.
- Although we do not use longer k-tuples, it is possible to generate the connecting sequences for two reasons. First, by using the large number of fragments in the ordering library very short segments are generated, thus reducing the chance of the multiple appearance of the k-1 sequences which are the reasons for generation of the subfragments. Second, longer segments, consisting of the various regions of the given fragment of the basic library, do not contain some of the repeated k-1 sequences. In every segment a connecting sequence (a connecting subfragment) is generated for a certain pair of the subfragments from the given fragment. The process of ordering consists of three steps: (1) generation of the k-tuple contents of each segment; (2) generation of subfragments in each segment; and (3) connection of the subfragments of the fragment of the basic library using subfragments of the segments. Primary segments are defined as significant intersections and differences of k-tuple contents of a given fragment of the basic library with the k-tuple contents of the pools of the ordering library. Secondary (shorter) segments are defined as intersections and differences of the k-tuple contents of the primary segments.
- There is a problem of accumulating both false positive and negative k-tuples in both the differences and intersections. The false negative k-tuples from starting sequences accumulate in the intersections (overlapping parts), as well as false positive k-tuples occurring randomly in both sequences, but not in the relevant overlapping region. On the other hand, the majority of false positives from either of the starting sequences is not taken up into intersections. This is an example of the reduction of experimental errors from individual fragments by using information from fragments overlapping with them. The false k-tuples accumulate in the differences for another reason. The set of false negatives from the original sequences are enlarged for false positives from intersections and the set of false positives for those k-tuples which are not included in the intersection by error, i.e. are false negatives in the intersection. If the starting sequences contains 10% false negative data, the primary and secondary intersections will contain 19% and 28% false negative k-tuples, respectively. On the other hand, a mathematical expectation of 77 false positives can be predicted if the basic fragment and the pools have lengths of 500 bp and 10,000 bp, respectively. However, there is a possibility of recovering most of the “lost” k-tuples and of eliminating most of the false positive k-tuples.
- First, one has to determine a basic content of She k-tuples for a given segment as the intersection of a given pair of the k-tuple contents. This is followed by including all k-tuples of the starting k-tuple contents in the intersection, which contain at one end k-1 and at the other end k-+ sequences which occur at the ends of two k-tuples of the basic set. This is done before generation of the differences thus preventing the accumulation of false positives in that process. Following that, the same type of enlargement of k-tuple set is applied to differences with the distinction that the borrowing is from the intersections. All borrowed k-tuples are eliminated from the intersection files as false positives.
- The intersection—a set of common k-tuples—is defined for each pair (a basic fragment)×(a pool of ordering library). If the number of k-tuples in the set is significant it is enlarged with the false negatives according to the described rule. The primary difference set is obtained by subtracting from a given basic fragment the obtained intersection set. The false negative k-tuples are appended to the difference set by borrowing from the intersection set according to the described rule and, at the same time, removed from the intersection set as false positive k-tuples. When the basic fragment is longer than the pooled fragments this difference can represent the two separate segments which somewhat reduces its utility in further steps. The primary segments are all generated intersections and differences of pairs (a basic fragment)×(a pool of ordering library) containing the-significant number of k-tuples. k-tuple sets of secondary segments are obtained by comparison of k-tuple sets of all possible pairs of primary segments. The two differences are defined for each pair which produces the intersection with the significant number of k-tuples. The majority of available information from overlapped fragments is recovered in this step so that there is little to be gained from the third round of forming intersections and differences.
- (2) Generation of the subfragments of the segments is performed identically as described for the fragments of the basic library.
- (3) The method of connection of subfragments consists of sequentially determining the correctly linked pairs of subfragments among the subfragments from a given basic library fragment which have some overlapped ends. In the case of 4 relevant subfragments, two of which contain the same beginning and two having the same end, there are 4 different pairs of subfragments that can be connected. In general 2 are correct and 2 are wrong. To find correct ones, the presence of the connecting sequence of each pair is tested in the subfragments generated from all primary and secondary segments for a given basic fragment. The length and the position of the connecting sequence are chosen to avoid interference with sequences which occur by chance. They are k+2 or longer, and include at least one element 2beside overlapping sequence in both subfragments of a given pair. The connection is permitted only if the two connecting sequences are found and the remaining two do not exist. The two linked subfragments replace former subfragments in the file and the process is cyclically repeated.
- As mentioned in the section “Generation of pSFs”, repeated sequences have to be generated in this step. This means that some subfragments have to be included in linked subfragments more than once. They will be recognized by finding the relevant connecting sequence which engages one subfragment in connection with two different subfragments.
- The recognition of misconnected subfragments generated in the processes of building pSFs and merging pSFs into longer subfragments is based on testing whether the sequences of subfragments from a given basic fragment exist in the sequences of subfragments generated in the segments for the fragment. The sequences from an incorrectly connected position will not be found indicating the misconnected subfragments.
- Beside the described three steps in ordering of subfragments some additional control steps or steps applicable to specific sequences will be necessary for the generation of more complete sequence without mistakes.
- The determination of which subfragment belongs to which segment is performed by comparison of contents of k-tuples in segments and subfragments. Because of the errors in the k-tuple contents (due to the primary error in pools and statistical errors due to the frequency of occurrences of k-tuples) the exact partitioning of subfragments is impossible. Thus, instead of “all or none” partition, the chance of coming from the given segment (P(sf,s)) is determined for each subfragment. This probability is the function of the lengths of k-tuples, the lengths of subfragments, the lengths of fragments of ordering library, the size of the pool, and of the percentage of false k-tuples in the file:
- P(sf,s)=(Ck−F)/Lsf,
- Where Lsf is the length of subfragment, Ck is the number of common k-tuples for a given subfragment/segment pair, and F is the parameter that includes relations between lengths of k-tuples, fragments of basic library, the size of the pool, and the error percentage.
- Subfragments attributed to a particular segment are treated as redundant short pSFs and are submitted to a process of unambiguous connection. The definition of unambiguous connection is slightly different in this case, since it is based on a probability that subfragments with overlapping end(s) belong to the segment considered. Besides, the accuracy of unambiguous connection is controlled by following the connection of these subfragments in other segments. After the connection in different segments, all of the obtained subfragments are merged together, shorter subfragments included within longer ones are eliminated, and the remaining ones are submitted to the ordinary connecting process. If the sequence is not regenerated completely, the process of partition and connection of subfragments is repeated with the same or less severe criterions of probability of belonging to the particular segment, followed by unambiguous connection.
- Using severe criteria for defining unambiguous overlap, some information is not used. Instead of a complete sequence, several subfragments that define a number of possibilities for a given fragment are obtained. Using less severe criteria an accurate and complete sequence is generated. In a certain number of situations, e.g. an erroneous connection, it is possible to generate a complete, but an incorrect sequence, or to generate “monster” subfragments with no possible connection among them. Thus, for each fragment of the basic library one obtains: a) several possible solutions where one is correct and b) the most probable correct solution. Also, in a very small number of cases, due to the mistake in the subfragment generation process or due to the specific ratio of the probabilities of belonging, it is impossible to generate any unambiguous solution or one, the most probable solution. These cases will remain as incomplete sequences, or the unambiguous solution will be obtained by comparing these data with other, overlapped fragments of basic library.
- 3. Implementation
- The described algorithm is tested on a randomly generated, 50 kb sequence, containing 40% GC to simulate the GC content of the human genome. In the middle part of this sequence were inserted various Alu, and some other repetitive sequences, of a total length of about 4 kb. To simulate an in vitro SBH experiment, the following operations were performed to prepare appropriate data.
- Positions of 605kb overlapping “clones” were randomly defined, to simulate preparation of basic library:
- Positions of 1000 500 bp “clones” were randomly determined, to simulate making the ordering library. These fragments were extracted from the sequence. Random pools of fragments were made, and k-tuple sets of pools were determined and stored on the hard disk. These data are used in the subfragment ordering phase: For the same density of
clones 4 million clones in basic library and 3 million clones in ordering library will be needed for the entire human genome. The total number of 7 million clones is several fold smaller than the number of clones a few kb long that are needed for random cloning of almost all of genomic DNA and its sequencing by a gel-based method. - From the data on the starts and ends of 5kb fragments, we determined 117 “informative fragments” in the Sequence. This was followed by determination of sets of overlapping k-tuples of which the single “informative fragment” consist. Only the subset of k-tuples matching our predetermined list were used. The list contains 65% 8-mers, 30% 9-mers and 5% 10-12-mers. We performed processes of generation and the ordering of subfragments on these data.
- The testing of the algorithm was performed on the simulated data in two experiments. The sequence of 50 informative fragments was regenerated with the 100% correct data set (over 20,000 bp), and 26 informative fragments about 10,000 bp) with 10% false k-tuples (5% positive and 5% negative ones).
- In the first experiment all subfragments were correct and in only one out of 50 informative fragments the sequence was not completely regenerated but remained in the form of 5 subfragments. The analysis of positions of overlapped fragments of ordering library has shown that they lack the information for the unique ordering of the 5 subfragments. The subfragments can be connected in two possible ways based on overlapping ends, 1-2-3-4-5 and 1-4-3-2-5. The only difference is the exchange of positions of
subfragments subfragments subfragment 3 region. - To simulate real sequencing, we included some false (“hybridization”) data as input in a number of experiments. In oligomer hybridization experiments, under proposed conditions the only situation possibly producing unreliable data is the end mismatch versus full match hybridization. Therefore, in simulation we considered as false positives only those k-tuples differing in a single element on either end from the real one. These “false” sets are made as follows. On the original set of a k-tuples of the informative fragment, a subset of 5% false positive k-tuples are added. False positive k-tuples are made by randomly picking a k-tuple from the set, copying it and altering a nucleotide on its beginning or end. This is followed by subtraction of a subset of 5% randomly chosen k-tuples. In this way the statistically expected number of the most complicated cases is-generated in which the correct k-tuple is replaced with a k-tuple with the wrong base on the end.
- Production of k-tuple sets as described leads to up to 10% of false data. This value varies from case to case, due to the randomness of choice of k-tuples to be copied-altered, and erased. Nevertheless, this percentage 3-4 times exceeds the amount of unreliable data in real hybridization experiments. The introduced error of 10% leads to the two fold increase in the number of subfragments both in fragments of basic library (basic library informative fragments) and in segments. About 10% of the final subfragments have a wrong base at the end as expected for the k-tuple set which contains false positives (see generation of primary subfragments). Neither the cases of misconnection of subfragments nor subfragments with wrong sequence were observed. In 4 informative fragments out of 26 examined in the ordering process the complete sequence was not regenerated. In all 4 cases the sequence was obtained in the form of several longer subfragments and several shorter subfragments contained in the same segment. This result shows that the algorithmic principles allow working with a large percentage of false data.
- 4. Discussion
- The success of the generation of the sequence from its k-tuple content can be described in terms of completeness and accuracy. In the process of generation, two particular situations can be defined: 1) Some part of the information is missing in the generated sequence, but one knows where the ambiguities are and to which type they belong, and 2) the regenerated sequence that is obtained does not match the sequence from which the k-tuple content is generated, but the mistake can not be detected. Assuming the algorithm is developed to its theoretical limits, as in the use of the exact k-tuple sets, only the first situation can take place. There the incompleteness results in a certain number of subfragments that can not be ordered unambiguously and the impossibility of determination of the exact length of monotonous sequences, i.e. number of perfect tandem repeats.
- With false k-tuples, there is also the possibility of generating incorrect sequence. The reason for mistakes does not lie in the shortcomings of the algorithm, but in the fact that a given content of k-tuples unambiguously represents the sequence that differs from the original one. One can define three classes of error, depending on the kind of the false k-tuples present in the file. False negative k-tuples (which are not accompanied with the false positives) are producing “deletions”. False positive k-tuples are producing “elongations (unequal crossing over)”. False positives accompanied with false negatives are the reason for generation of “insertions”, alone or combined with “deletions”. The deletions are produced when all of the k-tuples (or their majority) between two possible starts of the subfragments are false negatives. Since every position in the sequence is defined by k k-tuples, the occurrence of the deletions in a common case requires k consecutive false negatives. (With 10% of the false negatives and k=8, this situation takes place after every 108 elements). This situation will be extremely infreguent even in mammalian genome sequencing using random libraries containing ten genome equivalents.
- Elongation of the end of the sequence caused by false positive k-tuples is the special case of “insertions”, since the end of the sequence can be considered as the endless linear array of false negative k-tuples. One can consider a group of false positive k-tuples producing subfragments longer than one k-tuple. Situations of this kind can be detected if subfragments are generated in overlapped fragments, like random physical fragments of the ordering library. An insertion,-or insertion in place of a deletion, can arise as a result of specific combinations of false positive and false negative k-tuples. In the first case the number of consecutive false negatives has to be smaller than k. Both cases require several overlapping false positive k-tuples. The insertions and deletions are mostly theoretical possibilities without sizable practical repercussions since the requirements in the number and specificity of false k-tuples are simply too high.
- In every other situation of not meeting the theoretical requirement of the minimal number and the-kind of the false positives and/or negatives, mistakes in the k-tuples content can produce only the lesser completeness of a generated sequence.
Claims (6)
1. A method of analyzing a target nucleic acid sequence, comprising:
obtaining a plurality of H values, each H value corresponding to an extent of hybridization of an oligonucleotide probe with at least one nucleic acid sequence including said target sequence and each oligonucleotide probe differing from each other by at least a single base; and
comparing said H values to each other, thereby analyzing the target nucleic sequence.
2. A method of analyzing a target nucleic acid sequence, comprising:
obtaining a plurality of H values, each H value corresponding to an extent of hybridization of an oligonucleotide probe with said target nucleic acid sequence and each oligonucleotide probe differing from each other by at least a single base; and
calculating a D value from the ratio of a higher H value to a lower H value, thereby analyzing the target nucleic sequence.
3. A method of analyzing a target nucleic acid sequence, comprising:
obtaining a first set of H values, each H value in said first set corresponding to an extent of hybridization of an oligonucleotide probe with a control nucleic acid sequence and each probe differing from each other by at least a single base;
obtaining a second set of H values, each H value in said second set corresponding to an extent of hybridization of an oligonucleotide probe with said target nucleic acid sequence and each probe differing from each other by at least a single base; and
comparing at least one of said H values in said first set with at least one of said H values in said second set, thereby analyzing the target nucleic sequence.
4. A method of analyzing a target nucleic acid sequence, comprising:
obtaining a statistical data set corresponding to statistics about a plurality of experiments producing H values, each H value in said statistical data set corresponding to an extent of hybridization of an oligonucleotide probe with a control nucleic acid sequence and each oligonucleotide probe differing from each other by at least a single base;
obtaining a plurality of H values, each H value in said plurality corresponding to an extent of hybridization of an oligonucleotide probe with said target sequence and each probe differing from each other by at least a single base; and
comparing at least one of said H values in said plurality with said statistical data set, thereby analyzing the target nucleic sequence.
5. The method of any one of claims 1, 2, 3, or 4, wherein said method is performed on a computer.
6. The method according to claim 5 , wherein the oligonucleotide probes are in an array.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/954,426 US20020106673A1 (en) | 1987-04-01 | 2001-09-10 | Method of sequencing by hybridization of oligonucleotide probes |
US10/755,921 US7230093B2 (en) | 1987-04-01 | 2004-01-12 | Method of sequencing by hybridization of oligonucleotide probes |
Applications Claiming Priority (10)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
YU00570/87A YU57087A (en) | 1987-04-01 | 1987-04-01 | Process for obtaining genome by hebridization and oligonucleotidic tests |
YUP-570/87 | 1987-04-01 | ||
US17508888A | 1988-03-30 | 1988-03-30 | |
US57655990A | 1990-08-31 | 1990-08-31 | |
US4815293A | 1993-04-15 | 1993-04-15 | |
US08/203,502 US5525464A (en) | 1987-04-01 | 1994-02-28 | Method of sequencing by hybridization of oligonucleotide probes |
US08/460,853 US5695940A (en) | 1987-04-01 | 1995-06-05 | Method of sequencing by hybridization of oligonucleotide probes |
US82053497A | 1997-03-19 | 1997-03-19 | |
US09/205,865 US6316191B1 (en) | 1987-04-01 | 1998-12-04 | Computer-aided analysis system for sequencing by hybridization |
US09/954,426 US20020106673A1 (en) | 1987-04-01 | 2001-09-10 | Method of sequencing by hybridization of oligonucleotide probes |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US09/205,865 Continuation US6316191B1 (en) | 1987-04-01 | 1998-12-04 | Computer-aided analysis system for sequencing by hybridization |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/755,921 Division US7230093B2 (en) | 1987-04-01 | 2004-01-12 | Method of sequencing by hybridization of oligonucleotide probes |
Publications (1)
Publication Number | Publication Date |
---|---|
US20020106673A1 true US20020106673A1 (en) | 2002-08-08 |
Family
ID=46248968
Family Applications (5)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/203,502 Expired - Lifetime US5525464A (en) | 1987-04-01 | 1994-02-28 | Method of sequencing by hybridization of oligonucleotide probes |
US08/460,853 Expired - Fee Related US5695940A (en) | 1987-04-01 | 1995-06-05 | Method of sequencing by hybridization of oligonucleotide probes |
US09/176,819 Expired - Fee Related US5972619A (en) | 1987-04-01 | 1998-10-22 | Computer-aided analysis system for sequencing by hybridization |
US09/205,865 Expired - Fee Related US6316191B1 (en) | 1987-04-01 | 1998-12-04 | Computer-aided analysis system for sequencing by hybridization |
US09/954,426 Abandoned US20020106673A1 (en) | 1987-04-01 | 2001-09-10 | Method of sequencing by hybridization of oligonucleotide probes |
Family Applications Before (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/203,502 Expired - Lifetime US5525464A (en) | 1987-04-01 | 1994-02-28 | Method of sequencing by hybridization of oligonucleotide probes |
US08/460,853 Expired - Fee Related US5695940A (en) | 1987-04-01 | 1995-06-05 | Method of sequencing by hybridization of oligonucleotide probes |
US09/176,819 Expired - Fee Related US5972619A (en) | 1987-04-01 | 1998-10-22 | Computer-aided analysis system for sequencing by hybridization |
US09/205,865 Expired - Fee Related US6316191B1 (en) | 1987-04-01 | 1998-12-04 | Computer-aided analysis system for sequencing by hybridization |
Country Status (1)
Country | Link |
---|---|
US (5) | US5525464A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050266447A1 (en) * | 2004-04-19 | 2005-12-01 | Pioneer Hi-Bred International, Inc. | Method for identifying activators of gene transcription |
US20060008839A1 (en) * | 2001-08-29 | 2006-01-12 | Pioneer Hi-Bred International, Inc. | UTR tag assay for gene function discovery |
US7645596B2 (en) | 1998-05-01 | 2010-01-12 | Arizona Board Of Regents | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US7666593B2 (en) | 2005-08-26 | 2010-02-23 | Helicos Biosciences Corporation | Single molecule sequencing of captured nucleic acids |
US7897345B2 (en) | 2003-11-12 | 2011-03-01 | Helicos Biosciences Corporation | Short cycle methods for sequencing polynucleotides |
US7981604B2 (en) | 2004-02-19 | 2011-07-19 | California Institute Of Technology | Methods and kits for analyzing polynucleotide sequences |
US9096898B2 (en) | 1998-05-01 | 2015-08-04 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
Families Citing this family (556)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6270961B1 (en) * | 1987-04-01 | 2001-08-07 | Hyseq, Inc. | Methods and apparatus for DNA sequencing and DNA identification |
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 |
US5547839A (en) * | 1989-06-07 | 1996-08-20 | Affymax Technologies N.V. | Sequencing of surface immobilized polymers utilizing microflourescence detection |
US5424186A (en) | 1989-06-07 | 1995-06-13 | Affymax Technologies N.V. | Very large scale immobilized polymer synthesis |
US6406844B1 (en) | 1989-06-07 | 2002-06-18 | Affymetrix, Inc. | Very large scale immobilized polymer synthesis |
US6309822B1 (en) | 1989-06-07 | 2001-10-30 | Affymetrix, Inc. | Method for comparing copy number of nucleic acid sequences |
US5925525A (en) * | 1989-06-07 | 1999-07-20 | Affymetrix, Inc. | Method of identifying nucleotide differences |
US5800992A (en) | 1989-06-07 | 1998-09-01 | Fodor; Stephen P.A. | Method of detecting nucleic acids |
US6346413B1 (en) | 1989-06-07 | 2002-02-12 | Affymetrix, Inc. | Polymer arrays |
US5871928A (en) * | 1989-06-07 | 1999-02-16 | Fodor; Stephen P. A. | Methods for nucleic acid analysis |
US6040138A (en) | 1995-09-15 | 2000-03-21 | Affymetrix, Inc. | Expression monitoring by hybridization to high density oligonucleotide arrays |
US6551784B2 (en) | 1989-06-07 | 2003-04-22 | Affymetrix Inc | Method of comparing nucleic acid sequences |
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 |
US6379895B1 (en) | 1989-06-07 | 2002-04-30 | Affymetrix, Inc. | Photolithographic and other means for manufacturing arrays |
US5744101A (en) * | 1989-06-07 | 1998-04-28 | Affymax Technologies N.V. | Photolabile nucleoside protecting groups |
US6506558B1 (en) | 1990-03-07 | 2003-01-14 | Affymetrix Inc. | Very large scale immobilized polymer synthesis |
DE69132843T2 (en) | 1990-12-06 | 2002-09-12 | Affymetrix Inc N D Ges D Staat | Identification of nucleic acids in samples |
US6589726B1 (en) | 1991-09-04 | 2003-07-08 | Metrigen, Inc. | Method and apparatus for in situ synthesis on a solid support |
US6943034B1 (en) * | 1991-11-22 | 2005-09-13 | Affymetrix, Inc. | Combinatorial strategies for polymer synthesis |
US6468740B1 (en) | 1992-11-05 | 2002-10-22 | Affymetrix, Inc. | Cyclic and substituted immobilized molecular synthesis |
US5677195A (en) * | 1991-11-22 | 1997-10-14 | Affymax Technologies N.V. | Combinatorial strategies for polymer synthesis |
JPH0622798A (en) * | 1992-07-07 | 1994-02-01 | Hitachi Ltd | Method for determining base sequence |
AU7212494A (en) * | 1993-06-25 | 1995-01-17 | Affymax Technologies N.V. | Hybridization and sequencing of nucleic acids |
US6401267B1 (en) * | 1993-09-27 | 2002-06-11 | Radoje Drmanac | Methods and compositions for efficient nucleic acid sequencing |
US6331274B1 (en) | 1993-11-01 | 2001-12-18 | Nanogen, Inc. | Advanced active circuits and devices for molecular biological analysis and diagnostics |
US6225059B1 (en) | 1993-11-01 | 2001-05-01 | Nanogen, Inc. | Advanced active electronic devices including collection electrodes for molecular biological analysis and diagnostics |
US7625697B2 (en) | 1994-06-17 | 2009-12-01 | The Board Of Trustees Of The Leland Stanford Junior University | Methods for constructing subarrays and subarrays made thereby |
US5710000A (en) | 1994-09-16 | 1998-01-20 | Affymetrix, Inc. | Capturing sequences adjacent to Type-IIs restriction sites for genomic library mapping |
USRE43097E1 (en) | 1994-10-13 | 2012-01-10 | Illumina, Inc. | Massively parallel signature sequencing by ligation of encoded adaptors |
US8236493B2 (en) | 1994-10-21 | 2012-08-07 | Affymetrix, Inc. | Methods of enzymatic discrimination enhancement and surface-bound double-stranded DNA |
US6600996B2 (en) | 1994-10-21 | 2003-07-29 | Affymetrix, Inc. | Computer-aided techniques for analyzing biological sequences |
US5795716A (en) | 1994-10-21 | 1998-08-18 | Chee; Mark S. | Computer-aided visualization and analysis system for sequence evaluation |
US5763169A (en) * | 1995-01-13 | 1998-06-09 | Chiron Diagnostics Corporation | Nucleic acid probes for the detection and identification of fungi |
US6180339B1 (en) | 1995-01-13 | 2001-01-30 | Bayer Corporation | Nucleic acid probes for the detection and identification of fungi |
US5959098A (en) * | 1996-04-17 | 1999-09-28 | Affymetrix, Inc. | Substrate preparation process |
US5801155A (en) * | 1995-04-03 | 1998-09-01 | Epoch Pharmaceuticals, Inc. | Covalently linked oligonucleotide minor grove binder conjugates |
US6312894B1 (en) | 1995-04-03 | 2001-11-06 | Epoch Pharmaceuticals, Inc. | Hybridization and mismatch discrimination using oligonucleotides conjugated to minor groove binders |
US6017496A (en) | 1995-06-07 | 2000-01-25 | Irori | Matrices with memories and uses thereof |
US5925562A (en) * | 1995-04-25 | 1999-07-20 | Irori | Remotely programmable matrices with memories |
US6329139B1 (en) | 1995-04-25 | 2001-12-11 | Discovery Partners International | Automated sorting system for matrices with memory |
US5830655A (en) | 1995-05-22 | 1998-11-03 | Sri International | Oligonucleotide sizing using cleavable primers |
US20040175718A1 (en) * | 1995-10-16 | 2004-09-09 | Affymetrix, Inc. | Computer-aided visualization and analysis system for sequence evaluation |
US5736332A (en) * | 1995-11-30 | 1998-04-07 | Mandecki; Wlodek | Method of determining the sequence of nucleic acids employing solid-phase particles carrying transponders |
WO1997019958A1 (en) | 1995-11-30 | 1997-06-05 | Wlodek Mandecki | Screening of drugs from chemical combinatorial libraries employing transponders |
US5641634A (en) | 1995-11-30 | 1997-06-24 | Mandecki; Wlodek | Electronically-indexed solid-phase assay for biomolecules |
US6001571A (en) | 1995-11-30 | 1999-12-14 | Mandecki; Wlodek | Multiplex assay for nucleic acids employing transponders |
US6051377A (en) * | 1995-11-30 | 2000-04-18 | Pharmaseq, Inc. | Multiplex assay for nucleic acids employing transponders |
EP0880598A4 (en) | 1996-01-23 | 2005-02-23 | Affymetrix Inc | Nucleic acid analysis techniques |
US20020150921A1 (en) * | 1996-02-09 | 2002-10-17 | Francis Barany | Detection of nucleic acid sequence differences using the ligase detection reaction with addressable arrays |
US6852487B1 (en) | 1996-02-09 | 2005-02-08 | Cornell Research Foundation, Inc. | Detection of nucleic acid sequence differences using the ligase detection reaction with addressable arrays |
US6361937B1 (en) | 1996-03-19 | 2002-03-26 | Affymetrix, Incorporated | Computer-aided nucleic acid sequencing |
US6706875B1 (en) * | 1996-04-17 | 2004-03-16 | Affyemtrix, Inc. | Substrate preparation process |
EP0912761A4 (en) * | 1996-05-29 | 2004-06-09 | Cornell Res Foundation Inc | Detection of nucleic acid sequence differences using coupled ligase detection and polymerase chain reactions |
US5770370A (en) * | 1996-06-14 | 1998-06-23 | David Sarnoff Research Center, Inc. | Nuclease protection assays |
US6391550B1 (en) * | 1996-09-19 | 2002-05-21 | Affymetrix, Inc. | Identification of molecular sequence signatures and methods involving the same |
US6329140B1 (en) | 1996-09-19 | 2001-12-11 | Affymetrix, Inc. | Identification of molecular sequence signatures and methods involving the same |
US7094609B2 (en) | 1996-09-20 | 2006-08-22 | Burstein Technologies, Inc. | Spatially addressable combinatorial chemical arrays in encoded optical disk format |
US6096273A (en) * | 1996-11-05 | 2000-08-01 | Clinical Micro Sensors | Electrodes linked via conductive oligomers to nucleic acids |
US7381525B1 (en) * | 1997-03-07 | 2008-06-03 | Clinical Micro Sensors, Inc. | AC/DC voltage apparatus for detection of nucleic acids |
US6291164B1 (en) * | 1996-11-22 | 2001-09-18 | Invitrogen Corporation | Methods for preventing inhibition of nucleic acid synthesis by pyrophosphate |
US5905024A (en) * | 1996-12-17 | 1999-05-18 | University Of Chicago | Method for performing site-specific affinity fractionation for use in DNA sequencing |
US20020042048A1 (en) | 1997-01-16 | 2002-04-11 | Radoje Drmanac | Methods and compositions for detection or quantification of nucleic acid species |
US7144699B2 (en) | 1997-03-20 | 2006-12-05 | Affymetrix, Inc. | Iterative resequencing |
US6391622B1 (en) * | 1997-04-04 | 2002-05-21 | Caliper Technologies Corp. | Closed-loop biochemical analyzers |
US6235471B1 (en) * | 1997-04-04 | 2001-05-22 | Caliper Technologies Corp. | Closed-loop biochemical analyzers |
US6143496A (en) | 1997-04-17 | 2000-11-07 | Cytonix Corporation | Method of sampling, amplifying and quantifying segment of nucleic acid, polymerase chain reaction assembly having nanoliter-sized sample chambers, and method of filling assembly |
US5981166A (en) * | 1997-04-23 | 1999-11-09 | Pharmaseq, Inc. | Screening of soluble chemical compounds for their pharmacological properties utilizing transponders |
AU736321B2 (en) | 1997-05-23 | 2001-07-26 | Lynx Therapeutics, Inc. | System and apparatus for sequential processing of analytes |
CA2297681A1 (en) | 1997-07-22 | 1999-02-04 | Kristen Moynihan | Apparatus and methods for arraying solution onto a solid support |
EP1000176A1 (en) | 1997-07-22 | 2000-05-17 | Rapigene, Inc. | Computer method and system for correlating sequencing data by ms |
WO1999005321A1 (en) | 1997-07-22 | 1999-02-04 | Rapigene, Inc. | Amplification and other enzymatic reactions performed on nucleic acid arrays |
US6826296B2 (en) | 1997-07-25 | 2004-11-30 | Affymetrix, Inc. | Method and system for providing a probe array chip design database |
DE69823206T2 (en) * | 1997-07-25 | 2004-08-19 | Affymetrix, Inc. (a Delaware Corp.), Santa Clara | METHOD FOR PRODUCING A BIO-INFORMATICS DATABASE |
US6420108B2 (en) | 1998-02-09 | 2002-07-16 | Affymetrix, Inc. | Computer-aided display for comparative gene expression |
US7068830B2 (en) * | 1997-07-25 | 2006-06-27 | Affymetrix, Inc. | Method and system for providing a probe array chip design database |
AU9476298A (en) * | 1997-09-11 | 1999-03-29 | Seq, Ltd. | Method for dna sequencing analysis |
US7214298B2 (en) * | 1997-09-23 | 2007-05-08 | California Institute Of Technology | Microfabricated cell sorter |
US6607878B2 (en) | 1997-10-06 | 2003-08-19 | Stratagene | Collections of uniquely tagged molecules |
AU1118599A (en) * | 1997-10-24 | 1999-05-17 | Trustees Of Boston University | Inference sequencing by hybridization |
US6101946A (en) * | 1997-11-21 | 2000-08-15 | Telechem International Inc. | Microarray printing device including printing pins with flat tips and exterior channel and method of manufacture |
US6087102A (en) | 1998-01-07 | 2000-07-11 | Clontech Laboratories, Inc. | Polymeric arrays and methods for their use in binding assays |
US6893877B2 (en) | 1998-01-12 | 2005-05-17 | Massachusetts Institute Of Technology | Methods for screening substances in a microwell array |
US6428752B1 (en) | 1998-05-14 | 2002-08-06 | Affymetrix, Inc. | Cleaning deposit devices that form microarrays and the like |
US6407858B1 (en) | 1998-05-14 | 2002-06-18 | Genetic Microsystems, Inc | Focusing of microscopes and reading of microarrays |
US6722395B2 (en) | 1998-01-13 | 2004-04-20 | James W. Overbeck | Depositing fluid specimens on substrates, resulting ordered arrays, techniques for analysis of deposited arrays |
US6269846B1 (en) | 1998-01-13 | 2001-08-07 | Genetic Microsystems, Inc. | Depositing fluid specimens on substrates, resulting ordered arrays, techniques for deposition of arrays |
US20040035690A1 (en) * | 1998-02-11 | 2004-02-26 | The Regents Of The University Of Michigan | Method and apparatus for chemical and biochemical reactions using photo-generated reagents |
CA2319587C (en) * | 1998-02-11 | 2004-09-21 | University Of Houston | Method and apparatus for chemical and biochemical reactions using photo-generated reagents |
CA2321070C (en) * | 1998-02-23 | 2010-04-06 | Wisconsin Alumni Research Foundation | Method and apparatus for synthesis of arrays of dna probes |
US7095032B2 (en) * | 1998-03-20 | 2006-08-22 | Montagu Jean I | Focusing of microscopes and reading of microarrays |
US6683173B2 (en) | 1998-04-03 | 2004-01-27 | Epoch Biosciences, Inc. | Tm leveling methods |
US6426191B1 (en) | 1998-04-03 | 2002-07-30 | Hyseq, Inc. | Assays involving an IL-1 receptor antagonist |
US7715989B2 (en) * | 1998-04-03 | 2010-05-11 | Elitech Holding B.V. | Systems and methods for predicting oligonucleotide melting temperature (TmS) |
US6337072B1 (en) | 1998-04-03 | 2002-01-08 | Hyseq, Inc. | Interleukin-1 receptor antagonist and recombinant production thereof |
US7045610B2 (en) * | 1998-04-03 | 2006-05-16 | Epoch Biosciences, Inc. | Modified oligonucleotides for mismatch discrimination |
US6294655B1 (en) | 1998-04-03 | 2001-09-25 | Hyseq, Inc. | Anti-interleukin-1 receptor antagonist antibodies and uses thereof |
US6541623B1 (en) | 1998-04-03 | 2003-04-01 | Hyseq, Inc. | Interleukin—1 receptor antagonist and uses thereof |
US6127121A (en) * | 1998-04-03 | 2000-10-03 | Epoch Pharmaceuticals, Inc. | Oligonucleotides containing pyrazolo[3,4-D]pyrimidines for hybridization and mismatch discrimination |
US6949367B1 (en) | 1998-04-03 | 2005-09-27 | Epoch Pharmaceuticals, Inc. | Modified oligonucleotides for mismatch discrimination |
US6376619B1 (en) | 1998-04-13 | 2002-04-23 | 3M Innovative Properties Company | High density, miniaturized arrays and methods of manufacturing same |
US6355419B1 (en) | 1998-04-27 | 2002-03-12 | Hyseq, Inc. | Preparation of pools of nucleic acids based on representation in a sample |
US6268210B1 (en) | 1998-05-27 | 2001-07-31 | Hyseq, Inc. | Sandwich arrays of biological compounds |
US6271957B1 (en) * | 1998-05-29 | 2001-08-07 | Affymetrix, Inc. | Methods involving direct write optical lithography |
US6699659B2 (en) | 1998-06-01 | 2004-03-02 | Affymetrix, Inc. | Products and methods for analyzing nucleic acids including identification of substitutions, insertions and deletions |
US20050244954A1 (en) * | 1998-06-23 | 2005-11-03 | Blackburn Gary F | Binding acceleration techniques for the detection of analytes |
EP1090144A1 (en) * | 1998-06-24 | 2001-04-11 | Smithkline Beecham Corporation | Method for detecting, analyzing, and mapping rna transcripts |
US6593133B1 (en) | 1998-07-06 | 2003-07-15 | Nsgene A/S | Neurotrophic factors |
US6218118B1 (en) | 1998-07-09 | 2001-04-17 | Agilent Technologies, Inc. | Method and mixture reagents for analyzing the nucleotide sequence of nucleic acids by mass spectrometry |
US6476211B1 (en) | 1998-07-16 | 2002-11-05 | Hyseq, Inc. | Methods and materials relating to CD39-like polypeptides |
US6387645B1 (en) | 1998-07-16 | 2002-05-14 | Hyseq, Inc. | Methods and materials relating to novel CD39-like polypeptides |
US6447771B1 (en) * | 1999-03-19 | 2002-09-10 | Hyseq, Inc. | Methods and materials relating to novel CD39-like polypeptides |
AU762392B2 (en) * | 1998-07-16 | 2003-06-26 | Arca Biopharma, Inc. | Methods and materials relating to novel CD39-like polypeptides |
JP2002524091A (en) | 1998-08-21 | 2002-08-06 | ナックスコー・インコーポレイテッド | Assays using immobilized crosslinkable nucleic acids |
US6306643B1 (en) * | 1998-08-24 | 2001-10-23 | Affymetrix, Inc. | Methods of using an array of pooled probes in genetic analysis |
US6387234B1 (en) | 1998-08-31 | 2002-05-14 | Iowa State University Research Foundation, Inc. | Integrated multiplexed capillary electrophoresis system |
US6642062B2 (en) | 1998-09-03 | 2003-11-04 | Trellis Bioinformatics, Inc. | Multihued labels |
US6185561B1 (en) | 1998-09-17 | 2001-02-06 | Affymetrix, Inc. | Method and apparatus for providing and expression data mining database |
US20040199544A1 (en) * | 2000-11-02 | 2004-10-07 | Affymetrix, Inc. | Method and apparatus for providing an expression data mining database |
US6284466B1 (en) | 1998-09-18 | 2001-09-04 | The Board Of Regents Of University Of Nebraska | Method of detecting genetic polymorphisms using over represented sequences |
US20090119022A1 (en) * | 1998-09-25 | 2009-05-07 | Timberlake William E | Emericella Nidulans Genome Sequence On Computer Readable Medium and Uses Thereof |
WO2000018955A1 (en) * | 1998-09-28 | 2000-04-06 | MAX-PLANCK-Gesellschaft zur Förderung der Wissenschaften e.V. | Novel method for the preselection of shotgun clones of a genome or a portion thereof of an organism |
US7034143B1 (en) * | 1998-10-13 | 2006-04-25 | Brown University Research Foundation | Systems and methods for sequencing by hybridization |
US6489096B1 (en) | 1998-10-15 | 2002-12-03 | Princeton University | Quantitative analysis of hybridization patterns and intensities in oligonucleotide arrays |
US6221592B1 (en) | 1998-10-20 | 2001-04-24 | Wisconsin Alumi Research Foundation | Computer-based methods and systems for sequencing of individual nucleic acid molecules |
US6607888B2 (en) | 1998-10-20 | 2003-08-19 | Wisconsin Alumni Research Foundation | Method for analyzing nucleic acid reactions |
US6545264B1 (en) | 1998-10-30 | 2003-04-08 | Affymetrix, Inc. | Systems and methods for high performance scanning |
US6399299B1 (en) | 1998-11-02 | 2002-06-04 | Perkinelmer Life Sciences, Inc. | Amplified array analysis system |
US20110014706A2 (en) * | 1998-12-14 | 2011-01-20 | Monsanto Technology Llc | Arabidopsis thaliana Genome Sequence and Uses Thereof |
US6245518B1 (en) * | 1998-12-11 | 2001-06-12 | Hyseq, Inc. | Polynucleotide arrays and methods of making and using the same |
US6815163B1 (en) * | 1998-12-23 | 2004-11-09 | Rutgers, The State University | Methods and kits for screening nucleic acid duplex stability |
US7468250B2 (en) * | 1998-12-23 | 2008-12-23 | Rutgers, The State University Of New Jersey | Methods and kits for screening nucleic acid duplex stability |
US20030180789A1 (en) * | 1998-12-30 | 2003-09-25 | Dale Roderic M.K. | Arrays with modified oligonucleotide and polynucleotide compositions |
US6087112A (en) | 1998-12-30 | 2000-07-11 | Oligos Etc. Inc. | Arrays with modified oligonucleotide and polynucleotide compositions |
WO2000040758A2 (en) * | 1999-01-06 | 2000-07-13 | Hyseq Inc. | Enhanced sequencing by hybridization using pools of probes |
ATE343425T1 (en) | 1999-01-08 | 2006-11-15 | Applera Corp | FIBER MATRIX FOR BRINGING CHEMICAL SUBSTANCES TOGETHER, AND METHOD FOR THE PRODUCTION AND USE THEREOF |
US7595189B2 (en) * | 1999-01-08 | 2009-09-29 | Applied Biosystems, Llc | Integrated optics fiber array |
US20050026209A1 (en) * | 1999-01-08 | 2005-02-03 | Vann Charles S. | Optical fiber bundle for detecting binding of chemical species |
US6780977B1 (en) | 1999-01-29 | 2004-08-24 | Nuvelo, Inc. | Methods and compositions relating to CD39-like polypeptides and nucleic acids |
US6783959B1 (en) | 1999-01-29 | 2004-08-31 | Nuvelo, Inc. | Methods and compositions relating to CD39-like polypeptides and nucleic acids |
US6899875B1 (en) | 1999-01-29 | 2005-05-31 | Nuvelo, Inc. | Methods and compositions relating to CD39-like polypeptides and nucleic acids |
US6350447B1 (en) | 1999-01-29 | 2002-02-26 | Hyseq, Inc. | Methods and compositions relating to CD39-like polypeptides and nucleic acids |
US6500609B1 (en) | 1999-02-11 | 2002-12-31 | Scynexis Chemistry & Automation, Inc. | Method and apparatus for synthesizing characterizing and assaying combinatorial libraries |
US6506594B1 (en) * | 1999-03-19 | 2003-01-14 | Cornell Res Foundation Inc | Detection of nucleic acid sequence differences using the ligase detection reaction with addressable arrays |
WO2000056456A1 (en) | 1999-03-19 | 2000-09-28 | Genencor International, Inc. | Multi-through hole testing plate for high throughput screening |
US7014994B1 (en) | 1999-03-19 | 2006-03-21 | Cornell Research Foundation,Inc. | Coupled polymerase chain reaction-restriction-endonuclease digestion-ligase detection reaction process |
US6335013B1 (en) | 1999-03-19 | 2002-01-01 | Hyseq, Inc. | Methods and materials relating to CD39-like polypeptides |
US6537755B1 (en) | 1999-03-25 | 2003-03-25 | Radoje T. Drmanac | Solution-based methods and materials for sequence analysis by hybridization |
ATE413467T1 (en) | 1999-04-20 | 2008-11-15 | Illumina Inc | DETECTION OF NUCLEIC ACID REACTIONS ON BEAD ARRAYS |
US20060275782A1 (en) | 1999-04-20 | 2006-12-07 | Illumina, Inc. | Detection of nucleic acid reactions on bead arrays |
EP2360270B1 (en) | 1999-05-20 | 2016-11-09 | Illumina, Inc. | Combinatorial decoding of random nucleic acid arrays |
WO2000079007A1 (en) * | 1999-06-19 | 2000-12-28 | Hyseq Inc. | Improved methods of sequence assembly in sequencing by hybridization |
US6818395B1 (en) | 1999-06-28 | 2004-11-16 | California Institute Of Technology | Methods and apparatus for analyzing polynucleotide sequences |
ATE417127T1 (en) | 1999-07-26 | 2008-12-15 | Clinical Micro Sensors Inc | NUKELIC ACID SEQUENCE DETERMINATION USING ELECTRONIC DETECTION |
US6339147B1 (en) | 1999-07-29 | 2002-01-15 | Epoch Biosciences, Inc. | Attachment of oligonucleotides to solid supports through Schiff base type linkages for capture and detection of nucleic acids |
US6524863B1 (en) | 1999-08-04 | 2003-02-25 | Scynexis Chemistry & Automation, Inc. | High throughput HPLC method for determining Log P values |
US6413431B1 (en) | 1999-08-10 | 2002-07-02 | Scynexis Chemistry & Automation, Inc. | HPLC method for purifying organic compounds |
EP1212128A1 (en) | 1999-08-27 | 2002-06-12 | Scynexis Chemistry and Automation, Inc. | Sample preparation for high throughput purification |
US20050118618A1 (en) * | 1999-08-30 | 2005-06-02 | Dale Roderic M. | Method for detecting nucleic acid sequences |
US7244559B2 (en) * | 1999-09-16 | 2007-07-17 | 454 Life Sciences Corporation | Method of sequencing a nucleic acid |
US6274320B1 (en) | 1999-09-16 | 2001-08-14 | Curagen Corporation | Method of sequencing a nucleic acid |
US7211390B2 (en) * | 1999-09-16 | 2007-05-01 | 454 Life Sciences Corporation | Method of sequencing a nucleic acid |
US6844151B1 (en) | 1999-09-29 | 2005-01-18 | Oligos Etc. Inc. | Methods for production of arrays with modified oligonucleotide and polynucleotide compositions |
US6376177B1 (en) | 1999-10-06 | 2002-04-23 | Virtual Pro, Inc. | Apparatus and method for the analysis of nucleic acids hybridization on high density NA chips |
AU1075701A (en) * | 1999-10-08 | 2001-04-23 | Protogene Laboratories, Inc. | Method and apparatus for performing large numbers of reactions using array assembly |
US7668658B2 (en) | 1999-10-13 | 2010-02-23 | Sequenom, Inc. | Methods for generating databases and databases for identifying polymorphic genetic markers |
US8153367B2 (en) * | 1999-10-29 | 2012-04-10 | Perkinelmer Las, Inc. | Amplified array analysis system |
US6589778B1 (en) | 1999-12-15 | 2003-07-08 | Amersham Biosciences Ab | Method and apparatus for performing biological reactions on a substrate surface |
US6569674B1 (en) | 1999-12-15 | 2003-05-27 | Amersham Biosciences Ab | Method and apparatus for performing biological reactions on a substrate surface |
US6642046B1 (en) * | 1999-12-09 | 2003-11-04 | Motorola, Inc. | Method and apparatus for performing biological reactions on a substrate surface |
US7041445B2 (en) | 1999-11-15 | 2006-05-09 | Clontech Laboratories, Inc. | Long oligonucleotide arrays |
US6660845B1 (en) | 1999-11-23 | 2003-12-09 | Epoch Biosciences, Inc. | Non-aggregating, non-quenching oligomers comprising nucleotide analogues; methods of synthesis and use thereof |
US6727356B1 (en) * | 1999-12-08 | 2004-04-27 | Epoch Pharmaceuticals, Inc. | Fluorescent quenching detection reagents and methods |
US20040081959A9 (en) * | 1999-12-08 | 2004-04-29 | Epoch Biosciences, Inc. | Fluorescent quenching detection reagents and methods |
US7205105B2 (en) * | 1999-12-08 | 2007-04-17 | Epoch Biosciences, Inc. | Real-time linear detection probes: sensitive 5′-minor groove binder-containing probes for PCR analysis |
US6482638B1 (en) * | 1999-12-09 | 2002-11-19 | 3M Innovative Properties Company | Heat-relaxable substrates and arrays |
US6902891B2 (en) | 1999-12-17 | 2005-06-07 | Bio Merieux | Process for labeling a nucleic acid |
US6489114B2 (en) | 1999-12-17 | 2002-12-03 | Bio Merieux | Process for labeling a ribonucleic acid, and labeled RNA fragments which are obtained thereby |
US20030097222A1 (en) * | 2000-01-25 | 2003-05-22 | Craford David M. | Method, system, and computer software for providing a genomic web portal |
US7582420B2 (en) | 2001-07-12 | 2009-09-01 | Illumina, Inc. | Multiplex nucleic acid reactions |
US7611869B2 (en) | 2000-02-07 | 2009-11-03 | Illumina, Inc. | Multiplexed methylation detection methods |
US8076063B2 (en) | 2000-02-07 | 2011-12-13 | Illumina, Inc. | Multiplexed methylation detection methods |
US7955794B2 (en) | 2000-09-21 | 2011-06-07 | Illumina, Inc. | Multiplex nucleic acid reactions |
EP2261240B1 (en) | 2000-02-11 | 2015-09-02 | The Texas A & M University System | Biosensor compositions and methods of use |
WO2001061043A2 (en) * | 2000-02-16 | 2001-08-23 | Illumina, Inc. | Parallel genotyping of multiple patient samples |
US20020151040A1 (en) | 2000-02-18 | 2002-10-17 | Matthew O' Keefe | Apparatus and methods for parallel processing of microvolume liquid reactions |
CN1404415A (en) * | 2000-02-22 | 2003-03-19 | 基因谱公司 | Microarray fabrication techniques and apparatus |
CA2399189A1 (en) | 2000-02-22 | 2001-08-30 | Genospectra, Inc. | Microarray fabrication techniques and apparatus |
US20040014102A1 (en) * | 2000-02-22 | 2004-01-22 | Shiping Chen | High density parallel printing of microarrays |
CA2401781A1 (en) | 2000-03-01 | 2001-09-07 | Epoch Biosciences, Inc. | Modified oligonucleotides for mismatch discrimination |
EP1944310A3 (en) | 2000-03-01 | 2008-08-06 | Epoch Biosciences, Inc. | Modified oligonucleotides for mismatch discrimination |
AU2001240040A1 (en) * | 2000-03-03 | 2001-09-17 | California Institute Of Technology | Combinatorial array for nucleic acid analysis |
WO2001070974A1 (en) | 2000-03-24 | 2001-09-27 | Takeda Chemical Industries, Ltd. | Novel protein, process for producing the same and use thereof |
JP5508654B2 (en) | 2000-04-14 | 2014-06-04 | コーネル・リサーチ・ファンデーション・インコーポレイテッド | Method for designing positionable arrays for detection of nucleic acid sequence differences using ligase detection reactions |
ATE423221T1 (en) * | 2000-06-13 | 2009-03-15 | Univ Boston | USE OF MASS-MATCHED NUCLEOTIDES IN THE ANALYSIS OF OLIGONUCLEOTIDE MIXTURES AND IN HIGH-MULTIPLEX NUCLEIC ACID SEQUENCING |
WO2002000343A2 (en) * | 2000-06-27 | 2002-01-03 | Fluidigm Corporation | A microfluidic design automation method and system |
US6958214B2 (en) | 2000-07-10 | 2005-10-25 | Sequenom, Inc. | Polymorphic kinase anchor proteins and nucleic acids encoding the same |
US7062092B2 (en) | 2000-08-22 | 2006-06-13 | Affymetrix, Inc. | System, method, and computer software product for gain adjustment in biological microarray scanner |
US6865491B1 (en) | 2000-08-22 | 2005-03-08 | Ramot At Tel-Aviv University Ltd. | Method for sequencing polynucleotides |
US6713257B2 (en) | 2000-08-25 | 2004-03-30 | Rosetta Inpharmatics Llc | Gene discovery using microarrays |
AU1141602A (en) * | 2000-10-04 | 2002-04-15 | Celadon Lab Inc | Computer system for designing oligonucleotides used in biochemical methods |
EP1336097A4 (en) * | 2000-10-13 | 2006-02-01 | Fluidigm Corp | Microfluidic device based sample injection system for analytical devices |
WO2002042485A2 (en) * | 2000-11-21 | 2002-05-30 | Affymetrix, Inc. | Methods and computer software products for selecting nucleic acid probes |
US7442370B2 (en) | 2001-02-01 | 2008-10-28 | Biogen Idec Ma Inc. | Polymer conjugates of mutated neublastin |
US20050196785A1 (en) * | 2001-03-05 | 2005-09-08 | California Institute Of Technology | Combinational array for nucleic acid analysis |
US7297518B2 (en) | 2001-03-12 | 2007-11-20 | California Institute Of Technology | Methods and apparatus for analyzing polynucleotide sequences by asynchronous base extension |
US7276580B2 (en) * | 2001-03-12 | 2007-10-02 | Biogen Idec Ma Inc. | Neurotrophic factors |
US6949340B2 (en) | 2001-03-28 | 2005-09-27 | Creative Mines Llc | Optical phase modulator |
US20030143556A1 (en) * | 2001-04-03 | 2003-07-31 | Gary Blackburn | Nucleic acid reactions using labels with different redox potentials |
US6960437B2 (en) | 2001-04-06 | 2005-11-01 | California Institute Of Technology | Nucleic acid amplification utilizing microfluidic devices |
WO2002090599A1 (en) * | 2001-05-09 | 2002-11-14 | Genetic Id, Inc. | Universal microarray system |
US6855538B2 (en) * | 2001-06-27 | 2005-02-15 | The Regents Of The University Of California | High-efficiency microarray printing device |
US7402286B2 (en) * | 2001-06-27 | 2008-07-22 | The Regents Of The University Of California | Capillary pins for high-efficiency microarray printing device |
EP2246438B1 (en) | 2001-07-12 | 2019-11-27 | Illumina, Inc. | Multiplex nucleic acid reactions |
US6893822B2 (en) | 2001-07-19 | 2005-05-17 | Nanogen Recognomics Gmbh | Enzymatic modification of a nucleic acid-synthetic binding unit conjugate |
US20080026367A9 (en) * | 2001-08-17 | 2008-01-31 | Perlegen Sciences, Inc. | Methods for genomic analysis |
US6767731B2 (en) | 2001-08-27 | 2004-07-27 | Intel Corporation | Electron induced fluorescent method for nucleic acid sequencing |
US10539561B1 (en) | 2001-08-30 | 2020-01-21 | Customarray, Inc. | Enzyme-amplified redox microarray detection process |
AU2002364945A1 (en) * | 2001-10-25 | 2003-07-09 | Neurogenetics, Inc. | Genes and polymorphisms on chromosome 10 associated with alzheimer's disease and other neurodegenerative diseases |
US20030170678A1 (en) * | 2001-10-25 | 2003-09-11 | Neurogenetics, Inc. | Genetic markers for Alzheimer's disease and methods using the same |
US20030224380A1 (en) * | 2001-10-25 | 2003-12-04 | The General Hospital Corporation | Genes and polymorphisms on chromosome 10 associated with Alzheimer's disease and other neurodegenerative diseases |
US8440093B1 (en) | 2001-10-26 | 2013-05-14 | Fuidigm Corporation | Methods and devices for electronic and magnetic sensing of the contents of microfluidic flow channels |
US6821730B2 (en) | 2001-11-09 | 2004-11-23 | Intel Corporation | Carbon nanotube molecular labels |
WO2003042697A1 (en) * | 2001-11-14 | 2003-05-22 | Genospectra, Inc. | Biochemical analysis system with combinatorial chemistry applications |
US7691333B2 (en) * | 2001-11-30 | 2010-04-06 | Fluidigm Corporation | Microfluidic device and methods of using same |
US7118910B2 (en) | 2001-11-30 | 2006-10-10 | Fluidigm Corporation | Microfluidic device and methods of using same |
US20030177380A1 (en) * | 2002-03-18 | 2003-09-18 | Woods Stanley P. | Devices for storing array data and methods of using the same |
US7427480B2 (en) * | 2002-03-26 | 2008-09-23 | Perlegen Sciences, Inc. | Life sciences business systems and methods |
US7135286B2 (en) * | 2002-03-26 | 2006-11-14 | Perlegen Sciences, Inc. | Pharmaceutical and diagnostic business systems and methods |
US6916621B2 (en) | 2002-03-27 | 2005-07-12 | Spectral Genomics, Inc. | Methods for array-based comparitive binding assays |
JP2005521425A (en) | 2002-04-01 | 2005-07-21 | フルイディグム コーポレイション | Microfluidic particle analysis system |
CA2484676A1 (en) | 2002-05-03 | 2003-11-13 | Sequenom, Inc. | Kinase anchor protein muteins, peptides thereof, and related methods |
US7108976B2 (en) * | 2002-06-17 | 2006-09-19 | Affymetrix, Inc. | Complexity management of genomic DNA by locus specific amplification |
US9388459B2 (en) | 2002-06-17 | 2016-07-12 | Affymetrix, Inc. | Methods for genotyping |
US8277753B2 (en) | 2002-08-23 | 2012-10-02 | Life Technologies Corporation | Microfluidic transfer pin |
US7563600B2 (en) | 2002-09-12 | 2009-07-21 | Combimatrix Corporation | Microarray synthesis and assembly of gene-length polynucleotides |
JP4694201B2 (en) * | 2002-09-20 | 2011-06-08 | インテグレイテッド ディーエヌエイ テクノロジーズ インコーポレイテッド | Anthraquinone quenching dyes, their production and use |
US7143785B2 (en) | 2002-09-25 | 2006-12-05 | California Institute Of Technology | Microfluidic large scale integration |
EP1546412B1 (en) | 2002-10-02 | 2014-05-21 | California Institute Of Technology | Microfluidic nucleic acid analysis |
US7459273B2 (en) * | 2002-10-04 | 2008-12-02 | Affymetrix, Inc. | Methods for genotyping selected polymorphism |
US6742843B2 (en) * | 2002-10-25 | 2004-06-01 | Knoll, Inc. | Cable control with overload protection device |
JP4571776B2 (en) * | 2002-11-05 | 2010-10-27 | Jx日鉱日石エネルギー株式会社 | Lubricating oil composition |
EP2112229A3 (en) | 2002-11-25 | 2009-12-02 | Sequenom, Inc. | Methods for identifying risk of breast cancer and treatments thereof |
US8602774B2 (en) * | 2002-12-04 | 2013-12-10 | Bryan Wasylucha | Process of tooth whitening and apparatus therefor |
JP4395133B2 (en) * | 2002-12-20 | 2010-01-06 | カリパー・ライフ・サイエンシズ・インク. | Single molecule amplification and detection of DNA |
US20090246879A1 (en) * | 2002-12-20 | 2009-10-01 | Callida Genomics | Materials and Methods Relating to Nano-Tags and Nano-Brands |
WO2004074818A2 (en) | 2002-12-20 | 2004-09-02 | Biotrove, Inc. | Assay apparatus and method using microfluidic arrays |
WO2005003375A2 (en) * | 2003-01-29 | 2005-01-13 | 454 Corporation | Methods of amplifying and sequencing nucleic acids |
US7575865B2 (en) * | 2003-01-29 | 2009-08-18 | 454 Life Sciences Corporation | Methods of amplifying and sequencing nucleic acids |
CA2930400A1 (en) * | 2003-02-26 | 2004-09-10 | Callida Genomics, Inc. | Apparatus with pixel alignment system for analyzing nucleic acid |
US20070141570A1 (en) * | 2003-03-07 | 2007-06-21 | Sequenom, Inc. | Association of polymorphic kinase anchor proteins with cardiac phenotypes and related methods |
FR2852317B1 (en) | 2003-03-13 | 2006-08-04 | PROBE BIOPUCES AND METHODS OF USE | |
US7604965B2 (en) | 2003-04-03 | 2009-10-20 | Fluidigm Corporation | Thermal reaction device and method for using the same |
CA2521171C (en) * | 2003-04-03 | 2013-05-28 | Fluidigm Corp. | Microfluidic devices and methods of using same |
US20050145496A1 (en) | 2003-04-03 | 2005-07-07 | Federico Goodsaid | Thermal reaction device and method for using the same |
US8828663B2 (en) | 2005-03-18 | 2014-09-09 | Fluidigm Corporation | Thermal reaction device and method for using the same |
US7476363B2 (en) | 2003-04-03 | 2009-01-13 | Fluidigm Corporation | Microfluidic devices and methods of using same |
CA2522364C (en) | 2003-04-18 | 2014-12-09 | Biogen Idec Ma Inc. | Polymer-conjugated glycosylated neublastin |
US20050014173A1 (en) * | 2003-05-08 | 2005-01-20 | Farrer Matthew J. | Parkinson's disease markers |
EP1623020A4 (en) * | 2003-05-09 | 2007-05-09 | Sigma Aldrich Co | Genomic and proteomic approaches for the development of cell culture medium |
EP2107129B1 (en) | 2003-07-31 | 2012-08-22 | Sequenom, Inc. | Methods for high level multiplexed polymerase chain reactions and homogeneous mass extension reactions for genotyping of polymorhisms |
US8114978B2 (en) | 2003-08-05 | 2012-02-14 | Affymetrix, Inc. | Methods for genotyping selected polymorphism |
US7413712B2 (en) * | 2003-08-11 | 2008-08-19 | California Institute Of Technology | Microfluidic rotary flow reactor matrix |
US20050043245A1 (en) * | 2003-08-11 | 2005-02-24 | Siler -Khodr Theresa M. | Non-mammalian GnRH analogs and uses thereof in the immune system |
US7348146B2 (en) * | 2003-10-02 | 2008-03-25 | Epoch Biosciences, Inc. | Single nucleotide polymorphism analysis of highly polymorphic target sequences |
SI1687609T1 (en) * | 2003-10-28 | 2015-03-31 | Epoch Biosciences, Inc. | Fluorescent probes for dna detection by hybridization with improved sensitivity and low background |
US7439341B2 (en) * | 2003-11-14 | 2008-10-21 | Integrated Dna Technologies, Inc. | Fluorescence quenching azo dyes, their methods of preparation and use |
US7276338B2 (en) * | 2003-11-17 | 2007-10-02 | Jacobson Joseph M | Nucleotide sequencing via repetitive single molecule hybridization |
US20050164211A1 (en) * | 2004-01-22 | 2005-07-28 | Hannah Eric C. | Carbon nanotube molecular labels |
WO2005084134A2 (en) * | 2004-03-04 | 2005-09-15 | Dena Leshkowitz | Quantifying and profiling antibody and t cell receptor gene expression |
AU2005222618A1 (en) | 2004-03-12 | 2005-09-29 | Biotrove, Inc. | Nanoliter array loading |
US7476734B2 (en) | 2005-12-06 | 2009-01-13 | Helicos Biosciences Corporation | Nucleotide analogs |
JP2008512084A (en) | 2004-05-25 | 2008-04-24 | ヘリコス バイオサイエンシーズ コーポレイション | Methods and devices for nucleic acid sequencing |
EP2290071B1 (en) | 2004-05-28 | 2014-12-31 | Asuragen, Inc. | Methods and compositions involving microRNA |
US20060008833A1 (en) * | 2004-07-12 | 2006-01-12 | Jacobson Joseph M | Method for long, error-reduced DNA synthesis |
US20080199480A1 (en) * | 2004-07-22 | 2008-08-21 | Sequenom, Inc. | Methods for Identifying Risk of Type II Diabetes and Treatments Thereof |
EP2238983B1 (en) * | 2004-08-19 | 2014-07-23 | Biogen Idec MA Inc. | Neublastin variants |
CN101123978B (en) | 2004-08-19 | 2012-12-12 | 比奥根艾迪克Ma公司 | Neublastin variants |
US20110098554A1 (en) * | 2004-09-29 | 2011-04-28 | Tel Hashomer Medical Research Infrastructure And Services Ltd. | Monitoring of convection enhanced drug delivery |
US20070122817A1 (en) * | 2005-02-28 | 2007-05-31 | George Church | Methods for assembly of high fidelity synthetic polynucleotides |
JP2008523786A (en) * | 2004-10-18 | 2008-07-10 | コドン デバイシズ インコーポレイテッド | Method for assembling high fidelity synthetic polynucleotides |
EP2281887A1 (en) | 2004-11-12 | 2011-02-09 | Asuragen, Inc. | Methods and compositions involving miRNA and miRNA inhibitor molecules |
US20060102471A1 (en) | 2004-11-18 | 2006-05-18 | Karl Maurer | Electrode array device having an adsorbed porous reaction layer |
US7220549B2 (en) | 2004-12-30 | 2007-05-22 | Helicos Biosciences Corporation | Stabilizing a nucleic acid for nucleic acid sequencing |
US7482120B2 (en) | 2005-01-28 | 2009-01-27 | Helicos Biosciences Corporation | Methods and compositions for improving fidelity in a nucleic acid synthesis reaction |
US20060190184A1 (en) * | 2005-02-23 | 2006-08-24 | Incogen, Inc. | System and method using a visual or audio-visual programming environment to enable and optimize systems-level research in life sciences |
EP2166100B1 (en) | 2005-03-08 | 2012-07-18 | BASF Plant Science GmbH | Expression enhancing intron sequences |
WO2006095259A2 (en) | 2005-03-11 | 2006-09-14 | Novartis Ag | Biomarkers for cardiovascular side-effects induced by cox-2 inhibitory compounds |
JP2008537875A (en) | 2005-03-14 | 2008-10-02 | ザ ボード オブ トラスティーズ オブ ザ リーランド スタンフォード ジュニア ユニバーシティ | Methods and compositions for assessing graft survival in solid organ transplant recipients |
US20070034513A1 (en) | 2005-03-25 | 2007-02-15 | Combimatrix Corporation | Electrochemical deblocking solution for electrochemical oligomer synthesis on an electrode array |
US9394167B2 (en) | 2005-04-15 | 2016-07-19 | Customarray, Inc. | Neutralization and containment of redox species produced by circumferential electrodes |
US7452671B2 (en) * | 2005-04-29 | 2008-11-18 | Affymetrix, Inc. | Methods for genotyping with selective adaptor ligation |
EP1907560B1 (en) | 2005-05-20 | 2013-01-23 | Integrated DNA Technologies, Inc. | Compounds and methods for labeling oligonucleotides |
EP1731620A1 (en) | 2005-06-07 | 2006-12-13 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Method for the diagnosis of immune graft tolerance |
US20090264299A1 (en) * | 2006-02-24 | 2009-10-22 | Complete Genomics, Inc. | High throughput genome sequencing on DNA arrays |
EP3492602A1 (en) | 2005-06-15 | 2019-06-05 | Complete Genomics, Inc. | Single molecule arrays for genetic and chemical analysis |
US20070065877A1 (en) | 2005-09-19 | 2007-03-22 | Combimatrix Corporation | Microarray having a base cleavable succinate linker |
US20070065834A1 (en) * | 2005-09-19 | 2007-03-22 | Hillis William D | Method and sequences for determinate nucleic acid hybridization |
US7960104B2 (en) | 2005-10-07 | 2011-06-14 | Callida Genomics, Inc. | Self-assembled single molecule arrays and uses thereof |
WO2007047408A2 (en) * | 2005-10-12 | 2007-04-26 | Pathologica, Llc. | Promac signature application |
CN101310183A (en) | 2005-11-16 | 2008-11-19 | 诺瓦提斯公司 | Biomarkers for anti-nogo-a antibody treatment in spinal cord injury |
EP2206779B1 (en) * | 2005-12-06 | 2013-02-13 | Life Technologies Corporation | Reverse transcription primers and methods of design |
US20090305248A1 (en) * | 2005-12-15 | 2009-12-10 | Lander Eric G | Methods for increasing accuracy of nucleic acid sequencing |
US11306351B2 (en) | 2005-12-21 | 2022-04-19 | Affymetrix, Inc. | Methods for genotyping |
US20070168197A1 (en) * | 2006-01-18 | 2007-07-19 | Nokia Corporation | Audio coding |
EP2495337A1 (en) * | 2006-02-24 | 2012-09-05 | Callida Genomics, Inc. | High throughput genome sequencing on DNA arrays |
TWI501774B (en) * | 2006-02-27 | 2015-10-01 | Biogen Idec Inc | Treatments for neurological disorders |
ES2450065T3 (en) * | 2006-03-01 | 2014-03-21 | Biogen Idec Ma Inc. | Compositions and methods for the administration of GDNF ligand family proteins |
WO2007136834A2 (en) * | 2006-05-19 | 2007-11-29 | Codon Devices, Inc. | Combined extension and ligation for nucleic acid assembly |
EP3617321A3 (en) | 2006-05-31 | 2020-04-29 | Sequenom, Inc. | Methods and compositions for the extraction and amplification of nucleic acid from a sample |
AU2007257162A1 (en) | 2006-06-05 | 2007-12-13 | Cancer Care Ontario | Assessment of risk for colorectal cancer |
US7745129B1 (en) * | 2006-07-12 | 2010-06-29 | Kenneth David Schatz | Methods for sequencing of a necleic acid |
AU2007273055B2 (en) | 2006-07-14 | 2014-05-01 | The Regents Of The University Of California | Cancer biomarkers and methods of use thereof |
US8053191B2 (en) | 2006-08-31 | 2011-11-08 | Westend Asset Clearinghouse Company, Llc | Iterative nucleic acid assembly using activation of vector-encoded traits |
JP2010510964A (en) | 2006-09-19 | 2010-04-08 | アシュラジェン インコーポレイテッド | MiR-15, miR-26, miR-31, miR-145, miR-147, miR-188, miR-215, miR-216, miR-331, mmu-miR-292 as targets for therapeutic intervention Genes and pathways regulated by 3p |
AU2007299828C1 (en) | 2006-09-19 | 2014-07-17 | Interpace Diagnostics, Llc | MicroRNAs differentially expressed in pancreatic diseases and uses thereof |
GB0618514D0 (en) * | 2006-09-20 | 2006-11-01 | Univ Nottingham Trent | Method of detecting interactions on a microarray using nuclear magnetic resonance |
US20090105961A1 (en) * | 2006-11-09 | 2009-04-23 | Complete Genomics, Inc. | Methods of nucleic acid identification in large-scale sequencing |
US7902345B2 (en) | 2006-12-05 | 2011-03-08 | Sequenom, Inc. | Detection and quantification of biomolecules using mass spectrometry |
CN101657217B (en) | 2007-01-30 | 2013-11-20 | 药品循环公司 | Methods for determining cancer resistance to histone deacetylase inhibitors |
WO2008098142A2 (en) | 2007-02-08 | 2008-08-14 | Sequenom, Inc. | Nucleic acid-based tests for rhd typing, gender determination and nucleic acid quantification |
EP1961825A1 (en) | 2007-02-26 | 2008-08-27 | INSERM (Institut National de la Santé et de la Recherche Medicale) | Method for predicting the occurrence of metastasis in breast cancer patients |
CA2679954A1 (en) | 2007-03-05 | 2008-09-12 | Cancer Care Ontario | Assessment of risk for colorectal cancer |
WO2008109176A2 (en) | 2007-03-07 | 2008-09-12 | President And Fellows Of Harvard College | Assays and other reactions involving droplets |
US8652780B2 (en) | 2007-03-26 | 2014-02-18 | Sequenom, Inc. | Restriction endonuclease enhanced polymorphic sequence detection |
US7572990B2 (en) * | 2007-03-30 | 2009-08-11 | Intermec Ip Corp. | Keypad overlay membrane |
JP5583005B2 (en) * | 2007-05-01 | 2014-09-03 | バイオジェン・アイデック・エムエイ・インコーポレイテッド | Compositions and methods for increasing angiogenesis |
US20090068170A1 (en) * | 2007-07-13 | 2009-03-12 | President And Fellows Of Harvard College | Droplet-based selection |
EP2205634A2 (en) * | 2007-08-08 | 2010-07-14 | Biogen Idec MA, Inc. | Anti-neublastin antibodies and uses thereof |
WO2009021684A2 (en) | 2007-08-14 | 2009-02-19 | F. Hoffmann-La Roche Ag | Predictive marker for egfr inhibitor treatment |
AU2008286407A1 (en) | 2007-08-14 | 2009-02-19 | F. Hoffmann-La Roche Ag | Predictive markers for EGFR inhibitor treatment |
JP5276103B2 (en) | 2007-08-14 | 2013-08-28 | エフ.ホフマン−ラ ロシュ アーゲー | EGFR inhibitor therapeutic marker |
CN101784674A (en) | 2007-08-14 | 2010-07-21 | 霍夫曼-拉罗奇有限公司 | Predictive markers for EGFR inhibitors treatment |
US9404150B2 (en) | 2007-08-29 | 2016-08-02 | Sequenom, Inc. | Methods and compositions for universal size-specific PCR |
US9388457B2 (en) | 2007-09-14 | 2016-07-12 | Affymetrix, Inc. | Locus specific amplification using array probes |
EP2201136B1 (en) | 2007-10-01 | 2017-12-06 | Nabsys 2.0 LLC | Nanopore sequencing by hybridization of probes to form ternary complexes and variable range alignment |
US8951731B2 (en) * | 2007-10-15 | 2015-02-10 | Complete Genomics, Inc. | Sequence analysis using decorated nucleic acids |
US7897344B2 (en) * | 2007-11-06 | 2011-03-01 | Complete Genomics, Inc. | Methods and oligonucleotide designs for insertion of multiple adaptors into library constructs |
US8518640B2 (en) * | 2007-10-29 | 2013-08-27 | Complete Genomics, Inc. | Nucleic acid sequencing and process |
US8617811B2 (en) * | 2008-01-28 | 2013-12-31 | Complete Genomics, Inc. | Methods and compositions for efficient base calling in sequencing reactions |
US8415099B2 (en) | 2007-11-05 | 2013-04-09 | Complete Genomics, Inc. | Efficient base determination in sequencing reactions |
WO2009061840A1 (en) * | 2007-11-05 | 2009-05-14 | Complete Genomics, Inc. | Methods and oligonucleotide designs for insertion of multiple adaptors employing selective methylation |
WO2009073629A2 (en) * | 2007-11-29 | 2009-06-11 | Complete Genomics, Inc. | Efficient shotgun sequencing methods |
US8592150B2 (en) | 2007-12-05 | 2013-11-26 | Complete Genomics, Inc. | Methods and compositions for long fragment read sequencing |
WO2009085215A1 (en) | 2007-12-21 | 2009-07-09 | President And Fellows Of Harvard College | Systems and methods for nucleic acid sequencing |
WO2009094583A1 (en) * | 2008-01-23 | 2009-07-30 | Complete Genomics, Inc. | Methods and compositions for preventing bias in amplification and sequencing reactions |
US8211652B2 (en) * | 2008-01-29 | 2012-07-03 | University of Pittsburgh—of the Commonwealth System of Higher Education | FSTL-1 as a biomaker of inflammation |
US8658379B2 (en) | 2008-01-29 | 2014-02-25 | University of Pittsburgh—of the Commonwealth System of Higher Education | Follistatin-like protein-1 as a biomarker for sepsis |
ES2590177T3 (en) * | 2008-02-29 | 2016-11-18 | Monsanto Technology Llc | MON87460 corn plant event and compositions and procedures for its detection |
US9074244B2 (en) | 2008-03-11 | 2015-07-07 | Affymetrix, Inc. | Array-based translocation and rearrangement assays |
US8709726B2 (en) | 2008-03-11 | 2014-04-29 | Sequenom, Inc. | Nucleic acid-based tests for prenatal gender determination |
CA2718137A1 (en) | 2008-03-26 | 2009-10-01 | Sequenom, Inc. | Restriction endonuclease enhanced polymorphic sequence detection |
WO2009132028A1 (en) * | 2008-04-21 | 2009-10-29 | Complete Genomics, Inc. | Array structures for nucleic acid detection |
US8258111B2 (en) | 2008-05-08 | 2012-09-04 | The Johns Hopkins University | Compositions and methods related to miRNA modulation of neovascularization or angiogenesis |
WO2010096036A2 (en) | 2008-05-14 | 2010-08-26 | Millennium Pharmaceuticals, Inc. | Methods and kits for monitoring the effects of immunomodulators on adaptive immunity |
EP2324129A4 (en) | 2008-08-18 | 2012-06-20 | Univ Leland Stanford Junior | Methods and compositions for determining a graft tolerant phenotype in a subject |
ES2581774T3 (en) | 2008-08-29 | 2016-09-07 | Janssen Biotech, Inc. | Markers and methods of evaluation and treatment of ulcerative colitis and related disorders using a panel of 20 genes |
US8262879B2 (en) | 2008-09-03 | 2012-09-11 | Nabsys, Inc. | Devices and methods for determining the length of biopolymers and distances between probes bound thereto |
EP2342362B1 (en) | 2008-09-03 | 2017-03-01 | Nabsys 2.0 LLC | Use of longitudinally displaced nanoscale electrodes for voltage sensing of biomolecules and other analytes in fluidic channels |
US9650668B2 (en) | 2008-09-03 | 2017-05-16 | Nabsys 2.0 Llc | Use of longitudinally displaced nanoscale electrodes for voltage sensing of biomolecules and other analytes in fluidic channels |
US8476013B2 (en) * | 2008-09-16 | 2013-07-02 | Sequenom, Inc. | Processes and compositions for methylation-based acid enrichment of fetal nucleic acid from a maternal sample useful for non-invasive prenatal diagnoses |
US8962247B2 (en) * | 2008-09-16 | 2015-02-24 | Sequenom, Inc. | Processes and compositions for methylation-based enrichment of fetal nucleic acid from a maternal sample useful for non invasive prenatal diagnoses |
WO2010033200A2 (en) | 2008-09-19 | 2010-03-25 | President And Fellows Of Harvard College | Creation of libraries of droplets and related species |
US20120052066A1 (en) | 2008-11-07 | 2012-03-01 | Cesar Calderon | Markers and methods for assessing and treating lupus patients susceptible to photoprovocation |
EP2373812B1 (en) * | 2008-12-19 | 2016-11-09 | President and Fellows of Harvard College | Particle-assisted nucleic acid sequencing |
CN102428191A (en) * | 2009-03-18 | 2012-04-25 | 塞昆纳姆股份有限公司 | Use Of Thermostable Endonucleases For Generating Reporter Molecules |
US8455260B2 (en) | 2009-03-27 | 2013-06-04 | Massachusetts Institute Of Technology | Tagged-fragment map assembly |
EP3211095B1 (en) | 2009-04-03 | 2019-01-02 | Sequenom, Inc. | Nucleic acid preparation compositions and methods |
US8246799B2 (en) | 2009-05-28 | 2012-08-21 | Nabsys, Inc. | Devices and methods for analyzing biomolecules and probes bound thereto |
US9524369B2 (en) | 2009-06-15 | 2016-12-20 | Complete Genomics, Inc. | Processing and analysis of complex nucleic acid sequence data |
WO2011006169A1 (en) | 2009-07-10 | 2011-01-13 | Portola Pharmaceuticals, Inc. | Methods for diagnosis and treatment of thrombotic disorders mediated by cyp2c19*2 |
DK2430454T3 (en) | 2009-08-22 | 2013-02-11 | Univ Leland Stanford Junior | Imaging and assessment of embryos, oocytes and stem cells |
US9243290B2 (en) | 2009-10-09 | 2016-01-26 | Stc.Unm | Polony sequencing methods |
CN102648053B (en) | 2009-10-27 | 2016-04-27 | 哈佛学院院长等 | Drop formation technology |
WO2011056872A2 (en) | 2009-11-03 | 2011-05-12 | Gen9, Inc. | Methods and microfluidic devices for the manipulation of droplets in high fidelity polynucleotide assembly |
WO2011066185A1 (en) | 2009-11-25 | 2011-06-03 | Gen9, Inc. | Microfluidic devices and methods for gene synthesis |
EP3185013B1 (en) | 2009-12-02 | 2019-10-09 | The Board of Trustees of the Leland Stanford Junior University | Biomarkers for determining an allograft tolerant phenotype |
WO2011071935A2 (en) * | 2009-12-07 | 2011-06-16 | Brigham And Women's Hospital, Inc. | Methods for diagnosing and assessing risk of developing glomerulosclerosis |
GB0922085D0 (en) | 2009-12-17 | 2010-02-03 | Cambridge Entpr Ltd | Cancer diagnosis and treatment |
US9926593B2 (en) | 2009-12-22 | 2018-03-27 | Sequenom, Inc. | Processes and kits for identifying aneuploidy |
US9217144B2 (en) | 2010-01-07 | 2015-12-22 | Gen9, Inc. | Assembly of high fidelity polynucleotides |
WO2011090793A2 (en) | 2010-01-20 | 2011-07-28 | Customarray, Inc. | Multiplex microarray of serially deposited biomolecules on a microarray |
EP2542678B1 (en) | 2010-03-04 | 2017-04-12 | InteRNA Technologies B.V. | A MiRNA MOLECULE DEFINED BY ITS SOURCE AND ITS THERAPEUTIC USES IN CANCER ASSOCIATED WITH EMT |
AU2011230496B2 (en) | 2010-03-26 | 2015-09-17 | Integrated Dna Technologies, Inc. | Methods for enhancing nucleic acid hybridization |
US9506057B2 (en) | 2010-03-26 | 2016-11-29 | Integrated Dna Technologies, Inc. | Modifications for antisense compounds |
US20190300945A1 (en) | 2010-04-05 | 2019-10-03 | Prognosys Biosciences, Inc. | Spatially Encoded Biological Assays |
US10787701B2 (en) | 2010-04-05 | 2020-09-29 | Prognosys Biosciences, Inc. | Spatially encoded biological assays |
EP2560568B1 (en) | 2010-04-19 | 2016-09-21 | Syneron Medical Ltd. | Combined energy and topical composition application for regulating the condition of mammalian skin |
KR20120137407A (en) | 2010-04-19 | 2012-12-20 | 더 프록터 앤드 갬블 캄파니 | Genetic signatures and gene chips associated with administration of electrically conducted radio frequency current to skin and methods and treatments relating thereto |
EP3508854A1 (en) | 2010-04-27 | 2019-07-10 | The Regents of The University of California | Cancer biomarkers and methods of use thereof |
WO2011146725A1 (en) | 2010-05-19 | 2011-11-24 | Bayer Healthcare Llc | Biomarkers for a multikinase inhibitor |
WO2012005572A1 (en) | 2010-07-06 | 2012-01-12 | Interna Technologies Bv | Mirna and its diagnostic and therapeutic uses in diseases or conditions associated with melanoma, or in diseases or conditions associated with activated braf pathway |
WO2012019099A2 (en) | 2010-08-05 | 2012-02-09 | University Of Pittsburgh - Of The Commonwealth System Of Higher Education | Follistatin-like-protein-1 as a biomarker for inflammatory disorders |
CA2809457C (en) | 2010-09-07 | 2019-07-30 | Integrated Dna Technologies, Inc. | Modifications for antisense compounds |
KR20130115250A (en) | 2010-09-15 | 2013-10-21 | 알막 다이아그노스틱스 리미티드 | Molecular diagnostic test for cancer |
US8715933B2 (en) | 2010-09-27 | 2014-05-06 | Nabsys, Inc. | Assay methods using nicking endonucleases |
US20120095029A1 (en) | 2010-10-15 | 2012-04-19 | Hoffmann-La Roche Inc. | Ipp complex as marker for erlotinib treatment |
US10457935B2 (en) | 2010-11-12 | 2019-10-29 | Gen9, Inc. | Protein arrays and methods of using and making the same |
US9295965B2 (en) | 2010-11-12 | 2016-03-29 | Gen9, Inc. | Methods and devices for nucleic acid synthesis |
US8859201B2 (en) | 2010-11-16 | 2014-10-14 | Nabsys, Inc. | Methods for sequencing a biomolecule by detecting relative positions of hybridized probes |
BR112013012265A2 (en) | 2010-11-17 | 2016-08-02 | Asuragen Inc | mirnas as biomarkers to distinguish benign from malignant thyroid neoplasms |
EP2474617A1 (en) | 2011-01-11 | 2012-07-11 | InteRNA Technologies BV | Mir for treating neo-angiogenesis |
WO2012109574A2 (en) | 2011-02-11 | 2012-08-16 | Nabsys, Inc. | Assay methods using dna binding proteins |
EP2678675B1 (en) | 2011-02-23 | 2017-10-11 | The Board of Trustees of the Leland Stanford Junior University | Methods of detecting aneuploidy in human embryos |
WO2012158238A2 (en) | 2011-02-28 | 2012-11-22 | University Of Iowa Research Foundation | Anti-müllerian hormone changes in pregnancy and prediction ofadverse pregnancy outcomes and gender |
ES2611343T3 (en) | 2011-03-23 | 2017-05-08 | Elitech Holding B.V. | Functionalized 3-alkynyl pyrazolopyrimidine analogs as universal bases and methods of use |
US8969003B2 (en) | 2011-03-23 | 2015-03-03 | Elitech Holding B.V. | Functionalized 3-alkynyl pyrazolopyrimidine analogues as universal bases and methods of use |
GB201106254D0 (en) | 2011-04-13 | 2011-05-25 | Frisen Jonas | Method and product |
CA2834218C (en) | 2011-04-29 | 2021-02-16 | Sequenom, Inc. | Quantification of a minority nucleic acid species using inhibitory oligonucleotides |
EP2801626B1 (en) | 2011-05-24 | 2016-04-06 | ELITechGroup B.V. | Detection of methicillin-resistant staphylococcus aureus |
EA025926B1 (en) | 2011-06-02 | 2017-02-28 | Олмак Дайэгностикс Лимитед | Molecular diagnostic test for cancer |
AU2012300401B2 (en) | 2011-08-26 | 2018-02-08 | Gen9, Inc. | Compositions and methods for high fidelity assembly of nucleic acids |
US9644241B2 (en) | 2011-09-13 | 2017-05-09 | Interpace Diagnostics, Llc | Methods and compositions involving miR-135B for distinguishing pancreatic cancer from benign pancreatic disease |
WO2013055995A2 (en) | 2011-10-14 | 2013-04-18 | President And Fellows Of Harvard College | Sequencing by structure assembly |
WO2013063519A1 (en) | 2011-10-26 | 2013-05-02 | Asuragen, Inc. | Methods and compositions involving mirna expression levels for distinguishing pancreatic cysts |
WO2013063544A1 (en) | 2011-10-27 | 2013-05-02 | Asuragen, Inc. | Mirnas as diagnostic biomarkers to distinguish benign from malignant thyroid tumors |
CN106978487A (en) | 2011-10-28 | 2017-07-25 | 米伦纽姆医药公司 | Biomarkers of response to inhibitors of NEDD8 activating enzyme (NAE) |
US20150184246A1 (en) | 2011-11-11 | 2015-07-02 | Millennium Pharmaceuticals, Inc. | Biomarkers of response to proteasome inhibitors |
EP2776043B1 (en) | 2011-11-11 | 2018-02-21 | Millennium Pharmaceuticals, Inc. | Biomarkers of response to proteasome inhibitors |
WO2013087789A1 (en) | 2011-12-13 | 2013-06-20 | Glykos Finland Ltd. | Antibody isoform arrays and methods thereof |
EP4108782B1 (en) | 2011-12-22 | 2023-06-07 | President and Fellows of Harvard College | Compositions and methods for analyte detection |
US11021737B2 (en) | 2011-12-22 | 2021-06-01 | President And Fellows Of Harvard College | Compositions and methods for analyte detection |
EP3401399B1 (en) | 2012-03-02 | 2020-04-22 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
US9150853B2 (en) | 2012-03-21 | 2015-10-06 | Gen9, Inc. | Methods for screening proteins using DNA encoded chemical libraries as templates for enzyme catalysis |
CN111876409A (en) | 2012-04-24 | 2020-11-03 | Gen9股份有限公司 | Method for sorting nucleic acids and multiplex preparations in vitro cloning |
US9920361B2 (en) | 2012-05-21 | 2018-03-20 | Sequenom, Inc. | Methods and compositions for analyzing nucleic acid |
US9914967B2 (en) | 2012-06-05 | 2018-03-13 | President And Fellows Of Harvard College | Spatial sequencing of nucleic acids using DNA origami probes |
CN104685116A (en) | 2012-06-25 | 2015-06-03 | Gen9股份有限公司 | Methods for nucleic acid assembly and high throughput sequencing |
WO2014007623A1 (en) | 2012-07-03 | 2014-01-09 | Interna Technologies B.V. | Diagnostic portfolio and its uses |
JP2015521862A (en) | 2012-07-13 | 2015-08-03 | セクエノム, インコーポレイテッド | Process and composition for enrichment based on methylation of fetal nucleic acid from maternal samples useful for non-invasive prenatal diagnosis |
US9701998B2 (en) | 2012-12-14 | 2017-07-11 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10323279B2 (en) | 2012-08-14 | 2019-06-18 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10273541B2 (en) | 2012-08-14 | 2019-04-30 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10752949B2 (en) | 2012-08-14 | 2020-08-25 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US11591637B2 (en) | 2012-08-14 | 2023-02-28 | 10X Genomics, Inc. | Compositions and methods for sample processing |
US10584381B2 (en) | 2012-08-14 | 2020-03-10 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
MX364957B (en) | 2012-08-14 | 2019-05-15 | 10X Genomics Inc | Microcapsule compositions and methods. |
US9951386B2 (en) | 2014-06-26 | 2018-04-24 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10221442B2 (en) | 2012-08-14 | 2019-03-05 | 10X Genomics, Inc. | Compositions and methods for sample processing |
EP2904115B1 (en) | 2012-10-01 | 2018-08-08 | Millennium Pharmaceuticals, Inc. | Biomarkers and methods to predict response to inhibitors and uses thereof |
WO2014055117A1 (en) | 2012-10-04 | 2014-04-10 | Asuragen, Inc. | Diagnostic mirnas for differential diagnosis of incidental pancreatic cystic lesions |
GB201220924D0 (en) | 2012-11-21 | 2013-01-02 | Cancer Res Inst Royal | Materials and methods for determining susceptibility or predisposition to cancer |
BR112015012295A8 (en) | 2012-11-28 | 2023-03-14 | Merck Sharp & Dohme | USE OF A WEE1 INHIBITOR, E, KIT TO IDENTIFY A PATIENT WITH CANCER |
KR20150090246A (en) | 2012-12-03 | 2015-08-05 | 알막 다이아그노스틱스 리미티드 | Molecular diagnostic test for cancer |
CA2894694C (en) | 2012-12-14 | 2023-04-25 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10533221B2 (en) | 2012-12-14 | 2020-01-14 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US9914966B1 (en) | 2012-12-20 | 2018-03-13 | Nabsys 2.0 Llc | Apparatus and methods for analysis of biomolecules using high frequency alternating current excitation |
EP2956550B1 (en) | 2013-01-18 | 2020-04-08 | Nabsys 2.0 LLC | Enhanced probe binding |
US9896728B2 (en) | 2013-01-29 | 2018-02-20 | Arcticrx Ltd. | Method for determining a therapeutic approach for the treatment of age-related macular degeneration (AMD) |
BR112015019159A2 (en) | 2013-02-08 | 2017-07-18 | 10X Genomics Inc | polynucleotide barcode generation |
US10138509B2 (en) | 2013-03-12 | 2018-11-27 | President And Fellows Of Harvard College | Method for generating a three-dimensional nucleic acid containing matrix |
EP3597774A1 (en) | 2013-03-13 | 2020-01-22 | Sequenom, Inc. | Primers for dna methylation analysis |
WO2014143158A1 (en) | 2013-03-13 | 2014-09-18 | The Broad Institute, Inc. | Compositions and methods for labeling of agents |
WO2014159063A1 (en) | 2013-03-14 | 2014-10-02 | Elitech Holding B.V. | Functionalized 3-alkynyl pyrazolopyrimidine analogues as universal bases and methods of use |
EP2971132B1 (en) | 2013-03-15 | 2020-05-06 | Baylor Research Institute | Tissue and blood-based mirna biomarkers for the diagnosis, prognosis and metastasis-predictive potential in colorectal cancer |
EP3404116B1 (en) | 2013-03-15 | 2022-10-19 | The University of Chicago | Methods and compositions related to t-cell activity |
EP2971149B1 (en) | 2013-03-15 | 2018-05-09 | Baylor Research Institute | Ulcerative colitis (uc)-associated colorectal neoplasia markers |
EP3013983B1 (en) | 2013-06-25 | 2023-02-15 | Prognosys Biosciences, Inc. | Spatially encoded biological assays using a microfluidic device |
AU2014306867B2 (en) | 2013-08-12 | 2017-10-26 | Genentech, Inc. | Compositions and method for treating complement-associated conditions |
US11136625B2 (en) | 2013-08-28 | 2021-10-05 | Crown Bioscience, Inc. (Taicang) | Gene expression signatures predictive of subject response to a multi-kinase inhibitor and methods of using the same |
US10395758B2 (en) | 2013-08-30 | 2019-08-27 | 10X Genomics, Inc. | Sequencing methods |
US10655173B2 (en) | 2013-10-18 | 2020-05-19 | The Broad Institute, Inc. | Spatial and cellular mapping of biomolecules in situ by high-throughput sequencing |
WO2015069634A1 (en) | 2013-11-08 | 2015-05-14 | President And Fellows Of Harvard College | Microparticles, methods for their preparation and use |
US9824068B2 (en) | 2013-12-16 | 2017-11-21 | 10X Genomics, Inc. | Methods and apparatus for sorting data |
US20160304969A1 (en) | 2013-12-17 | 2016-10-20 | Merck Sharp & Dohme Corp. | Ifn-gamma gene signature biomarkers of tumor response to pd-1 antagonists |
GB201409479D0 (en) | 2014-05-28 | 2014-07-09 | Almac Diagnostics Ltd | Molecular diagnostic test for cancer |
CA2938731A1 (en) | 2014-02-08 | 2015-08-13 | Genentech, Inc. | Methods of treating alzheimer's disease |
US11365447B2 (en) | 2014-03-13 | 2022-06-21 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
CN114534806B (en) | 2014-04-10 | 2024-03-29 | 10X基因组学有限公司 | Fluidic devices, systems and methods for packaging and partitioning reagents and uses thereof |
SG11201610691QA (en) | 2014-06-26 | 2017-01-27 | 10X Genomics Inc | Processes and systems for nucleic acid sequence assembly |
WO2015200893A2 (en) | 2014-06-26 | 2015-12-30 | 10X Genomics, Inc. | Methods of analyzing nucleic acids from individual cells or cell populations |
US20180363070A9 (en) | 2014-08-12 | 2018-12-20 | Univ Wayne State | Systems and methods to detect stem cell stress and uses thereof |
CN107002128A (en) | 2014-10-29 | 2017-08-01 | 10X 基因组学有限公司 | The method and composition being sequenced for target nucleic acid |
US9975122B2 (en) | 2014-11-05 | 2018-05-22 | 10X Genomics, Inc. | Instrument systems for integrated sample processing |
US10507210B2 (en) | 2014-12-03 | 2019-12-17 | Auckland Uniservices Limited | Kinase inhibitor prodrug for the treatment of cancer |
US11543411B2 (en) | 2014-12-05 | 2023-01-03 | Prelude Corporation | DCIS recurrence and invasive breast cancer |
RU2017123117A (en) | 2014-12-09 | 2019-01-10 | Мерк Шарп И Доум Корп. | SYSTEM AND METHODS FOR PRODUCING BIOMARKERS OF GENE SIGNATURES OF RESPONSE TO PD-1 ANTAGONISTS |
SG11201705615UA (en) | 2015-01-12 | 2017-08-30 | 10X Genomics Inc | Processes and systems for preparing nucleic acid sequencing libraries and libraries prepared using same |
CA2968417A1 (en) | 2015-01-13 | 2016-07-21 | 10X Genomics, Inc. | Systems and methods for visualizing structural variation and phasing information |
CN107208156B (en) | 2015-02-09 | 2021-10-08 | 10X基因组学有限公司 | System and method for determining structural variation and phasing using variation recognition data |
US10697000B2 (en) | 2015-02-24 | 2020-06-30 | 10X Genomics, Inc. | Partition processing methods and systems |
CA2975958A1 (en) | 2015-02-24 | 2016-09-01 | 10X Genomics, Inc. | Methods for targeted nucleic acid sequence coverage |
FI3901281T3 (en) | 2015-04-10 | 2023-01-31 | Spatially distinguished, multiplex nucleic acid analysis of biological specimens | |
EP3283882B1 (en) | 2015-04-17 | 2020-12-16 | Merck Sharp & Dohme Corp. | Blood-based biomarkers of tumor sensitivity to pd-1 antagonists |
GB201512869D0 (en) | 2015-07-21 | 2015-09-02 | Almac Diagnostics Ltd | Gene signature for minute therapies |
US11123297B2 (en) | 2015-10-13 | 2021-09-21 | President And Fellows Of Harvard College | Systems and methods for making and using gel microspheres |
WO2017075294A1 (en) | 2015-10-28 | 2017-05-04 | The Board Institute Inc. | Assays for massively combinatorial perturbation profiling and cellular circuit reconstruction |
MX2018005611A (en) | 2015-11-03 | 2018-11-09 | Harvard College | Method and apparatus for volumetric imaging of a three-dimensional nucleic acid containing matrix. |
US10522243B2 (en) * | 2015-11-20 | 2019-12-31 | Bio-Rad Laboratories, Inc. | Sparse identity spaces in droplet sequencing |
EP4144861A1 (en) | 2015-12-04 | 2023-03-08 | 10X Genomics, Inc. | Methods and compositions for nucleic acid analysis |
CN108779491B (en) | 2016-02-11 | 2021-03-09 | 10X基因组学有限公司 | Systems, methods, and media for de novo assembly of whole genome sequence data |
JP6707181B2 (en) | 2016-04-20 | 2020-06-10 | エルディエックス・プログノスティクス・リミテッド・カンパニーLdx Prognostics Limited Co. | Kits or packages for identifying pregnant women at risk of preterm birth and use of such kits or packages |
WO2017189525A1 (en) | 2016-04-25 | 2017-11-02 | President And Fellows Of Harvard College | Hybridization chain reaction methods for in situ molecular detection |
WO2017197338A1 (en) | 2016-05-13 | 2017-11-16 | 10X Genomics, Inc. | Microfluidic systems and methods of use |
WO2017197573A1 (en) | 2016-05-17 | 2017-11-23 | Ldx Prognostics Limited Co. | Methods and compositions for providing preeclampsia assessment |
US11352667B2 (en) | 2016-06-21 | 2022-06-07 | 10X Genomics, Inc. | Nucleic acid sequencing |
GB2569252A (en) | 2016-08-31 | 2019-06-12 | Harvard College | Methods of combining the detection of biomolecules into a single assay using fluorescent in situ sequencing |
GB2570412A (en) | 2016-08-31 | 2019-07-24 | Harvard College | Methods of generating libraries of nucleic acid sequences for detection via fluorescent in situ sequencing |
US10815525B2 (en) | 2016-12-22 | 2020-10-27 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10550429B2 (en) | 2016-12-22 | 2020-02-04 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
US10011872B1 (en) | 2016-12-22 | 2018-07-03 | 10X Genomics, Inc. | Methods and systems for processing polynucleotides |
WO2018140966A1 (en) | 2017-01-30 | 2018-08-02 | 10X Genomics, Inc. | Methods and systems for droplet-based single cell barcoding |
JP2020517715A (en) | 2017-04-28 | 2020-06-18 | メルク・シャープ・エンド・ドーム・コーポレイション | Biomarkers for cancer therapy |
CN110870018A (en) | 2017-05-19 | 2020-03-06 | 10X基因组学有限公司 | System and method for analyzing a data set |
US10844372B2 (en) | 2017-05-26 | 2020-11-24 | 10X Genomics, Inc. | Single cell analysis of transposase accessible chromatin |
EP4230746A3 (en) | 2017-05-26 | 2023-11-01 | 10X Genomics, Inc. | Single cell analysis of transposase accessible chromatin |
US11597965B2 (en) | 2017-10-06 | 2023-03-07 | 10X Genomics, Inc. | RNA templated ligation |
US11497762B2 (en) | 2017-11-03 | 2022-11-15 | Interna Technologies B.V. | MiRNA molecule, equivalent, antagomir, or source thereof for treating and/or diagnosing a condition and/or a disease associated with neuronal deficiency or for neuronal (re)generation |
SG11201913654QA (en) | 2017-11-15 | 2020-01-30 | 10X Genomics Inc | Functionalized gel beads |
US10829815B2 (en) | 2017-11-17 | 2020-11-10 | 10X Genomics, Inc. | Methods and systems for associating physical and genetic properties of biological particles |
US20200377958A1 (en) | 2017-12-01 | 2020-12-03 | Millennium Pharmaceuticals, Inc. | Biomarkers and methods for treatment with nae inhibitors |
US11332736B2 (en) | 2017-12-07 | 2022-05-17 | The Broad Institute, Inc. | Methods and compositions for multiplexing single cell and single nuclei sequencing |
EP3728321A1 (en) | 2017-12-22 | 2020-10-28 | F. Hoffmann-La Roche AG | Use of pilra binding agents for treatment of a disease |
CN112262218A (en) | 2018-04-06 | 2021-01-22 | 10X基因组学有限公司 | System and method for quality control in single cell processing |
WO2019213660A2 (en) | 2018-05-04 | 2019-11-07 | The Broad Institute, Inc. | Compositions and methods for modulating cgrp signaling to regulate innate lymphoid cell inflammatory responses |
CN112770776A (en) | 2018-07-30 | 2021-05-07 | 瑞德库尔有限责任公司 | Method and system for sample processing or analysis |
EP3837380A4 (en) | 2018-08-17 | 2022-05-11 | Cellecta, Inc. | Multiplex preparation of barcoded gene specific dna fragments |
AU2019339508A1 (en) | 2018-09-14 | 2021-04-15 | Prelude Corporation | Method of selection for treatment of subjects at risk of invasive breast cancer |
US20220401460A1 (en) | 2018-10-10 | 2022-12-22 | Dana-Farber Cancer Institute, Inc. | Modulating resistance to bcl-2 inhibitors |
EP3894590A2 (en) | 2018-12-10 | 2021-10-20 | 10X Genomics, Inc. | Methods of using master / copy arrays for spatial detection |
US11926867B2 (en) | 2019-01-06 | 2024-03-12 | 10X Genomics, Inc. | Generating capture probes for spatial analysis |
US11649485B2 (en) | 2019-01-06 | 2023-05-16 | 10X Genomics, Inc. | Generating capture probes for spatial analysis |
CA3131740A1 (en) | 2019-02-27 | 2020-09-03 | Madera Therapeutics, LLC | Use of caseinolytic protease p function as a biomarker of drug response to imipridone-like agents |
US20220154282A1 (en) | 2019-03-12 | 2022-05-19 | The Broad Institute, Inc. | Detection means, compositions and methods for modulating synovial sarcoma cells |
US20220143148A1 (en) | 2019-03-14 | 2022-05-12 | The Broad Institute, Inc. | Compositions and methods for modulating cgrp signaling to regulate intestinal innate lymphoid cells |
US20220152148A1 (en) | 2019-03-18 | 2022-05-19 | The Broad Institute, Inc. | Modulation of type 2 immunity by targeting clec-2 signaling |
US20220142948A1 (en) | 2019-03-18 | 2022-05-12 | The Broad Institute, Inc. | Compositions and methods for modulating metabolic regulators of t cell pathogenicity |
US20220162198A1 (en) | 2019-04-12 | 2022-05-26 | The Regents Of The University Of California | Compositions and methods for increasing muscle mass and oxidative metabolism |
CA3139791A1 (en) | 2019-05-31 | 2020-12-03 | 10X Genomics, Inc. | Method of detecting target nucleic acid molecules |
CN111406117B (en) * | 2019-08-09 | 2023-06-30 | 北京大北农生物技术有限公司 | Nucleic acid sequence for detecting soybean plant DBN8002 and detection method thereof |
WO2021030627A1 (en) | 2019-08-13 | 2021-02-18 | The General Hospital Corporation | Methods for predicting outcomes of checkpoint inhibition and treatment thereof |
US11793787B2 (en) | 2019-10-07 | 2023-10-24 | The Broad Institute, Inc. | Methods and compositions for enhancing anti-tumor immunity by targeting steroidogenesis |
US11844800B2 (en) | 2019-10-30 | 2023-12-19 | Massachusetts Institute Of Technology | Methods and compositions for predicting and preventing relapse of acute lymphoblastic leukemia |
WO2021092433A2 (en) | 2019-11-08 | 2021-05-14 | 10X Genomics, Inc. | Enhancing specificity of analyte binding |
US20230044602A1 (en) | 2019-12-20 | 2023-02-09 | EDWARD Via COLLEGE OF OSTEOPATHIC MEDICINE | Cancer signatures, methods of generating cancer signatures, and uses thereof |
CN114885610A (en) | 2019-12-23 | 2022-08-09 | 10X基因组学有限公司 | Methods for spatial analysis using RNA templated ligation |
US11702693B2 (en) | 2020-01-21 | 2023-07-18 | 10X Genomics, Inc. | Methods for printing cells and generating arrays of barcoded cells |
US11732299B2 (en) | 2020-01-21 | 2023-08-22 | 10X Genomics, Inc. | Spatial assays with perturbed cells |
US11898205B2 (en) | 2020-02-03 | 2024-02-13 | 10X Genomics, Inc. | Increasing capture efficiency of spatial assays |
US11732300B2 (en) | 2020-02-05 | 2023-08-22 | 10X Genomics, Inc. | Increasing efficiency of spatial analysis in a biological sample |
US11891654B2 (en) | 2020-02-24 | 2024-02-06 | 10X Genomics, Inc. | Methods of making gene expression libraries |
EP4242325A3 (en) | 2020-04-22 | 2023-10-04 | 10X Genomics, Inc. | Methods for spatial analysis using targeted rna depletion |
WO2021237087A1 (en) | 2020-05-22 | 2021-11-25 | 10X Genomics, Inc. | Spatial analysis to detect sequence variants |
EP4153775A1 (en) | 2020-05-22 | 2023-03-29 | 10X Genomics, Inc. | Simultaneous spatio-temporal measurement of gene expression and cellular activity |
EP4162074A1 (en) | 2020-06-08 | 2023-04-12 | 10X Genomics, Inc. | Methods of determining a surgical margin and methods of use thereof |
EP4172362A1 (en) | 2020-06-25 | 2023-05-03 | 10X Genomics, Inc. | Spatial analysis of dna methylation |
US11761038B1 (en) | 2020-07-06 | 2023-09-19 | 10X Genomics, Inc. | Methods for identifying a location of an RNA in a biological sample |
US11926822B1 (en) | 2020-09-23 | 2024-03-12 | 10X Genomics, Inc. | Three-dimensional spatial analysis |
US11827935B1 (en) | 2020-11-19 | 2023-11-28 | 10X Genomics, Inc. | Methods for spatial analysis using rolling circle amplification and detection probes |
WO2022140028A1 (en) | 2020-12-21 | 2022-06-30 | 10X Genomics, Inc. | Methods, compositions, and systems for capturing probes and/or barcodes |
WO2022192419A2 (en) | 2021-03-09 | 2022-09-15 | Massachusetts Institute Of Technology | Methods of treating inflammatory bowel disease (ibd) with anti- tnf-blockade |
WO2022197933A1 (en) | 2021-03-18 | 2022-09-22 | The Broad Institute, Inc. | Compositions and methods for characterizing lymphoma and related conditions |
WO2023034489A1 (en) | 2021-09-01 | 2023-03-09 | 10X Genomics, Inc. | Methods, compositions, and kits for blocking a capture probe on a spatial array |
WO2023122723A1 (en) | 2021-12-23 | 2023-06-29 | The Broad Institute, Inc. | Panels and methods for diagnosing and treating lung cancer |
WO2023200843A1 (en) | 2022-04-12 | 2023-10-19 | The Broad Institute, Inc. | Compositions and methods for screening cis regulatory elements |
WO2024028794A1 (en) | 2022-08-02 | 2024-02-08 | Temple Therapeutics BV | Methods for treating endometrial and ovarian hyperproliferative disorders |
Family Cites Families (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4562159A (en) * | 1981-03-31 | 1985-12-31 | Albert Einstein College Of Medicine, A Division Of Yeshiva Univ. | Diagnostic test for hepatitis B virus |
US4591567A (en) * | 1982-04-21 | 1986-05-27 | California Institute Of Technology | Recombinant DNA screening system including fixed array replicator and support |
US5260190A (en) * | 1982-11-05 | 1993-11-09 | Fuji Photo Film Co., Ltd. | Autoradiographic process |
US4888695A (en) * | 1983-01-08 | 1989-12-19 | Fuji Photo Film Co., Ltd. | Signal processing method in autoradiography |
EP0113677A3 (en) * | 1983-01-08 | 1989-02-01 | Fuji Photo Film Co., Ltd. | Signal processing method in autoradiography |
EP0113672A3 (en) * | 1983-01-08 | 1988-03-02 | Fuji Photo Film Co., Ltd. | Signal processing method in autoradiography |
US4672040A (en) * | 1983-05-12 | 1987-06-09 | Advanced Magnetics, Inc. | Magnetic particles for use in separations |
JPS6010174A (en) * | 1983-06-29 | 1985-01-19 | Fuji Photo Film Co Ltd | Screening method of gene by auto radiography |
US5270162A (en) * | 1983-09-19 | 1993-12-14 | Fuji Photo Film Co., Ltd. | Autoradiographic gene-screening method |
US4613566A (en) * | 1984-01-23 | 1986-09-23 | President And Fellows Of Harvard College | Hybridization assay and kit therefor |
CA1223222A (en) * | 1984-02-22 | 1987-06-23 | Nanibhushan Dattagupta | Immobilized nucleic acid-containing probes |
US4766062A (en) * | 1984-05-07 | 1988-08-23 | Allied Corporation | Displacement polynucleotide assay method and polynucleotide complex reagent therefor |
US4675283A (en) * | 1984-07-19 | 1987-06-23 | Massachusetts Institute Of Technology | Detection and isolation of homologous, repeated and amplified nucleic acid sequences |
DE3506703C1 (en) * | 1985-02-26 | 1986-04-30 | Sagax Instrument AB, Sundbyberg | Process for sequence analysis of nucleic acids, in particular deoxyribonucleic acid (DNA) and ribonucleic acid (RNA), as well as carriers for carrying out the method and process for producing the carrier |
US4683195A (en) * | 1986-01-30 | 1987-07-28 | Cetus Corporation | Process for amplifying, detecting, and/or-cloning nucleic acid sequences |
US4683202A (en) * | 1985-03-28 | 1987-07-28 | Cetus Corporation | Process for amplifying nucleic acid sequences |
US4865968A (en) * | 1985-04-01 | 1989-09-12 | The Salk Institute For Biological Studies | DNA sequencing |
DE3684030D1 (en) * | 1985-04-19 | 1992-04-09 | Fuji Photo Film Co Ltd | SIGNAL PROCESSING METHOD FOR DETECTING THE ORDER OF THE BASES OF NUCLEIC ACID. |
US4941092A (en) * | 1985-05-23 | 1990-07-10 | Fuji Photo Film Co., Ltd. | Signal processing method for determining base sequence of nucleic acid |
US4794073A (en) * | 1985-07-10 | 1988-12-27 | Molecular Diagnostics, Inc. | Detection of nucleic acid hybrids by prolonged chemiluminescence |
JPS6242057A (en) * | 1985-08-19 | 1987-02-24 | Fuji Photo Film Co Ltd | Signal processing method for determining base sequence of nucleic acid |
JPS6285861A (en) * | 1985-10-11 | 1987-04-20 | Fuji Photo Film Co Ltd | Signal processing method for determining base sequence of nucleic acid |
US4741043B1 (en) * | 1985-11-04 | 1994-08-09 | Cell Analysis Systems Inc | Method of and apparatus for image analyses of biological specimens |
US4770992A (en) * | 1985-11-27 | 1988-09-13 | Den Engh Gerrit J Van | Detection of specific DNA sequences by flow cytometry |
EP0240729A3 (en) * | 1986-03-05 | 1988-08-24 | Fuji Photo Film Co., Ltd. | Method of analyzing autoradiograph for determining base sequence of nucleic acid |
US4885696A (en) * | 1986-03-26 | 1989-12-05 | Fuji Photo Film Co., Ltd. | Signal processing method for determining base sequence of nucleic acid |
JPS62228165A (en) * | 1986-03-29 | 1987-10-07 | Fuji Photo Film Co Ltd | Signal processing method for determining base sequence of nucleic acid |
US4811218A (en) * | 1986-06-02 | 1989-03-07 | Applied Biosystems, Inc. | Real time scanning electrophoresis apparatus for DNA sequencing |
CA1340806C (en) * | 1986-07-02 | 1999-11-02 | James Merrill Prober | Method, system and reagents for dna sequencing |
ATE88762T1 (en) * | 1986-08-11 | 1993-05-15 | Siska Diagnostics Inc | METHODS AND COMPOSITIONS FOR TESTING WITH NUCLEIC ACID PROBES. |
JPS63167290A (en) * | 1986-12-27 | 1988-07-11 | Fuji Photo Film Co Ltd | Signal processing method for autoradiographic analysis |
JPH0664057B2 (en) * | 1987-01-06 | 1994-08-22 | 富士写真フイルム株式会社 | Signal processing method for autoradiographic analysis |
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 |
US4849334A (en) * | 1987-06-09 | 1989-07-18 | Life Technologies, Inc. | Human papillomavirus 43 nucleic acid hybridization probes and methods for employing the same |
US4942124A (en) * | 1987-08-11 | 1990-07-17 | President And Fellows Of Harvard College | Multiplex sequencing |
US5149625A (en) * | 1987-08-11 | 1992-09-22 | President And Fellows Of Harvard College | Multiplex analysis of DNA |
GB8810400D0 (en) * | 1988-05-03 | 1988-06-08 | Southern E | Analysing polynucleotide sequences |
US4965725B1 (en) * | 1988-04-08 | 1996-05-07 | Neuromedical 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 |
GB8822228D0 (en) * | 1988-09-21 | 1988-10-26 | Southern E M | Support-bound oligonucleotides |
US5297288A (en) * | 1989-11-28 | 1994-03-22 | United States Biochemical Corporation | System for use with a high resolution scanner for scheduling a sequence of software tools for determining the presence of bands in DNA sequencing samples |
-
1994
- 1994-02-28 US US08/203,502 patent/US5525464A/en not_active Expired - Lifetime
-
1995
- 1995-06-05 US US08/460,853 patent/US5695940A/en not_active Expired - Fee Related
-
1998
- 1998-10-22 US US09/176,819 patent/US5972619A/en not_active Expired - Fee Related
- 1998-12-04 US US09/205,865 patent/US6316191B1/en not_active Expired - Fee Related
-
2001
- 2001-09-10 US US09/954,426 patent/US20020106673A1/en not_active Abandoned
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9540689B2 (en) | 1998-05-01 | 2017-01-10 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US9096898B2 (en) | 1998-05-01 | 2015-08-04 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US7645596B2 (en) | 1998-05-01 | 2010-01-12 | Arizona Board Of Regents | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US10208341B2 (en) | 1998-05-01 | 2019-02-19 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US9957561B2 (en) | 1998-05-01 | 2018-05-01 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US9725764B2 (en) | 1998-05-01 | 2017-08-08 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US10214774B2 (en) | 1998-05-01 | 2019-02-26 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US9458500B2 (en) | 1998-05-01 | 2016-10-04 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US9212393B2 (en) | 1998-05-01 | 2015-12-15 | Life Technologies Corporation | Method of determining the nucleotide sequence of oligonucleotides and DNA molecules |
US20060008839A1 (en) * | 2001-08-29 | 2006-01-12 | Pioneer Hi-Bred International, Inc. | UTR tag assay for gene function discovery |
US7405282B2 (en) | 2001-08-29 | 2008-07-29 | Pioneer Hi-Bred International, Inc. | UTR tag assay for gene function discovery |
US7026123B1 (en) | 2001-08-29 | 2006-04-11 | Pioneer Hi-Bred International, Inc. | UTR tag assay for gene function discovery |
US9657344B2 (en) | 2003-11-12 | 2017-05-23 | Fluidigm Corporation | Short cycle methods for sequencing polynucleotides |
US7897345B2 (en) | 2003-11-12 | 2011-03-01 | Helicos Biosciences Corporation | Short cycle methods for sequencing polynucleotides |
US9012144B2 (en) | 2003-11-12 | 2015-04-21 | Fluidigm Corporation | Short cycle methods for sequencing polynucleotides |
US7981604B2 (en) | 2004-02-19 | 2011-07-19 | California Institute Of Technology | Methods and kits for analyzing polynucleotide sequences |
US20050266447A1 (en) * | 2004-04-19 | 2005-12-01 | Pioneer Hi-Bred International, Inc. | Method for identifying activators of gene transcription |
US9868978B2 (en) | 2005-08-26 | 2018-01-16 | Fluidigm Corporation | Single molecule sequencing of captured nucleic acids |
US7666593B2 (en) | 2005-08-26 | 2010-02-23 | Helicos Biosciences Corporation | Single molecule sequencing of captured nucleic acids |
Also Published As
Publication number | Publication date |
---|---|
US6316191B1 (en) | 2001-11-13 |
US5972619A (en) | 1999-10-26 |
US5525464A (en) | 1996-06-11 |
US5695940A (en) | 1997-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6316191B1 (en) | Computer-aided analysis system for sequencing by hybridization | |
US7230093B2 (en) | Method of sequencing by hybridization of oligonucleotide probes | |
US6403315B1 (en) | Method and apparatus for DNA sequencing and DNA identification | |
Dramanac et al. | Sequencing of megabase plus DNA by hybridization: theory of the method | |
US6291182B1 (en) | Methods, software and apparati for identifying genomic regions harboring a gene associated with a detectable trait | |
US5492806A (en) | Method of determining an ordered sequence of subfragments of a nucleic acid fragment by hybridization of oligonucleotide probes | |
CN106591441B (en) | Alpha and/or beta-thalassemia mutation detection probe, method and chip based on whole gene capture sequencing and application | |
Dover et al. | [38] Detection and quantification of concerted evolution and molecular drive | |
EP1056889B1 (en) | Methods related to genotyping and dna analysis | |
US10648030B2 (en) | Methods of determining the presence or absence of a plurality of target polynucleotides in a sample | |
AU746682B2 (en) | Biallelic markers for use in constructing a high density disequilibrium map of the human genome | |
Arcot et al. | High-resolution cartography of recently integrated human chromosome 19-specific Alu fossils | |
Hong et al. | Genomic organization of Trypanosoma brucei kinetoplast DNA minicircles | |
CN112280863A (en) | Method and kit for effectiveness of targeted drug apatinib | |
Yager et al. | The human genome project | |
Weir | Statistical analysis of DNA sequences | |
WO1999058721A1 (en) | Multiplex dna amplification using chimeric primers | |
CN109825552B (en) | Primer and method for enriching target region | |
Ruchi | Bioinformatics: genomics and proteomics | |
Zischler | Mitochondrial DNA: Diversity analysis and possible pitfalls | |
Singh et al. | Sequence-based markers | |
Consevage et al. | Basic elements of gene mapping and identification | |
Chauhan et al. | Genomics: The Rise of Sequencing Technology | |
Barnes | 8 G aining Insight into Human Population-Specific Selection Pressure | |
Dizikes | Update on the human genome project |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |