Patent Grant
8050516
The present invention relates generally to microfluidic techniques. In particular, the invention provides a method and system for imaging one or more entities suspended in a volume of fluid in a chamber of a microfluidic device. More particularly, the present method and system for imaging uses indications from a fluorescence signal associated with the one or more entities in the microfluidic device. Merely by way of example, the techniques for microfluidic methods and systems are applied using fluorescent, chemiluminescent, and bioluminescent readers coupled to the microfluidic device, but it would be recognized that the invention has a much broader range of applicability.
Concerted efforts to develop and manufacture microfluidic systems to perform various chemical and biochemical analyses and syntheses have occurred. Such systems have been developed for preparative and analytical applications. A goal to make such micro-sized devices arises from significant benefits achieved from miniaturization of conventional macro scale analyses and syntheses, which are often cumbersome and less efficient. A substantial reduction in time, lower costs, and more efficient space allocation are achieved as benefits using these microfluidic systems. Additional benefits may include a reduction in human operator involvement with automated systems using these microfluidic devices. Automated systems also decrease operator errors and other operator type limitations. Microfluidic devices have been proposed for use in a variety of applications, including, for instance, capillary electrophoresis, gas chromatography and cell separations.
Microfluidic devices adapted to conduct nucleic acid amplification processes are potentially useful in a wide variety of applications. For example, such devices could be used to determine the presence or absence of a particular target nucleic acid in a sample, as an analytical tool. Examples of utilizing microfluidic device as an analytical tool include:
Alternatively, the devices can be used in a preparative fashion to amplify nucleic acids, producing an amplified product at sufficient levels needed for further analysis. Examples of these analysis processes include sequencing of the amplified product, cell-typing, DNA fingerprinting, and the like. Amplified products can also be used in various genetic engineering applications. These genetic engineering applications include (but are not limited to) the production of a desired protein product, accomplished by insertion of the amplified product into a vector that is then used to transform cells into the desired protein product.
Despite these potential applications, imaging systems (also referred to as readers) adapted to collect and process imaging data, for example, fluorescence data, from such microfluidic devices have various shortcomings related to the collection and processing of images. Therefore, there is a need in the art for improved methods and systems for imaging one or more entities suspended in a volume of fluid in a chamber of a microfluidic device.
According to the present invention, techniques for microfluidic systems are provided. In particular, the invention provides a method and system for imaging one or more entities suspended in a volume of fluid in a chamber of a microfluidic device. More particularly, the present method and system for imaging uses indications from a fluorescence signal associated with the one or more entities in the microfluidic device. Merely by way of example, the techniques for microfluidic methods and systems are applied using fluorescent, chemiluminescent, and bioluminescent readers coupled to the microfluidic device, but it would be recognized that the invention has a much broader range of applicability.
According to a specific embodiment according to the present invention, a method of processing data associated with fluorescent emissions from a microfluidic device is provided. The method includes performing an auto-focus process associated with a first image of the microfluidic device and performing an auto-exposure process associated with the first image of the microfluidic device. The method also includes capturing a plurality of images of the microfluidic device. The plurality of images are associated with a plurality of thermal cycles. The method further includes performing image analysis of the plurality of captured images to determine a series of optical intensities and performing data analysis of the series of optical intensities to provide a series of change in threshold values.
According to another embodiment of the present invention, a method of determining a focal plane for imaging a microfluidic device comprising a plurality of reaction chambers is provided. The method includes acquiring a stack of images. Each image in the stack of images is associated with a predetermined focal distance. The method also includes determining, for each image in the stack of images, a contrast score associated with each image and selecting a baseline image associated with a contrast score selected from the contrast scores associated with each image. The method further includes storing a focal distance associated with the baseline image and determining a plurality of positions associated with the plurality of reaction chambers. Each of the plurality of positions includes a number of chamber pixels. Moreover, the method includes determining, for each image in the stack of images, a contrast score associated with the plurality of positions associated with the plurality of reaction chambers and selecting an updated image associated with a contrast score selected from the contrast scores associated with the plurality of positions.
According to an alternative embodiment of the present invention, a method of performing an auto-focus process for a microfluidic device including a plurality of chambers is provided. The method includes acquiring an image stack comprising a plurality of images, determining a focal plane characterized by a predetermined contrast score, and determining a plurality of positions associated with the plurality of chambers. The method also includes computing a contrast score for each of the plurality of images in the image stack, determining an updated focal plane characterized by a second predetermined contrast score, and storing the updated focal plane in one or more memories.
According to another alternative embodiment of the present invention, a method of determining a baseline value associated with fluorescent emissions from a microfluidic device is provided. The method includes computing an initial baseline value, providing a first baseline value having a different value than the initial baseline value, and providing a series of emission data points as a function of PCR cycle number. The method also includes subtracting the initial baseline value from the series of emission data points to provide an initial baseline removed curve and computing an initial quality factor associated with initial baseline removed curve. The method further includes subtracting the first baseline value from the series of emission data points to provide a first baseline removed curve and computing a first quality factor associated with the first baseline removed curve. Moreover, the method includes determining that the first quality factor is greater than the second quality factor and determining that the baseline value is the first baseline value.
According to a specific embodiment of the present invention, a method of performing an auto-exposure process for an imaging system adapted to characterize a microfluidic device is provided. The method includes providing a histogram of pixel intensities, computing an average intensity for a range of pixels, and computing an exposure time. The method also includes computing a number of saturated pixels, and determining that the exposure time is less than a predetermined time and that the number of saturated pixels is less than a predetermined value.
According to another specific embodiment of the present invention, a method of determining a plurality of positions associated with a plurality of reaction chambers of a microfluidic device is provided. The method includes a) providing a baseline image, b) providing a template image of a reaction chamber, and c) selecting a region of the baseline image. The method also includes d) performing a matching process including matching the template image to one or more portions of the region of the baseline image, e) determining a position of a first chamber, and f) predicting a position of a second chamber. The method further includes g) repeating steps c) through f) for subsequent chambers.
According to yet another embodiment of the present invention, a method of removing a baseline from PCR ratio data is provided. The method includes obtaining a plurality of images. Each of the plurality of images is associated with a cycle of a PCR process. The method also includes determining a ratio of reporter dye to passive reference dye for each of the plurality of images, thereby providing a ratio data curve. The method further includes determining a corner point of the ratio data curve, determining a baseline value associated with the ratio data curve, and subtracting the baseline value from each value in the ratio data curve, thereby providing a baseline removed ratio data curve.
According to yet another specific embodiment of the present invention, a method for determining a ratio data value associated with a microfluidic device is provided. The method includes measuring a passive reference image for the microfluidic device. The passive reference image includes an array of reference intensity values. The method also includes measuring a reporter dye image for the microfluidic device. The reporter dye image includes an array of reporter intensity values. The method further includes determining a position of a reaction chamber of the microfluidic device, determining a low intensity value for the passive reference image in a vicinity of the position of the reaction chamber of the microfluidic device, and determining a low intensity value for the reporter dye image in a vicinity of the position of the reaction chamber of the microfluidic device. Additionally, the method includes subtracting the low intensity value for the passive reference image from the array of reference intensity values to provide a denominator value, subtracting the low intensity value for the reporter dye image from the array of reporter intensity values to provide a numerator value, and computing a ratio data value equal to the numerator value divided by the denominator value.
According to yet another alternative embodiment of the present invention, a method of determining a change in threshold value for a plurality of ratio data curves is provided. The method includes determining a first quality factor for a first ratio data curve of the plurality of ratio data curves, determining that the first quality factor is greater than a predetermined value, and determining a second quality factor for a second ratio data curve of the plurality of ratio data curves. The method also includes determining that the second quality factor is greater than the predetermined value and determining the change in threshold value based on the first ratio data curve and the second ratio data curve.
According to a particular embodiment of the present invention, a method of baseline correcting a ratio data curve associated with a PCR process for a microfluidic device is provided. The method includes providing a plurality of baseline values for the ratio data curve, computing a plurality of baseline removed curves, and computing a plurality of fitting errors for the plurality of baseline removed curves. The method also includes determining the minimum fitting error of the plurality of fitting errors, thereby determining a baseline value and computing the baseline corrected ratio data curve by removing the baseline value from the ratio data curve.
According to another specific embodiment of the present invention, a method of acquiring reporter images during a PCR process is provided. The method includes acquiring a passive reference image and determining an exposure time for the passive reference image. The method also includes determining a focal plane for the passive reference image and acquiring a plurality of reporter images at the determined focal plane using the determined exposure time.
According to yet another specific embodiment of the present invention, a method of determining a corner point for a PCR ratio data curve is provided. The method includes computing an estimate of a plurality of amplification factors. Each of the plurality of amplification factors is associated with a cycle of the PCR ratio data curve. The method also includes determining a maximum value of the plurality of amplification factors and identifying the cycle of the PCR ratio data curve associated with the maximum value as the corner point.
In an alternative embodiment of the present invention a method of determining a corner point for a PCR ratio data curve is provided. The method includes computing a plurality of derivative values. Each of the derivative values is associated with a cycle of the PCR ratio data curve. The method also includes determining a maximum value of the plurality of derivative values and identifying the cycle of the PCR ratio data curve associated with the maximum value as the corner point.
In yet another alternative embodiment of the present invention, a method of determining a corner point for a PCR ratio data curve is provided. The method includes a) computing an estimate of an amplification factor for a predetermined cycle number; b) determining a fitting window for the PCR ratio data curve, wherein the fitting window is associated with the predetermined cycle number; and c) fitting a first function to the PCR ratio data curve over the fitting window. The method also includes d) determining a fitting error for the first function; e) fitting a second function to the PCR ratio data curve over the fitting window; and f) determining a fitting error for the second function. The method further includes g) computing a difference between the fitting error for the first function and the fitting error for the second function, h) repeating steps a) through g) for subsequent cycle numbers; and i) identifying a maximum of the computed difference as the corner point.
Numerous benefits are achieved using the present invention over conventional techniques. Some embodiments provide optical imaging systems to generate and detect fluorescence from a microfluidic device. Additionally, embodiments of the present invention provide image processing systems that increase the detection sensitivity for fluorescent processes. Depending upon the embodiment, one or more of these benefits may exist. These and other benefits have been described throughout the present specification and more particularly below.
Various additional objects, features and advantages of the present invention can be more fully appreciated with reference to the detailed description and accompanying drawings that follow.
According to the present invention, techniques for microfluidic systems are provided. In particular, the invention provides a method and system for imaging one or more entities suspended in a volume of fluid in a chamber of a microfluidic device. More particularly, the present method and system for imaging uses indications from a fluorescence signal associated with the one or more entities in the microfluidic device. Merely by way of example, the techniques for microfluidic methods and systems are applied using fluorescent, chemiluminescent, and bioluminescent readers coupled to the microfluidic device, but it would be recognized that the invention has a much broader range of applicability.
In some embodiments, a variety of devices and methods for conducting microfluidic analyses are utilized herein, including devices that can be utilized to conduct thermal cycling reactions such as nucleic acid amplification reactions. The devices differ from conventional microfluidic devices in that they include elastomeric components; in some instances, much or all of the device is composed of elastomeric material. For example, amplification reactions can be linear amplifications, (amplifications with a single primer), as well as exponential amplifications (i.e., amplifications conducted with a forward and reverse primer set).
The methods and systems provided by some embodiments of the present invention utilize blind channel type devices in performing nucleic acid amplification reactions. In these devices, the reagents that are typically deposited within the reaction sites are those reagents necessary to perform the desired type of amplification reaction. Usually this means that some or all of the following are deposited: primers, polymerase, nucleotides, metal ions, buffer, and cofactors, for example. The sample introduced into the reaction site in such cases is the nucleic acid template. Alternatively, however, the template can be deposited and the amplification reagents flowed into the reaction sites. As discussed in more detail throughout the present specification, when a matrix device is utilized to conduct an amplification reaction, samples containing nucleic acid template are flowed through the vertical flow channels and the amplification reagents through the horizontal flow channels or vice versa.
PCR is perhaps the best known amplification technique. The devices utilized in embodiments of the present invention are not limited to conducting PCR amplifications. Other types of amplification reactions that can be conducted include, but are not limited to, (i) ligase chain reaction (LCR) (see Wu and Wallace, Genomics 4:560 (1989) and Landegren et al., Science 241:1077 (1988)); (ii) transcription amplification (see Kwoh et al., Proc. Natl. Acad. Sci. USA 86:1173 (1989)); (iii) self-sustained sequence replication (see Guatelli et al., Proc. Nat. Acad. Sci. USA, 87:1874 (1990)); and (iv) nucleic acid based sequence amplification (NASBA) (see, Sooknanan, R. and Malek, L., BioTechnology 13: 563-65 (1995)). Each of the foregoing references are incorporated herein by reference in their entirety for all purposes.
Moreover, certain devices are designed to conduct thermal cycling reactions (e.g., PCR) with devices that include one or more elastomeric valves to regulate solution flow through the device. Thus, methods for conducting amplification reactions with devices of this design are also provided.
Amplicons can be detected and distinguished (whether isolated in a reaction chamber or at any subsequent time) using routine methods for detecting nucleic acids. Amplicons comprising double-stranded DNA can be detected using intercalation dyes such as SYBR™, Pico Green (Molecular Probes, Inc., Eugene, Oreg.), ethidium bromide and the like (see Zhu et al., 1994, Anal. Chem. 66:1941-48) and/or gel electrophoresis. More often, sequence-specific detection methods are used (i.e., amplicons are detected based on their nucleotide sequence). Examples of detection methods include hybridization to arrays of immobilized oligo or polynucleotides, and use of differentially labeled molecular beacons or other “fluorescence resonance energy transfer” (FRET)-based detection systems. FRET-based detection is a preferred method for detection according to some embodiments of the present invention. In FRET-based assays a change in fluorescence from a donor (reporter) and/or acceptor (quencher) fluorophore in a donor/acceptor fluorophore pair is detected. The donor and acceptor fluorophore pair are selected such that the emission spectrum of the donor overlaps the excitation spectrum of the acceptor. Thus, when the pair of fluorophores are brought within sufficiently close proximity to one another, energy transfer from the donor to the acceptor can occur and can be detected. A variety of assays are known including, for example and not limitation, template extension reactions, quantitative RT-PCR, Molecular Beacons, and Invader assays, these are described briefly below.
FRET and template extension reactions utilize a primer labeled with one member of a donor/acceptor pair and a nucleotide labeled with the other member of the donor/acceptor pair. Prior to incorporation of the labeled nucleotide into the primer during an template-dependent extension reaction, the donor and acceptor are spaced far enough apart that energy transfer cannot occur. However, if the labeled nucleotide is incorporated into the primer and the spacing is sufficiently close, then energy transfer occurs and can be detected. These methods are particularly useful in conducting single base pair extension reactions in the detection of single nucleotide polymorphisms and are described in U.S. Pat. No. 5,945,283 and PCT Publication WO 97/22719. The reactions can optionally be thermocycled to increase signal using the temperature control methods and apparatus described throughout the present specification.
A variety of so-called “real time amplification” methods or “real time quantitative PCR” methods can also be used to determine the quantity of a target nucleic acid present in a sample by measuring the amount of amplification product formed during or after the amplification process itself. Fluorogenic nuclease assays are one specific example of a real time quantitation method which can be used successfully with the devices described herein. This method of monitoring the formation of amplification product involves the continuous measurement of PCR product accumulation using a dual-labeled fluorogenic oligonucleotide probe—an approach frequently referred to in the literature as the “TaqMan” method. See, for example, U.S. Pat. No. 5,723,591.
With molecular beacons, a change in conformation of the probe as it hybridizes to a complementary region of the amplified product results in the formation of a detectable signal. The probe itself includes two sections: one section at the 5′ end and the other section at the 3′ end. These sections flank the section of the probe that anneals to the probe binding site and are complementary to one another. One end section is typically attached to a reporter dye and the other end section is usually attached to a quencher dye. In solution, the two end sections can hybridize with each other to form a hairpin loop. In this conformation, the reporter and quencher dye are in sufficiently close proximity that fluorescence from the reporter dye is effectively quenched by the quencher dye. Hybridized probe, in contrast, results in a linearized conformation in which the extent of quenching is decreased. Thus, by monitoring emission changes for the two dyes, it is possible to indirectly monitor the formation of amplification product. Probes of this type and methods of their use are described further, for example, by Piatek et al., 1998, Nat. Biotechnol. 16:359-63; Tyagi, and Kramer, 1996, Nat. Biotechnology 14:303-308; and Tyagi, et al., 1998, Nat. Biotechnol. 16:49-53 (1998).
The Scorpion detection method is described, for example, by Thelwell et al. 2000, Nucleic Acids Research, 28:3752-3761 and Solinas et al., 2001, “Duplex Scorpion primers in SNP analysis and FRET applications” Nucleic Acids Research 29:20. Scorpion primers are fluorogenic PCR primers with a probe element attached at the 5′-end via a PCR stopper. They are used in real-time amplicon-specific detection of PCR products in homogeneous solution. Two different formats are possible, the ‘stem-loop’ format and the ‘duplex’ format. In both cases the probing mechanism is intramolecular. The basic elements of Scorpions in all formats are: (i) a PCR primer; (ii) a PCR stopper to prevent PCR read-through of the probe element; (iii) a specific probe sequence; and (iv) a fluorescence detection system containing at least one fluorophore and quencher. After PCR extension of the Scorpion primer, the resultant amplicon contains a sequence that is complementary to the probe, which is rendered single-stranded during the denaturation stage of each PCR cycle. On cooling, the probe is free to bind to this complementary sequence, producing an increase in fluorescence, as the quencher is no longer in the vicinity of the fluorophore. The PCR stopper prevents undesirable read-through of the probe by Taq DNA polymerase.
Invader assays (Third Wave Technologies, Madison, Wis.) are used particularly for SNP genotyping and utilize an oligonucleotide, designated the signal probe, that is complementary to the target nucleic acid (DNA or RNA) or polymorphism site. A second oligonucleotide, designated the Invader Oligo, contains the same 5′ nucleotide sequence, but the 3′ nucleotide sequence contains a nucleotide polymorphism. The Invader Oligo interferes with the binding of the signal probe to the target nucleic acid such that the 5′ end of the signal probe forms a “flap” at the nucleotide containing the polymorphism. This complex is recognized by a structure specific endonuclease, called the Cleavase enzyme. Cleavase cleaves the 5′ flap of the nucleotides. The released flap binds with a third probe bearing FRET labels, thereby forming another duplex structure recognized by the Cleavase enzyme. This time the Cleavase enzyme cleaves a fluorophore away from a quencher and produces a fluorescent signal. For SNP genotyping, the signal probe will be designed to hybridize with either the reference (wild type) allele or the variant (mutant) allele. Unlike PCR, there is a linear amplification of signal with no amplification of the nucleic acid. Further details sufficient to guide one of ordinary skill in the art are provided by, for example, Neri, B. P., et al., Advances in Nucleic Acid and Protein Analysis 3826:117-125, 2000) and U.S. Pat. No. 6,706,471.
As described above, a variety of multiplex amplification systems can be used in conjunction with the present invention. In one type, several different targets can be detected simultaneously by using multiple differently labeled probes each of which is designed to hybridize only to a particular target. Since each probe has a different label, binding to each target to be detected based on the fluorescence signals. By judicious choice of the different labels that are utilized, analyses can be conducted in which the different labels are excited and/or detected at different wavelengths in a single reaction. See, e.g., Fluorescence Spectroscopy (Pesce et al., Eds.) Marcel Dekker, New York, (1971); White et al., Fluorescence Analysis: A Practical Approach, Marcel Dekker, New York, (1970); Berlman, Handbook of Fluorescence Spectra of Aromatic Molecules, 2nd ed., Academic Press, New York, (1971); Griffiths, Colour and Constitution of Organic Molecules, Academic Press, New York, (1976); Indicators (Bishop, Ed.). Pergamon Press, Oxford, 19723; and Haugland, Handbook of Fluorescent Probes and Research Chemicals, Molecular Probes, Eugene (1992).
As used herein, a “detectable label” has the ordinary meaning in the art and refers to an atom (e.g., radionuclide), molecule (e.g., fluorescein), or complex, that is or can be used to detect (e.g., due to a physical or chemical property), indicate the presence of a molecule or to enable binding of another molecule to which it is covalently bound or otherwise associated. The term “label” also refers to covalently bound or otherwise associated molecules (e.g., a biomolecule such as an enzyme) that act on a substrate to produce a detectable atom, molecule or complex. Detectable labels suitable for use in the present invention include any composition detectable by spectroscopic, photochemical, biochemical, immunochemical, electrical, optical or chemical means. Labels useful in the present invention include biotin for staining with labeled streptavidin conjugate, magnetic beads (e.g., Dynabeads.#.), fluorescent dyes (e.g., fluorescein, Texas red, rhodamine, green fluorescent protein, enhanced green fluorescent protein, lissamine, phycoerythrin, Cy2, Cy3, Cy3.5, Cy5, Cy5.5, Cy7, FluorX [Amersham], SyBR Green I & II [Molecular Probes], and the like), radiolabels (e.g., 3H, 125I, 35S, 14C, or 32P), enzymes (e.g., hydrolases, particularly phosphatases such as alkaline phosphatase, esterases and glycosidases, or oxidoreductases, particularly peroxidases such as horse radish peroxidase, and others commonly used in ELISAs), substrates, cofactors, inhibitors, chemiluminescent groups, chromogenic agents, and colorimetric labels such as colloidal gold or colored glass or plastic (e.g., polystyrene, polypropylene, latex, etc.) beads. Patents teaching the use of such labels include U.S. Pat. Nos. 3,817,837; 3,850,752; 3,939,350; 3,996,345; 4,277,437; 4,275,149; and 4,366,241. Means of detecting such labels are well known to those of skill in the art. Thus, for example, radiolabels and chemiluminescent labels may be detected using photographic film or scintillation counters, fluorescent markers may be detected using a photodetector to detect emitted light (e.g., as in fluorescence-activated cell sorting). Enzymatic labels are typically detected by providing the enzyme with a substrate and detecting the reaction product produced by the action of the enzyme on the substrate, and colorimetric labels are detected by simply visualizing the colored label. Thus, a label is any composition detectable by spectroscopic, photochemical, biochemical, immunochemical, electrical, optical or chemical means. The label may be coupled directly or indirectly to the desired component of the assay according to methods well known in the art. Non-radioactive labels are often attached by indirect means. Generally, a ligand molecule (e.g., biotin) is covalently bound to the molecule. The ligand then binds to an anti-ligand (e.g., streptavidin) molecule which is either inherently detectable or covalently bound to a signal generating system, such as a detectable enzyme, a fluorescent compound, or a chemiluminescent compound. A number of ligands and anti-ligands can be used. Where a ligand has a natural anti-ligand, for example, biotin, thyroxine, and cortisol, it can be used in conjunction with the labeled, naturally occurring anti-ligands. Alternatively, any haptenic or antigenic compound can be used in combination with an antibody. The molecules can also be conjugated directly to signal generating compounds, e.g., by conjugation with an enzyme or fluorophore. Means of detecting labels are well known to those of skill in the art. Thus, for example, where the label is a radioactive label, means for detection include a scintillation counter, photographic film as in autoradiography, or storage phosphor imaging. Where the label is a fluorescent label, it may be detected by exciting the fluorochrome with the appropriate wavelength of light and detecting the resulting fluorescence. The fluorescence may be detected visually, by means of photographic film, by the use of electronic detectors such as charge coupled devices (CCDs) or photomultipliers and the like. Similarly, enzymatic labels may be detected by providing the appropriate substrates for the enzyme and detecting the resulting reaction product. Also, simple calorimetric labels may be detected by observing the color associated with the label. It will be appreciated that when pairs of fluorophores are used in an assay, it is often preferred that the they have distinct emission patterns (wavelengths) so that they can be easily distinguished.
Many diseases linked to genome modifications, either of the host organism or of infectious organisms, are the consequence of a change in a small number of nucleotides, frequently involving a change in a single nucleotide. Such single nucleotide changes are referred to as single nucleotide polymorphisms or simply SNPs, and the site at which the SNP occurs is typically referred to as a polymorphic site. The devices described herein can be utilized to determine the identify of a nucleotide present at such polymorphic sites. As an extension of this capability, the devices can be utilized in genotyping analyses. Genotyping involves the determination of whether a diploid organism (i.e., an organism with two copies of each gene) contains two copies of a reference allele (a reference-type homozygote), one copy each of the reference and a variant allele (i.e., a heterozygote), or contains two copies of the variant allele (i.e., a variant-type homozygote). When conducting a genotyping analysis, the methods of the invention can be utilized to interrogate a single variant site. However, as described further below in the section on multiplexing, the methods can also be used to determine the genotype of an individual in many different DNA loci, either on the same gene, different genes or combinations thereof.
Devices to be utilized for conducting genotyping analyses are designed to utilize reaction sites of appropriate size to ensure from a statistical standpoint that a copy of each of the two alleles for a diploid subject are present in the reaction site at a workable DNA concentrations. Otherwise, an analysis could yield results suggesting that a heterozygote is a homozygote simply because a copy of the second allele is not present at the reaction site. Table 1 below indicates the number of copies of the genome present in a 1 nl reaction volume at various exemplary DNA concentrations that can be utilized with the devices described herein.
As a general matter, due to stochastic proportioning of the sample, the copy number present before an amplification reaction is commenced determines the likely error in the measurement. Genotyping analyses using certain devices are typically conducted with samples having a DNA concentration of approximately 0.10 μg/μL, although the current inventors have run successful TaqMan reactions at concentrations in which there is a single genome per reaction site.
Genotyping analyses can be conducted using a variety of different approaches. In these methods, it is generally sufficient to obtain a “yes” or “no” result, i.e., detection need only be able to answer the question whether a given allele is present. Thus, analyses can be conducted only with the primers or nucleotides necessary to detect the presence of one allele potentially at a polymorphic site. However, more typically, primers and nucleotides to detect the presence of each allele potentially at the polymorphic site are included.
Single Base Pair Extension (SBPE) reactions are one technique specifically developed for conducting genotyping analyses. Although a number of SPBE assays have been developed, the general approach is quite similar. Typically, these assays involve hybridizing a primer that is complementary to a target nucleic acid such that the 3′ end of the primer is immediately 5′ of the variant site or is adjacent thereto. Extension is conducted in the presence of one or more labeled non-extendible nucleotides that are complementary to the nucleotide(s) that occupy the variant site and a polymerase. The non-extendible nucleotide is a nucleotide analog that prevents further extension by the polymerase once incorporated into the primer. If the added non-extendible nucleotide(s) is (are) complementary to the nucleotide at the variant site, then a labeled non-extendible nucleotide is incorporated onto the 3′ end of the primer to generate a labeled extension product. Hence, extended primers provide an indication of which nucleotide is present at the variant site of a target nucleic acid. Such methods and related methods are discussed, for example, in U.S. Pat. Nos. 5,846,710; 6,004,744; 5,888,819; 5,856,092; and 5,710,028; and in WO 92/16657.
Genotyping analyses can also be conducted using quantitative PCR methods. In this case, differentially labeled probes complementary to each of the allelic forms are included as reagents, together with primers, nucleotides and polymerase. However, reactions can be conducted with only a single probe, although this can create ambiguity as to whether lack of signal is due to absence of a particular allele or simply a failed reaction. For the typical biallelic case in which two alleles are possible for a polymorphic site, two differentially labeled probes, each perfectly complementary to one of the alleles are usually included in the reagent mixture, together with amplification primers, nucleotides and polymerase. Sample containing the target DNA is introduced into the reaction site. If the allele to which a probe is complementary is present in the target DNA, then amplification occurs, thereby resulting in a detectable signal as described in the detection above. Based upon which of the differential signal is obtained, the identity of the nucleotide at the polymorphic site can be determined. If both signals are detected, then both alleles are present. Thermocycling during the reaction is performed as described in the temperature control section supra.
Gene expression analysis involves determining the level at which one or more genes is expressed in a particular cell. The determination can be qualitative, but generally is quantitative. In a differential gene expression analysis, the levels of the gene(s) in one cell (e.g., a test cell) are compared to the expression levels of the same genes in another cell (control cell). A wide variety of such comparisons can be made. Examples include, but are not limited to, a comparison between healthy and diseased cells, between cells from an individual treated with one drug and cells from another untreated individual, between cells exposed to a particular toxicant and cells not exposed, and so on. Genes whose expression levels vary between the test and control cells can serve as markers and/or targets for therapy. For example, if a certain group of genes is found to be up-regulated in diseased cells rather than healthy cells, such genes can serve as markers of the disease and can potentially be utilized as the basis for diagnostic tests. These genes could also be targets. A strategy for treating the disease might include procedures that result in a reduction of expression of the up-regulated genes.
The design of the microfluidic devices utilized in embodiments of the present invention is helpful in facilitating a variety of gene expression analyses. Because the devices contain a large number of reaction sites, a large number of genes and/or samples can be tested at the same time. Using the blind flow channel devices, for instance, the expression levels of hundreds or thousands of genes can be determined at the same time. The devices also facilitate differential gene expression analyses. With the matrix design, for example, a sample obtained from a healthy cell can be tested in one flow channel, with a sample from a diseased cell run in an immediately adjacent channel. This feature enhances the ease of detection and the accuracy of the results because the two samples are run on the same device at the same time and under the same conditions.
A variety of matrix or array-based devices are also utilized according to embodiments of the present invention. Certain of these devices include: (i) a first plurality of flow channels formed in an elastomeric substrate, (ii) a second plurality of flow channels formed in the elastomeric substrate that intersect the first plurality of flow channels to define an array of reaction sites, (iii) a plurality of isolation valves disposed within the first and second plurality of flow channels that can be actuated to isolate solution within each of the reaction sites from solution at other reaction sites, and (iv) a plurality of guard channels overlaying one or more of the flow channels and/or one or more of the reaction sites to prevent evaporation of solution therefrom. The foregoing devices can be utilized to conduct a number of different types of reactions, including those involving temperature regulation (e.g., thermocycling of nucleic acid analyses).
Some of the microfluidic devices utilize a design typically referred to herein as “blind channel” or “blind fill” and are characterized in part by having a plurality of blind channels, which are flow channels having a dead end or isolated end such that solution can only enter and exit the blind channel at one end (i.e., there is not a separate inlet and outlet for the blind channel). These devices require only a single valve for each blind channel to isolate a region of the blind channel to form an isolated reaction site. During manufacture of this type of device, one or more reagents for conducting an analysis are deposited at the reaction sites, thereby resulting in a significant reduction in the number of input and outputs. Additionally, the blind channels are connected to an interconnected network of channels such that all the reaction sites can be filled from a single, or limited number, of sample inputs. Because of the reduction in complexity in inputs and outputs and the use of only a single valve to isolate each reaction site, the space available for reaction sites is increased. Thus, the features of these devices means that each device can include a large number of reaction sites (e.g., up to tens of thousands) and can achieve high reaction site densities (e.g., over 1,000-4,000 reaction sites/cm2). Individually and collectively, these features also directly translate into a significant reduction in the size of these devices compared to traditional microfluidic devices.
Other microfluidic devices that are disclosed herein utilize a matrix design. In general, microfluidic devices of this type utilize a plurality of intersecting horizontal and vertical flow channels to define an array of reaction sites at the points of intersection. Thus, devices of this design also have an array or reaction sites; however, there is a larger number of sample inputs and corresponding outputs to accommodate the larger number of samples with this design. A valve system referred to as a switchable flow array architecture enables solution be flowed selectively through just the horizontal or flow channels, thus allowing for switchable isolation of various flow channels in the matrix. Hence, whereas the blind channel devices are designed to conduct a large number of analyses under different conditions with a limited number of samples, the matrix devices are constructed to analyze a large number of samples under a limited number of conditions. Still other devices are hybrids of these two general design types.
The microfluidic devices that are described herein are further characterized in part by utilizing various components such as flow channels, control channels, valves and/or pumps fabricated from elastomeric materials. In some instances, essentially the entire device is made of elastomeric materials. Consequently, such devices differ significantly in form and function from the majority of conventional microfluidic devices that are formed from plastics or silicon-based materials. The number of reaction chambers provided according to embodiments of the present invention
The design of the devices enables them to be utilized in combination with a number of different heating systems. Thus, the devices are useful in conducting diverse analyses that require temperature control. Additionally, those microfluidic devices adapted for use in heating applications can incorporate a further design feature to minimize evaporation of sample from the reaction sites. Devices of this type in general include a number of guard channels and/or reservoirs or chambers formed within the elastomeric device through which water can be flowed to increase the water vapor pressure within the elastomeric material from which the device is formed, thereby reducing evaporation of sample material from the reaction sites.
In another embodiment, a temperature cycling device may be used to control the temperature of the microfluidic devices. Preferably, the microfluidic device would be adapted to make thermal contact with the microfluidic device. Where the microfluidic device is supported by a substrate material, such as a glass slide or the bottom of a carrier plate, such as a plastic carrier, a window may be formed in a region of the carrier or slide such that the microfluidic device, preferably a device having an elastomeric block, may directly contact the heating/cooling block of the temperature cycling device. In a preferred embodiment, the heating/cooling block has grooves therein in communication with a vacuum source for applying a suction force to the microfluidic device, preferably a portion adjacent to where the reactions are taking place. Alternatively, a rigid thermally conductive plate may be bonded to the microfluidic device that then mates with the heating and cooling block for efficient thermal conduction resulting.
The array format of certain of the devices means the devices can achieve high throughput. Collectively, the high throughput and temperature control capabilities make the devices useful for performing large numbers of nucleic acid amplifications (e.g., polymerase chain reaction (PCR)). Such reactions will be discussed at length herein as illustrative of the utility of the devices, especially of their use in any reaction requiring temperature control. However, it should be understood that the devices are not limited to these particular applications. The devices can be utilized in a wide variety of other types of analyses or reactions. Examples include analyses of protein-ligand interactions and interactions between cells and various compounds. Further examples are provided throughout the present specification.
The microfluidic devices disclosed herein are typically constructed at least in part from elastomeric materials and constructed by single and multilayer soft lithography (MSL) techniques and/or sacrificial-layer encapsulation methods (see, e.g., Unger et al. (2000) Science 288:113-116, and PCT Publication WO 01/01025, both of which are incorporated by reference herein in their entirety for all purposes). Utilizing such methods, microfluidic devices can be designed in which solution flow through flow channels of the device is controlled, at least in part, with one or more control channels that are separated from the flow channel by an elastomeric membrane or segment. This membrane or segment can be deflected into or retracted from the flow channel with which a control channel is associated by applying an actuation force to the control channels. By controlling the degree to which the membrane is deflected into or retracted out from the flow channel, solution flow can be slowed or entirely blocked through the flow channel. Using combinations of control and flow channels of this type, one can prepare a variety of different types of valves and pumps for regulating solution flow as described in extensive detail in Unger et al. (2000) Science 288:113-116, and PCT Publication WO/02/43615 and WO 01/01025.
The devices provided herein incorporate such pumps and/or valves to isolate selectively a reaction site at which reagents are allowed to react. Alternatively, devices without pumps and/or valves are utilized that use pressure driven flow or polymerization processes to close appropriate channels and thereby selectively isolate reaction sites. The reaction sites can be located at any of a number of different locations within the device. For example, in some matrix-type devices, the reaction site is located at the intersection of a set of flow channels. In blind channel devices, the reaction site is located at the end of the blind channel.
If the device is to be utilized in temperature control reactions (e.g., thermocycling reactions), then, as described in greater detail infra, the elastomeric device is typically fixed to a support (e.g., a glass slide). The resulting structure can then be placed on a temperature control plate, for example, to control the temperature at the various reaction sites. In the case of thermocycling reactions, the device can be placed on any of a number of thermocycling plates.
Because the devices are made of elastomeric materials that are relatively optically transparent, reactions can be readily monitored using a variety of different detection systems at essentially any location on the microfluidic device. Most typically, however, detection occurs at the reaction site itself (e.g., within a region that includes an intersection of flow channels or at the blind end of a flow channel). The fact that the device is manufactured from substantially transparent materials also means that certain detection systems can be utilized with the current devices that are not usable with traditional silicon-based microfluidic devices. Detection can be achieved using detectors that are incorporated into the device or that are separate from the device but aligned with the region of the device to be detected.
Devices utilizing the matrix design generally have a plurality of vertical and horizontal flow channels that intersect to form an array of junctions. Because a different sample and reagent (or set of reagents) can be introduced into each of the flow channels, a large number of samples can be tested against a relatively large number of reaction conditions in a high throughput format. Thus, for example, if a different sample is introduced into each of M different vertical flow channels and a different reagent (or set of reagents) is introduced into each of N horizontal flow channels, then M×N different reactions can be conducted at the same time. Typically, matrix devices include valves that allow for switchable isolation of the vertical and horizontal flow channels. Said differently, the valves are positioned to allow selective flow just through the vertical flow channels or just through the horizontal flow channels. Because devices of this type allow flexibility with respect to the selection of the type and number of samples, as well as the number and type of reagents, these devices are well-suited for conducting analyses in which one wants to screen a large number of samples against a relatively large number of reaction conditions. The matrix devices can optionally incorporate guard channels to help prevent evaporation of sample and reactants.
Some high-density matrix designs utilize fluid communication vias between layers of the microfluidic device to weave control lines and fluid lines through the device. For example, by having a fluid line in each layer of a two layer elastomeric block, higher density reaction cell arrangements are possible. As will be evident to one of skill in the art, multi-layer devices allow fluid lines to cross over or under each other without being in fluid communication. For example, in a particular design, a reagent fluid channel in a first layer is connected to a reagent fluid channel in a second layer through a via, while the second layer also has sample channels therein, the sample channels and the reagent channels terminating in sample and reagent chambers, respectively. The sample and reagent chambers are in fluid communication with each other through an interface channel that has an interface valve associated therewith to control fluid communication between each of the chambers of a reaction cell. In use, the interface is first closed, then reagent is introduced into the reagent channel from the reagent inlet and sample is introduced into the sample channel through the sample inlet. Containment valves are then closed to isolate each reaction cell from other reaction cells. Once the reaction cells are isolated, the interface valve is opened to cause the sample chamber and the reagent chamber to be in fluid communication with each other so that a desired reaction may take place. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
Accordingly, a particular design for a microfluidic device provides for a microfluidic device adapted to react M number of different samples with N number of different reagents comprising: a plurality of reaction cells, each reaction cell comprising a sample chamber and a reagent chamber, the sample chamber and the reagent chamber being in fluid communication through an interface channel having an interface valve associated therewith for controlling fluid communication between the sample chamber and the reagent chamber; a plurality of sample inlets each in fluid communication with the sample chambers; a plurality of reagent inlets each in fluid communication with the reagent chambers; wherein one of the sample inlets or reagent inlets is in fluid communication with one of the sample chambers or one of the reagent chambers, respectively, through a via. Certain embodiments include having the reaction cells be formed within an elastomeric block formed from a plurality of layers bonded together and the interface valve is a deflectable membrane; having the sample inlets be in communication with the sample chamber through a sample channel and the reagent inlet in fluid communication with the reagent chamber through a reagent channel, a portion of the sample channel and a portion of the reagent channel being oriented about parallel to each other and each having a containment valve associated therewith for controlling fluid communication therethrough; having the valve associated with the sample channel and the valve associated with the reagent channel in operable communication with each other through a common containment control channel; having the containment common control channel located along a line about normal to one of the sample channel or the reagent channel.
The microfluidic devices utilized in embodiments of the present invention may be further integrated into the carrier devices described in co-pending and commonly owned U.S. Patent Application No. 60/557,715 by Unger filed on Mar. 29, 2004, which is incorporated herein for all purposes. The carrier of Unger provides on-board continuous fluid pressure to maintain valve closure away from a source of fluid pressure, e.g., house air pressure. Unger further provides for an automated system for charging and actuating the valves of the present invention as described therein. An another preferred embodiment, the automated system for charging accumulators and actuating valves employs a device having a platen that mates against one or more surfaces of the microfluidic device, wherein the platen has at least two or more ports in fluid communication with a controlled vacuum or pressure source, and may include mechanical portions for manipulating portions of the microfluidic device, for example, but not limited to, check valves.
Another device utilized in embodiments of the present invention provides a carrier used as a substrate for stabilizing an elastomeric block. Preferably the carrier has one or more of the following features; a well or reservoir in fluid communication with the elastomeric block through at least one channel formed in or with the carrier; an accumulator in fluid communication with the elastomeric block through at least one channel formed in or with the carrier; and, a fluid port in fluid communication with the elastomeric block, wherein the fluid port is preferably accessible to an automated source of vacuum or pressure, such as the automated system described above, wherein the automated source further comprises a platen having a port that mates with the fluid port to form an isolated fluid connection between the automated system for applying fluid pressure or vacuum to the elastomeric block. In devices utilized in certain embodiments, the automated source can also make fluid communication with one or more accumulators associated with the carrier for charging and discharging pressure maintained in an accumulator. In certain embodiments, the carrier may further comprise a region located in an area of the carrier that contacts the microfluidic device, wherein the region is made from a material different from another portion of the carrier, the material of the region being selected for improved thermal conduction and distribution properties that are different from the other portion of the carrier. Preferred materials for improved thermal conduction and distribution include, but are not limited to silicon, preferably silicon that is highly polished, such as the type of silicon available in the semiconductor field as a polished wafer or a portion cut from the wafer, e.g., chip.
As described more fully below, embodiments of the present invention utilize a thermal source, for example, but not limited to a PCR thermocycler, which may have been modified from its original manufactured state. Generally the thermal source has a thermally regulated portion that can mate with a portion of the carrier, preferably the thermal conduction and distribution portion of the carrier, for providing thermal control to the elastomeric block through the thermal conduction and distribution portion of the carrier. In a preferred embodiment, thermal contact is improved by applying a source of vacuum to a one or more channels formed within the thermally regulated portion of the thermal source, wherein the channels are formed to contact a surface of the thermal conduction and distribution portion of the carrier to apply suction to and maintain the position of the thermal conduction and distribution portion of the carrier. In a preferred embodiment, the thermal conduction and distribution portion of the carrier is not in physical contact with the remainder of the carrier, but is associated with the remainder of the carrier and the elastomeric block by affixing the thermal conduction and distribution portion to the elastomeric block only and leaving a gap surrounding the edges of the thermal conduction and distribution portion to reduce parasitic thermal effects caused by the carrier. It should be understood that in many aspects of the invention described herein, the preferred elastomeric block could be replaced with any of the known microfluidic devices in the art not described herein, for example devices produced such as the GeneChip® by Affymetrix® of Santa Clara, Calif., USA, or by Caliper of Mountain View, Calif., USA. U.S. patents issued to Soane, Parce, Fodor, Wilding, Ekstrom, Quake, or Unger, describe microfluidic or mesoscale fluidic devices that can be substituted for the elastomeric block of the present invention to take advantage of the thermal advantages and improvements, e.g., suction positioning, reducing parasitic thermal transfer to other regions of the fluidic device, which are described above in the context of using an elastomeric block.
Utilizing systems and methods provided according to embodiments of the present invention, throughput increases are provided over 384 well systems. As an example, throughput increases of a factor of 4, 6, 12, and 24 and greater are provided in some embodiments. These throughput increases are provided while reducing the logistical friction of operations. Moreover the systems and methods of embodiments of the present invention enable multiple assays for multiple samples. For example, in a specific embodiment 96 samples and 96 assays are utilized to provide a total of 9,216 data points. In a particular example, the 96 assays are components of a TaqMan 5′ Nuclease Assay.
Furthermore, embodiments of the present invention provide reduced reaction volumes. In embodiments of the present invention, reaction volumes ranging from 10 picoliters to 100 nanoliters are utilized. In some embodiments, reaction volumes greater than 100 nanoliters are utilized. Merely by way of example, in an embodiment, the methods and systems of the present invention are utilized with reaction volumes of 10 picoliters, 50 picoliters, 100 picoliters, 250 picoliters, 500 picoliters, and 1 nanoliter. In alternative embodiments, reaction volumes of 2 nanoliters, 5 nanoliters, 10 nanoliters, 20 nanoliters, 30 nanoliters, 40 nanoliters, 50 nanoliters, 75 nanoliters, and 100 nanoliters are utilized.
Depending on the geometry of the particular microfluidic device and the size of the microfluidic device and the arrangement of the fluid communication paths and processing site, embodiments of the present invention provide for a range of processing site (or reaction chamber) densities. In some embodiments, the methods and systems of the present invention are utilized with chamber densities ranging from about 100 chambers per cm2 to about 1 million chambers per cm2. Merely by way of example, microfluidic devices with chamber densities of 250, 1,000, 2,500, 10,000, 25,000, 100,000, and 250,000 chambers per cm2 are utilized according to embodiments of the present invention. In some embodiments, chamber densities in excess of 1,000,000 chambers per cm2 are utilized, although this is not required by the present invention.
Operating microfluidic devices with such small reaction volumes reduces reagent usage as well as sample usage. Moreover, some embodiments of the present invention provide methods and systems adapted to perform real-time detection, when used in combination with real-time quantitative PCR. Utilizing these systems and methods, six orders of linear dynamic range are provided for some applications as well as quantitative resolution high enough to allow for the detection of sub-nanoMolar fluorophore concentrations in 10 nanoliter volumes. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
Methods conducted with certain blind channel type devices involve providing a microfluidic device that comprises a flow channel formed within an elastomeric material; and a plurality of blind flow channels in fluid communication with the flow channel, with an end region of each blind flow channel defining a reaction site. At least one reagent is introduced into each of the reaction sites, and then a reaction is detected at one or more of the reaction sites. The method can optionally include heating the at least one reagent within the reaction site. Thus, for example, a method can involve introducing the components for a nucleic acid amplification reaction and then thermocycling the components to form amplified product. As more fully described below, an optical imaging system adapted to characterize reactions occurring in certain microfluidic devices is provided according to embodiments of the present invention.
As illustrated in
Excitation filter wheel 244 is illustrated in
Filter 1: A filter with a center wavelength of 485 nm and a spectral bandwidth of 20 nm;
Filter 2: A filter with a center wavelength of 530 nm and a spectral bandwidth of 20 nm; and
Filter 3: A filter with a center wavelength of 580 nm and a spectral bandwidth of 20 nm.
As will be evident to one of skill in the art, embodiments of the present invention are not limited to these particular spectral filters, but will utilize spectral filters adapted for fluorescence processes for particular samples. Moreover, although the previous discussion related to the use of a filter wheel, this is not required by the present invention. In alternative embodiments, spectral filters are provided in geometries other than a wheel. For example, spectral filters that drop into a filter holder, electro-optic filters, filters placed into the optical path by actuators, and the like are included according to embodiments of the present invention. Moreover, in other embodiments, the optical source is a tunable laser adapted to emit radiation at predetermined wavelengths suitable for excitation of fluorescence. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
As illustrated in
When the excitation shutter is operated in an open position, the optical excitation signal passes through a fiber bundle 248 and is directed so as to impinge on a microfluidic device 205 provided in chip carrier to a seven. Other embodiments of the present invention utilize quartz light guides, liquid light guides, other scrambling systems, and the like to increase illumination homogeneity. As illustrated in
In some embodiments, the illumination of the microfluidic device with light produced by the excitation source is provided over a two-dimensional area of the sample. In these embodiments, a large field of view is provided, which enables the performance of fluorescence applications that involve imaging of time resolved chemical processes and reactions. As an example, fluorescent imaging of protein calorimetry and nucleic acid amplification processes are time resolved processes that benefit from embodiments of the present invention. In some of these processes, simultaneously excitation of the fluorescent samples provided in a number of reaction chambers and simultaneous collection of the fluorescent signals produced by the reactions occurring in the number of reaction chambers is desirable. In other processes, for instance, fluorescence lifetime imaging, a brief excitation pulse is followed by detection (and analysis) of the fluorescent signal as it decays in time from an initial level. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
As an example, nucleic acid amplification processes typically include the target DNA, a thermostable DNA polymerase, two oligonucleotide primers, deoxynucleotide triphosphates (dNTPs), a reaction buffer, and magnesium. Once assembled, the reaction is placed in a thermal cycler, an instrument that subjects the reaction to a series of different temperatures for varying amounts of time. This series of temperature and time adjustments is referred to as one cycle of amplification. Each cycle theoretically doubles the amount of targeted sequence (amplicon) in the reaction. Ten cycles theoretically multiply the amplicon by a factor of about one thousand; 20 cycles, by a factor of more than a million in a matter of hours. In some applications, it is desirable to acquire fluorescent imaging data from a large area (e.g., on the order of several cm2) in a time period ranging from seconds to minutes.
In some embodiments of the present invention, the methods and systems provided by embodiments of the present invention facilitate image capture processes that are performed in a predetermined time period. Merely by way of example, in an embodiment of the present invention a method of imaging microfluidic devices is provided. The method includes capturing an image of a spatial region associated with at least a determined number of chambers of a microfluidic device using an image detection spatial region during a time frame of less than one minute, whereupon the capturing of the image of the spatial region is substantially free from a stitching and/or scanning process.
Embodiments of the present invention provide a variety of time frames for image capture, ranging from 1 millisecond to 1 minute. In some embodiments, time frames for image capture are greater than one minute. Depending on the emission properties associated with the processes performed in the chambers of the microfluidic device, the time frame for image capture will vary. For example, in an embodiment, the time frame is 10 ms, 50 ms, 100 ms, 250 ms, 500 ms, 750 ms, or 1 second. In other embodiments, the time frame is 2 seconds, 5 seconds, 10 seconds, 15 seconds, 20 seconds, 30 seconds, 40 seconds, 50 seconds, or 1 minute. Of course, the time frame will depend on the particular applications.
In some embodiments, the image capture process is performed in a synchronous manner, capturing an image of a determined number of chambers simultaneously. As an example, in an exemplary PCR process, the microfluidic device is maintained at a temperature of 90° C. for a time period of 15 seconds. Subsequently, the microfluidic device is maintained at a temperature of 60° C. for 45 seconds. The heating and cooling cycle is repeated at a one minute cycle period for a number of cycles. Utilizing embodiments of the present invention, images of a determined number of chambers present in the microfluidic device are acquired synchronously, while the chambers are maintained at a uniform temperate as a function of position. For example, a two-dimensional image of an entire microfluidic device may be acquired utilizing a 30 second exposure while the microfluidic device is maintained at the temperature of 60° C. One of skill in the art will appreciate the benefits provided by the present invention over raster scanning or stitching systems, in which images of chambers in a first portion (e.g., an upper left quadrant) of the microfluidic device are acquired prior to images of chambers in a second portion (e.g., a lower right quadrant) of the microfluidic device.
In other embodiments, multiple images are acquired of the determined number of chambers during a time frame of less than one minute. As an example of these embodiments, multiple images associated with multiple fluorophores are acquired in a particular embodiment. During the 45 second time period during which the microfluidic device is maintained at the temperature of 60° C., three consecutive images utilizing exposures of 15 seconds may be acquired for three different fluorophores, for example, ROX™, VIC®, and FAM™. Utilizing these multiple images, differential fluorescence ratios can be calculated and analyzed. Of course, depending on the strength of the fluorescent emissions, the exposure times for the various fluorophores may be modified as appropriate the particular application. In this way, embodiments of the present invention provide for imaging of a microfluidic device in multiple spectral bands while the microfluidic device is maintained a constant temperature. The constant temperature, as illustrated by the previous example, may be a portion of a PCR process including cyclical temperature processes.
Embodiments of the present invention provide methods and systems are also adapted to perform and analyze chemiluminescence processes. In some of these processes, reactions occur on a first time scale and an image of the chemiluminescence process is acquired on a second time scale. In a particular process, the second time scale is less than the first time scale. Thus, embodiments of the present invention are adapted to capture synchronous images of chemiluminescence processes when the samples in the reaction chambers of interest have been reacting for an equal amount of time. In some of these processes, temperature control, including temperature cycling of the samples is provided, whereas in other embodiments, the reaction chambers are maintained at a constant temperature.
As illustrated in
Generally, the devices are placed on a thermal cycling plate to thermal cycle the device. A variety of such plates are readily available from commercial sources, including for example the ThermoHybaid Px2 (Franklin, Mass.), MJ Research PTC-200 (South San Francisco, Calif.), Eppendorf Part# E5331 (Westbury, N.Y.), Techne Part# 205330 (Princeton, N.J.).
In some embodiments, the microfluidic device is contacted with a thermal control device such that the thermal control device is in thermal communication with the thermal control source so that a temperature of the reaction in at least one of the reaction chamber is changed as a result of a change in temperature of the thermal control source. In different embodiments, the thermal transfer device may comprise a semiconductor, such as silicon, may comprise a reflective material, and/or may comprise a metal.
The thermal control device may be adapted to apply a force to the thermal transfer device to urge the thermal transfer device towards the thermal control source. The force may comprise a mechanical pressure, a magnetic force, an electrostatic force, or a vacuum force in different embodiments. For example, in one embodiment, the force comprises a vacuum force applied towards the thermal transfer device through channels formed in a surface of the thermal control device or the thermal transfer device. A level of vacuum achieved between the surface of the thermal control device and a surface (or a portion of a surface) of the thermal transfer device may be detected. Such detection may be performed with a vacuum level detector located at a position along the channel or channels distal from a location of a source of vacuum. When the vacuum does not exceed a preset level, an alert may be manifested or a realignment protocol may be engaged.
The array device may be contacted with the thermal control device by employment of one or more mechanical or electromechanical positioning devices. Carrying out of the method may be automatically controlled and monitored. For example, such automatic control and monitoring may be performed with an automatic control system in operable communication with a robotic control system for introducing and removing the array device from the thermal control device. The progress of the reactions may also be monitored.
A unit may be provided comprising the thermal control device. A system may be provided comprising the array device and the thermal control device. To ensure the accuracy of thermal cycling steps, in certain devices it is useful to incorporate sensors detecting temperature at various regions of the device. One structure for detecting temperature is a thermocouple. Such a thermocouple could be created as thin film wires patterned on the underlying substrate material, or as wires incorporated directly into the microfabricated elastomer material itself.
Temperature can also be sensed through a change in electrical resistance. For example, change in resistance of a thermistor fabricated on an underlying semiconductor substrate utilizing conventional techniques can be calibrated to a given temperature change. Alternatively, a thermistor could be inserted directly into the microfabricated elastomer material. Still another approach to detection of temperature by resistance is described in Wu et al. in “MEMS Flow Sensors for Nano-fluidic Applications”, Sensors and Actuators A 89 152-158 (2001), which is hereby incorporated by reference in its entirety. This paper describes the use of doped polysilicon structures to both control and sense temperature. For polysilicon and other semiconductor materials, the temperature coefficient of resistance can be precisely controlled by the identity and amount of dopant, thereby optimizing performance of the sensor for a given application.
Thermo-chromatic materials are another type of structure available to detect temperature on regions of an amplification device. Specifically, certain materials dramatically and reproducibly change color as they pass through different temperatures. Such a material could be added to the solution as they pass through different temperatures. Thermo-chromatic materials could be formed on the underlying substrate or incorporated within the elastomer material. Alternatively, thermo-chromatic materials could be added to the sample solution in the form of particles.
Another approach to detecting temperature is through the use of an infrared camera. An infrared camera in conjunction with a microscope could be utilized to determine the temperature profile of the entire amplification structure. Permeability of the elastomer material to radiation of appropriate wavelengths (e.g. thermal, infrared, and the like) would facilitate this analysis.
Yet another approach to temperature detection is through the use of pyroelectric sensors. Specifically, some crystalline materials, particularly those materials also exhibiting piezoelectric behavior, exhibit the pyroelectric effect. This effect describes the phenomena by which the polarization of the material's crystal lattice, and hence the voltage across the material, is highly dependent upon temperature. Such materials could be incorporated onto the substrate or elastomer and utilized to detect temperature. Other electrical phenomena, such as capacitance and inductance, can be exploited to detect temperature in accordance with embodiments of the present invention. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
Imaging system 200 operates, in one embodiment, in the following manner. First, microfluidic device 205 is securely placed on carrier 207. Based on a fixed feature of the microfluidic device 205, for example, an edge of the base support of microfluidic device, computer 270 then causes and x,y drive (not shown) to move the carrier 207 to align the microfluidic device in a first x,y position. In some embodiments, one or more fiducial markings are utilized during the alignment and positioning process. In a specific embodiment, a user of the system then registers the precise coordinate of one or more fiducial marks with the imaging system. In other embodiments, this process is performed automatically as the centroids of the fiducials can be calculated precisely by locating a symmetric XY fiducial object and removing any non-symmetric components. In some embodiments, features of the fiducials, such as edges and corners are utilized during alignment processes. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
Under the control of computer 270, either adjustments of the carrier 207 to position it in the focal plane of the optical elements 210 and 212 or adjustments of the optical elements 210 and 212 to position the focal plane of the optical elements 210 and 212 to the carrier 207 are performed. In preferred embodiments, the field of view can embrace an entire microfluidic device, including the number of reaction chambers present on the microfluidic device.
A fluorescent, chemiluminescent, or optical signal emitted by the chemical processes occurring in the reaction chambers of the microfluidic device is collected by a first lens system 210. In some embodiments of the present invention, the first lens system is a multi-element optical train including one or more lenses and one or more apertures. As illustrated in
Referring once again to
In some embodiments, the optical filter device 213 and the emission shutter 215 are located between the first lens system and the second lens system. In some of these embodiments, light rays passing through the optical filter device propagate at small angles with respect to the optic axis. As will be evident to one of skill in the art, spectral filters (e.g., interference filters) placed in regions with small incident ray angle are simpler to design and can potentially provide narrower total spectral bandwidth, through such narrow spectral bandwidth characteristics and/or filter positioning are required by the present invention. As illustrated in
In other embodiments, certain intercalation dyes that have dramatic fluorescent enhancement upon binding to double-stranded DNA, and/or show strong chemical affinity for double-stranded DNA, can be used to detect double-stranded amplified DNA. Examples of suitable dyes include, but are not limited to, SYBR™ and Pico Green (from Molecular Probes, Inc. of Eugene, Oreg.), ethidium bromide, propidium iodide, chromomycin, acridine orange, Hoechst 33258, Toto-1, Yoyo-1, and DAPI (4′,6-diamidino-2-phenylindole hydrochloride). Additional discussion regarding the use of intercalation dyes is provided by Zhu et al., Anal. Chem. 66:1941-1948 (1994), which is incorporated by reference in its entirety.
An second lens system 212 is also illustrated in
A number of different detection strategies can be utilized with the microfluidic devices that are provided herein. Selection of the appropriate system is informed in part on the type of event and/or agent being detected. The detectors can be designed to detect a number of different signal types including, but not limited to, signals from radioisotopes, fluorophores, chromophores, electron dense particles, magnetic particles, spin labels, molecules that emit chemiluminescence, electrochemically active molecules, enzymes, cofactors, enzymes linked to nucleic acid probes and enzyme substrates.
Illustrative detection methodologies suitable for use with the present microfluidic devices include, but are not limited to, light scattering, multichannel fluorescence detection, UV and visible wavelength absorption, luminescence, differential reflectivity, and confocal laser scanning. Additional detection methods that can be used in certain application include scintillation proximity assay techniques, radiochemical detection, fluorescence polarization anisotropy, fluorescence lifetime, fluorescence correlation spectroscopy (FCS), time-resolved energy transfer (TRET), fluorescence resonance energy transfer (FRET) and variations such as bioluminescence resonance energy transfer (BRET). Additional detection options include electrical resistance, resistivity, impedance, and voltage sensing.
In some embodiments, detection occurs at a “detection section”, or “detection region”. These terms and other related terms refer to the portion of the microfluidic device at which detection occurs. In some microfluidic devices, the detection section is generally the reaction chambers present in the microfluidic device. The detection section for matrix-based devices is usually within regions of flow channels that are adjacent an intersection, the intersection itself, or a region that encompasses the intersection and a surrounding region.
The detection section can be in communication with one or more microscopes, diodes, light stimulating devices (e.g., lasers), photomultiplier tubes, processors and combinations of the foregoing, which cooperate to detect a signal associated with a particular event and/or agent. Often the signal being detected is an optical signal that is detected in the detection section by one or more optical detectors. The optical detector can include one or more photodiodes (e.g., avalanche photodiodes), a fiber-optic light guide leading, for example, to a photomultiplier tube or tubes, a microscope, and/or a video camera (e.g., a CCD camera).
Detectors can be microfabricated within the microfluidic device, or can be a separate element. If the detector exists as a separate element and the microfluidic device includes a plurality of detection sections, detection can occur within a single detection section at any given moment. As a specific illustrative example, the microfluidic device can be attached to a translatable stage and scanned under a microscope objective. A signal so acquired is then routed to a processor for signal interpretation and processing. Arrays of photomultiplier tubes can also be utilized. Additionally, optical systems that have the capability of collecting signals from all the different detection sections simultaneously while determining the signal from each section can be utilized.
External detectors are usable because the devices that are provided are completely or largely manufactured of materials that are optically transparent at the wavelength being monitored. This feature enables the devices described herein to utilize a number of optical detection systems that are not possible with conventional silicon-based microfluidic devices.
A particular embodiment of the present invention utilizes a detector in the form of a CCD camera and an optical path that provides for a large field of view and a high numerical aperture to maximize the amount of light collected from each reaction chamber, thereby increasing detection sensitivity. In this embodiment, the CCD is used as an array of photodetectors wherein each pixel or group of pixels corresponds to a reaction chamber rather than being used to produce an image of the array. Thus, the optics may be designed or altered such that image quality is reduced or the image is blurred at the detector in order to increase the useable depth of field of the optical system to collect more light from each reaction chamber. Particularly because the assays contemplated in some embodiments of the present invention include biological assays using fluorescent dyes, which dyes photobleach due to exposure to excitation light hence limiting the total number of signal photons obtainable from a given sample, efficient collection of the limited signal photons can be of importance in instruments such as that discussed. Etendue considerations relate the object and image NA (numerical aperture) and total system magnification for any optical system; since image-side NA can be limited (e.g. by reflection losses at the CCD surface for high-incident-angle rays), in general, arbitrarily high object (sample)-side NA is not achievable simultaneously with arbitrary system magnification. In fact, a larger system magnification can allow a higher object-side NA without requiring a simultaneous (and potentially deleterious for reasons described above) rise in image-side NA. Consequently, in the system described, a large CCD (e.g., 30.7 mm×30.7 mm) focal-plane array has been used to allow for a 1:1 optical system (i.e., a system magnification of 1). This allows a collection NA of 0.36 simultaneous with an image-side NA of 0.36 onto the CCD, which provides reasonable performance with respect to surface reflection losses.
In some embodiments, larger object-side NAs result in reduced object-side depth-of-focus, and hence larger blurring at the detector (assuming blur due to depth of focus greater than or equal to blur due to lens aberrations and other issues) for a given depth of reaction chamber in the sample, limiting the allowable minimum spacing between reaction chambers at the sample if low crosstalk in signal between chambers is to be achieved. In conjunction with a 1:1 optical system, this object-side NA consideration is in good keeping with the ˜0.5 NA maximum generally desirable NA onto a CCD (or silicon detector) if one is to avoid reflection losses at the surface thereof. The 1:1 imaging lens system is furthermore inherently free of most odd-order aberrations, increasing the advantage of this particular magnification (M=1). The use of a 1:1 optical system with a detector as large or larger than the microfluidic system to be imaged is thus provided by some embodiments of the present invention as a design for the detailed system.
In other embodiments, there may be a cost constraint related to the size of the detector (e.g. a CCD focal-plane array). For example, some current high quantum-efficiency, full-frame CCD arrays have dimensions of 27.6 mm×27.6 mm. This detector is slightly smaller than a microfluidic device with dimensions of 30.7 mm×30.7 mm, resulting in a system magnification of 0.88 as a design for the system described. Being near system magnification M=1, constraints related to the detector (image-side) incident NA described above are satisfied for such a magnification.
In other embodiments, a given XY-plane (perpendicular to the optical axis) spacing and size of the reaction chambers may be specified (e.g. to achieve a desired density of sample-chambers in the XY-plane), while constraints on the minimum total volume of the chambers remain (e.g. to achieve minimum required chemical volumes, for instance to avoid over-large statistical fluctuations due to small numbers of reagent or target molecules, or simply to achieve a required minimum number of fluorescent or otherwise optically-emitting molecules or objects). In such a case, it may be necessary to extend the chambers parallel to the Z (optical)-axis such that the total volume of each chamber remains equal to or greater than some minimum figure. Greater extension along Z (creating high-aspect ratio, or columnar chambers which concentrate the sample to be interrogated along the Z-axis) will generally result in a larger blur of the chamber image at the detector for given object-side NA, due to depth-of-focus considerations, assuming blur due to depth of focus is greater than or equal to blur due to lens aberrations and other issues. In some situations, this will lead to the user of a lower object-side NA. Use of a lower NA lens system allows for greater depth of focus and hence light collection from a chambers extended parallel to the optic axis without generally incurring inordinate crosstalk in the optical signal between adjacent or nearby chambers. In this way, a greater density of chambers in the X-Y plane (the place perpendicular to the optic axis) may be used without inordinate crosstalk, while the total chamber volume may be kept large by extension of the chambers in Z (parallel to the optic axis). In this case, or other cases where a lower object-side NA is acceptable (e.g., cases where a larger XY spacing of reaction chambers allows for more chamber-image blur at the detector without undue crosstalk; in non-light-limited applications, where higher NA is not essential; where there is sufficient sample that photobleaching is not an issue; non-photobleaching samples, circumstances such as lower acceptable system sensitivity), a lower system magnification (M<1) may be suitable, particularly if M≧NA0/0.5, or more preferably M≧NA0/0.36, where NA0=object side NA, or more generally M≧NA0/NAdet where NAdet=maximum NA allowable onto the detector face without overlarge reflection/insertion losses to the detector (NAdet=0.36 to 0.5 for a typical CCD).
In cases where object-side depth-of-focus and/or blur requirements do not necessitate an object-side NA≦0.36, or possibly 0.5, or more generally NA0≦NAdet, a larger detector is desirable since due to Etendue considerations (as discussed above), since a larger M (generally requiring a larger detector for a given sample size) will allow a smaller NAi (image-side NA) for a given NA0. Hence where light-collection requirements (e.g. to achieve a certain assay sensitivity) call for a large NA0 (defined by NA0>NAdet) and depth-of-focus and other design considerations (e.g. cost) allow for a large NA0, a larger M is desirable such that losses are minimized at the detector. In such embodiments it can be useful to use a detector device, for example, one or more CCD devices, having a size of, or larger than, the area of the microfluidic device to be imaged. Use of such a large detector allows an increase in the magnification of the optical system, and hence (via etendue considerations) higher NA light collection from the sample for a fixed incident NA onto the detector (the latter set, e.g., by reflection losses at the CCD surface at high incoming ray incident angles).
A particularly preferred detector uses a CCD camera and an optical path that provides for a large field of view and a high numerical aperture to maximize the amount of light collected from each reaction chamber, thereby increasing detection sensitivity. In this regard, the CCD is used as an array of photodetectors wherein each pixel or group of pixels corresponds to a reaction chamber rather than being used to produce an image of the array. Thus, the optics may be altered such that image quality is reduced or defocused to increase the depth of field of the optical system to collect more light from each reaction chamber. In some embodiments, it is useful to employ high aspect ratio, or columnar chambers, to concentrate the sample to be interrogated by the detector along the optical axis of the optical system, and preferably by defocussing the image to increase the depth of field. Use of a low NA lens system, preferably a bilaterally symmetrical lens system is used. It is also useful to use a detector device, for example, one or more CCD devices, having a size of, or larger than, the area of the microfluidic device to be imaged. Used in conjunction with the low NA optics, improved detection sensitivity can be realized.
A detector system can include a light source for stimulating a reporter that generates a detectable signal. The type of light source utilized depends in part on the nature of the reporter being activated. Suitable light sources include, but are not limited to, lasers, laser diodes, white light sources, and high intensity lamps. If a laser is utilized, the laser can be utilized to scan across a set of detection sections or a single detection section. Laser diodes can be microfabricated into the microfluidic device itself. Alternatively, laser diodes can be fabricated into another device that is placed adjacent to the microfluidic device being utilized to conduct a thermal cycling reaction such that the laser light from the diode is directed into the detection section.
Detection can involve a number of non-optical approaches as well. For example, the detector can also include, for example, a temperature sensor, a conductivity sensor, a potentiometric sensor (e.g., pH electrode) and/or an amperometric sensor (e.g., to monitor oxidation and reduction reactions).
Certain intercalation dyes that that have dramatic fluorescent enhancement upon binding to double-stranded DNA, and/or show strong chemical affinity for double-stranded DNA, can be used to detect double-stranded amplified DNA. Examples of suitable dyes include, but are not limited to, SYBR™ and Pico Green (from Molecular Probes, Inc. of Eugene, Oreg.), ethidium bromide, propidium iodide, chromomycin, acridine orange, Hoechst 33258, Toto-1, Yoyo-1, and DAPI (4′,6-diamidino-2-phenylindole hydrochloride). Additional discussion regarding the use of intercalation dyes is provided by Zhu et al., Anal. Chem. 66:1941-1948 (1994), which is incorporated by reference in its entirety.
As illustrated in
Additionally, embodiments of the present invention provide for imaging systems that are slightly reducing, forming, for example, an image that ranges from about the same size as the object to about half the object size. For example, in an embodiment, an image of a spatial region of a microfluidic device is transmitted and captured, the spatial region being associated with more than 96 chambers. An image detecting device is used to capture the image of the spatial region using an image detection spatial region that is about equal to or slightly less in size than the spatial region of the microfluidic device. Merely by way of example, the ratio of the area of the spatial region of the microfluidic device to the area of the image of the spatial region can be 1:1, 1:0.99, 1:0.95, 1:0.9, 1:0.88, and 1:0.85. These particular ratios are merely exemplary, as the ratio selected for the imaging system will depend on the particular application.
In some embodiments, the optical imaging system includes a field of view of about 3 cm×3 cm. In other embodiments, the optical imaging system includes a field of view that ranges from about 1 cm×1 cm to about 5 cm×5 cm. In particular embodiments, an object field of view of 2 cm×2 cm, 2.5 cm×2.5 cm, 2.76 cm×2.76 cm, 3.07 cm×3.07 cm, 3.5 cm×3.5 cm, and 4 cm×4 cm, is provided. In general, the field of view of the optical imaging system is selected to correspond to the spatial region of the microfluidic device, for example, an area including a number of reaction chambers of interest.
Moreover, embodiments of the present invention provide optical imaging systems with a range of numerical apertures. As an example, an NA ranging from 0.1 to 0.5 is provided according to various embodiments. In a particular embodiment, NAs of 0.15, 0.18, 0.2, 0.23, 0.25, 0.3, 0.36, and 0.4 are provided.
The spatial resolution of the optical imaging system will generally be a function of the size of the pixels in the image detecting device. In some embodiments of the present invention, the magnification (equal to one for some embodiments) and the size of the pixels present in the detector will determine the number of pixels associated with each reaction chamber. Generally, it is preferable to have multiple detector pixels associated with each reaction chamber. For example, if a reaction chamber is 45 μm on a side, up to nine square pixels having a side dimension equal to 15 μm will overlap with the reaction chamber in the 1:1 imaging system. Thus, according to embodiments of the present invention, the number of pixels associated with each reaction chamber ranges from 1 to 100. For example, 4 pixel regions, 9 pixel regions, 16 pixel regions, 25 pixel regions, 36 pixel regions, 49 pixel regions, 64 pixel regions, and 81 pixel regions are associated with each reaction chamber according to some embodiments of the present invention.
In embodiments of the present invention, a range of pixel sizes from 1 μm2 to 900 μm2 are utilized. For example, square pixels 1 μm on a side, 2 μm on a side, 3 μm on a side, 4 μm on a side, 5 μm on a side, 10 μm on a side, 13.5 μm on a side, 15 μm on a side, 20 μm on a side, 25 μm on a side, and 30 μm on a side are utilized in various embodiments of the present invention. As will be evident to one of skill in the art, the pixel size, the detector array dimensions, and the number of pixels per array are related. In alternative embodiments, rectangular pixels with pixel areas ranging from 1 μm2 to 900 μm2 are utilized.
Moreover, detector arrays, also referred to as image detecting devices, including a range of pixel counts are utilized according to various embodiments of the present invention. Array dimensions range from 512×512 pixel regions to 3,000×3,000 pixel regions. Depending on the availability of detector arrays, greater numbers of pixels per array may be provided in some embodiments. In particular embodiments, array dimensions of 1,024×1,024 pixel regions and 2,048 by 2,048 pixel regions are utilized.
Embodiments of the present invention provide an optical imaging system characterized by several system parameters. For example, a working distance of greater than 35 mm, for instance, 45.92 mm is available through embodiments of the present invention. In another embodiment, a Root-Mean-Square (RMS) spot diameter averaging 13.44 μm with a maximum value of 17.85 μm is provided. Moreover, through embodiments of the present invention, an illumination variation of about ±5% is achieved. In some embodiments, the overall length of the optical imaging system is 542.1 mm with a maximum filter AOI of 12.56 degrees, a maximum beam diameter at the filter of 76 mm, a distortion of <0.10%, and a maximum lens diameter of 5.512 inches.
Imaging system (10) operates, in one embodiment, in the following manner. First, microfluidic device (30) is securely placed on stage (20). Based on a fixed feature of the microfluidic device (30), for example, an edge of the base support of microfluidic device (30), computer (110) then causes x,y drive (25) to move stage (20) about to align microfluidic device (30) in a first x,y position with a first of a plurality of fiducial markings, wherein the fiducial markings are embedded within the microfluidic device at a known z dimension distance from a chamber center point, comes into focus by imaging device (60) based on dead reckoning from the fixed feature. A user of the system then registers the precise coordinate of the fiducial with the imaging system. Two or more additional fiducial marks are then likewise mapped with the assistance of a user. In other embodiments, this process is automatic as the centroids of the fiducials can be calculated precisely by locating the symmetric XY fiducial object and removing any non-symmetric components. Imaging device (60), under the control of computer (110) then adjusts the z dimension location of focal plane (100) to focus upon the fiducial marking. For example, once focused upon the first fiducial marking, the imaging system then obtains a first x,y coordinate image of microfluidic device (30) looking for additional fiducial markings within the field of view of imaging device (60). In preferred embodiments, the field of view can embrace an entire metering cell. The computer then analyzes the first x,y coordinate image to determine whether the microfluidic device has skew and stretch, and if skew or stretch are determined, transforms the first x,y image to align the image and coordinate map of the microfluidic device to an idealized coordinate map. The idealized coordinate map is used later during image subtraction and masking steps.
In preferred embodiments, with the microfluidic device x,y coordinate image aligned against the ideal coordinate map, the system then determines whether the stretch, distortion or lack of co-registration between the various microfluidic layers is present in the microfluidic device by comparing the location of the fiducial markings in the x,y coordinate image with the fiducial markings locations in the x,y coordinate image of the ideal stored image map. If differences are present between the actual fiducial locations and the imaged fiducial locations, a matrix transformation, preferably an Affine transformation, is performed to transform the imaged shape of the metering cell into a virtual shape of the ideal metering cell shape. By converting the actual image to a known and fixed ideal image using the matrix transformation computed from the differences between the measured actual fiducial locations and the stored ideal fiducial locations, image subtraction and other image analysis are made possible.
By computing the differences between the coordinate maps through matrix analysis, a matrix transformation may be developed to reform the actual image into an ideal image for use in further image processing described herein. By causing the imaged microfluidic device to conform to a standard shape, image subtraction and masking is possible to maximize the viewable area of a metering cell chamber. Moreover, if defects or debris are present within the chamber at time zero in a series of time based images, such defects or debris can be masked out of subsequent images to avoid false signals when applying automated analysis. In addition to masking off areas of the chambers which contain defects or debris, the walls of the chambers may be subtracted from subsequent images, again so as to not cause false readings in the subsequent analysis. The discrepancy between various layers, such as between the control layer and the channel layer, can also be calculated based on the position of a found object in the control layer, such as the control lines themselves. In another example, this correction is determined based on the control layer fiducials themselves. For certain embodiments, this extra transformation is important since the control layer partitions the protein chamber from the rest of the control line.
Images are captured for multiple thermal cycles (1014). As an example, 40 PCR thermal cycles are performed in some applications and images are captured for three dyes (one passive reference and two reporters). Thus, 120 images are captured during the 40 thermal cycles. In this image capture process, the computed auto-exposure and auto-focus parameters determined in steps 1010 and 1012 are utilized. Image analysis is performed (1016) on the sequence of images to extract optical signal intensities in the PCR reaction chambers as a function of cycle number. In an embodiment utilizing a 48×48 microfluidic device, each of the 2,304 chambers are analyzed for each image associated with the series of 40 cycles for each of the dyes. Additional discussion of a method of extracting optical signal data from chambers in the images as a function of the PCR cycles is provided with reference to
Data analysis is performed (1018) of the optical intensities extracted in step 1016. For gene expression applications, the 40 point DNA amplification curve is analyzed to determine the change in threshold (Ct) value associated with each reaction chamber. Accordingly, for a 48×48 microfluidic device, 2,304 Ct values are provided by the method illustrated in
Chamber positions are then determined (1114) for passive reference images at subsequent cycles (e.g., 2, 3, 4, . . . 40). The chamber positions determined in step 1110 and 1112 are utilized as initial estimates of the chamber positions for this process. After initial chamber positions are determined, fine tuning may be performed as appropriate to the particular applications. After the completion of step 1114, chamber positions for the passive reference images at all thermal cycles are available for subsequent processing.
Chamber positions are determined (1116) for reporter dye images for all PCR cycles. The chamber positions determined for the passive reference images are utilized as initial estimates in determining these chamber positions. In an embodiment, 40 images are utilized for each reporter dye. After initial chamber positions are determined, fine tuning may be performed as appropriate to the particular applications.
One or more statistical parameters are extracted from the intensity data associated with the chambers in the captured images (1118). In an embodiment, an average value of the intensity for each chamber is determined by considering pixel chambers as well as a number of pixels from an area adjacent the particular chamber. As an example, the area analyzed in slightly larger (e.g., by including a one pixel wide boundary around the perimeter of the chamber) than the chamber dimensions so that noise in the extracted intensity value as a result of shift errors in chamber locations is reduced.
The method also includes determining a z-plane with a predetermined contrast score (284). An edge detection operator is used to analyze the images from the stack of images and to determine contrast scores associated with each image. Chambers provided on microfluidic devices analyzed using embodiments of the present invention typically have a chamber thickness on the order of 300 μm. Thus, many of the images in the stack of images will be focused on planes other than a plane containing reaction chambers. For some images in the stack of images, control lines and other features present on the microfluidic device will provide a certain contrast score. Preferably, the highest contrast scores will be associated with the planes including cross-sections of the reaction chambers. In a specific embodiment, several of the images will include cross-sections of the reaction chambers and will be characterized by contrast scores higher than a threshold value. Generally, the threshold value is selected to be greater than a contrast value associated with background features, for example, control lines.
In a particular embodiment, the image that is characterized by the maximum contrast score is selected and referred to as a baseline image. Of course, it is not required that the baseline image be the image with the maximum contrast score selected from all the images in the stack. It will be appreciated that other images will also be suitable for additional processing as described below. In some embodiments, a sample of the images in the image stack are utilized for processing by the edge detection operator. Thus, sampling of the images increases stack scanning speed, with additional processing then performed on images in the vicinity of images with higher contrast scores. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
Although step 284 has been described with respect to an edge detection operator, this particular operator is not required by embodiments of the present invention. Other image processing operators that provide space-frequency analysis are also included within the scope of embodiments of the present invention. Merely by way of example, wavelet operators may also be utilized.
The chamber positions are determined (286) using an algorithm that performs image processing of the baseline image. For example, an algorithm that considers individual pixels or groups of pixels at a predetermined position on the image is used in an embodiment to determine the chamber positions. Generally, data related to the design of the microfluidic device, which has been previously stored in one or more memories, is utilized during the chamber position determination process. The chambers or sample wells may contain reagents and samples during the chamber position determining process or may be empty. According to embodiments of the present invention, the chambers are also referred to as sample wells and reaction chambers. Such references are not intended to limit the scope of embodiments as described throughout the present specification.
In a particular embodiment, several locations with match scores greater than the predetermined match score may be utilized, for example, using an averaging or weighting technique to determine the position of the first chamber. As illustrated in
The locations of additional chambers are predicted (518) and located (520). In an embodiment, prediction of additional chamber location is performed by moving radially outwards from the location of the first chamber. The prediction is based, in an embodiment, on design information previously provided. Finding the locations of additional chambers is performed by using the template matching process described with reference to step 514. Appropriate regions are selected for use in the template matching process. In a specific embodiment, a linear model of the grid pattern present in the microfluidic device is generated and refined as the number of additional chambers located increases. Refined versions of the linear model may be used as part of the prediction model, thereby increasing the accuracy of the prediction model as more additional chambers are located.
If the fit between the predicted locations and the determined locations for a certain number of chambers is less than a predetermined value, then portions of process 286 are repeated as illustrated by branch 524 in
In a particular embodiment, the intensity distribution of the pixels in the baseline image is analyzed to determine the position of a first chamber. Then, the determination of a first chamber position is utilized in some applications, along with design data, to predict chamber position for subsequently determined chambers. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
The process of determining the chamber positions is not limited to the use of the baseline image. Utilizing embodiments of the present invention, multiple images, including, but not limited to the baseline image, are candidates for use with the chamber positioning algorithms.
Because the images in the stack of images are registered with some precision in the x-y plane, once the chamber positions have been determined for the baseline image, the same chamber position locations may be applied to the other images in the stack. It is possible to perform other chamber positioning processes, including processing of images other than the baseline image, as appropriate to a particular application.
The contrast scores are recomputed for the images in the stack, now considering the chamber pixels (288). In an embodiment, only the chamber pixels are utilized in recomputing the contrast scores. Generally, the edge detection operator utilized in step 284 is used to recompute the contrast scores. The use of the chamber pixels enables the computation process to effectively remove the portion of the contrast score resulting from features such as control lines and the like. Accordingly, the contrast score recomputed in step 288 is largely determined by the presence of chamber features. The image acquired in the image plane characterized by the maximum recomputed contrast score is labeled as the updated baseline image.
Generally, the z-position of the updated baseline image will be in the vicinity of the baseline image. In a specific embodiment, the updated baseline image will be the same as the baseline image, although this result is not required by embodiments of the present invention. The reduction of background noise provided by the use of chamber pixels typically provides the selection of an updated baseline image more closely aligned with the reaction chambers.
An additional image stack is optionally acquired (290) at a series of planes at z-positions on either side of the updated baseline image. The spacing between the images in the additional image stack is typically less than the spacing between the images acquired in step 282. In an embodiment, 16 images are acquired at a spacing of 20 μm between z-positions. Of course, alternative embodiments acquire a larger or smaller number of images at larger or smaller spacings. One of ordinary skill in the art would recognize many variations, modifications, and alternatives.
The images in the additional image stack are optionally analyzed (292) to compute the contrast scores for these images. In an embodiment, only chamber pixels are utilized during this analysis to reduce the impact of background features on the computed contrast scores. The process of acquisition and computations in optional steps 290 and 292 provides an enhanced resolution updated baseline image located at an enhanced z-plane. The processes performed in steps 290 and 292 may be repeated a number of times (294) depending on the desired z-position resolution for the final baseline image. In a particular embodiment, steps 290 and 292 are performed once, improving the resolution associated with the enhanced resolution updated baseline image.
The z-plane associated with the updated baseline image determined in step 284 or the enhanced z-plane associated with the enhanced resolution updated baseline image (if optional steps 290 and 292 are utilized) is stored in a memory (296). This stored position, associated with the z-position of the reaction chambers is utilized in subsequent processes as described more fully throughout the present specification.
The passive reference image is analyzed to determine an exposure time for acquisition of a passive reference image (312).
Referring to
A new exposure time (t1) is computed (414) utilizing the formula in equation (1):
where J1 is the target intensity. The target intensity is a predetermined value dependent, in part, on the type of dye utilized. The computation in equation (1) includes the assumption that intensity is linearly proportional to exposure time, an assumption that is valid for many applications provided according to embodiments of the present invention. The new exposure time is compared to a set of predetermined threshold values to determine if the new exposure time is outside of a desired range. In an embodiment, the minimum predetermined threshold values is 100 ms and the maximum predetermined threshold value is 2,000 ms. If the new exposure time is outside the range defined by the set of predetermined threshold values, then the new exposure time is determined to be out of range.
The number of saturated pixels is computed (416) and compared to a predetermined threshold. If the number of saturated pixels exceeds the predetermined threshold, then the image acquired in step 310 is considered to be saturated. Merely by way of example, the predetermined threshold is 2% according to an embodiment of the present invention. In other embodiments, the predetermined threshold is larger or smaller depending on the particular application. If the image is saturated or if the new exposure time is out of range (418), then the process returns to step 310, in which a new passive reference image is acquired. It should be noted that when the process returns to step 310, the new exposure time (rather than the original default exposure time) is utilized during the image acquisition process. If the image is not saturated and the new exposure time is in the range defined by the set of predetermined threshold values, then the process continues to step 314 as illustrated in
A new exposure time (t1) is computed (466) utilizing the formula in equation (2):
where J1 is the target intensity. The target intensity is a predetermined value dependent, in part, on the type of dye utilized. The computation in equation (2) includes the assumption that intensity is linearly proportional to exposure time, an assumption that is valid for many applications provided according to embodiments of the present invention. The new exposure time is compared to a set of predetermined threshold values to determine if the new exposure time is outside of a desired range. In an embodiment, the minimum predetermined threshold values is 100 ms and the maximum predetermined threshold value is 2,000 ms. If the new exposure time is outside the range defined by the set of predetermined threshold values, then the new exposure time is determined to be out of range.
The number of saturated pixels is computed (468) and compared to a predetermined threshold. If the number of saturated pixels exceeds the predetermined threshold, then the image acquired in step 310 is considered to be saturated. Merely by way of example, the predetermined threshold is 2% according to an embodiment of the present invention. In other embodiments, the predetermined threshold is larger or smaller depending on the particular application. If the image is saturated or if the new exposure time is out of range (470), then the process returns to step 310, in which a new passive reference image is acquired. It should be noted that when the process returns to step 310, the new exposure time (rather than the original default exposure time) is utilized during the image acquisition process. If the image is not saturated and the new exposure time is in the range defined by the set of predetermined threshold values, then the process continues to step 314 as illustrated in
Returning to
An updated exposure time is optionally determined (316) in some embodiments. For purposes of determining the updated exposure time, the focal plane for the passive reference image is utilized and chamber pixels are considered during the image analysis process. The processes described in steps 460-470 are utilized to determine the updated exposure time, with the process proceeding to step 320 if the image is not saturated and the new exposure time is within the predetermined range. If the image is saturated or the new exposure time is out of range, a new image is acquired by returning the process to step 314 as indicated by process flow 318. In an embodiment, if the results provided by the iterative process flow 318 are not acceptable to the user after a predetermined number of cycles, a manual exposure mode is utilized.
An updated focal plane is optionally determined (320) as illustrated in
The reporter images are acquired (324). These reporter images may include images including various reporter dyes such as FAM and VIC. These reporter images are acquired at either the focal plane determined in step 314 or the updated focal plane determined in step 320. Respective default exposure times are utilized in acquiring the images in step 324. In order to improve the image quality, the method of determining an exposure time for acquisition of a reporter image described in
RatioData[c]=K·Bc, (3)
where K is the original amount of DNA, B is the amplification factor, and c is the PCR cycle number. In an ideal experiment, the value of B=2. In practice, the measured value of the amplification factor B is less than two and represented by an approximate value (A), for example, A=1.9. For each cycle number c, the approximate value A of the amplification factor B is computed as the minimum of equation (4), (5) or (6):
According to embodiments of the present invention, several methods are provided to determine the Corner Point given the RatioData plot.
As an alternative to the method of determining the Corner Point described with respect to
As another alternative to the methods of determining the Corner Point described with respect to
In another embodiment, an alternative method to determine the Corner Point is provided. For purposes of clarity, this alternative method is referred to herein as the slope method. According to this alternative method, for each RatioData curve, the linear portion of the curve in log domain is determined. This is generally performed using a regression (line fitting) analysis. For each cycle number c, the longest possible straight line with a goodness-of-fit greater than 0.95 is fit to the RatioData curve in the log domain (log2). The Corner Point is determined by the cycle number which provides the maximum product of the length of the fitted line and the amplification factor at the given cycle number. In this computation, the amplification factor is the slope of the RatioData curve in the log domain. The amplification region is defined as the span of the fitted line associated with the Corner Point cycle number.
A measure of the quality of the curve fit is provided by QualityG, which is the goodness-of-fit, a value that ranges between zero and one. In other embodiments, other ranges are provided as appropriate to the particular applications. Referring to the slope method described above, the slope of the fitted straight line associated with the Corner Point is assigned the value QualityS. An improved quality measure is provided by finding the minimum of the two values: QualityG and QualityS.
As illustrated in
Examining the group of curves in succession, the thresholds for the various curves are considered. The overall average (or other statistical parameter) for the group of curves is used to determine the auto-threshold for the entire collection of amplification curves. The intersection of the amplification curve with the auto-threshold value is computed using the log graph and the following equation:
Embodiments of the present invention provide methods and systems for baseline removal. In order to provide the final amplification curve, the baseline of the RatioData curve is subtracted from the ratio data.
It is also understood that the examples and embodiments described herein are for illustrative purposes only and that various modifications or changes in light thereof will be suggested to persons skilled in the art and are to be included within the spirit and purview of this application and scope of the appended claims.
The present application claims benefit under 35 U.S.C. §119(e) of U.S. Provisional Patent Application No. 60/844,465, filed Sep. 13, 2006, entitled “Methods and Systems for Image Processing of Microfluidic Devices”, which is incorporated herein by reference in its entirety. The following two regular U.S. patent applications (including this one) are being filed concurrently, and the entire disclosure of the other application is incorporated by reference into this application for all purposes. U.S. patent application Ser. No. 11/900,927; filed Sep. 13, 2007, entitled “Methods and Systems for Image Processing of Microfluidic Devices”; andU.S. patent application Ser. No. 11/900,926; filed Sep. 13, 2007, entitled “Methods and Systems for Determining a Baseline During Image Processing”.
| Number | Name | Date | Kind |
|---|---|---|---|
| 3817837 | Rubenstein et al. | Jun 1974 | A |
| 3850752 | Schuurs et al. | Nov 1974 | A |
| 3939350 | Kronick et al. | Feb 1976 | A |
| 3996345 | Ullman et al. | Dec 1976 | A |
| 4275149 | Litman et al. | Jun 1981 | A |
| 4277437 | Maggio | Jul 1981 | A |
| 4366241 | Tom et al. | Dec 1982 | A |
| 5710028 | Eyal et al. | Jan 1998 | A |
| 5723591 | Livak et al. | Mar 1998 | A |
| 5846710 | Bajaj | Dec 1998 | A |
| 5856092 | Dale et al. | Jan 1999 | A |
| 5888819 | Goelet et al. | Mar 1999 | A |
| 5945283 | Kwok et al. | Aug 1999 | A |
| 6004744 | Goelet et al. | Dec 1999 | A |
| 6706471 | Brow et al. | Mar 2004 | B1 |
| 6716620 | Bashir et al. | Apr 2004 | B2 |
| 6727071 | Dunlay et al. | Apr 2004 | B1 |
| 6902883 | Dunlay et al. | Jun 2005 | B2 |
| 6939452 | Foret et al. | Sep 2005 | B2 |
| 7244396 | Chait et al. | Jul 2007 | B2 |
| 7247274 | Chow | Jul 2007 | B1 |
| 7351376 | Quake et al. | Apr 2008 | B1 |
| 7394943 | Kinney et al. | Jul 2008 | B2 |
| 7403647 | Deutsch et al. | Jul 2008 | B2 |
| 7583853 | Taylor et al. | Sep 2009 | B2 |
| 7616791 | Soderman et al. | Nov 2009 | B2 |
| 7695683 | Quan et al. | Apr 2010 | B2 |
| 7738695 | Shorte et al. | Jun 2010 | B2 |
| 7792345 | Taylor et al. | Sep 2010 | B2 |
| 7883669 | Sun et al. | Feb 2011 | B2 |
| 7911617 | Padmanabhan et al. | Mar 2011 | B2 |
| 20050205005 | Hansen et al. | Sep 2005 | A1 |
| 20050282175 | Taylor et al. | Dec 2005 | A1 |
| 20060281183 | Sun et al. | Dec 2006 | A1 |
| 20080063251 | Deutsch | Mar 2008 | A1 |
| 20080176230 | Owen et al. | Jul 2008 | A1 |
| 20080243523 | Beilis | Oct 2008 | A1 |
| 20090262994 | Haussecker et al. | Oct 2009 | A1 |
| 20100027857 | Wang | Feb 2010 | A1 |
| Number | Date | Country |
|---|---|---|
| WO 9216657 | Oct 1992 | WO |
| WO 9722719 | Jun 1997 | WO |
| WO 0101025 | Jan 2001 | WO |
| WO 0243615 | Jun 2002 | WO |
| Number | Date | Country | |
|---|---|---|---|
| 20080075322 A1 | Mar 2008 | US |
| Number | Date | Country | |
|---|---|---|---|
| 60844465 | Sep 2006 | US |
