Methods and apparatus for measuring analytes using large scale FET arrays
11530444 · 2022-12-20
Assignee
Inventors
- Jonathan M. Rothberg (Guilford, CT)
- Wolfgang Hinz (Kilingworth, CT, US)
- Kim L. Johnson (Carlsbad, CA, US)
- James Bustillo (Castro Valley, CA, US)
- John Leamon (Stonington, CT, US)
- Jonathan Schultz (Guilford, CT, US)
Cpc classification
G01N27/4148
PHYSICS
G01N27/4145
PHYSICS
B01L2300/046
PERFORMING OPERATIONS; TRANSPORTING
C12Q1/6874
CHEMISTRY; METALLURGY
B01L2400/086
PERFORMING OPERATIONS; TRANSPORTING
B01L3/502761
PERFORMING OPERATIONS; TRANSPORTING
International classification
C12Q1/6874
CHEMISTRY; METALLURGY
B01L3/00
PERFORMING OPERATIONS; TRANSPORTING
Abstract
Methods and apparatus relating to very large scale FET arrays for analyte measurements. ChemFET (e.g., ISFET) arrays may be fabricated using conventional CMOS processing techniques based on improved FET pixel and array designs that increase measurement sensitivity and accuracy, and at the same time facilitate significantly small pixel sizes and dense arrays. Improved array control techniques provide for rapid data acquisition from large and dense arrays. Such arrays may be employed to detect a presence and/or concentration changes of various analyte types in a wide variety of chemical and/or biological processes. In one example, chemFET arrays facilitate DNA sequencing techniques based on monitoring changes in the concentration of inorganic pyrophosphate (PPi), hydrogen ions, and nucleotide triphosphates.
Claims
1. A method of using a semiconductor device, the method comprising: flowing a reagent solution uniformly through a flow cell of an apparatus, the apparatus comprising: a substrate including an array of sensors; an array of wells disposed over the substrate, each well capacitively coupled to at least one sensor of the array of sensors, each well of array of wells forming a reaction chamber for a chemical reaction; and a flow cell cover disposed over the array of wells to cover the flow cell and defining a flow chamber between the array of wells and the flow cell cover, the flow chamber comprising a flow cross section including a curved wall providing a uniform fluid flow front over the semiconductor device; and measuring a signal from each sensor of the array of sensors in response to a reaction in the reaction chamber.
2. The method of claim 1, wherein a nucleic acid is disposed within at least one well in the array of wells and the reagent solution includes a nucleotide, wherein measuring the signal includes measuring a byproduct released as a result of nucleotide incorporation during an extension reaction.
3. The method of claim 2, where the array of sensors is an array of ISFET sensors, and where measuring the byproduct released comprises measuring hydrogen ion released as a result of nucleotide incorporation during an extension reaction.
4. The method of claim 2, where the array of sensors is an array of chemFET sensors, and where measuring the byproduct released comprises measuring pyrophosphate released as a result of nucleotide incorporation during an extension reaction.
5. The method of claim 1, wherein measuring the signal from each sensor comprises measuring the signal at a frame rate of at least 20 frames per second.
6. The method of claim 1, wherein measuring a signal from each sensor includes measuring a signal as a result of a change of surface potential of an analyte-sensitive passivation layer of a sensor.
7. The method of claim 6, wherein the analyte-sensitive passivation layer is sensitive to hydrogen ion and measuring a signal from each sensor includes measuring the change of surface potential as a result of the change in hydrogen ion concentration.
8. The method of claim 7, wherein measuring the change in hydrogen ion concentration includes measuring a change in a pH range of about pH 7 to about pH 9.
9. A semiconductor device comprising: an array of chemFET sensors formed on a substrate; wherein an array of wells is disposed over the substrate, each well capacitively coupled to at least one chemFET sensor of the array of sensors; a flow cell housing the array of chemFET sensors, the flow cell comprising a cover disposed over the array of wells and defining a flow chamber between the array of wells and the cover; wherein a flow cross section of the flow chamber comprises a curved wall enabling uniform transit times of fluid across the array of chemFET sensors; and a fluid inlet and a fluid outlet in formed in the cover, the fluid inlet in fluid communication with a fluid source, and the fluid outlet in fluid communication with waste.
10. The semiconductor device of claim 9, wherein the curved wall is a top wall.
11. The semiconductor device of claim 9, wherein the curved wall is a bottom wall.
12. The semiconductor device of claim 9, wherein each chemFET sensor in the array of chemFET sensors includes an analyte-sensitive passivation layer disposed upon a floating gate.
13. The semiconductor device of claim 12, wherein the array of chemFET sensors is an array of ISFET sensors.
14. The semiconductor device of claim 13, wherein the analyte-sensitive passivation layer is sensitive to hydrogen ion.
15. The semiconductor device of claim 12, wherein the analyte-sensitive passivation layer is sensitive to pyrophosphate.
16. The semiconductor device of claim 9, wherein the fluid inlet is in fluid communication with a fluid delivery system and the fluid outlet is in fluid communication with a waste container.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) In the drawings, like reference characters generally refer to the same parts throughout the different views. Also, the drawings are not necessarily to scale, emphasis instead being placed upon generally illustrating the various concepts discussed herein.
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
(25)
(26)
(27)
(28)
(29)
(30)
(31)
(32)
(33)
(34)
(35)
(36)
(37)
(38)
(39)
(40)
(41)
(42)
(43)
(44)
(45)
(46)
(47)
(48)
(49)
(50)
(51)
(52)
(53)
(54)
(55)
(56)
(57)
(58)
(59) FIG. 42F1 is a diagrammatic illustration of an example of a ceiling baffle arrangement for a flow cell in which fluid is introduced at one corner of the chip and exits at a diagonal corner, the baffle arrangement facilitating a desired fluid flow across the array.
(60) FIGS. 42F2-42F8 comprise a set of illustrations of an exemplary flow cell member that may be manufactured by injection molding and may incorporate baffles to facilitate fluid flow, as well as a metalized surface for serving as a reference electrode, including an illustration of said member mounted to a sensor array package over a sensor array, to form a flow chamber thereover.
(61)
(62)
(63)
(64)
(65)
(66)
(67)
(68)
(69)
(70)
(71)
(72)
(73)
(74)
(75)
(76)
DETAILED DESCRIPTION
(77) Following below are more detailed descriptions of various concepts related to, and embodiments of, inventive methods and apparatus relating to large scale chemFET arrays for analyte detection and/or measurement. It should be appreciated that various concepts introduced above and discussed in greater detail below may be implemented in any of numerous ways, as the disclosed concepts are not limited to any particular manner of implementation. Examples of specific implementations and applications are provided primarily for illustrative purposes.
(78) Various inventive embodiments according to the present disclosure are directed at least in part to a semiconductor-based/microfluidic hybrid system that combines the power of microelectronics with the biocompatibility of a microfluidic system. In some examples below, the microelectronics portion of the hybrid system is implemented in CMOS technology for purposes of illustration. It should be appreciated, however, that the disclosure is not intended to be limiting in this respect, as other semiconductor-based technologies may be utilized to implement various aspects of the microelectronics portion of the systems discussed herein.
(79) One embodiment disclosed herein is directed to a large sensor array (e.g., a two-dimensional array) of chemically-sensitive field effect transistors (chemFETs), wherein the individual chemFET sensor elements or “pixels” of the array are configured to detect analyte presence (or absence), analyte levels (or amounts), and/or analyte concentration in an unmanipulated sample, or as a result of chemical and/or biological processes (chemical reactions, cell cultures, neural activity, nucleic acid sequencing processes, etc.) occurring in proximity to the array. Examples of chemFETs contemplated by various embodiments discussed in greater detail below include, but are not limited to, ion-sensitive field effect transistors (ISFETs) and enzyme-sensitive field effect transistors (ENFETs). In one exemplary implementation, one or more microfluidic structures is/are fabricated above the chemFET sensor array to provide for containment and/or confinement of a biological or chemical reaction in which an analyte of interest may be produced or consumed, as the case may be. For example, in one implementation, the microfluidic structure(s) may be configured as one or more “wells” (e.g., small reaction chambers or “reaction wells”) disposed above one or more sensors of the array, such that the one or more sensors over which a given well is disposed detect and measure analyte presence, level, and/or concentration in the given well.
(80) In another exemplary implementation, the invention encompasses a system for high-throughput sequencing comprising at least one two-dimensional array of reaction chambers, wherein each reaction chamber is coupled to a chemically-sensitive field effect transistor (“chemFET”) and each reaction chamber is no greater than 10 μm.sup.3 (i.e., 1 pL) in volume. Preferably, each reaction chamber is no greater than 0.34 pL, and more preferably no greater than 0.096 pL or even 0.012 pL in volume. A reaction chamber can optionally be 2.sup.2, 3.sup.2, 4.sup.2, 5.sup.2, 6.sup.2, 7.sup.2, 8.sup.2, 9.sup.2 or 10.sup.2 square microns in cross-sectional area at the top. Preferably, the array has at least 100, 1,000, 10,000, 100,000, or 1,000,000 reaction chambers. The reaction chambers may be capacitively coupled to the chemFETs, and preferably are capacitively coupled to the chemFETs.
(81) In some embodiments, such a chemFET array/microfluidics hybrid structure may be used to analyze solution(s)/material(s) of interest containing nucleic acids. For example, such structures may be employed to monitor sequencing of nucleic acids. Sequencing of nucleic acids may be performed to determine partial or complete nucleotide sequence of a nucleic acid, to detect the presence and in some instances nature of a single nucleotide polymorphism in a nucleic acid, to determine what therapeutic regimen will be most effective to treat a subject having a particular condition as can be determined by the subject's genetic make-up, to determine and compare nucleic acid expression profiles of two or more states (e.g., comparing expression profiles of diseased and normal tissue, or comparing expression profiles of untreated tissue and tissue treated with drug, enzymes, radiation or chemical treatment), to haplotype a sample (e.g., comparing genes or variations in genes on each of the two alleles present in a human subject), to karyotype a sample (e.g., analyzing chromosomal make-up of a cell or a tissue such as an embryo, to detect gross chromosomal or other genomic abnormalities), and to genotype (e.g., analyzing one or more genetic loci to determine for example carrier status and/or species-genus relationships).
(82) The systems described herein can be utilized to sequence the nucleic acids of an entire genome, or any portion thereof. Genomes that can be sequenced include mammalian genomes, and preferably human genomes. Thus, in one exemplary embodiment, the invention encompasses a method for sequencing a genome or part thereof comprising: delivering fragmented nucleic acids from the genome or part thereof to a system for high-throughput sequencing comprising at least one two-dimensional array of reaction chambers, wherein each reaction chamber is coupled to a chemically-sensitive field effect transistor (“chemFET”) and each reaction chamber is no greater than 1 pL in volume; and detecting a sequencing reaction in at least one of the reaction chambers via a signal from the chemFET of the reaction chamber.
(83) In an alternative exemplary embodiment, the method comprises: delivering fragmented nucleic acids from the genome or part thereof to a sequencing apparatus comprising a two-dimensional array of reaction chambers, wherein each of the reaction chambers is disposed in a sensing relationship with an associated chemFET; and detecting a sequencing reaction in at least one of the reaction chambers via a signal from the associated chemFET.
(84) Preferably, the sequencing reaction is performed by delivering a first deoxyribonucleotide triphosphate (“dNTP”) to each of the reaction chambers. Preferably, the delivering step comprises delivering the first dNTP at substantially the same time to each of the reaction chambers. Once the first dNTP is delivered to the reaction chambers, an enzyme is typically delivered to the reaction chambers to degrade any unused dNTP, followed by washing to remove substantially all of the enzyme from the reaction chambers. Preferably, the enzyme also degrades PPi. The washing can also remove any degraded dNTP or degraded PPi.
(85) The nucleic acids to be sequenced can be naturally or non-naturally occurring nucleic acids, and preferably are naturally occurring nucleic acids. The nucleic acids can be obtained from any of several sources including, but not limited to, deoxyribonucleic acid (“DNA”) (e.g., messenger DNA, complementary DNA, or nuclear DNA), ribonucleic acid (“RNA”) (e.g., micro RNA, transfer RNA, messenger RNA, or small interfering RNA), or peptides. When the source is DNA, the DNA can be obtained from any bodily fluid or tissue that contains DNA, including, but not limited to, blood, saliva, cerebrospinal fluid (“CSF”), skin, hair, urine, stool, and mucus. The starting amounts of nucleic acids to be sequenced determine the minimum sample requirements. Considering the following bead sizes, with an average of 450 bases in the single strand, with an average molecular weight of 325 g/mol per base we have the following:
(86) TABLE-US-00002 Bead Size (um) femto gram of DNA 0.2 0.124 0.3 0.279 0.7 1.52 1.05 3.42 2.8 24.3 5.9 108
(87) Given the number of beads and microwells contemplated for use in an array, it will thus be apparent that a sample taken from a subject to be tested need only be on the order of 3 μg.
(88) Preferably, at least 10.sup.6 base pairs are sequenced per hour, more preferably at least 10.sup.7 base pairs are sequenced per hour, even more preferably at least 10.sup.8 base pairs are sequenced per hour, even more preferably at least 10.sup.9 base pairs are sequenced per hour, and most preferably at least 10.sup.10 base pairs are sequenced per hour using the above-described method. Thus, the method may be used to sequence an entire human genome within about 24 hours, more preferably within about 20 hours, even more preferably within about 15 hours, even more preferably within about 10 hours, even more preferably within about 5 hours, and most preferably within about 1 hour.
(89) The systems described herein can be utilized to sequence an entire genome of an organism from about 3 μg of DNA or less. In one embodiment, the invention encompasses a method for sequencing an entire genome of an organism comprising delivering about 3 μg of DNA or less from the organism to an array of chemically sensitive field effect transistors and determining a sequence of the genome. In another embodiment, the invention encompasses a device adapted for sequencing a complete human genome from about 1 ng of DNA or less without the use of optics or labels. Typically, the device comprises an array of chemFETs and/or an array of microfluidic reaction chambers and/or a semiconductor material coupled to a dielectric material.
(90) The above-described sequencing can be performed in a reaction mixture having an ionic strength of up to 500 μM, 400 μM, or 300 μM. The ionic strength, I, of a solution is a function of the concentration of all ions present in a solution.
(91)
(92) where c.sub.B is the molar concentration of ion B (mol dm.sup.−3), z.sub.B is the charge number of that ion, and the sum is taken over all ions in the solution. The Mg.sup.+2 or Mn.sup.2+ concentration may be 1000 μM, 500 μM, 200 μM, 100 μM, 50 μM, 10 μM, 5 μM, or even 1 μM.
(93) The above-described method may be automated such that the sequencing reactions are performed via robotics. In addition, the sequencing information obtained via the signal from the chemFET of the reaction chamber may be provided to a personal computer, a personal digital assistant, a cellular phone, a video game system, or a television so that a user can monitor the progress of the sequencing reactions remotely. This process is illustrated, for example, in
(94) The systems described herein can also be used to aid in the identification and treatment of disease. For example, the system can be used for identifying a sequence associated with a particular disease or for identifying a sequence associated with a positive response to a particular active ingredient.
(95) In one embodiment, the invention encompasses a method for identifying a sequence associated with a condition comprising: delivering nucleic acids from a plurality of subjects having the condition to a sequencing apparatus comprising a two-dimensional array of reaction chambers, wherein each of the reaction chambers is capacitively coupled to a chemFET; determining sequences of the nucleic acids from signal from said chemFETs; and identifying a common sequence between the DNA from the plurality of subjects. Preferably, the subject is a mammal, and more preferably a human. Preferably, the condition is cancer, an immunosuppressant condition, a neurological condition, or a viral infection.
(96) In another embodiment, the invention encompasses a method for identifying a sequence associated with a positive response to a particular active agent, comprising: sequencing DNA from a plurality of subjects that have exhibited a positive response and from a plurality of subjects having a negative response to an active agent using one or more sequencing apparatuses, wherein each sequencing apparatus comprises an array of chemFETs; and identifying a common DNA sequence in the plurality of subjects that have exhibited a positive response or from the subjects that have exhibited a negative response that is not present in the other plurality of subjects. Preferably, the subject is a mammal, and more preferably a human.
(97) It should be appreciated, however, that while some illustrative examples of the concepts disclosed herein focus on nucleic acid processing, the invention contemplates a broader application of these concepts and is not intended to be limited to these examples.
(98)
(99) The system 1000 includes a semiconductor/microfluidics hybrid structure 300 comprising an ISFET sensor array 100 and a microfluidics flow cell 200. In one aspect, the flow cell 200 may comprise a number of wells (not shown in
(100) As illustrated in
(101) As will be discussed in greater detail herein, various embodiments of the sequencing methods provided in accordance with the invention use beads having attached to their surface multiple identical copies of a template nucleic acid. Such beads are generally referred to herein as “loaded” with nucleic acid. Preferably, each reaction well comprises only a single bead. The nucleic acid loaded beads, of which there may be tens, hundreds, thousands, or more, first enter the flow cell and then individual beads enter individual wells. The beads may enter the wells passively or otherwise. For example, the beads may enter the wells through gravity without any applied external force. The beads may enter the wells through an applied external force including but not limited to a magnetic force or a centrifugal force. In some embodiments, if an external force is applied, it is applied in a direction that is parallel to the well height/depth rather than transverse to the well height/depth, with the aim being to “capture” as many beads as possible. Preferably, the wells (or well arrays) are not agitated, as for example may occur through an applied external force that is perpendicular to the well height/depth. Moreover, once the wells are so loaded, they are not subjected to any other force that could dislodge the beads from the wells.
(102) Furthermore, as will be discussed herein, in addition to the nucleic acid loaded beads, each well may also comprise a plurality of smaller beads, referred to herein as “packing beads”. The packing beads may be composed of any inert material that does not interact or interfere with analytes, reagents, reaction parameters, and the like, present in the wells. Packing beads may be positioned between the chemFET surface and the nucleic acid loaded bead, in which case they may be introduced into the wells before the nucleic acid loaded beads. Alternatively, they may be positioned all around the nucleic acid loaded beads, in which case they may be added to the wells before, during and/or after the nucleic acid loaded beads. In still other embodiments, the majority of the packing beads may be positioned on top of the nucleic acid loaded beads, in which case they may be added to the wells after the nucleic acid loaded beads.
(103) Packing beads may serve one or more various functions. For example, the packing beads may act to minimize or prevent altogether dislodgement of nucleic acid loaded beads from wells, particularly if some fraction of them are positioned on top of the nucleic acid loaded beads. The packing beads may act to decrease the diffusion rate of one or more component in the solution in the well, including but not limited to PPi and unincorporated nucleotides. In some embodiments, the packing beads are paramagnetic beads. Packing beads are commercially available from sources such as Bangs Laboratories.
(104) Diffusion may also be impacted by including in the reaction chambers viscosity increasing agents. An example of such an agent is a polymer that is not a nucleic acid (i.e., a non-nucleic acid polymer). The polymer may be naturally or non-naturally occurring, and it may be of any nature provided it does not interfere with nucleotide incorporation and detection thereof except for slowing the diffusion of PPi, unincorporated nucleotides, and/or other reaction byproducts or reagents, towards the reaction chamber bottom. An example of a suitable polymer is polyethylene glycol (PEG). Other examples include PEA, dextrans, acrylamides, celluloses (e.g. methyl cellulose), and the like. The polymer may be free in solution or it may be immobilized (covalently or non-covalently) to one or more sides of the reaction chamber, to the capture bead, and/or to any packing beads that may be present. Non-covalent attachment may be accomplished via a biotin-avidin interaction.
(105) In still other sequencing embodiments, instead of employing capture beads, the wells can be coated with one or more nucleic acids, including for example a pair of primer nucleic acids, and the template nucleic acid having adaptor nucleotide sequences complementary to the primer nucleotide sequence may be introduced into the wells. These and other agents useful in immobilizing analytes or template nucleic acids may be provided to the sensor array 100, to individual dies as part of the chip packaging, or to wells immediately before the processing of a sample). Other methods involving solgels may be used to immobilize agents such as nucleic acids near the surface of the ISFET array 100.
(106) As will be discussed in greater detail herein, in some aspects of the sequencing methods contemplated by the invention, the template nucleic acid may be amplified prior to or after placement in the well. Various methods exist to amplify nucleic acids. Thus, in one aspect, once a template nucleic acid is loaded into a well of the flow cell 200, amplification may be performed in the well, the resulting amplified product denatured, and sequencing-by-synthesis or then performed. Amplification methods include but are not limited to bridge amplification, rolling circle amplification, or other strategies using isothermal or non-isothermal amplification techniques.
(107) In sum, the flow cell 200 in the system of
(108) In the system 1000 of
(109) As will be discussed in greater detail herein, various embodiments of the present invention may relate to monitoring/measurement techniques that involve the static and/or dynamic responses of an ISFET. In one embodiment relating to detection of nucleotide incorporation during a nucleic acid synthesis or sequencing reaction, detection/measurement techniques particularly rely on the transient or dynamic response of an ISFET (ion-step response, or “ion pulse” output), as discussed above in connection with
(110) In one exemplary implementation, beyond the step-wise or essentially instantaneous pH changes in the analyte solution contemplated by prior research efforts, detection/measurement techniques relying on the dynamic response of an ISFET according to some embodiments of the present invention are based at least in part on the differential diffusion of various ionic species proximate to the analyte/passivation layer interface of the ISFET(s) (e.g., at the bottom of a reaction well over an ISFET). In particular, Applicants have recognized and appreciated that if a given stimulus constituted by a change in ionic strength proximate to the analyte/passivation layer interface, due to the appropriate diffusion of respective species of interest, occurs at a rate that is significantly faster than the ability of the passivation layer to adjust its surface charge density in response to the stimulus of the concentration change (e.g., faster than a characteristic response time constant r associated with the passivation layer surface), a step-wise or essentially instantaneous change in ionic strength is not necessarily required to observe an ion pulse output from the ISFET. This principle is applicable not only to the example of DNA sequencing, but also to other types of chemical and chemical reaction sensing, as well/
(111) For purposes of the present disclosure, the pulse response of an ISFET to a stimulus constituted by any significant change in ionic strength in the analyte solution that occurs on a sufficiently fast time scale is referred to as an “ion pulse,” wherein the amplitude of the pulse is described by Eq. (17) above, and wherein the width and shape of the pulse (rise and decay of the pulse over time) relate to the diffusion parameters associated with the ion species giving rise to the change in ionic strength and the equilibrium reaction kinetics at the analyte/passivation layer interface (as reflected in one or more characteristic time constants r, as discussed above in connection with Eq. (18)). In various embodiments discussed herein, one or more ion pulses due to concentration changes of ionic species proximate to the analyte/passivation layer surface may be observed in an ISFET output signal if such concentration changes occur at a rate that is faster than the ability of the passivation layer to adjust its surface charge density in response to the stimulus of the concentration change (e.g., ionic strength changes occur in a time t<<τ). Again, from the foregoing, it should be appreciated that step-wise or instantaneous changes in ionic strength (e.g., a high flow rate directed normal to the ISFET passivation layer surface) is not required to observe an ion pulse output from the ISFET.
(112) More specifically, with reference again to Eqs. (10) and (11) (assuming the applicability of the Debye theory in a regime of relatively lower ionic strengths, e.g., <1 mole/liter), consider that an ISFET re-equilibrates to the same surface potential after a change in ionic strength. This suggests that if the Debye screening length changes to λ′, the surface charge density changes by the opposite amount to keep the surface potential constant, i.e.:
(113)
As discussed above in connection with Eq. (18), the surface charge density σ may be described approximately as an exponential function of time with time constant τ. Accordingly, if there is a time-varying Debye screening length k(t) due to a time varying ionic strength I(t) in the analyte solution (see Eqs. (12) and (13)), the surface charge density obeys the differential equation:
(114)
From Eq. (21), it may be observed that the surface charge density is continuously trying to get back to a value that gives the initial surface potential as the Debye screening length changes with time. The surface potential as a function of time may be expressed in terms of the Debye screening length. from Eqs. (10) and (11), as:
(115)
Introducing a new variable g(t) that includes σ(t) and referring to Eq. (12), Eq. (22) may be rewritten in terms of a time varying ionic strength I(t) as:
(116)
from which is obtained:
(117)
Eqs. (23) and (24) show that the generation and overall profile of an ion pulse in an ISFET output signal depends both on the time constant t associated with the surface kinetics of the analyte/passivation layer interface and the ionic strength as a function of time I(t), which is given by the concentrations of respective ionic species in the reaction well as they change over time due to reaction and diffusion. As noted above, from the foregoing, Applicants have recognized and appreciated that as long as I(t) varies at a rate that is faster than the ability of the surface charge density to re-equilibrate (as characterized by the time constant τ), an ion pulse output may be observed, wherein the profile of the onset of the pulse is determined primarily by I(t) and the decay profile of the pulse is determined primarily by the time constant τ. Again, the foregoing analysis applies generally to virtually any ion species present in the analyte solution and, accordingly, the transient or dynamic response of an ISFET may be exploited to monitor a variety of chemical and/or biological activity.
(118) With respect to the specific example of a nucleic acid synthesis or sequencing reaction, in one embodiment an exemplary ISFET of the array 100 may be employed to provide an output signal that includes one or more ion pulses in response to addition of synthesis or sequencing reagents to a corresponding reaction well associated with the ISFET and containing a previously loaded template. With reference to
(119) More specifically, in a method according to one embodiment for detecting nucleotide incorporation, if one or more of the nucleotides introduced into a reaction well are incorporated upon their introduction, PPi is generated.
(120)
(121) In
(122) More specifically, in an exemplary method to detect nucleotide incorporation, if PPi is generated pursuant to one or more incorporation events, at some point incorporation ceases, PPi is no longer generated, and the concentration of unincorporated nucleotide (dNTP) increases, as shown in
(123) More generally, Applicants have recognized and appreciated that for any of the ISFET output signals in
(124) In another aspect of the foregoing method for detecting nucleotide incorporation, if no nucleotides are incorporated upon introduction of a given nucleotide dNTP into a well, no PPi is generated, and the first and only ion stimulus to the ISFET associated with the well is the rise in concentration of dNTP as it is admitted to the well. As a result, only a single ion pulse is generated, again assuming appropriate diffusion parameters for the dNTP. This situation is illustrated by the output signal labeled as N.sub.poly=0 in
(125) Accordingly, the above-described method provides for significant characterization of the nucleic acid synthesis or sequencing reaction as follows: 1) If only a single ion pulse is present in the ISFET output signal, there is no incorporation of the nucleotide introduced to a reaction well (i.e., only dNTP stimulus); 2) If two successive ion pulses are present in the ISFET output signal, there is incorporation of the nucleotide introduced to the reaction well (i.e., PPi stimulus followed by dNTP stimulus); and 3) The time interval between two successive ion pulses is a function of the number of incorporation events.
(126) In various exemplary implementations of methods pursuant to the various concepts discussed immediately above for nucleotide incorporation detection, it should be appreciated that merely determining whether one or two pulses are present in an ISFET output signal provides useful information regarding nucleotide incorporation. Hence, not all methods according to the present invention need also analyze the time interval between multiple pulses, if present. Again, however, this time interval provides additional useful information relating to the number of incorporation events, and indeed some exemplary implementations analyze not only the number of pulses present in an ISFET output signal but the timing between pulses, in some instances taking into consideration the various diffusion parameters involved for the respective species giving rise to one or more ion pulses.
(127) Regarding the method for nucleotide incorporation detection discussed above, and with reference again to Eq. (17) above, Applicants have recognized and appreciated that various parameters of a nucleic acid synthesis or sequencing reaction may be improved and/or optimized to increase the signal-to-noise ratio (SNR) associated with a given ion pulse in an ISFET output signal (e.g., increase the amplitude ΔΨ.sub.0 of the ion pulse).
(128) First, it may be appreciated from Eq. (17) that a higher initial equilibrium surface potential Ψ.sub.1, (e.g., the steady state surface potential in the presence of a wash buffer before introduction of a given nucleotide dNTP to a reaction well), generally results in a higher amplitude ΔΨ.sub.0 for an ion pulse. In one aspect, with reference again to Eq. (16), Applicants have appreciated that a higher Ψ.sub.1 is facilitated by the selection of a wash buffer having a pH that is significantly different than the pH.sub.pzc of the ISFET's passivation layer (i.e., the bulk analyte pH that results in zero charge density at the passivation layer surface). Hence, in various implementations, wash buffers may be appropriately selected based on the type of material employed for the passivation layer.
(129) From Eq. (15), Applicants also have appreciated that a relatively high pH sensitivity (ΔΨ.sub.0/ΔpH) of the ISFET also may contribute to a higher Ψ.sub.1; in this manner, as discussed above in connection with Eq. (15), a passivation layer material with a large intrinsic buffering capacity β.sub.int may be desirable, as the dimensionless sensitivity factor a given in Eq. (15) is increased by a large β.sub.int. It should also be appreciated in this regard, however, that a high sensitivity to pH is not necessarily a requirement in the nucleotide incorporation detection method according to the specific embodiment described above, as the method is predicated primarily on detection of PPi and dNTP concentration changes, and not measurement of hydrogen ion concentration per se. This situation underscores the applicability to the method described above of various ISFET array optimization techniques discussed herein based on a consideration of a more limited pH measurement range and reduced linearity requirements (e.g., in consideration of the MOSFET body effect). Of course, methods according to other embodiments of the invention discussed herein indeed may rely on hydrogen ion concentration and concentration changes as reaction indicators. In any case, as noted above, Applicants have appreciated that a relatively high pH sensitivity of the ISFET, even over a significantly limited pH range, may in some instances facilitate an increased amplitude ΔΨ.sub.0 for an ion pulse and hence an increased SNR of the ISFET output signal.
(130) Additionally, Applicants have appreciated from Eq. (17) that reducing or minimizing the initial double layer capacitance C.sub.dl,1 due to the initial ionic strength in the analyte solution (i.e., arising from the wash buffer, prior to introduction of dNTP), and/or increasing or maximizing the subsequent double layer capacitance C.sub.dl,2 arising from the concentration change stimulus (i.e., PPi generation or dNTP increase), also results in a higher amplitude ΔΨ.sub.0 for an ion pulse. Accordingly, in various implementations, low ionic strength buffers may be appropriately selected to facilitate an increased amplitude ΔΨ.sub.0 for an ion pulse and hence an increased SNR of the ISFET output signal.
(131) With respect to time resolution, as alluded to above, some aspect of pulse shape as well as the time interval between any two pulses is a function of the time and rate at which various ionic species of interest are generated proximate to the ISFET's analyte/passivation layer interface. Examples of significant parameters in these respects include, but are not limited to, the flow rate at which various constituents of the analyte solution (e.g., various reagents) are admitted to the reaction well, respective and/or collective concentration of various species, amount of template material present in the well, kinetics of the nucleic acid sequencing or synthesis reactions that may occur in the well, and various factors that affect diffusion of species in the reaction well (e.g., diffusion coefficients of various species, geometry/dimensions of the reaction well, presence and size of packing beads, additives to the analyte solution that affect density). Parameters generally relating to time resolution of one or more ion pulses in an ISFET output signal are discussed in greater detail below.
(132) As noted above, the ISFET may be employed to measure steady state pH values, since in some embodiments pH change is proportional to the number of nucleotides incorporated into the newly synthesized nucleic acid strand. In other embodiments discussed in greater detail below, the FET sensor array may be particularly configured for sensitivity to other analytes that may provide relevant information about the chemical reactions of interest. An example of such a modification or configuration is the use of analyte-specific receptors to bind the analytes of interest, as discussed in greater detail herein.
(133) Via an array controller 250 (also under operation of the computer 260), the ISFET array may be controlled so as to acquire data (e.g., output signals of respective ISFETs of the array) relating to analyte detection and/or measurements, and collected data may be processed by the computer 260 to yield meaningful information associated with the processing (including sequencing) of the template nucleic acid.
(134) With respect to the ISFET array 100 of the system 1000 shown in
(135) The ISFET array 100 is not limited to any particular size, as one- or two-dimensional arrays, including but not limited to as few as two to 256 pixels (e.g., 16 by 16 pixels in a two-dimensional implementation) or as many as 54 mega-pixels (e.g., 7400 by 7400 pixels in a two-dimensional implementation) or even greater may be fabricated and employed for various chemical/biological analysis purposes pursuant to the concepts disclosed herein. In one embodiment of the exemplary system shown in
(136) More generally, a chemFET array according to various embodiments of the present disclosure may be configured for sensitivity to any one or more of a variety of analytes. In one embodiment, one or more chemFETs of an array may be particularly configured for sensitivity to one or more analytes and/or one or more binding events, and in other embodiments different chemFETs of a given array may be configured for sensitivity to different analytes. For example, in one embodiment, one or more sensors (pixels) of the array may include a first type of chemFET configured to be sensitive to a first analyte, and one or more other sensors of the array may include a second type of chemFET configured to be sensitive to a second analyte different from the first analyte. In one exemplary implementation, both a first and a second analyte may indicate a particular reaction such as for example nucleotide incorporation in a sequencing-by-synthesis method. Of course, it should be appreciated that more than two different types of chemFETs may be employed in any given array to detect and/or measure different types of analytes and/or other reactions. In general, it should be appreciated in any of the embodiments of sensor arrays discussed herein that a given sensor array may be “homogeneous” and include chemFETs of substantially similar or identical types to detect and/or measure a same type of analyte (e.g., pH or other ion concentration), or a sensor array may be “heterogeneous” and include chemFETs of different types to detect and/or measure different analytes. For simplicity of discussion, again the example of an ISFET is discussed below in various embodiments of sensor arrays, but the present disclosure is not limited in this respect, and several other options for analyte sensitivity are discussed in further detail below (e.g., in connection with
(137) The chemFET arrays configured for sensitivity to any one or more of a variety of analytes may be disposed in electronic chips, and each chip may be configured to perform one or more different biological reactions. The electronic chips can be connected to the portions of the above-described system which read the array output by means of pins coded in a manner such that the pins convey information to the system as to characteristics of the array and/or what kind of biological reaction(s) is(are) to be performed on the particular chip.
(138) In one embodiment, the invention encompasses an electronic chip configured for conducting biological reactions thereon, comprising one or more pins for delivering information to a circuitry identifying a characteristic of the chip and/or a type of reaction to be performed on the chip. Such t may include, but are not limited to, a short nucleotide polymorphism detection, short tandem repeat detection, or sequencing.
(139) In another embodiment, the invention encompasses a system adapted to performing more than one biological reaction on a chip comprising: a chip receiving module adapted for receiving the chip; and a receiver for detecting information from the electronic chip, wherein the information determines a biological reaction to be performed on the chip. Typically, the system further comprises one or more reagents to perform the selected biological reaction.
(140) In another embodiment, the invention encompasses an apparatus for sequencing a polymer template comprising: at least one integrated circuit that is configured to relay information about spatial location of a reaction chamber, type of monomer added to the spatial location, time required to complete reaction of a reagent comprising a plurality of the monomers with an elongating polymer.
(141) In exemplary implementations based on 0.35 micrometer CMOS processing techniques (or CMOS processing techniques capable of smaller feature sizes), each pixel of the ISFET array 100 may include an ISFET and accompanying enable/select components, and may occupy an area on a surface of the array of approximately ten micrometers by ten micrometers (i.e., 100 micrometers.sup.2) or less; stated differently, arrays having a pitch (center of pixel-to-center of pixel spacing) on the order of 10 micrometers or less may be realized. An array pitch on the order of 10 micrometers or less using a 0.35 micrometer CMOS processing technique constitutes a significant improvement in terms of size reduction with respect to prior attempts to fabricate ISFET arrays, which resulted in pixel sizes on the order of at least 12 micrometers or greater.
(142) More specifically, in some embodiments discussed further below based on the inventive concepts disclosed herein, an array pitch of approximately nine (9) micrometers allows an ISFET array including over 256,000 pixels (e.g., a 512 by 512 array), together with associated row and column select and bias/readout electronics, to be fabricated on a 7 millimeter by 7 millimeter semiconductor die, and a similar sensor array including over four million pixels (e.g., a 2048 by 2048 array yielding over 4 Mega-pixels) to be fabricated on a 21 millimeter by 21 millimeter die. In other examples, an array pitch of approximately 5 micrometers allows an ISFET array including approximately 1.55 Mega-pixels (e.g., a 1348 by 1152 array) and associated electronics to be fabricated on a 9 millimeter by 9 millimeter die, and an ISFET sensor array including over 14 Mega-pixels and associated electronics on a 22 millimeter by 20 millimeter die. In yet other implementations, using a CMOS fabrication process in which feature sizes of less than 0.35 micrometers are possible (e.g., 0.18 micrometer CMOS processing techniques), ISFET sensor arrays with a pitch significantly below 5 micrometers may be fabricated (e.g., array pitch of 2.6 micrometers or pixel area of less than 8 or 9 micrometers.sup.2), providing for significantly dense ISFET arrays. Of course, it should be appreciated that pixel sizes greater than 10 micrometers (e.g., on the order of approximately 20, 50, 100 micrometers or greater) may be implemented in various embodiments of chemFET arrays according to the present disclosure also.
(143) As will be understood by those of skill in the art, the ability to miniaturize sequencing reactions reduces the time, cost and labor involved in sequencing of large genomes (such as the human genome).
(144) In other aspects of the system shown in
(145) As used herein, an array is planar arrangement of elements such as sensors or wells. The array may be one or two dimensional. A one dimensional array is an array having one column (or row) of elements in the first dimension and a plurality of columns (or rows) in the second dimension. An example of a one dimensional array is a 1×5 array. A two dimensional array is an array having a plurality of columns (or rows) in both the first and the second dimension. The number of columns (or rows) in the first and second dimensions may or may not be the same. An example of a two dimensional array is a 5×10 array.
(146) Having provided a general overview of the role of a chemFET (e.g., ISFET) array 100 in an exemplary system 1000 for measuring one or more analytes, following below are more detailed descriptions of exemplary chemFET arrays according to various inventive embodiments of the present disclosure that may be employed in a variety of applications. Again, for purposes of illustration, chemFET arrays according to the present disclosure are discussed below using the particular example of an ISFET array, but other types of chemFETs may be employed in alternative embodiments. Also, again, for purposes of illustration, chemFET arrays are discussed in the context of nucleic acid sequencing applications, however, the invention is not so limited and rather contemplates a variety of applications for the chemFET arrays described herein.
(147) As noted above, various inventive embodiments disclosed herein specifically improve upon the ISFET array design of Milgrew et al. discussed above in connection with
(148) To this end,
(149) In one aspect of the embodiment shown in
(150) As illustrated in
(151) By employing the diode-connected MOSFET Q6 in the bias/readout circuitry 110.sub.j of
(152) In
(153)
(154) In another aspect of the embodiment shown in
(155) In yet another aspect of the embodiment shown in
(156) By not tying the body connection of each ISFET to its source, the possibility of some non-zero source-to-body voltage V.sub.SB may give rise to the “body effect,” as discussed above in connection with
(157)
(158) In the top view of
(159) With reference now to the cross-sectional view of
(160) In the composite cross-sectional view of
(161) Above the substrate, gate oxide, and polysilicon layers shown in
(162) As indicated above,
(163) Applicants have recognized and appreciated that, at least in some applications, pixel capacitance may be a salient parameter for some type of analyte measurements. Accordingly, in another embodiment related to pixel layout and design, various via and metal layers may be reconfigured so as to at least partially mitigate the potential for parasitic capacitances to arise during pixel operation. For example, in one such
(164) embodiment, pixels are designed such that there is a greater vertical distance between the signal lines 112.sub.1, 114.sub.1, 116.sub.1 and 118.sub.1, and the topmost metal layer 304 constituting the floating gate structure 170.
(165) In the embodiment described immediately above, with reference again to
(166) To this end, in another embodiment some via and metal layers are reconfigured such that the signal lines 112.sub.1, 114.sub.1, 116.sub.1 and 118.sub.1 are implemented in the Metal1 and Metal2 layers, and the Metal3 layer is used only as a jumper between the Metal2 layer component of the floating gate structure 170 and the topmost metal layer 304, thereby ensuring a greater distance between the signal lines and the metal layer 304.
(167) In
(168) With reference now to the cross-sectional view of
(169) More specifically, as in the embodiment of
(170) which the correspondence between the lettered top views of respective layers and the cross-sectional view of
(171) Accordingly, by consolidating the signal lines 112.sub.1, 114.sub.1, 116.sub.1 and 118.sub.1 to the Metal1 and Metal2 layers and thereby increasing the distance between these signal lines and the topmost layer 304 of the floating gate structure 170 in the Metal4 layer, parasitic capacitances in the ISFET may be at least partially mitigated. It should be appreciated that this general concept (e.g., including one or more intervening metal layers between signal lines and topmost layer of the floating gate structure) may be implemented in other fabrication processes involving greater numbers of metal layers. For example, distance between pixel signal lines and the topmost metal layer may be increased by adding additional metal layers (more than four total metal layers) in which only jumpers to the topmost metal layer are formed in the additional metal layers. In particular, a six-metal layer fabrication process may be employed, in which the signal lines are fabricated using the Metal1 and Metal2 layers, the topmost metal layer of the floating gate structure is formed in the Metal6 layer, and jumpers to the topmost metal layer are formed in the Metal3, Metal4 and Metal5 layers, respectively (with associated vias between the metal layers). In another exemplary implementation based on a six-metal-layer fabrication process, the general pixel configuration shown in
(172) In yet another aspect relating to reduced capacitance, a dimension “f” of the topmost metal layer 304 (and thus the ISFET sensitive area 178) may be reduced so as to reduce cross-capacitance between neighboring pixels. As may be observed in
(173) dimension “g” so as to provide for additional space between the top metal layers of neighboring pixels. In some illustrative non-limiting implementations, for pixels having a dimension “e” on the order of 9 micrometers the dimension “f” may be on the order of 6 micrometers (as opposed to 7 micrometers, as discussed above), and for pixels having a dimension “e” on the order of 5 micrometers the dimension “f” may be on the order of 3.5 micrometers.
(174) Thus, the pixel chip layout designs respectively shown in
(175) In one exemplary implementation, the gate oxide 165 for the ISFET may be fabricated to have a thickness on the order of approximately 75 Angstroms, giving rise to a gate oxide capacitance per unit area C.sub.ox of 4.5 fF/μm.sup.2. Additionally, the polysilicon gate 164 may be fabricated with dimensions corresponding to a channel width W of 1.2 μm and a channel length L of from 0.35 to 0.6 μm (i.e., W/L ranging from approximately 2 to 3.5), and the doping of the region 160 may be selected such that the carrier mobility for the p-channel is 190 cm.sup.2/V.Math.s (i.e., 1.9E10 μm.sup.2/V.Math.s). From Eq. (2) above, this results in an ISFET transconductance parameter R on the order of approximately 170 to 300 μA/V.sup.2. In other aspects of this exemplary implementation, the analog supply voltage VDDA is 3.3 Volts, and VB1 and VB2 are biased so as to provide a constant ISFET drain current I.sub.Dj on the order of 5 μA (in some implementations, VB1 and VB2 may be adjusted to provide drain currents from approximately 1 μA to 20 μA). Additionally, the MOSFET Q6 (see bias/readout circuitry 110j in
(176) With respect to the analyte-sensitive passivation layer 1720 shown in
(177) For CMOS processes involving aluminum as the metal (which has a melting point of approximately 650 degrees Celsius), a silicon nitride and/or silicon oxynitride passivation layer generally is formed via plasma-enhanced chemical vapor deposition (PECVD), in which a glow discharge at 250-350 degrees Celsius ionizes the constituent gases that form silicon nitride or silicon oxynitride, creating active species that react at the wafer surface to form a laminate of the respective materials. In one exemplary process, a passivation layer having a thickness on the order of approximately 1.0 to 1.5 μm may be formed by an initial deposition of a thin layer of silicon oxynitride (on the order of 0.2 to 0.4 μm) followed by a slighting thicker deposition of silicon oxynitride (on the order of 0.5 μm) and a final deposition of silicon nitride (on the order of 0.5 μm). Because of the low deposition temperature involved in the PECVD process, the aluminum metallization is not adversely affected.
(178) However, Applicants have recognized and appreciated that while a low-temperature PECVD process provides adequate passivation for conventional CMOS devices, the low-temperature process results in a generally low-density and somewhat porous passivation layer, which in some cases may adversely affect ISFET threshold voltage stability. In particular, during ISFET device operation, a low-density porous passivation layer over time may absorb and become saturated with ions from the solution, which may in turn cause an undesirable time-varying drift in the ISFETs threshold voltage V.sub.TH, making accurate measurements challenging.
(179) In view of the foregoing, in one embodiment a CMOS process that uses tungsten metal instead of aluminum may be employed to fabricate ISFET arrays according to the present disclosure. The high melting temperature of Tungsten (above 3400 degrees Celsius) permits the use of a higher temperature low pressure chemical vapor deposition (LPCVD) process (e.g., approximately 700 to 800 degrees Celsius) for a silicon nitride or silicon oxynitride passivation layer. The LPCVD process typically results in significantly more dense and less porous films for the passivation layer, thereby mitigating the potentially adverse effects of ion absorption from the analyte solution leading to ISFET threshold voltage drift.
(180) In yet another embodiment in which an aluminum-based CMOS process is employed to fabricate ISFET arrays according to the present disclosure, the passivation layer 172 shown in
(181) Examples of materials suitable for the second portion 172B (or other additional portions) of the passivation layer 172 include, but are not limited to, silicon nitride, silicon oxynitride, aluminum oxide (Al.sub.2O.sub.3), tantalum oxide (Ta.sub.3O.sub.5), tin oxide (SnO.sub.2) and silicon dioxide (SiO.sub.2). In one aspect, the second portion 172B (or other additional portions) may be deposited via a variety of relatively low-temperature processes including, but not limited to, RF sputtering, DC magnetron sputtering, thermal or e-beam evaporation, and ion-assisted depositions. In another aspect, a pre-sputtering etch process may be employed, prior to deposition of the second portion 172B, to remove any native oxide residing on the first portion 172A (alternatively, a reducing environment, such as an elevated temperature hydrogen environment, may be employed to remove native oxide residing on the first portion 172A). In yet another aspect, a thickness of the second portion 172B may be on the order of approximately 0.04 μm to 0.06 μm (400 to 600 Angstroms) and a thickness of the first portion may be on the order of 1.0 to 1.5 μm, as discussed above. In some exemplary implementations, the first portion 172A may include multiple layers of silicon oxynitride and silicon nitride having a combined thickness of 1.0 to 1.5 μm, and the second portion 172B may include a single layer of either aluminum oxide or tantalum oxide having a thickness of approximately 400 to 600 Angstroms. Again, it should be appreciated that the foregoing exemplary thicknesses are provided primarily for purposes of illustration, and that the disclosure is not limited in these respects.
(182) It has been found according to the invention that hydrogen ion sensitive passivation layers are also sensitive to other analytes including but not limited to PPi and unincorporated nucleotide triphosphates. As an example, a silicon nitride passivation layer is able to detect changes in the concentration of PPi and nucleotide triphosphates. The ability to measure the concentration change of these analytes using the same chemFET greatly facilitates the ability to sequence a nucleic acid using a single array, thereby simplifying the sequencing method.
(183) Thus it is to be understood that the chemFET arrays described herein may be used to detect and/or measure various analytes and, by doing so, may monitor a variety of reactions and/or interactions. It is also to be understood that the discussion herein relating to hydrogen ion detection (in the form of a pH change) is for the sake of convenience and brevity and that static or dynamic levels/concentrations of other analytes (including other ions) can be substituted for hydrogen in these descriptions. In particular, sufficiently fast concentration changes of any one or more of various ion species present in the analyte may be detected via the transient or dynamic response of a chemFET, as discussed above in connection with
(184) The chemFETs, including ISFETs, described herein are capable of detecting any analyte that is itself capable of inducing a change in electric field when in contact with or otherwise sensed or detected by the chemFET surface. The analyte need not be charged in order to be detected by the sensor. For example, depending on the embodiment, the analyte may be positively charged (i.e., a cation), negatively charged (i.e., an anion), zwitterionic (i.e., capable of having two equal and opposite charges but being neutral overall), and polar yet neutral. This list is not intended as exhaustive as other analyte classes as well as species within each class will be readily contemplated by those of ordinary skill in the art based on the disclosure provided herein.
(185) In the broadest sense of the invention, the passivation layer may or may not be coated and the analyte may or may not interact directly with the passivation layer. As an example, the passivation layer may be comprised of silicon nitride and the analyte may be something other than hydrogen ions. As a specific example, the passivation layer may be comprised of silicon nitride and the analyte may be PPi. In these instances, PPi is detected directly (i.e., in the absence of PPi receptors attached to the passivation layer either directly or indirectly).
(186) If the analyte being detected is hydrogen (or alternatively hydroxide), then it is preferable to use weak buffers so that changes in either ionic species can be detected at the passivation layer. If the analyte being detected is something other than hydrogen (or hydroxide) but there is some possibility of a pH change in the solution during the reaction or detection step, then it is preferable to use a strong buffer so that changes in pH do not interfere with the detection of the analyte. A buffer is an ionic molecule (or a solution comprising an ionic molecule) that resists to varying extents changes in pH. Some buffers are able to neutralize acids or bases added to or generated in a solution, resulting in no effective pH change in the solution. It is to be understood that any buffer is suitable provided it has a pKa in the desired range. For some embodiments, a suitable buffer is one that functions in about the pH range of 6 to 9, and more preferably 6.5 to 8.5. In other embodiments, a suitable buffer is one that functions in about the pH range of 7-10, including 8.5-9.5.
(187) The strength of a buffer is a relative term since it depends on the nature, strength and concentration of the acid or base added to or generated in the solution of interest. A weak buffer is a buffer that allows detection (and therefore is not able to otherwise control) pH changes of about at least +/−0.005, about at least +/−0.01, about at least +/−0.015, about at least +/−0.02, about at least +/−0.03, about at least +/−0.04, about at least +/−0.05, about at least +/−0.10, about at least +/−0.15, about at least +/−0.20, about at least +/−0.25, about at least +/−0.30, about at least +/−0.35, about at least +/−0.45, about at least +/−0.50, or more. In some embodiments, the pH change is on the order of about 0.005 (e.g., per nucleotide incorporation) and is preferably an increase in pH. A strong buffer is a buffer that controls pH changes of about at least +/−0.005, about at least +/−0.01, about at least +/−0.015, about at least +/−0.02, about at least +/−0.03, about at least +/−0.04, about at least +/−0.05, about at least +/−0.10, about at least +/−0.15, about at least +/−0.20, about at least +/−0.25, about at least +/−0.30, about at least +/−0.35, about at least +/−0.45, about at least +/−0.50, or more. Buffer strength can be varied by varying the concentration of the buffer species itself. Thus low concentration buffers can be low strength buffers. Examples include those having less than 1 mM (e.g., 50-100 μM) buffer species. A non-limiting example of a weak buffer suitable for the sequencing reactions described herein wherein pH change is the readout is 0.1 mM Tris or Tricine. Examples of suitable weak buffers are provided in the Examples and are also known in the art. Higher concentration buffers can be stronger buffers. Examples include those having 1-25 mM buffer species. A non-limiting example of a strong buffer suitable for the sequencing reactions described herein wherein PPi and/or nucleotide triphosphates are read directly is 1, 5 or 25 mM (or more) Tris or Tricine. One of ordinary skill in the art will be able to determine the optimal buffer for use in the reactions and detection methods encompassed by the invention.
(188) In some embodiments, the passivation layer and/or the layers and/or molecules coated thereon dictate the analyte specificity of the array readout.
(189) Detection of hydrogen ions (in the form of pH), and other analytes as determined by the invention, can be carried out using a passivation layer made of silicon nitride (Si.sub.3N.sub.4), silicon oxynitride (Si.sub.2N.sub.2O), silicon oxide (SiO.sub.2), aluminum oxide (Al.sub.2O.sub.3), tantalum pentoxide (Ta.sub.2O.sub.5), tin oxide or stannic oxide (SnO.sub.2), and the like.
(190) The passivation layer can also detect other ion species directly including but not limited to calcium, potassium, sodium, iodide, magnesium, chloride, lithium, lead, silver, cadmium, nitrate, phosphate, dihydrogen phosphate, and the like.
(191) In some embodiments, the passivation layer is coated with a receptor for the analyte of interest. Preferably, the receptor binds selectively to the analyte of interest or in some instances to a class of agents to which the analyte belongs. As used herein, a receptor that binds selectively to an analyte is a molecule that binds preferentially to that analyte (i.e., its binding affinity for that analyte is greater than its binding affinity for any other analyte). Its binding affinity for the analyte of interest may be 2-fold, 3-fold, 4-fold, 5-fold, 6-fold, 7-fold, 8-fold, 9-fold, 10-fold, 15-fold, 20-fold, 25-fold, 30-fold, 40-fold, 50-fold, 100-fold or more than its binding affinity for any other analyte. In addition to its relative binding affinity, the receptor must also have an absolute binding affinity that is sufficiently high to efficiently bind the analyte of interest (i.e., it must have a sufficient sensitivity). Receptors having binding affinities in the picomolar to micromolar range are suitable. Preferably such interactions are reversible.
(192) The receptor may be of any nature (e.g., chemical, nucleic acid, peptide, lipid, combinations thereof and the like). In such embodiments, the analyte too may be of any nature provided there exists a receptor that binds to it selectively and in some instances specifically. It is to be understood however that the invention further contemplates detection of analytes in the absence of a receptor. An example of this is the detection of PPi and Pi by the passivation layer in the absence of PPi or Pi receptors.
(193) In one aspect, the invention contemplates receptors that are ionophores. As used herein, an ionophore is a molecule that binds selectively to an ionic species, whether anion or cation. In the context of the invention, the ionophore is the receptor and the ion to which it binds is the analyte. Ionophores of the invention include art-recognized carrier ionophores (i.e., small lipid-soluble molecules that bind to a particular ion) derived from microorganisms. Various ionophores are commercially available from sources such as Calbiochem.
(194) Detection of some ions can be accomplished through the use of the passivation layer itself or through the use of receptors coated onto the passivation layer. For example, potassium can be detected selectively using polysiloxane, valinomycin, or salinomycin; sodium can be detected selectively using monensin, nystatin, or SQI-Pr; calcium can be detected selectively using ionomycin, calcimycine (A23187), or CA 1001 (ETH 1001).
(195) Receptors able to bind more than one ion can also be used in some instances. For example, beauvericin can be used to detect calcium and/or barium ions, nigericin can be used to detect potassium, hydrogen and/or lead ions, and gramicidin can be used to detect hydrogen, sodium and/or potassium ions. One of ordinary skill in the art will recognize that these compounds can be used in applications in which single ion specificity is not required or in which it is unlikely (or impossible) that other ions which the compounds bind will be present or generated. Similarly, receptors that bind multiple species of a particular genus may also be useful in some embodiments including those in which only one species within the genus will be present or in which the method does not require distinction between species.
(196) As another example, receptors for neurotoxins are described in Simonian Electroanalysis 2004, 16: 1896-1906.
(197) In other embodiments, including but not limited to nucleic acid sequencing applications, receptors that bind selectively to PPi can be used. Examples of PPi receptors include those compounds shown in
(198) Receptors may be attached to the passivation layer covalently or non-covalently. Covalent attachment of a receptor to the passivation layer may be direct or indirect (e.g., through a linker).
(199) A bifunctional linker is a compound having at least two reactive groups to which two entities may be bound. In some instances, the reactive groups are located at opposite ends of the linker. In some embodiments, the bifunctional linker is a universal bifunctional linker such as that shown in
(200) The bifunctional linker may be a homo-bifunctional linker or a hetero-bifunctional linker, depending upon the nature of the molecules to be conjugated. Homo-bifunctional linkers have two identical reactive groups. Hetero-bifunctional linkers are have two different reactive groups. Various types of commercially available linkers are reactive with one or more of the following groups: primary amines, secondary amines, sulphydryls, carboxyls, carbonyls and carbohydrates. Examples of amine-specific linkers are bis(sulfosuccinimidyl) suberate, bis[2-(succinimidooxycarbonyloxy)ethyl] sulfone, disuccinimidyl suberate, disuccinimidyl tartarate, dimethyl adipimate.2 HCl, dimethyl pimelimidate.2 HCl, dimethyl suberimidate.2 HCl, and ethylene glycolbis-[succinimidyl-[succinate]]. Linkers reactive with sulfhydryl groups include bismaleimidohexane, 1,4-di-[3′-(2′-pyridyldithio)-propionamido)] butane, 1-[p-azidosalicylamido]-4-[iodoacetamido] butane, and N-[4-(p-azidosalicylamido) butyl]-3′-[2′-pyridyldithio] propionamide. Linkers preferentially reactive with carbohydrates include azidobenzoyl hydrazine. Linkers preferentially reactive with carboxyl groups include 4-[p-azidosalicylamido] butylamine.
(201) Heterobifunctional linkers that react with amines and sulfhydryls include N-succinimidyl-3-[2-pyridyldithio] propionate, succinimidyl [4-iodoacetyl]aminobenzoate, succinimidyl 4-[N-maleimidomethyl] cyclohexane-1-carboxylate, m-maleimidobenzoyl-N-hydroxysuccinimide ester, sulfosuccinimidyl 6-[3-[2-pyridyldithio]propionamido]hexanoate, and sulfosuccinimidyl 4-[N-maleimidomethyl]cyclohexane-1-carboxylate. Heterobifunctional linkers that react with carboxyl and amine groups include 1-ethyl-3-[3-dimethylaminopropyl]-carbodiimide hydrochloride. Heterobifunctional linkers that react with carbohydrates and sulfhydryls include 4-[N-maleimidomethyl]-cyclohexane-1-carboxylhydrazide.2 HCl, 4-(4-N-maleimidophenyl)-butyric acid hydrazide.2 HCl, and 3-[2-pyridyldithio] propionyl hydrazide.
(202) Alternatively, receptors may be non-covalently coated onto the passivation layer. Non-covalent deposition of the receptor onto the passivation layer may involve the use of a polymer matrix. The polymer may be naturally occurring or non-naturally occurring and may be of any type including but not limited to nucleic acid (e.g., DNA, RNA, PNA, LNA, and the like, or mimics, derivatives, or combinations thereof), amino acid (e.g., peptides, proteins (native or denatured), and the like, or mimics, derivatives, or combinations thereof, lipids, polysaccharides, and functionalized block copolymers. The receptor may be adsorbed onto and/or entrapped within the polymer matrix. The nature of the polymer will depend on the nature of the receptor being used and/or analyte being detected.
(203) Alternatively, the receptor may be covalently conjugated or crosslinked to the polymer (e.g., it may be “grafted” onto a functionalized polymer).
(204) An example of a suitable peptide polymer is poly-lysine (e.g., poly-L-lysine). Examples of other polymers include block copolymers that comprise polyethylene glycol (PEG), polyamides, polycarbonates, polyalkylenes, polyalkylene glycols, polyalkylene oxides, polyalkylene terepthalates, polyvinyl alcohols, polyvinyl ethers, polyvinyl esters, polyvinyl halides, polyvinylpyrrolidone, polyglycolides, polysiloxanes, polyurethanes, alkyl cellulose, hydroxyalkyl celluloses, cellulose ethers, cellulose esters, nitrocelluloses, polymers of acrylic and methacrylic esters, methyl cellulose, ethyl cellulose, hydroxypropyl cellulose, hydroxypropyl methyl cellulose, hydroxybutyl methyl cellulose, cellulose acetate, cellulose propionate, cellulose acetate butyrate, cellulose acetate phthalate, carboxylethyl cellulose, cellulose triacetate, cellulose sulphate sodium salt, poly(methyl methacrylate), poly(ethyl methacrylate), poly(butylmethacrylate), poly(isobutyl methacrylate), poly(hexylmethacrylate), poly(isodecyl methacrylate), poly(lauryl methacrylate), poly(phenyl methacrylate), poly(methyl acrylate), poly(isopropyl acrylate), poly(isobutyl acrylate), poly(octadecyl acrylate), polyethylene, polypropylene, poly(ethylene glycol), poly(ethylene oxide), poly(ethylene terephthalate), poly(vinyl alcohols), polyvinyl acetate, polyvinyl chloride, polystyrene, polyhyaluronic acids, casein, gelatin, glutin, polyanhydrides, polyacrylic acid, alginate, chitosan, poly(methyl methacrylates), poly(ethyl methacrylates), poly(butylmethacrylate), poly(isobutyl methacrylate), poly(hexylmethacrylate), poly(isodecyl methacrylate), poly(lauryl methacrylate), poly(phenyl methacrylate), poly(methyl acrylate), poly(isopropyl acrylate), poly(isobutyl acrylate), and poly(octadecyl acrylate), poly(lactide-glycolide), copolyoxalates, polycaprolactones, polyesteramides, polyorthoesters, polyhydroxybutyric acid, polyanhydrides, poly (styrene-b-isobutylene-b-styrene) (SIBS) block copolymer, ethylene vinyl acetate, poly(meth)acrylic acid, polymers of lactic acid and glycolic acid, polyanhydrides, poly(ortho)esters, polyurethanes, poly(butic acid), poly(valeric acid), and poly(lactide-cocaprolactone), and natural polymers such as alginate and other polysaccharides including dextran and cellulose, collagen, albumin and other hydrophilic proteins, zein and other prolamines and hydrophobic proteins, copolymers and mixtures thereof, and chemical derivatives thereof including substitutions and/or additions of chemical groups, for example, alkyl, alkylene, hydroxylations, oxidations, and other modifications routinely made by those skilled in the art.
(205) Another issue that relates to ISFET threshold voltage stability and/or predictability involves trapped charge that may accumulate (especially) on metal layers of CMOS-fabricated devices as a result of various processing activities during or following array fabrication (e.g., back-end-of-line processing such as plasma metal etching, wafer cleaning, dicing, packaging, handling, etc.). In particular, with reference to
(206) One opportunity for trapped charge to accumulate includes plasma etching of the topmost metal layer 304. Applicants have recognized and appreciated that other opportunities for charge to accumulate on one or more conductors of the floating gate structure or other portions of the FETs includes wafer dicing, during which the abrasive process of a dicing saw cutting through a wafer generates static electricity, and/or various post-processing wafer handling/packaging steps, such as die-to-package wire bonding, where in some cases automated machinery that handles/transports wafers may be sources of electrostatic discharge (ESD) to conductors of the floating gate structure. If there is no connection to the silicon substrate (or other semi-conductor substrate) to provide an electrical path to bleed off such charge accumulation, charge may build up to the point of causing undesirable changes or damage to the gate oxide 165 (e.g., charge injection into the oxide, or low-level oxide breakdown to the underlying substrate). Trapped charge in the gate oxide or at the gate oxide-semiconductor interface in turn can cause undesirable and/or unpredictable variations in ISFET operation and performance, such as fluctuations in threshold voltage.
(207) In view of the foregoing, other inventive embodiments of the present disclosure are directed to methods and apparatus for improving ISFET performance by reducing trapped charge or mitigating the antenna effect. In one embodiment, trapped charge may be reduced after a sensor array has been fabricated, while in other embodiments the fabrication process itself may be modified to reduce trapped charge that could be induced by some conventional process steps. In yet other embodiments, both “during fabrication” and “post fabrication” techniques may be employed in combination to reduce trapped charge and thereby improve ISFET performance.
(208) With respect to alterations to the fabrication process itself to reduce trapped charge, in one embodiment the thickness of the gate oxide 165 shown in
(209) In another embodiment, the topmost metal layer 304 of the ISFETs floating gate structure 170 shown in
(210) In yet another embodiment, the metal etch process for the topmost metal layer 304 may be modified to include wet chemistry or ion-beam milling rather than plasma etching. For example, the metal layer 304 could be etched using an aqueous chemistry selective to the underlying dielectric (e.g., see website for Transene relating to aluminum, which is hereby incorporated herein by reference). Another alternative approach employs ion-milling rather than plasma etching for the metal layer 304. Ion-milling is commonly used to etch materials that cannot be readily removed using conventional plasma or wet chemistries. The ion-milling process does not employ an oscillating electric field as does a plasma, so that charge build-up does not occur in the metal layer(s). Yet another metal etch alternative involves optimizing the plasma conditions so as to reduce the etch rate (i.e. less power density).
(211) In yet another embodiment, architecture changes may be made to the metal layer to facilitate complete electrical isolation during definition of the floating gate. In one aspect, designing the metal stack-up so that the large area ISFET floating gate is not connected to anything during its final definition may require a subsequent metal layer serving as a “jumper” to realize the electrical connection to the floating gate of the transistor. This “jumper” connection scheme prevents charge flow from the large floating gate to the transistor. This method may be implemented as follows (M=metal layer): i) M1 contacting Poly gate electrode; ii) M2 contacting M1; iii) M3 defines floating gate and separately connects to M2 with isolated island; iv) M4 jumper, having very small area being etched over the isolated islands and connections to floating gate M3, connects the M3 floating gate to the M1/M2/M3 stack connected to the Poly gate immediately over the transistor active area; and v) M3 to M4 interlayer dielectric is removed only over the floating gate so as to expose the bare M3 floating gate. In the method outlined immediately above, step v) need not be done, as the ISFET architecture according to some embodiments discussed above leaves the M4 passivation in place over the M4 floating gate. In one aspect, removal may nonetheless improve ISFET performance in other ways (i.e. sensitivity). In any case, the final sensitive passivation layer may be a thin sputter-deposited ion-sensitive metal-oxide layer. It should be appreciated that the over-layer jumpered architecture discussed above may be implemented in the standard CMOS fabrication flow to allow any of the first three metal layers to be used as the floating gates (i.e. M1, M2 or M3).
(212) With respect to post-fabrication processes to reduce trapped charge, in one embodiment a “forming gas anneal” may be employed as a post-fabrication process to mitigate potentially adverse effects of trapped charge. In a forming gas anneal, CMOS-fabricated ISFET devices are heated in a hydrogen and nitrogen gas mixture. The hydrogen gas in the mixture diffuses into the gate oxide 165 and neutralizes certain forms of trapped charges. In one aspect, the forming gas anneal need not necessarily remove all gate oxide damage that may result from trapped charges; rather, in some cases, a partial neutralization of some trapped charge is sufficient to significantly improve ISFET performance. In exemplary annealing processes according to the present disclosure, ISFETs may be heated for approximately 30 to 60 minutes at approximately 400 to 425 degrees Celsius in a hydrogen/nitrogen mixture that includes 10% to 15% hydrogen. In one particular implementation, annealing at 425 degrees Celsius at 30 minutes in a hydrogen/nitrogen mixture that includes 10% hydrogen is observed to be particularly effective at improving ISFET performance. For aluminum CMOS processes, the temperature of the anneal should be kept at or below 450 degrees Celsius to avoid damaging the aluminum metallurgy. In another aspect of an annealing process according to the present disclosure, the forming gas anneal is performed after wafers of fabricated ISFET arrays are diced, so as to ensure that damage due to trapped charge induced by the dicing process itself, and/or other pre-dicing processing steps (e.g., plasma etching of metals) may be effectively ameliorated. In yet another aspect, the forming gas anneal may be performed after die-to-package wirebonding to similarly ameliorate damage due to trapped charge. At this point in the assembly process, a diced array chip is typically in a heat and chemical resistant ceramic package, and low-tolerance wirebonding procedures as well as heat-resistant die-to-package adhesives may be employed to withstand the annealing procedure. Thus, in one exemplary embodiment, the invention encompasses a method for manufacturing an array of FETs, each having or coupled to a floating gate having a trapped charge of zero or substantially zero comprising: fabricating a plurality of FETs in a common semiconductor substrate, each of a plurality of which is coupled to a floating gate; applying a forming gas anneal to the semiconductor prior to a dicing step; dicing the semiconductor; and applying a forming gas anneal to the semiconductor after the dicing step. Preferably, the semiconductor substrate comprises at least 100,000 FETs. Preferably, the plurality of FETs are chemFETs. The method may further comprise depositing a passivation layer on the semiconductor, depositing a polymeric, glass, ion-reactively etchable or photodefineable material layer on the passivation layer and etching the polymeric, glass ion-reactively etchable or photodefineable material to form an array of reaction chambers in the glass layer.
(213) In yet other processes for mitigating potentially adverse effects of trapped charge according to embodiments of the present disclosure, a variety of “electrostatic discharge (ESD)-sensitive protocols” may be adopted during any of a variety of wafer post-fabrication handling/packaging steps. For example, in one exemplary process, anti-static dicing tape may be employed to hold wafer substrates in place (e.g., during the dicing process). Also, although high-resistivity (e.g., 10 MΩ) deionized water conventionally is employed in connection with cooling of dicing saws, according to one embodiment of the present disclosure less resistive/more conductive water may be employed for this purpose to facilitate charge conduction via the water; for example, deionized water may be treated with carbon dioxide to lower resistivity and improve conduction of charge arising from the dicing process. Furthermore, conductive and grounded die-ejection tools may be used during various wafer dicing/handling/packaging steps, again to provide effective conduction paths for charge generated during any of these steps, and thereby reduce opportunities for charge to accumulate on one or more conductors of the floating gate structure of respective ISFETs of an array.
(214) In yet another embodiment involving a post-fabrication process to reduce trapped charge, the gate oxide region of an ISFET may be irradiated with UV radiation. With reference again to
(215) To facilitate a UV irradiation process to reduce trapped charge, Applicants have recognized and appreciated that materials other than silicon nitride and silicon oxynitride generally need to be employed in the passivation layer 172 shown in
(216) In another aspect of an embodiment involving UV irradiation, each ISFET of a sensor array must be appropriately biased during a UV irradiation process to facilitate reduction of trapped charge. In particular, high energy photons from the UV irradiation, impinging upon the bulk silicon region 160 in which the ISFET conducting channel is formed, create electron-hole pairs which facilitate neutralization of trapped charge in the gate oxide as current flows through the ISFETs conducting channel. To this end, an array controller, discussed further below in connection with
(217) Utilizing at least one of the above-described techniques for reducing trapped charge, we have been able to fabricate FETs floating gates having a trapped charge of zero or substantially zero. Thus, in some embodiments, an aspect of the invention encompasses a floating gate having a surface area of about 4 μm.sup.2 to about 50 μm.sup.2 having baseline threshold voltage and preferably a trapped charge of zero or substantially zero. Preferably the FETs are chemFETs. The trapped charge should be kept to a level that does not cause appreciable variations from FET to FET across the array, as that would limit the dynamic range of the devices, consistency of measurements, and otherwise adversely affect performance.
(218)
(219) Also, as discussed above, it should be appreciated that arrays according to various embodiments of the present invention may be fabricated according to conventional CMOS fabrications techniques, as well as modified CMOS fabrication techniques (e.g., to facilitate realization of various functional aspects of the chemFET arrays discussed herein, such as additional deposition of passivation materials, process steps to mitigate trapped charge, etc.) and other semiconductor fabrication techniques beyond those conventionally employed in CMOS fabrication. Additionally, various lithography techniques may be employed as part of an array fabrication process. For example, in one exemplary implementation, a lithography technique may be employed in which appropriately designed blocks are “stitched” together by overlapping the edges of a step and repeat lithography exposures on a wafer substrate by approximately 0.2 micrometers. In a single exposure, the maximum die size typically is approximately 21 millimeters by 21 millimeters. By selectively exposing different blocks (sides, top & bottoms, core, etc.) very large chips can be defined on a wafer (up to a maximum, in the extreme, of one chip per wafer, commonly referred to as “wafer scale integration”).
(220) In one aspect of the array 100 shown in
(221) In yet another implementation of an array similar to that shown in
(222) In
(223)
(224) Regarding the column select shift registers 194.sub.1 and 194.sub.2, these are implemented in a manner similar to that of the row select shift registers, with each column select shift register comprising 256 series-connected flip-flops and responsible for enabling readout from either the odd columns of the array or the even columns of the array. For example,
(225) With reference again for the moment to
(226) In the embodiment of
(227) In one exemplary implementation, the switches of both the even and odd output drivers 198.sub.1 and 198.sub.2 (e.g., the switches 191.sub.2, 191.sub.4, . . . 191.sub.512 shown in
(228) The ability of the bus 175 to settle quickly following enabling of successive switches in turn facilitates rapid data acquisition from the array. To this end, in some embodiments the switches 191 of the output drivers 198.sub.1 and 198.sub.2 are particularly configured to significantly reduce the settling time of the bus 175. Both the n-channel and the p-channel MOSFETs of a given switch add to the capacitance of the bus 175; however, n-channel MOSFETs generally have better frequency response and current drive capabilities than their p-channel counterparts. In view of the foregoing, Applicants have recognized and appreciated that some of the superior characteristics of n-channel MOSFETs may be exploited to improve settling time of the bus 175 by implementing “asymmetric” switches in which respective sizes for the n-channel MOSFET and p-channel MOSFET of a given switch are different.
(229) For example, in one embodiment, with reference to
(230) While the example above describes asymmetric switches 191 for the output drivers 198.sub.1 and 198.sub.2 in which the n-channel MOSFET is larger than the p-channel MOSFET, it should be appreciated that in another embodiment, the converse may be implemented, namely, asymmetric switches in which the p-channel MOSFET is larger than the n-channel MOSFET. In one aspect of this embodiment, with reference again to
(231) In yet another embodiment directed to facilitating rapid settling of the bus 175 shown in
(232) For purposes of illustration, the bus 175 may have a capacitance in the range of approximately 5 pF to 20 pF in any of the embodiments discussed immediately above (e.g. symmetric switches, asymmetric switches, greater numbers of output drivers, etc.). Of course, it should be appreciated that the capacitance of the bus 175 is not limited to these exemplary values, and that other capacitance values are possible in different implementations of an array according to the present disclosure.
(233) In one aspect of the array design discussed above in connection with
(234)
(235) Generally, the array controller 250 provides various supply voltages and bias voltages to the array 100, as well as various signals relating to row and column selection, sampling of pixel outputs and data acquisition. In particular, the array controller 250 reads one or more analog output signals (e.g., Vout1 and Vout2) including multiplexed respective pixel voltage signals from the array 100 and then digitizes these respective pixel signals to provide measurement data to the computer 260, which in turn may store and/or process the data. In some implementations, the array controller 250 also may be configured to perform or facilitate various array calibration and diagnostic functions, and an optional array UV irradiation treatment as discussed above in connection with
(236) As illustrated in
(237) In another aspect, the power supply 258 includes one or more digital-to-analog converters (DACs) that may be controlled by the computer 260 to allow any or all of the bias voltages, reference voltage, and supply voltages to be changed under software control (i.e., programmable bias settings). For example, a power supply 258 responsive to computer control (e.g., via software execution) may facilitate adjustment of one or more of the supply voltages (e.g., switching between 3.3 Volts and 1.8 Volts depending on chip type as represented by an identification code), and or adjustment of one or more of the bias voltages VB1 and VB2 for pixel drain current, VB3 for column bus drive, VB4 for column amplifier bandwidth, and VBO0 for column output buffer current drive. In some aspects, one or more bias voltages may be adjusted to optimize settling times of signals from enabled pixels. Additionally, the common body voltage V.sub.BODY for all ISFETs of the array may be grounded during an optional post-fabrication UV irradiation treatment to reduce trapped charge, and then coupled to a higher voltage (e.g., VDDA) during diagnostic analysis, calibration, and normal operation of the array for measurement/data acquisition. Likewise, the reference voltage VREF may be varied to facilitate a variety of diagnostic and calibration functions.
(238) As also shown in
(239) Regarding data acquisition from the array 100, in one embodiment the array controller 250 of
(240) The array controller 250 of
(241) In the embodiment of
(242)
(243) In the example of
(244) For example, with respect to the method for detecting nucleotide incorporation, discussed above in connection with
(245) In one implementation, the array controller 250 controls the array 100 to enable rows successively, one at a time. For example, with reference again for the moment to
(246)
(247) In
(248)
(249) As discussed above in connection with
(250) In one embodiment, once pixel values are sampled and digitized by the ADC(s) 254, the computer 260 may be programmed to process pixel data obtained from the array 100 and the array controller 250 so as to facilitate high data acquisition rates that in some cases may exceed a sufficient settling time for pixel voltages represented in a given array output signal. A flow chart illustrating an exemplary method according to one embodiment of the present invention that may be implemented by the computer 260 for processing and correction of array data acquired at high acquisition rates is illustrated in
(251) Regarding pixel settling time, with reference again to
(252)
where A is the difference (V.sub.COLj−V.sub.COLj-1) between two pixel voltage values and k is a time constant associated with a capacitance of the bus 175.
(253) For purposes of the present discussion, pixel “settling time” t.sub.settle defined as the time t is at which ΔV.sub.PIX(t) attains a value that differs from it's final value by an amount that is equal to the peak noise level of the array output signal. If the peak noise level of the array output signal is denoted as n.sub.p, then the voltage at the settling time t.sub.settle is given by ΔV.sub.PIX (t.sub.settle)=A[1−(n.sub.p/A)]. Substituting in Eq. (PP) and solving for t.sub.settle yields
(254)
(255) As indicated above, in one embodiment pixel data may be acquired from the array at data rates that exceed those dictated by the pixel settling time.
(256) Subsequently, in block 506 of
(257) In block 510 of
(258) In addition to controlling the sensor array and ADCs, the timing generator 256 may be configured to facilitate various array calibration and diagnostic functions, as well as an optional UV irradiation treatment. To this end, the timing generator may utilize the signal LSTV indicating the selection of the last row of the array and the signal LSTH to indicate the selection of the last column of the array. The timing generator 256 also may be responsible for generating the CAL signal which applies the reference voltage VREF to the column buffer amplifiers, and generating the UV signal which grounds the drains of all ISFETs in the array during a UV irradiation process (see
(259) With respect to the computer interface 252 of the array controller 250, in one exemplary implementation the interface is configured to facilitate a data rate of approximately 200 MB/sec to the computer 260, and may include local storage of up to 400 MB or greater. The computer 260 is configured to accept data at a rate of 200 MB/sec, and process the data so as to reconstruct an image of the pixels (e.g., which may be displayed in false-color on a monitor). For example, the computer may be configured to execute a general-purpose program with routines written in C++ or Visual Basic to manipulate the data and display is as desired.
(260) The systems described herein, when used for sequencing, typically involve a chemFET array supporting reaction chambers, the chemFETs being coupled to an interface capable of executing logic that converts the signals from the chemFETs into sequencing information.
(261) In some embodiments, the invention encompasses logic (preferably computer executable logic) for polymer sequencing, comprising logic for determining ion pulses associated with an ionic interaction with a PPi or a dNTP or both. Typically, the logic converts characteristic(s) of the ion pulses into polymer sequencing information.
(262) In some embodiments, the invention encompasses logic (preferably computer executable logic) comprising logic for determining a sequence of a nucleic acid template based on time between ion pulses or a characteristic of a single ion pulse. The logic may optionally further comprise logic for determining spatial location of the ion pulse on an array of chemFETs.
(263) In some embodiments, the invention encompasses logic (preferably computer executable logic) comprising logic for determining a sequence of a nucleic acid template based on a duration of time it takes for a particular dNTP to be utilized in a sequencing reaction. Typically, the logic receives signal from one or more chemFETs. Preferably, the sequence is displayed in substantially real time.
(264) In some embodiments, the invention encompasses logic (preferably computer executable logic) for processing ion pulses from an array of chemFETs to determine the sequence of a polymer of interest. The logic may optionally further comprise logic for file management, file storage, and visualization. The logic may also optionally further comprise logic for converting the ion pulses into nucleotide sequences. Preferably, the sequence is displayed in substantially real time.
(265) The sequencing information obtained from the system may be delivered to a handheld computing device, such as a personal digital assistant. Thus, in one embodiment, the invention encompasses logic for displaying a complete genome of an organism on a handheld computing device. The invention also encompasses logic adapted for sending data from a chemFET array to a handheld computing device. Any of such logic may be computer-implemented.
(266) Having discussed several aspects of an exemplary ISFET array and an array controller according to the present disclosure,
(267)
(268) In one aspect of the embodiment shown in
(269) In particular,
(270) For each of the first and second groups of rows, the array 100A of
(271) In one exemplary implementation of the array 100A of
(272) Like the array 100 of
(273)
(274) As noted in
(275)
(276) While the exemplary arrays discussed above in connection with
(277) The array 100D of
(278) The array 100E of
(279) Thus, in various examples of ISFET arrays based on the inventive concepts disclosed herein, an array pitch of approximately nine (9) micrometers (e.g., a sensor surface area of less than ten micrometers by ten micrometers) allows an ISFET array including over 256,000 pixels (i.e., a 512 by 512 array), together with associated row and column select and bias/readout electronics, to be fabricated on a 7 millimeter by 7 millimeter semiconductor die, and a similar sensor array including over four million pixels (i.e., a 2048 by 2048 array, over 4 Mega-pixels) to be fabricated on a 21 millimeter by 21 millimeter die. In other examples, an array pitch of approximately 5 micrometers allows an ISFET array including approximately 1.55 Mega-pixels (i.e., a 1348 by 1152 array) and associated electronics to be fabricated on a 9 millimeter by 9 millimeter die, and an ISFET sensor array including over 14 Mega-pixels and associated electronics on a 22 millimeter by 20 millimeter die. In yet other implementations, using a CMOS fabrication process in which feature sizes of less than 0.35 micrometers are possible (e.g., 0.18 micrometer CMOS processing techniques), ISFET sensor arrays with a pixel size/pitch significantly below 5 micrometers may be fabricated (e.g., array pitch of 2.6 micrometers or pixel/sensor area of less than 8 or 9 micrometers.sup.2), providing for significantly dense ISFET arrays.
(280) In the embodiments of ISFET arrays discussed above, array pixels employ a p-channel ISFET, as discussed above in connection with
(281) For example,
(282) One of the primary differences between the n-channel ISFET pixel design of
(283) In addition to the pixel designs shown in
(284)
(285) In
(286) In
(287)
(288)
(289) Turning from the sensor discussion, we will now be addressing the combining of the ISFET array with a microwell array and the attendant fluidics. As most of the drawings of the microwell array structure are presented only in cross-section or showing that array as only a block in a simplified diagram,
(290) For many uses, to complete a system for sensing chemical reactions or chemical agents using the above-explained high density electronic arrays, techniques and apparatus are required for delivery to the array elements (called “pixels”) fluids containing chemical or biochemical components for sensing. In this section, exemplary techniques and methods will be illustrated, which are useful for such purposes, with desirable characteristics.
(291) As high speed operation of the system may be desired, it is preferred that the fluid delivery system, insofar as possible, not limit the speed of operation of the overall system.
(292) Accordingly, needs exist not only for high-speed, high-density arrays of ISFETs or other elements sensitive to ion concentrations or other chemical attributes, or changes in chemical attributes, but also for related mechanisms and techniques for supplying to the array elements the samples to be evaluated, in sufficiently small reaction volumes as to substantially advance the speed and quality of detection of the variable to be sensed.
(293) There are two and sometimes three components or subsystems, and related methods, involved in delivery of the subject chemical samples to the array elements: (1) macrofluidic system of reagent and wash fluid supplies and appropriate valving and ancillary apparatus, (2) a flow cell and (3) in many applications, a microwell array. Each of these subsystems will be discussed, though in reverse order.
(294) As discussed elsewhere, for many uses, such as in DNA sequencing, it is desirable to provide over the array of semiconductor sensors a corresponding array of microwells, each microwell being small enough preferably to receive only one DNA-loaded bead, in connection with which an underlying pixel in the array will provide a corresponding output signal.
(295) The use of such a microwell array involves three stages of fabrication and preparation, each of which is discussed separately: (1) creating the array of microwells to result in a chip having a coat comprising a microwell array layer; (2) mounting of the coated chip to a fluidic interface; and in the case of DNA sequencing, (3) loading DNA-loaded bead or beads into the wells. It will be understood, of course, that in other applications, beads may be unnecessary or beads having different characteristics may be employed.
(296) The systems described herein can include an array of microfluidic reaction chambers integrated with a semiconductor comprising an array of chemFETs. In some embodiments, the invention encompasses such an array. The reaction chambers may, for example, be formed in a glass, dielectric, photodefineable or etchable material. The glass material may be silicon dioxide.
(297) Preferably, the array comprises at least 100,000 chambers. Preferably, each reaction chamber has a horizontal width and a vertical depth that has an aspect ratio of about 1:1 or less. Preferably, the pitch between the reaction chambers is no more than about 10 microns.
(298) The above-described array can also be provided in a kit for sequencing. Thus, in some embodiments, the invention encompasses a kit comprising an array of microfluidic reaction chambers integrated with an array of chemFETs, and one or more amplification reagents.
(299) In some embodiments, the invention encompasses a sequencing apparatus comprising a dielectric layer overlying a chemFET, the dielectric layer having a recess laterally centered atop the chemFET. Preferably, the dielectric layer is formed of silicon dioxide.
(300) Microwell fabrication may be accomplished in a number of ways. The actual details of fabrication may require some experimentation and vary with the processing capabilities that are available.
(301) In general, fabrication of a high density array of microwells involves photo-lithographically patterning the well array configuration on a layer or layers of material such as photoresist (organic or inorganic), a dielectric, using an etching process. The patterning may be done with the material on the sensor array or it may be done separately and then transferred onto the sensor array chip, of some combination of the two. However, techniques other than photolithography are not to be excluded if they provide acceptable results.
(302) One example of a method for forming a microwell array is now discussed, starting with reference to
(303) After the semiconductor structures, as shown, are formed, the microwell structure is applied to the die. That is, the microwell structure can be formed right on the die or it may be formed separately and then mounted onto the die, either approach being acceptable. To form the microwell structure on the die, various processes may be used. For example, the entire die may be spin-coated with, for example, a negative photoresist such as Microchem's SU-8 2015 or a positive resist/polyimide such as HD Microsystems HD8820, to the desired height of the microwells. The desired height of the wells (e.g., about 4-12 μm in the example of one pixel per well, though not so limited as a general matter) in the photoresist layer(s) can be achieved by spinning the appropriate resist at predetermined rates (which can be found by reference to the literature and manufacturer specifications, or empirically), in one or more layers. (Well height typically may be selected in correspondence with the lateral dimension of the sensor pixel, preferably for a nominal 1:1-1.5:1 aspect ratio, height:width or diameter. Based on signal-to-noise considerations, there is a relationship between dimensions and the required data sampling rates to achieve a desired level of performance. Thus there are a number of factors that will go into selecting optimum parameters for a given application.) Alternatively, multiple layers of different photoresists may be applied or another form of dielectric material may be deposited. Various types of chemical vapor deposition may also be used to build up a layer of materials suitable for microwell formation therein.
(304) Once the photoresist layer (the singular form “layer” is used to encompass multiple layers in the aggregate, as well) is in place, the individual wells (typically mapped to have either one or four ISFET sensors per well) may be generated by placing a mask (e.g., of chromium) over the resist-coated die and exposing the resist to crosslinking (typically UV) radiation. All resist exposed to the radiation (i.e., where the mask does not block the radiation) becomes cross-linked and as a result will form a permanent plastic layer bonded to the surface of the chip (die). Unreacted resist (i.e., resist in areas which are not exposed, due to the mask blocking the light from reaching the resist and preventing crosslinking) is removed by washing the chip in a suitable solvent (i.e., developer) such as propyleneglycolmethylethylacetate (PGMEA) or other appropriate solvent. The resultant structure defines the walls of the microwell array.
(305)
(306)
(307) After exposure of the die/resist to the UV radiation, a second layer of resist may be coated on the surface of the chip. This layer of resist may be relatively thick, such as about 400-450 μm thick, typically. A second mask 3210 (
(308) Other photolithographic approaches may be used for formation of the microwell array, of course, the foregoing being only one example.
(309) For example, contact lithography of various resolutions and with various etchants and developers may be employed. Both organic and inorganic materials may be used for the layer(s) in which the microwells are formed. The layer(s) may be etched on a chip having a dielectric layer over the pixel structures in the sensor array, such as a passivation layer, or the layer(s) may be formed separately and then applied over the sensor array. The specific choice or processes will depend on factors such as array size, well size, the fabrication facility that is available, acceptable costs, and the like.
(310) Among the various organic materials which may be used in some embodiments to form the microwell layer(s) are the above-mentioned SU-8 type of negative-acting photoresist, a conventional positive-acting photoresist and a positive-acting photodefineable polyimide. Each has its virtues and its drawbacks, well known to those familiar with the photolithographic art.
(311) Naturally, in a production environment, modifications will be appropriate.
(312) Contact lithography has its limitations and it may not be the production method of choice to produce the highest densities of wells—i.e., it may impose a higher than desired minimum pitch limit in the lateral directions. Other techniques, such as a deep UV step-and-repeat process, are capable of providing higher resolution lithography and can be used to produce small pitches and possibly smaller well diameters. Of course, for different desired specifications (e.g., numbers of sensors and wells per chip), different techniques may prove optimal. And pragmatic factors, such as the fabrication processes available to a manufacturer, may motivate the use of a specific fabrication method. While novel methods are discussed, various aspects of the invention are limited to use of these novel methods.
(313) Preferably the CMOS wafer with the ISFET array will be planarized after the final metallization process. A chemical mechanical dielectric planarization prior to the silicon nitride passivation is suitable. This will allow subsequent lithographic steps to be done on very flat surfaces which are free of back-end CMOS topography.
(314) By utilizing deep-UV step-and-repeat lithography systems, it is possible to resolve small features with superior resolution, registration, and repeatability. However, the high resolution and large numerical aperture (NA) of these systems precludes their having a large depth of focus. As such, it may be necessary, when using such a fabrication system, to use thinner photodefinable spin-on layers (i.e., resists on the order of 1-2 μm rather than the thicker layers used in contact lithography) to pattern transfer and then etch microwell features to underlying layer or layers. High resolution lithography can then be used to pattern the microwell features and conventional SiO.sub.2 etch chemistries can be used—one each for the bondpad areas and then the microwell areas—having selective etch stops; the etch stops then can be on aluminum bondpads and silicon nitride passivation (or the like), respectively. Alternatively, other suitable substitute pattern transfer and etch processes can be employed to render microwells of inorganic materials.
(315) Another approach is to form the microwell structure in an organic material. For example, a dual-resist “soft-mask” process may be employed, whereby a thin high-resolution deep-UV resist is used on top of a thicker organic material (e.g., cured polyimide or opposite-acting resist). The top resist layer is patterned. The pattern can be transferred using an oxygen plasma reactive ion etch process. This process sequence is sometimes referred to as the “portable conformable mask” (PCM) technique. See B. J. Lin et al., “Practicing the Novolac deep-UV portable conformable masking technique”, Journal of Vacuum Science and Technology 19, No. 4, 1313-1319 (1981); and A. Cooper et al, “Optimization of a photosensitive spin-on dielectric process for copper inductor coil and interconnect protection in RF SoC devices.”
(316) Alternatively a “drill-focusing” technique may be employed, whereby several sequential step-and-repeat exposures are done at different focal depths to compensate for the limited depth of focus (DOF) of high-resolution steppers when patterning thick resist layers. This technique depends on the stepper NA and DOF as well as the contrast properties of the resist material.
(317) Another PCM technique may be adapted to these purposes, such as that shown in U.S. patent application publication no. 2006/0073422 by Edwards et al. This is a three-layer PCM process and it is illustrated in
(318) In a first step, 3320, a layer of high contrast negative-acting photoresist such as type Shipley InterVia Photodielectric Material 8021 (IV8021) 3322 is spun on the surface of a wafer, which we shall assume to be the wafer providing the substrate 3312 of
(319) Although as shown above, the wells bottom out (i.e. terminate) on the top passivation layer of the ISFETs, it is believed that an improvement in ISFET sensor performance (i.e. such as signal-to-noise ratio) can be obtained if the active bead(s) is(are) kept slightly elevated from the ISFET passivation layer. One way to do so is to place a spacer “bump” within the boundary of the pixel microwell. An example of how this could be rendered would be not etching away a portion of the layer-or-layers used to form the microwell structure (i.e. two lithographic steps to form the microwells—one to etch part way done, the other to pattern the bump and finish the etch to bottom out), by depositing and lithographically defining and etching a separate layer to form the “bump”, by using a permanent photo-definable material for the bump once the microwells are complete, or by forming the bump prior to forming the microwell. The bump feature is shown as 3350 in
(320) Using a 6 um (micron) thick layer of tetra-methyl-ortho-silicate (TEOS) as a SiO.sub.2-like layer for microwell formation,
(321) Microwells 3302A are formed in the TEOS layer 3304A. The wells extend about 4 um into the 6 um thick layer. Typically, the etched well bottoms on an etch-stop material which may be, for example, an oxide, an organic material or other suitable material known in semiconductor processing for etch-stopping use. A thin layer of etch stop material may be formed on top of a thicker layer of polyimide or other suitable dielectric, such that there is about 2 um of etch stop+polyimide between the well bottom and the Metal4 (M4) layer of the chip in which the extended gate electrode 3308A is formed for each underlying ISFET in the array. As labeled on the side, the CMOS metallization layers M3, M2 and M1, which form lower level interconnects and structures, are shown, with the ISFET channels being formed in the areas indicated by arrows 3310A.
(322) In the orthogonal cross-sectional view (i.e., looking down from the top), the wells may be formed in either round or square shape. Round wells may improve bead capture and may obviate the need for packing beads at the bottom or top of the wells.
(323) The tapered slopes to the sides of the microwells also may be used to advantage. Referring to
(324) Thus, microwells can be fabricated by any high aspect ratio photo-definable or etchable thin-film process, that can provide requisite thickness (e.g., about 4-10 um). Among the materials believed to be suitable are photosensitive polymers, deposited silicon dioxide, non-photosensitive polymer which can be etched using, for example, plasma etching processes, etc. In the silicon dioxide family, TEOS and silane nitrous oxide (SILOX) appear suitable. The final structures are similar but the various materials present differing surface compositions that may cause the target biology or chemistry to react differently.
(325) When the microwell layer is formed, it may be necessary to provide an etch stop layer so that the etching process does not proceed further than desired. For example, there may be an underlying layer to be preserved, such as a low-K dielectric. The etch stop material should be selected according to the application. SiC and SiN materials may be suitable, but that is not meant to indicate that other materials may not be employed, instead. These etch-stop materials can also serve to enhance the surface chemistry which drives the ISFET sensor sensitivity, by choosing the etch-stop material to have an appropriate point of zero charge (PZC). Various metal oxides may be suitable addition to silicon dioxide and silicon nitride.
(326) The PZCs for various metal oxides may be found in various texts, such as “Metal Oxides—Chemistry and Applications” by J. Fierro. We have learned that Ta.sub.2O.sub.5 may be preferred as an etch stop over Al.sub.2O.sub.3 because the PZC of Al.sub.2O.sub.3 is right at the pH being used (i.e., about 8.8) and, hence, right at the point of zero charge. In addition Ta.sub.2O.sub.5 has a higher sensitivity to pH (i.e., mV/pH), another important factor in the sensor performance. Optimizing these parameters may require judicious selection of passivation surface materials.
(327) Using thin metal oxide materials for this purpose (i.e., as an etch stop layer) is difficult due to the fact of their being so thinly deposited (typically 200-500 A). A post-microwell fabrication metal oxide deposition technique may allow placement of appropriate PZC metal oxide films at the bottom of the high aspect ratio microwells.
(328) Electron-beam depositions of (a) reactively sputtered tantalum oxide, (b) non-reactive stoichiometric tantalum oxide, (c) tungsten oxide, or (d) Vanadium oxide may prove to have superior “down-in-well” coverage due to the superior directionality of the deposition process.
(329) The array typically comprises at least 100 microfluidic wells, each of which is coupled to one or more chemFET sensors. Preferably, the wells are formed in at least one of a glass (e.g., SiO.sub.2), a polymeric material, a photodefinable material or a reactively ion etchable thin film material. Preferably, the wells have a width to height ratio less than about 1:1. Preferably the sensor is a field effect transistor, and more preferably a chemFET. The chemFET may optionally be coupled to a PPi receptor. Preferably, each of the chemFETs occupies an area of the array that is 10.sup.2 microns or less.
(330) In some embodiments, the invention encompasses a sequencing device comprising a semiconductor wafer device coupled to a dielectric layer such as a glass (e.g., SiO.sub.2), polymeric, photodefinable or reactive ion etchable material in which reaction chambers are formed. Typically, the glass, dielectric, polymeric, photodefinable or reactive ion etchable material is integrated with the semiconductor wafer layer. In some instances, the glass, polymeric, photodefinable or reactive ion etchable layer is non-crystalline. In some instances, the glass may be SiO.sub.2. The device can optionally further comprise a fluid delivery module of a suitable material such as a polymeric material, preferably an injection moldable material. More preferably, the polymeric layer is polycarbonate.
(331) In some embodiments, the invention encompasses a method for manufacturing a sequencing device comprising: using photolithography, generating wells in a glass, dielectric, photodefinable or reactively ion etchable material on top of an array of transistors.
(332) The process of using the assembly of an array of sensors on a chip combined with an array of microwells to sequence the DNA in a sample is referred to as an “experiment.” Executing an experiment requires loading the wells with the DNA-bound beads and the flowing of several different fluid solutions (i.e., reagents and washes) across the wells. A fluid delivery system (e.g., valves, conduits, pressure source(s), etc.) coupled with a fluidic interface is needed which flows the various solutions across the wells in a controlled even flow with acceptably small dead volumes and small cross contamination between sequential solutions. Ideally, the fluidic interface to the chip (sometimes referred to as a “flow cell”) would cause the fluid to reach all microwells at the same time. To maximize array speed, it is necessary that the array outputs be available at as close to the same time as possible. The ideal clearly is not possible, but it is desirable to minimize the differentials, or skews, of the arrival times of an introduced fluid, at the various wells, in order to maximize the overall speed of acquisition of all the signals from the array.
(333) Flow cell designs of many configurations are possible; thus the system and methods presented herein are not dependent on use of a specific flow cell configuration. It is desirable, though, that a suitable flow cell substantially conform to the following set of objectives:
(334) have connections suitable for interconnecting with a fluidics delivery system—e.g., via appropriately-sized tubing;
(335) have appropriate head space above wells;
(336) minimize dead volumes encountered by fluids;
(337) minimize small spaces in contact with liquid but not quickly swept clean by flow of a wash fluid through the flow cell (to minimize cross contamination);
(338) be configured to achieve uniform transit time of the flow over the array;
(339) generate or propagate minimal bubbles in the flow over the wells;
(340) be adaptable to placement of a removable reference electrode inside or as close to the flow chamber as possible;
(341) facilitate easy loading of beads;
(342) be manufacturable at acceptable cost; and
(343) be easily assembled and attached to the chip package.
(344) Satisfaction of these criteria so far as possible will contribute to system performance positively. For example, minimization of bubbles is important so that signals from the array truly indicate the reaction in a well rather than being spurious noise.
(345) Each of several example designs will be discussed, meeting these criteria in differing ways and degrees. In each instance, one typically may choose to implement the design in one of two ways: either by attaching the flow cell to a frame and gluing the frame (or otherwise attaching it) to the chip or by integrating the frame into the flow cell structure and attaching this unified assembly to the chip. Further, designs may be categorized by the way the reference electrode is integrated into the arrangement. Depending on the design, the reference electrode may be integrated into the flow cell (e.g., form part of the ceiling of the flow chamber) or be in the flow path (typically to the outlet or downstream side of the flow path, after the sensor array).
(346) A first example of a suitable experiment apparatus 3410 incorporating such a fluidic interface is shown in
(347) The apparatus comprises a semiconductor chip 3412 (indicated generally, though hidden) on or in which the arrays of wells and sensors are formed, and a fluidics assembly 3414 on top of the chip and delivering the sample to the chip for reading. The fluidics assembly includes a portion 3416 for introducing fluid containing the sample, a portion 3418 for allowing the fluid to be piped out, and a flow chamber portion 3420 for allowing the fluid to flow from inlet to outlet and along the way interact with the material in the wells. Those three portions are unified by an interface comprising a glass slide 3422 (e.g., Erie Microarray Cat #C22-5128-M20 from Erie Scientific Company, Portsmouth, N. H., cut in thirds, each to be of size about 25 mm×25 mm).
(348) Mounted on the top face of the glass slide are two fittings, 3424 and 3426, such as nanoport fittings Part #N-333 from Upchurch Scientific of Oak Harbor, Wash. One port (e.g., 3424) serves as an inlet delivering liquids from the pumping/valving system described below but not shown here. The second port (e.g., 3426) is the outlet which pipes the liquids to waste. Each port connects to a conduit 3428, 3432 such as flexible tubing of appropriate inner diameter. The nanoports are mounted such that the tubing can penetrate corresponding holes in the glass slide. The tube apertures should be flush with the bottom surface of the slide.
(349) On the bottom of the glass slide, flow chamber 3420 may comprise various structures for promoting a substantially laminar flow across the microwell array. For example, a series of microfluidic channels fanning out from the inlet pipe to the edge of the flow chamber may be patterned by contact lithography using positive photoresists such as SU-8 photoresist from MicroChem Corp. of Newton, Mass. Other structures will be discussed below.
(350) The chip 3412 will in turn be mounted to a carrier 3430, for packaging and connection to connector pins 3432.
(351) For ease of description, to discuss fabrication starting with
(352) A layer of photoresist 3810 is applied to the “top” of the slide (which will become the “bottom” side when the slide and its additional layers is turned over and mounted to the sensor assembly of ISFET array with microwell array on it). Layer 3810 may be about 150 μm thick in this example, and it will form the primary fluid carrying layer from the end of the tubing in the nanoports to the edge of the sensor array chip. Layer 3810 is patterned using a mask such as the mask 3910 of
(353) A second layer of photoresist is formed quite separately, not on the resist 3810 or slide 3422. Preferably it is formed on a flat, flexible surface (not shown), to create a peel-off, patterned plastic layer. As shown in
(354) The other alignment mark or set of marks produced by pattern 4022 is used for alignment with a subsequent layer to be discussed.
(355) The second layer is preferably about 150 μm deep and it will cover the fluid-carrying channel with the exception of a slit about 150 μm long at each respective edge of the sensor array chip, under slit-forming regions 4014 and 4016.
(356) Once the second layer of photoresist is disposed on the first layer, a third patterned layer of photoresist is formed over the second layer, using a mask such as mask 4110, shown in
(357)
(358) The fluidics assembly may be secured to the sensor array chip assembly by applying an adhesive to parts of mating surfaces of those two assemblies, and pressing them together, in alignment.
(359) Though not illustrated in
(360) Another way to introduce the reference electrode is shown in
(361) Achieving a uniform flow front and eliminating problematic flow path areas is desirable for a number of reasons. One reason is that very fast transition of fluid interfaces within the system's flow cell is desired for many applications, particularly gene sequencing. In other words, an incoming fluid must completely displace the previous fluid in a short period of time. Uneven fluid velocities and diffusion within the flow cell, as well as problematic flow paths, can compete with this requirement. Simple flow through a conduit of rectangular cross section can exhibit considerable disparity of fluid velocity from regions near the center of the flow volume to those adjacent the sidewalls, one sidewall being the top surface of the microwell layer and the fluid in the wells. Such disparity leads to spatially and temporally large concentration gradients between the two traveling fluids. Further, bubbles are likely to be trapped or created in stagnant areas like sharp corners interior the flow cell. (The surface energy (hydrophilic vs. hydrophobic) can significantly affect bubble retention. Avoidance of surface contamination during processing and use of a surface treatment to create a more hydrophilic surface should be considered if the as-molded surface is too hydrophobic.) Of course, the physical arrangement of the flow chamber is probably the factor which most influences the degree of uniformity achievable for the flow front.
(362) One approach is to configure the flow cross section of the flow chamber to achieve flow rates that vary across the array width so that the transit times are uniform across the array. For example, the cross section of the diffuser (i.e., flow expansion chamber) section 3416, 3610 may be made as shown at 4204A in
(363) Another configuration, shown in
(364) FIGS. 42F2-42F8 illustrate an example of a single-piece, injection-molded (preferably of polycarbonate) flow cell member 42F200 which may be used to provide baffles 4220F, a ceiling to the flow chamber, fluid inlet and outlet ports and even the reference electrode. FIG. 42F7 shows an enlarged view of the baffles on the bottom of member 42F200 and the baffles are shown as part of the underside of member 42F200 in FIG. 42F6. As it is difficult to form rectangular features in small dimensions by injection molding, the particular instance of these baffles, shown as 4220F′, are triangular in cross section.
(365) In FIG. 42F2, there is a top, isometric view of member 42F200 mounted onto a sensor array package 42F300, with a seal 42F202 formed between them and contact pins 42F204 depending from the sensor array chip package. FIGS. 42F3 and 42F4 show sections, respectively, through section lines H-H and I-I of FIG. 42F5, permitting one to see in relationship the sensor array chip 42F250, the baffles 4220F′ and fluid flow paths via inlet 42F260 and outlet 42F270 ports.
(366) By applying a metallization to bottom 42F280 of member 42F200, the reference electrode may be formed.
(367) Various other locations and approaches may be used for introducing fluid flow into the flow chamber, as well. In addition to embodiments in which fluid may be introduced across the width of an edge of the chip assembly 42F1, as in
(368)
(369) A variation on this idea is depicted in
(370) In all cases, attention should be given to assuring a thorough washing of the entire flow chamber, along with the microwells, between reagent cycles. Flow disturbances may exacerbate the challenge of fully cleaning out the flow chamber.
(371) Flow disturbances may also induce or multiply bubbles in the fluid. A bubble may prevent the fluid from reaching a microwell, or delay its introduction to the microwell, introducing error into the microwell reading or making the output from that microwell useless in the processing of outputs from the array. Thus, care should be taken in selecting configurations and dimensions for the flow disruptor elements to manage these potential adverse factors. For example, a tradeoff may be made between the heights of the disruptor elements and the velocity profile change that is desired.
(372)
(373) In the illustrated embodiment, the reference electrode is introduced to the top of the flow chamber via a bore 4325 in the member 4320. The placement of the removable reference electrode is facilitated by a silicone sleeve 4360 and an epoxy stop ring 4370 (see the blow-up of
(374)
(375)
(376) Yet another alternative for a fluidics assembly, as shown in
(377) Some of the foregoing alternative embodiments also may be implemented in a hybrid plastic/PDMS configuration. For example, as shown in
(378) The fluidic structure may also be made from glass as discussed above, such as photo-definable (PD) glass. Such a glass may have an enhanced etch rate in hydrofluoric acid once selectively exposed to UV light and features may thereby be micromachined on the top-side and back-side, which when glued together can form a three-dimensional low aspect ratio fluidic cell.
(379) An example is shown in
(380) Nanoports may be secured over the nanoport fluidic holes to facilitate connection of input and output tubing.
(381) A central bore 5550 may be etched through the glass layers for receiving a reference electrode, 5560. The electrode may be secured and sealed in place with a silicone collar 5570 or like structure; or the electrode may be equipped integrally with a suitable washer for effecting the same purpose.
(382) By using glass materials for the two-layer fluidic cell, the reference electrode may also be a conductive layer or pattern deposited on the bottom surface of the second glass layer (not shown). Or, as shown in
(383) Another alternative is to integrate the reference electrode to the sequencing chip/flow cell by using a metalized surface on the ceiling of the flow chamber—i.e., on the underside of the member forming the ceiling of the fluidic cell. An electrical connection to the metalized surface may be made in any of a variety of ways, including, but not limited to, by means of applying a conductive epoxy to the ceramic package seal ring that, in turn, may be electrically connected through a via in the ceramic substrate to a spare pin at the bottom of the chip package. Doing this would allow system-level control of the reference potential in the fluid cell by means of inputs through the chip socket mount to the chip's control electronics.
(384) By contrast, an externally inserted electrode requires extra fluid tubing to the inlet port, which requires additional fluid flow between cycles.
(385) Ceramic pin grid array (PGA) packaging may be used for the ISFET array, allowing customized electrical connections between various surfaces on the front face with pins on the back.
(386) The flow cell can be thought of as a “lid” to the ISFET chip and its PGA. The flow cell, as stated elsewhere, may be fabricated of many different materials. Injection molded polycarbonate appears to be quite suitable. A conductive metal (e.g., gold) may be deposited using an adhesion layer (e.g., chrome) to the underside of the glow cell roof (the ceiling of the flow chamber). Appropriate low-temperature thin-film deposition techniques preferably are employed in the deposition of the metal reference electrode due to the materials (e.g., polycarbonate) and large step coverage topography at the bottom-side of the fluidic cell (i.e., the frame surround of ISFET array). One possible approach would be to use electron-beam evaporation in a planetary system.
(387) The active electrode area is confined to the central flow chamber inside the frame surround of the ISFET array, as that is the only metalized surface that would be in contact with the ionic fluid during sequencing.
(388) Once assembly is complete—conductive epoxy (e.g., Epo-Tek H20E or similar) may be dispensed on the seal ring with the flow cell aligned, placed, pressed and cured—the ISFET flow cell is ready for operation with the reference potential being applied to the assigned pin of the package.
(389) The resulting fluidic system connections to the ISFET device thus incorporate shortened input and output fluidic lines, which is desirable.
(390) Still another example embodiment for a fluidic assembly is shown in
(391) Still further examples of flow cell structures are shown in
(392) Whether glass or plastic or other material is used to form the flow cell, it may be desirable, especially with larger arrays, to include in the inlet chamber of the flow cell, between the inlet conduit and the front edge of the array, not just a gradually expanding (fanning out) space, but also some structure to facilitate the flow across the array being suitably laminar. Using the bottom layer 5990 of an injection molded flow cell as an example, one example type of structure for this purpose, shown in
(393) The above-described systems for sequencing typically utilize a laminar fluid flow system to sequence a biological polymer. In part, the fluid flow system preferably includes a flow chamber formed by the sensor chip and a single piece, injection molded member comprising inlet and outlet ports and mountable over the chip to establish the flow chamber. The surface of such member interior to the chamber is preferably formed to facilitate a desired expedient fluid flow, as described herein.
(394) In some embodiments, the invention encompasses an apparatus for detection of ion pulses comprising a laminar fluid flow system. Preferably, the apparatus is used for sequencing a plurality of nucleic acid templates, wherein the nucleic acid templates are optionally deposited on an array.
(395) The apparatus typically includes a fluidics assembly comprising a member comprising one or more apertures for non-mechanically directing a fluid to flow to an array of at least 100 K, 500 K, or 1 M microfluidic reaction chambers such that the fluid reaches all of the microfluidic reaction chambers at the same time or substantially the same time. Typically, the fluid flow is parallel to the sensor surface. Typically, the assembly has a Reynolds number of less than 1000, 500, 200, 100, 50, 20, or 10. Preferably, the member further comprises a first aperture for directing fluid towards the sensor array and a second aperture for directing fluid away from the sensor array.
(396) In some embodiments, the invention encompasses a method for directing a fluid to a sensor array comprising: providing a fluidics assembly comprising an aperture fluidly coupling a fluid source to the sensor array; and non-mechanically directing a fluid to the sensor array. By “non-mechanically” it is meant that the fluid is moved under pressure from a gaseous pressure source, as opposed to a mechanical pump.
(397) In some embodiments, the invention encompasses an array of wells, each of which is coupled to a lid having an inlet port and an outlet port and a fluid delivery system for delivering and removing fluid from said inlet and outlet ports non-mechanically.
(398) In some embodiments, the invention encompasses a method for sequencing a biological polymer utilizing the above-described apparatus, comprising: directing a fluid comprising a monomer to an array of reaction chambers wherein the fluid has a fluid flow Reynolds number of at most 2000, 1000, 200, 100, 50, or 20. The method may optionally further comprise detecting an ion pulse from each said reaction chamber. The ion pulse is typically detected by ion diffusion to the sensor surface. There are various other ways of providing a fluidics assembly for delivering an appropriate fluid flow across the microwell and sensor array assembly, and the forgoing examples are thus not intended to be exhaustive.
(399) Commercial flow-type fluidic electrodes, such as silver chloride proton-permeable electrodes, may be inserted in series in a fluidic line and are generally designed to provide a stable electrical potential along the fluidic line for various electrochemical purposes. In the above-discussed system, however, such a potential must be maintained at the fluidic volume in contact with the microwell ISFET chip. With conventional silver chloride electrodes, it has been found difficult, due to an electrically long fluidic path between the chip surface and the electrode (through small channels in the flow cell), to achieve a stable potential. This led to reception of noise in the chip's electronics. Additionally, the large volume within the flow cavity of the electrode tended to trap and accumulate gas bubbles that degraded the electrical connection to the fluid. With reference to
(400) A complete system for using the sensor array will include suitable fluid sources, valving and a controller for operating the valving to low reagents and washes over the microarray or sensor array, depending on the application. These elements are readily assembled from off-the-shelf components, with and the controller may readily be programmed to perform a desired experiment.
(401) It should be understood that the readout at the chemFET may be current or voltage (and change thereof) and that any particular reference to either readout is intended for simplicity and not to the exclusion of the other readout. Therefore any reference in the following text to either current or voltage detection at the chemFET should be understood to contemplate and apply equally to the other readout as well. In important embodiments, the readout reflects a rapid, transient change in concentration of an analyte. The concentration of more than one analyte may be detected at different times. Such measurements are to be contrasted with prior art methods which focused on steady state concentration measurements.
(402) As already discussed, the apparatus and systems of the invention can be used to detect and/or monitor interactions between various entities. These interactions include biological and chemical reactions and may involve enzymatic reactions and/or non-enzymatic interactions such as but not limited to binding events. As an example, the invention contemplates monitoring enzymatic reactions in which substrates and/or reagents are consumed and/or reaction intermediates, byproducts and/or products are generated. An example of a reaction that can be monitored according to the invention is a nucleic acid synthesis method such as one that provides information regarding nucleic acid sequence. This reaction will be discussed in greater detail herein.
(403) In the context of a sequencing reaction, the apparatus and system provided herein is able to detect nucleotide incorporation based on changes in the chemFET current and/or voltage, as those latter parameters are interrelated. Current changes may be the result of one or more of the following events either singly or some combination thereof: generation of PPi, generation of Pi (e.g., in the presence of pyrophosphatase), generation of hydrogen (and concomitant changes in pH for example in the presence of low strength buffer), reduced concentration of unincorporated dNTP at the chemFET surface, delayed arrival of unincorporated dNTP at the chemFET surface, and the like. The methods described herein are able to detect changes in analyte concentration at the chemFET surface, and such changes may result from one or more of the afore-mentioned events. The invention contemplates the use of a chemFET such as an ISFET in the sequencing methods described herein, even if the readout is independent of (or insensitive to) pH. In other words, the invention contemplates the use of an ISFET for the detection of analytes such as PPi and unincorporated nucleotides. The methods provided herein in regards to sequencing can be contrasted to those in the literature including Pourmand et al. PNAS 2006 103(17):6466-6470. As discussed herein, the invention contemplates methods for determining the nucleotide sequence (i.e., the “sequence”) of a nucleic acid. Such methods involve the synthesis of a new nucleic acid (primed by a pre-existing nucleic acid, as will be appreciated by those of ordinary skill), based on the sequence of a template nucleic acid. That is, the sequence of the newly synthesized nucleic acid is complimentary to the sequence of the template nucleic acid and therefore knowledge of sequence of the newly synthesized nucleic acid yields information about the sequence of the template nucleic acid. Knowledge of the sequence of the newly synthesized nucleic acid is derived by determining whether a known nucleotide has been incorporated into the newly synthesized nucleic acid and, if so, how many of such known nucleotides have been incorporated. Nucleotide incorporation can be monitored in a number of ways, including the production of products such as PPi, Pi and/or H.sup.+.
(404)
(405)
(406) Thus, some embodiments the invention contemplate that upon introduction of nucleotides that are complementary to the next position to be sequenced on the template nucleic acid, such nucleotides will be effectively consumed yielding PPi which is then detected at the chemFET. Only once the synthesis reactions have gone to completion (i.e., the appropriate number of nucleotides have been incorporated into all newly synthesized nucleic acids) is the population of unincorporated nucleotides able to diffuse to and be detected by the chemFET surface. The delay in the detection of the unincorporated nucleotides is indicative of how many such nucleotides have been incorporated. These delays are modeled in
(407) The incorporation of a dNTP into the nucleic acid strand releases PPi which can then be hydrolyzed to two orthophosphates (Pi) and one hydrogen ion. The generation of the hydrogen ion therefore can be detected as an indicator of nucleotide incorporation. Alternatively, Pi may be detected directly or indirectly. Any and all of these events (and more as described herein) may be detected at the chemFET thereby causing a current change that correlates with nucleotide incorporation.
(408) The limit of this dual pulse detection occurs when the two pulses cannot be resolved in time (i.e., the two peaks may merge into one, where to and t.sub.1 peaks are not distinguishable). This limit, for the case N.sub.poly=0 vs N.sub.poly=1 is given by the following: Consider a first peak N.sub.poly=0 which has a width=W.sub.0 and standard deviation σ.sub.W0 and a second peak N.sub.poly=1 which has a width=W.sub.1 and standard deviation σ.sub.W1. For a P.sub.e of 1% [where P.sub.e=erfc(SNR/2*(sqr(2)))], the signal-to-noise ratio (SNR)=5.15. Therefore the limit of detection is the difference in peak widths, ΔW, is given by ΔW=(W.sub.1−W.sub.0)>5.15*σ.sub.W0. Similar calculations may be done for other SNRs.
(409) In certain embodiments the chemFET (or ISFET) readouts are independent of (or insensitive to) changes in the concentration of Pi and/or H.sup.+. In these embodiments, the rate of PPi hydrolysis to Pi with the concomitant release of H.sup.+ is reduced, sometimes significantly. The PPi to Pi conversion may be catalyzed by an enzyme such as a pyrophosphatase and/or by certain ions. In the absence of the both (and in some cases, either) the enzyme and the particular ion, the rate of conversion is sufficiently slow to be negligible. (See for example J. Chem, Soc. Farady Trans. 2007, 93:4295 which reports the rate constant for the conversion in the absence of enzyme or ion to be on the order of about 3.sup.−1 years.sup.−1, as compared to Biochemistry, 2002, 41:12025 which reports the rate in the presence of enzyme and Mg.sup.2+ to be on the order of about 10.sup.3 s.sup.−1.) A pH insensitive environment may be one that is buffered sufficiently to obscure any changes in pH resulting from release of hydrogen. The reaction may be pH insensitive or independent as a result of there being few if any hydrogen ions released as a result of the reaction.
(410) The nucleic acid being sequenced is referred to herein as the target nucleic acid. Target nucleic acids include but are not limited to DNA such as but not limited to genomic DNA, mitochondrial DNA, cDNA and the like, and RNA such as but not limited to mRNA, miRNA, and the like. The nucleic acid may be from any source including naturally occurring sources or synthetic sources. The nucleic acids may be PCR products, cosmids, plasmids, naturally occurring or synthetic libraries, and the like. The invention is not intended to be limited in this regard. The methods provided herein can be used to sequence nucleic acids of any length. To be clear, the Examples provide a proof of principle demonstration of the sequencing of four templates of known sequence. This artificial model is intended to show that embodiments of the apparatuses and systems described herein are able to readout nucleotide incorporation that correlates to the known sequence of the templates. This is not intended to represent typical use of the method or system in the field. The following is a brief description of these methods.
(411) Target nucleic acids are prepared using any manner known in the art. As an example, genomic DNA may be harvested from a sample according to techniques known in the art (see for example Sambrook et al. “Maniatis”). Following harvest, the DNA may be fragmented to yield nucleic acids of smaller length. The resulting fragments may be on the order of hundreds, thousands, or tens of thousands nucleotides in length. In some embodiments, the fragments are 200-1000 base pairs in size, or 300-800 base pairs in size, about 200, about 300, about 400, about 500, about 600, about 700, about 800, about 900, or about 1000 base pairs in length, although they are not so limited. Nucleic acids may be fragmented by any means including but not limited to mechanical, enzymatic or chemical means. Examples include shearing, sonication, nebulization, endonuclease (e.g., DNase I) digestion, amplification such as PCR amplification, or any other technique known in the art to produce nucleic acid fragments, preferably of a desired length. As used herein, fragmentation also embraces the use of amplification to generate a population of smaller sized fragments of the target nucleic acid. That is, the target nucleic acids may be melted and then annealed to two (and preferably more) amplification primers and then amplified using for example a thermostable polymerase (such as Taq). An example is a massively parallel PCR-based amplification. Fragmentation can be followed by size selection techniques to enrich or isolate fragments of a particular length or size. Such techniques are also known in the art and include but are not limited to gel electrophoresis or SPRI.
(412) Alternatively, target nucleic acids that are already of sufficient small size (or length) may be used. Such target nucleic acids include those derived from an exon enrichment process. Thus, rather than fragmenting (randomly or non-randomly) longer target nucleic acids, the targets may be nucleic acids that naturally exist or can be isolated in shorter, useable lengths such as mRNAs, cDNAs, exons, PCR products (as described above), and the like. See Albert et al. Nature Methods 2007 4(11):903-905 (microarray hybridization of exons and locus-specific regions), Porreca et al. Nature Methods 2007 4(11):931-936, and Okou et al. Nature Methods 2007 4(11):907-909 for methods of isolating and/or enriching sequences such as exons prior to sequencing.
(413) In some embodiments, the size selected target nucleic acids are ligated to adaptor sequences on both the 5′ and 3′ ends. These adaptor sequences comprise sequences complementary to amplification primer sequences, to be used in amplifying the target nucleic acids. One adaptor sequence may also comprise a sequence complementary to the sequencing primer. The opposite adaptor sequence may comprise a moiety that facilitates binding of the nucleic acid to a solid support such as but not limited to a bead. An example of such a moiety is a biotin molecule (or a double biotin moiety, as described by Diehl et al. Nature Methods, 2006, 3(7):551-559) and such a labeled nucleic acid can therefore be bound to a solid support having avidin or streptavidin groups. Another moiety that can be used is the NHS-ester and amine affinity pair. It is to be understood that the invention is not limited in this regard and one of ordinary skill is able to substitute these affinity pairs with other binding pairs. The resulting nucleic acid is referred to herein as a template nucleic acid. The template nucleic acid comprises at least the target nucleic acid and usually comprises nucleotide sequences in addition to the target at both the 5′ and 3′ ends.
(414) In some embodiments, the template nucleic acid is able to self-anneal thereby creating a 3′ end from which to incorporate nucleotide triphosphates. In such instances, there is no need for a separate sequencing primer since the template acts as both template and primer. See Eriksson et al. Electrophoresis 2004 25:20-27 for a discussion of the use of self-annealing template in a pyrosequencing reaction.
(415) In still other embodiments, the template is used in a double stranded (or partially double stranded) from which is nicked preferably within an adaptor sequence (at the free end of the bead bound template). The opening in the double stranded template allows a polymerase such as DNA polymerase to enter at the nicked site and begin nucleotide incorporation. These openings can be introduced into the template in a controlled manner particularly since the sequence acted upon by DNA nickase (i.e., the nicking enzyme in this instance) is know (i.e., 5′GAGTC3′). See Zyrina et al. for a discussion of these consensus nicking target sequences.
(416) In some instances a spacer is used to distance the template nucleic acid (and in particular the target nucleic acid sequence comprised therein) from a solid support such as a bead. This facilitates sequencing of the end of the target closest to the bead, for instance. Examples of suitable linkers are known in the art (see Diehl et al. Nature Methods, 2006, 3(7):551-559) and include but are not limited to carbon-carbon linkers such as but not limited to iSp18.
(417) The solid support to which the template nucleic acids are bound is referred to herein as the “capture solid support”. If the solid support is a bead, then such bead is referred to herein as a “capture bead”. The beads can be made of any material including but not limited to cellulose, cellulose derivatives, gelatin, acrylic resins, glass, silica gels, polyvinyl pyrrolidine (PVP), copolymers of vinyl and acrylamide, polystyrene, polystyrene cross-linked with divinylbenzene or the like (see, Merrifield Biochemistry 1964, 3, 1385-1390), polyacrylamides, latex gels, dextran, crosslinked dextrans (e.g., Sephadex™), rubber, silicon, plastics, nitrocellulose, natural sponges, metal, and agarose gel (Sepharose™). In one embodiment, the beads are streptavidin-coated beads. The bead diameter will depend on the density of the chemFET and microwell array used with larger arrays (and thus smaller sized wells) requiring smaller beads. Generally the bead size may be about 1-10 μM, and more preferably 2-6 μM. In some embodiments, the beads are about 5.91 μM while in other embodiments the beads are about 2.8 μM. In still other embodiments, the beads are about 1.5 μm, or about 1 μm in diameter. It is to be understood that the beads may or may not be perfectly spherical in shape. It is also to be understood that other beads may be used and other mechanisms for attaching the nucleic acid to the beads may be used. In some instances the capture beads (i.e., the beads on which the sequencing reaction occurs) are the same as the template preparation beads including the amplification beads.
(418) Important aspects of the invention contemplate sequencing a plurality of different template nucleic acids simultaneously. This may be accomplished using the sensor arrays described herein. In one embodiment, the sensor arrays are overlayed (and/or integral with) an array of microwells (or reaction chambers or wells, as those terms are used interchangeably herein), with the proviso that there be at least one sensor per microwell. Present in a plurality of microwells is a population of identical copies of a template nucleic acid. There is no requirement that any two microwells carry identical template nucleic acids, although in some instances such templates may share overlapping sequence. Thus, each microwell comprises a plurality of identical copies of a template nucleic acid, and the templates between microwells may be different.
(419) The microwells may vary in size between arrays. The microwell size may be described in terms of cross section. The cross section may refer to a “slice” parallel to the depth (or height) of the well, or it may be a slice perpendicular to the depth (or height) of the well.
(420) The size of these microwells may be described in terms of a width (or diameter) to height ratio. In some embodiments, this ratio is 1:1 to 1:1.5. The bead to well size (e.g., the bead diameter to well width, diameter, or height) is preferably in the range of 0.6-0.8.
(421) The microwells may be square in cross-section, but they are not so limited. The dimensions at the bottom of a microwell (i.e., in a cross section that is perpendicular to the depth of the well) may be 1.5 μm by 1.5 μm, or it may be 1.5 μm by 2 μm. Various diameters are shown in the Examples and include but are not limited to diameters at or about 100 μm, 95 μm, 90 μm, 85 μm, 80 μm, 75 μm, 70 μm, 65 μm, 60 μm, 55 μm, 50 μm, 45 μm, 40 μm, 35 μm, 30 μm, 25 μm, 20 μm, 15 μm, 10 μm, 9 μm, 8 μm, 7 μm, 6 μm, 5 μm, 4 μm, 3 μm, 2 μm, 1 μm or less. In some particular embodiments, the diameters may be at or about 44 μm, 32 μm, 8 μm, 4 μm, or 1.5 μm. Various heights are shown in the Examples and include but are not limited to heights at or about 100 μm, 95 μm, 90 μm, 85 μm, 80 μm, 75 μm, 70 μm, 65 μm, 60 μm, 55 μm, 50 μm, 45 μm, 40 μm, 35 μm, 30 μm, 25 μm, 20 μm, 15 μm, 10 μm, 9 μm, 8 μm, 7 μm, 6 μm, 5 μm, 4 μm, 3 μm, 2 μm, 1 μm or less. In some particular embodiments, the heights may be at or about 55 μm, 48 μm, 32 μm, 12 μm, 8 μm, 6 μm, 4 μm, 2.25 μm, 1.5 μm, or less. Various embodiments of the invention contemplate the combination of any of these diameters with any of these heights. In still other embodiments, the reaction well dimensions may be (diameter in μm by height in μm) 44 by 55, 32 by 32, 32 by 48, 8 by 8, 8 by 12, 4 by 4, 4 by 6, 1.5 by 1.5, or 1.5 by 2.25.
(422) The reaction well volume may range (between arrays, and preferably not within a single array) based on the well dimensions. This volume may be at or about 100 picoliter (pL), 90, 80, 70, 60, 50, 40, 30, 20, 10, or fewer pL. In important embodiments, the well volume is less than 1 pL, including equal to or less than 0.5 pL, equal to or less than 0.1 pL, equal to or less than 0.05 pL, equal to or less than 0.01 pL, equal to or less than 0.005 pL, or equal to or less than 0.001 pL. The volume may be 0.001 to 0.9 pL, 0.001 to 0.5 pL, 0.001 to 0.1 pL, 0.001 to 0.05 pL, or 0.005 to 0.05 pL. In particular embodiments, the well volume is 75 pL, 34 pL, 23 pL, 0.54 pL, 0.36 pL, 0.07 pL, 0.045 pL, 0.0024 pL, or 0.004 pL. The plurality of templates in each microwell may be introduced into the microwells (e.g., via a nucleic acid loaded bead), or it may be generated in the microwell itself. A plurality is defined herein as at least two, and in the context of template nucleic acids in a microwell or on a nucleic acid loaded bead includes tens, hundreds, thousands, ten thousands, hundred thousands, millions, or more copies of the template nucleic acid. The limit on the number of copies will depend on a number of variables including the number of binding sites for template nucleic acids (e.g., on the beads or on the walls of the microwells), the size of the beads, the length of the template nucleic acid, the extent of the amplification reaction used to generate the plurality, and the like. It is generally preferred to have as many copies of a given template per well in order to increase signal to noise ratio as much as possible. Amplification and conjugation of nucleic acids to solid supports such as beads may be accomplished in a number of ways, including but not limited to emulsion PCR (i.e., water in oil emulsion amplification) as described by Margulies et al. Nature 2005 437(15):376-380 and accompanying supplemental materials. In some embodiments, the amplification is a representative amplification. A representative amplification is an amplification that does not alter the relative representation of any nucleic acid species. The wells generally also include sequencing primers, polymerases and other substrates or catalysts necessary for the synthesis reaction.
(423) The degree of saturation of any capture (i.e., sequencing) bead with template nucleic acid to be sequenced may not be 100%. In some embodiments, a saturation level of 10%-100% exists. As used herein, the degree of saturation of a capture bead with a template refers to the proportion of sites on the bead that are conjugated to template. In some instances this may be at least 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, or it may be 100%.
(424) It will be understood that the amount of sequencing primers and polymerases may be saturating, above saturating level, or in some instances below saturating levels. As used herein, a saturating level of a sequencing primer or a polymerase is a level at which every template nucleic acid is hybridized to a sequencing primer or bound by a polymerase, respectively. Thus the saturating amount is the number of polymerases or primers that is equal to the number of templates on a single bead. In some embodiments, the level is at greater than this, including at least 2 fold, 3 fold, 4 fold, 5 fold, 10 fold, or more over the level of the template nucleic acid. In other embodiments, the number of polymerases and/or primers may be 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or up to 100% of the number of templates on a single bead in a single well.
(425) Thus, for example, before and/or while in the wells, the template nucleic acids are incubated with a sequencing primer that binds to its complementary sequence located on the 3′ end of the template nucleic acid (i.e., either in the amplification primer sequence or in another adaptor sequence ligated to the 3′ end of the target nucleic acid) and with a polymerase for a time and under conditions that promote hybridization of the primer to its complementary sequence and that promote binding of the polymerase to the template nucleic acid. The primer can be of virtually any sequence provided it is long enough to be unique. The hybridization conditions are such that the primer will hybridize to only its true complement on the 3′ end of the template. Suitable conditions are disclosed in Margulies et al. Nature 2005 437(15):376-380 and accompanying supplemental materials.
(426) As described herein, the template nucleic acids may be engineered such that different templates have identical 5′ ends and identical 3′ ends. In some embodiments, however, the invention contemplates the use of a plurality of template populations, wherein each member of a given plurality shares the same 3′ end but different template populations differ from each other based on their 3′ end sequences. As an example, the invention contemplates in some instances sequencing nucleic acids from more than one subject or source. Nucleic acids from first source may have a first 3′ sequence, nucleic acids from a second source may have a second 3′ sequence, and so on, provided that the first and second 3′ sequences are different. In this respect, the 3′ end, which is typically a unique sequence, can be used as a barcode or identifier to label (or identify) the source of the particular nucleic acid in a given well. Reference can be made to Meyer et al. Nucleic Acids Research 2007 35(15):e97 for a discussion of labeling nucleic acid with barcodes followed by sequencing. In some instances, the sequencing primers (if used) may be hybridized (or annealed, as the terms are used interchangeably herein) to the templates prior to loading (or introducing) the beads to the wells or after such loading.
(427) The 5′ and 3′ ends on every individual template however are preferably different in sequence. In particular, the templates share identical primer binding sequences. This facilitates the use of an identical primer across microwells and also ensures that a similar (or identical) degree of primer hybridization occurs across microwells. Once annealed to complementary primers such as sequencing primers, the templates are in a complex referred to herein as a template/primer hybrid. In this hybrid, one region of the template is double stranded (i.e., where it is bound to its complementary primer) and one region is single stranded. It is this single stranded region that acts as the template for the incorporation of nucleotides to the end of the primer and thus it is also this single stranded region which is ultimately sequenced according to the invention.
(428) Suitable polymerases include but are not limited to DNA polymerase, RNA polymerase, or a subunit thereof, provided it is capable of synthesizing a new nucleic acid strand based on the template and starting from the hybridized primer. An example of a suitable polymerase subunit is the exo-version of the Klenow fragment of E. coli DNA polymerase I which lacks 3′ to 5′ exonuclease activity. Other polymerases include T4 exo-, Therminator, and Bst polymerases. The polymerase may be free in solution (and may be present in wash and dNTP solutions) or it may be bound for example to the beads (or corresponding solid support) or to the walls of the chemFET but preferably not to the ISFET surface itself. In another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of template nucleic acids into a plurality of reaction chambers and each reaction chamber is in contact with a chemFET. Sequencing includes synthesizing a new nucleic acid strand by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting incorporation of the one or more known nucleotide triphosphates. The incorporation is carried out by a DNA polymerase. In some embodiments, the polymerase has high affinity for the template, and has a processivity of 100 bp, 250 bp, 500 bp, 750 bp or even 1000 bp. In some embodiments, the polymerase has high activity and a high catalytic rate in high pH conditions. In some embodiments, the pH of the reaction mixture is about 7 to 10, and preferably about 9. In some embodiments, the enzyme has high activity and a high catalytic rate in low concentrations of dNTPs. In some embodiments the dNTP concentration is 50 μM, 40 μM, 30 μM, 20 μM, 10 μM, 5 μM, and preferably 20 μM or less. In some embodiments the enzyme has high activity a catalytic rate where the reaction mixture is of low ionic strength. In some embodiments the concentration of MgCl.sub.2 or MnCl.sub.2 is 100 μM or less. In some embodiments the concentration of Tris is 0.5 mM or less.
(429) Various aspects of the invention relate to nucleotide incorporation into a pre-existing nucleic acid (such as for example a sequencing primer) in the presence of low ionic strength solutions or environments. As used herein, a low ionic strength solution is a solution that has an ion concentration equal to or less than 3 mM total ion concentration. An example of a low ionic strength solution is a 0.5 mM Tris-HCl solution having less than 100 μM MgCl.sub.2. The MgCl.sub.2 concentration may be equal to or less than 90 μM, equal to or less than 80 μM, equal to or less than 70 μM, equal to or less than 60 μM, equal to or less than 50 μM, equal to or less than 40 μM, equal to or less than 30 μM, equal to or less than 20 μM, equal to or less than 10 μM, or less. The polymerase in some embodiments therefore must be capable of functioning at low ionic strength. As described in the Examples, polymerases preferably will be able to recognize and bind to the primer/template hybrid, have sufficient affinity for the hybrid resulting for example in increased processivity of the polymerase, and extend the primer, in low ionic strength conditions. Example 2 demonstrates these various properties for T4 exo-, Therminator, Bst, and Klenow exo-polymerases. The Example further demonstrates that these polymerases are able to incorporate nucleotides into and thereby extend a primer in the context of a silicon well in contact with a chemFET.
(430) In some embodiments, the polymerase must also be capable of incorporating nucleotides in low concentrations of dNTPs. In some embodiments, a low concentration of dNTPs is a concentration that is less 20 μM total dNTP. It will be understood by those of ordinary skill in the art that this concentration relates to the concentration of each dNTP added to the reaction wells, in those embodiments contemplating sequential rather than concurrent introduction of dNTP to a reaction well. In other embodiments, the dNTP concentration (either per dNTP or in some instances total dNTP present in a reaction well at any time during a reaction) is equal to or less than 19 μM, equal to or less than 18 μM, equal to or less than 17 μM, equal to or less than 16 μM, equal to or less than 15 μM, equal to or less than 14 μM, equal to or less than 13 μM, equal to or less than 12 μM, equal to or less than 11 μM, equal to or less than 10 μM, equal to or less than 9 μM, equal to or less than 8 μM, equal to or less than 7 μM, equal to or less than 6 μM, equal to or less than 5 μM, equal to or less than 4 μM, equal to or less than 3 μM, equal to or less than 2 μM, equal to or less than 1 μM, or less. In some embodiments, the dNTP concentration is between 0-20 μM, between 1-19 μM, or between 5-15 μM.
(431) Some embodiments of the invention require that the polymerase have sufficient processivity. As used herein, processivity is the ability of a polymerase to remain bound to a single primer/template hybrid. As used herein, it is measured by the number of nucleotides that a polymerase incorporates into a nucleic acid (such as a sequencing primer) prior to disassociation of the polymerase from the primer/template hybrid. In some embodiments, the polymerase has a processivity of at least 100 nucleotides, although in other embodiments it has a processivity of at least 200 nucleotides, at least 300 nucleotides, at least 400 nucleotides, or at least 500 nucleotides. It will be understood those of ordinary skill in the art that the higher the processivity of the polymerase, the more nucleotides that can be incorporated prior to disassociation, and therefore the longer the nucleic acid that can be sequenced. In other words, polymerases having low processivity (e.g., a polymerase that dissociates from the hybrid after five incorporations will only provide sequence of 5 nucleotides in length, according to some aspects of the invention.
(432) Still other aspects and embodiments of the invention contemplate performing nucleotide incorporation at a high pH (i.e., a pH that is greater than 7). Thus, some reactions are carried out at a pH equal to or greater than 7.5, equal to or greater than 8, equal to or greater than 8.5, equal to or greater than 9, equal to or greater than 9.5, equal to or greater than 10, or equal to or greater than 11. The polymerase may be one that incorporates nucleotides into a nucleic acid (such as a sequencing primer) at a pH of 7-11, 7.5-10.5, 8-10, 8.5-9.5, or at about 9.
(433) It will be understood based on the teachings herein that some aspects and embodiments of the invention will employ polymerases demonstrating high activity in low ionic strength solutions, and/or low dNTP concentrations, and/or high pH. Suitable rates of incorporation will vary depending on the embodiment. In some embodiments, nucleotide incorporation occurs at a rate faster than the diffusion of reagents within a reaction chamber. For example, in one embodiment, nucleotide incorporation occurs at a rate that exceeds the rate of diffusion of dNTP in the reaction well and preferably of dNTP towards the bottom or the reaction well (i.e., to the sensor). In these instances, the newly introduced dNTP do not diffuse towards the sensor until the incorporation reactions are complete. At that time, the remaining dNTP (i.e., the unincorporated dNTP) diffuse past the hybrids (or the beads, as the case may be) and towards the sensor bottom. In these instances, until the incorporation reactions are complete, the dNTP diffusing in the direction of the bead and the well bottom are used up before even reaching the well bottom.
(434) The rate at which a polymerase incorporates nucleotides will vary depending on the particular application, although generally faster rates of incorporation are preferable. The rate of incorporation of a single polymerase to a single nucleic acid (such as a sequencing primer) in some instances is faster than the rate of diffusion of the unincorporated nucleotide triphosphates in the well, and particularly their diffusion towards the bottom of the well. The rate of “sequencing” will depend on the number of arrays on chip, the size of the wells, the temperature and conditions at which the reactions are run, etc. As discussed elsewhere herein, the diffusion rate can also be manipulated by effectively impeding the movement of unincorporated nucleotide triphosphates. This can be done for example by increasing the viscosity in the reaction chamber. This can be accomplished by adding to the reaction chamber a viscosity increasing agent such as but not limited to polyethylene glycol (PEG) or PEA. Diffusion can also be delayed by the use of packing beads in the reaction wells. In some embodiments of the invention, the time for a 4 nucleotide cycle may be 50-100 seconds, 60-90 seconds, or about 70 seconds. In other embodiments, this cycle time can be equal to or less than 70 seconds, including equal to or less than 60 seconds, equal to or less than 50 seconds, equal to or less than 40 seconds, or equal to or less than 30 seconds. A read length of about 400 bases may take on the order of 30 minutes, 60 minutes, 1.5 hours, 2 hours, 2.5 hours, 3 hours, 3.5 hours, 4 hours, 4.5 hours, or in some instance 5 or more hours. These times are sufficient for the sequencing of megabases, and more preferably gigabases of sequence, with greater amounts of sequence being attainable through the use of denser arrays (i.e., arrays with greater numbers of reaction wells and FETs) and/or the simultaneous use of multiple arrays.
(435) Table 2 provides estimates for the rates sequencing based on various array, chip and system configurations contemplated herein. It is to be understood that the invention contemplates even denser arrays than those shown in Table 2. These denser arrays can be characterized as 90 nm CMOS with a pitch of 1.4 μm and a well size of 1 μm which may be used with 0.7 μm beads, or 65 nm CMOS with a pitch of 1 μm and a well size of 0.5 μm which may be used with 0.3 μm beads, or 45 μm CMOS with a pitch of 0.7 μm and a well size of 0.3 μm which can be used with 0.2 μm beads.
(436) TABLE-US-00003 TABLE 2 Reaction Parameters and Read Rates. chip type A B C D E pixel/CMOS 2.8 um/0.18 um 5.1 um/0.35 um 5.1 um/0.35 um 9 um/3.5 um 9 um/0.35 um chip size 17.5 × 17.5 17.5 × 17.5 12 × 12 17.5 × 17.5 12 × 12 # possible reads 27,800,000 7,220,000 2,950,000 2,320,000 1,060,000 read length 400 400 400 400 400 (assumption) # chips/board 4 4 4 4 4 bead load 0.80 0.80 0.80 0.80 0.80 efficiency # yielded Gbp* 35.6 9.2 3.8 3.0 1.4 per run # times HG** 11.86 3.08 1.26 0.99 0.45 (3 Gbp/HG) *Gbp is gigabases **HG is human genome
(437) The template nucleic acid is also contacted with other reagents and/or cofactors including but not limited to buffer, detergent, reducing agents such as dithiothrietol (DTT, Cleland's reagent), single stranded binding proteins, and the like before and/or while in the well. In one embodiment, the polymerase comprises one or more single stranded binding proteins (e.g., the polymerase may be one that is engineered to include one or more single stranded binding proteins). In one embodiment, the template nucleic acid is contacted with the primer and the polymerase prior to its introduction into the flow chamber and wells thereof.
(438) In one embodiment, the nucleic acid loaded beads are introduced into the flow chamber and ultimately the wells situated above the chemFET array. The method requires that each well in the flow chamber contain only one nucleic acid loaded bead since the presence of two beads per well will yield unusable sequencing information derived from two different template nucleic acids. The Examples provides a brief description of an exemplary bead loading protocol in the context of magnetic beads. It is to be understood that a similar approach could be used to load other bead types. The protocol has been demonstrated to reduce the likelihood and incidence of trapped air in the wells of the flow chamber, uniformly distribute nucleic acid loaded beads in the totality of wells of the flow chamber, and avoid the presence and/or accumulation of excess beads in the flow chamber. In some embodiments, the microwell array may be analyzed to determine the degree of loading of beads into the microwells, and in some instances to identify those microwells having beads and those lacking beads. The ability to know which microwells lack beads provides another internal control for the sequencing reaction. The presence or absence of a bead in a well can be determined by standard microscopy or by the sensor itself.
(439) The percentage of occupied wells in the well array may vary depending on the methods being performed. If the method is aimed at extracting maximum sequence data in the shortest time possible, then higher occupancy is desirable. If speed and throughout is not as critical, then lower occupancy may be tolerated. Therefore depending on the embodiment, suitable occupancy percentages may be at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, at least 99%, or 100% of the wells. As used herein, occupancy refers to the presence of one nucleic acid loaded bead in a well and the percentage occupancy refers to the proportion of total wells in an array that are occupied by a single bead. Wells that are occupied by more than one bead cannot be used in the analyses contemplated by the invention.
(440) Ultimately a homogeneous population of template nucleic acids is placed into each of a plurality of wells, each well situated over and thus corresponding to at least one sensor. As discussed above, preferably the well contains at least 10, at least 100, at least 1000, at least 10.sup.4, at least 10.sup.5, at least 10.sup.6, or more copies of an identical template nucleic acid. Identical template nucleic acids means that the templates are identical in sequence. Most and preferably all the template nucleic acids within a well are uniformly hybridized to a primer. Uniform hybridization of the template nucleic acids to the primers means that the primer hybridizes to the template at the same location (i.e., the sequence along the template that is complementary to the primer) as every other template/primer hybrid in the well. The uniform positioning of the primer on every template allows the co-ordinated synthesis of all new nucleic acid strands within a well, thereby resulting in a greater signal-to-noise ratio.
(441) In some embodiments, nucleotides are then added in flow, or by any other suitable method, in sequential order to the flow chamber and thus the wells. The nucleotides can be added in any order provided it is known and for the sake of simplicity kept constant throughout a run.
(442) As mentioned herein, in some embodiments, the method relies on detection of PPi and unincorporated nucleotides. In other embodiments, unincorporated nucleotides are not detected. This can be accomplished for example by adding ATP to the wash buffer so that dNTPs flowing into a well displace ATP from the well. The ATP matches the ionic strength of the dNTPs entering the wells and it also has a similar diffusion profile as those dNTPs. In this way, influx and efflux of dNTPs during the sequencing reaction do not interfere with measurements at the chemFET. The concentration of ATP used is on the order of the concentration of dNTP used.
(443) In some embodiments, particularly those in which a low ionic strength environment is used, the dNTP and/or the polymerase will be pre-incubated with divalent cation such as but not limited to Mg.sup.2+ (for example in the form of MgCl.sub.2) or Mn.sup.2+ (for example in the form of MnCl.sub.2). Other divalent cations can also be used including but not limited to Ca.sup.2+, Co.sup.2+. This pre-incubation (and thus “pre-loading” of the dNTP and/or the polymerase can ensure that the polymerase is exposed to a sufficient amount of divalent cation for proper and necessary functioning even if it is present in a low ionic strength environment. Pre-incubation may occur for 1-60 minutes, 5-45 minutes, or 10-30 minutes, depending on the embodiment, although the invention is not limited to these time ranges.
(444) A sequencing cycle may therefore proceed as follows washing of the flow chamber (and wells) with wash buffer (optionally containing ATP), introduction of a first dNTP species (e.g., dATP) into the flow chamber (and wells), release and detection of PPi and then unincorporated nucleotides (if incorporation occurred) or detection of solely unincorporated nucleotides (if incorporation did not occur) (by any of the mechanisms described herein), washing of the flow chamber (and wells) with wash buffer, washing of the flow chamber (and wells) with wash buffer containing apyrase (to remove as many of the unincorporated nucleotides as possible prior to the flow through of the next dNTP, washing of the flow chamber (and wells) with wash buffer, and introduction of a second dNTP species. This process is continued until all 4 dNTP (i.e., dATP, dCTP, dGTP and dTTP) have been flowed through the chamber and allowed to incorporate into the newly synthesized strands. This 4-nucleotide cycle may be repeated any number of times including but not limited to 10, 25, 50, 100, 200 or more times. The number of cycles will be governed by the length of the template being sequenced and the need to replenish reaction reagents, in particular the dNTP stocks and wash buffers.
(445) As part of the sequencing reaction, a dNTP will be ligated to (or “incorporated into” as used herein) the 3′ of the newly synthesized strand (or the 3′ end of the sequencing primer in the case of the first incorporated dNTP) if its complementary nucleotide is present at that same location on the template nucleic acid. Incorporation of the introduced dNTP (and concomitant release of PPi) therefore indicates the identity of the corresponding nucleotide in the template nucleic acid. If no dNTP has been incorporated, then only a single ion pulse is detected corresponding to unincorporated nucleotides at the surface of the chemFET. One can therefore conclude that the complementary nucleotide was not present in the template at that location. If the introduced dNTP has been incorporated into the newly synthesized strand, then the chemFET will detect a first ion pulse corresponding to the PPi and a second pulse corresponding to the unincorporated nucleotides. It is further possible to quantitate the number of dNTP incorporated by for example measuring the time difference between the time for of the first (PPi) ion pulse (i.e., to) and the second (unincorporated nucleotides) ion pulse (i.e., t.sub.1). The greater the time difference, the greater the number of nucleotides incorporated. The result is that no sequence information is lost through the sequencing of a homopolymer stretch (e.g., poly A, poly T, poly C, or poly G) in the template.
(446) Apyrase is an enzyme that degrades residual unincorporated nucleotides converting them into monophosphate and releasing inorganic phosphate in the process. It is useful for degrading dNTPs that are not incorporated and/or that are in excess. It is important that excess and/or unincorporated dNTP be washed away from all wells after measurements are complete and before introduction of the subsequent dNTP. Accordingly, addition of apyrase between the introduction of different dNTPs is useful to remove unincorporated dNTPs that would otherwise obscure the sequencing data.
(447) Additional sequencing reaction reagents such as those described above may be introduced throughout the reaction, although in some cases this may not be necessary. For example additional polymerase, DTT, SBB and the like may be added if necessary.
(448) The invention therefore contemplates performing a plurality of different sequencing reactions simultaneously. A plurality of identical sequencing reactions is occurring in each occupied well simultaneously. It is this simultaneous and identical incorporation of dNTP within each well that increases the signal to noise ratio. By performing sequencing reactions in a plurality of wells simultaneously, a plurality of different nucleic acids are simultaneously sequenced.
(449) The methods aim to maximize complete incorporation across all microwells for any given dNTP, reduce or decrease the number of unincorporated dNTPs that remain in the wells after signal detection is complete, and achieve as a high a signal to noise ratio as possible.
(450) The sequencing reaction can be run at a range of temperatures. Typically, the reaction is run in the range of 30-60° C., 35-55° C., or 40-45° C. It is preferable to run the reaction at temperatures that prevent formation of secondary structure in the nucleic acid. However this must be balanced with the binding of the primer (and the newly synthesized strand) to the template nucleic acid and the reduced half-life of apyrase at higher temperatures. A suitable temperature is about 41° C. The solutions including the wash buffers and the dNTP solutions are generally warmed to these temperatures in order not to alter the temperature in the wells. The wash buffer containing apyrase however is preferably maintained at a lower temperature in order to extend the half-life of the enzyme. Typically, this solution is maintained at about 4-15° C., and more preferably 4-10° C. It will be appreciated that sequencing in a low ionic concentration can lower the Tm of the sequencing primer, possibly below that optimal for the polymerase. The invention therefore contemplates that sequencing primers may comprise nucleotide modifications such as modified bases, linked nucleic acid (LNA) monomers, peptide nucleic acid (PNA) monomers, or alternatively minor groove binders could be used to elevate the primer Tm into a suitable range. The nucleotide incorporation reaction can occur very rapidly. As a result, it may be desirable in some instances to slow the reaction down or to slow the diffusion of analytes in the well in order to ensure maximal data capture during the reaction. The diffusion of reagents and/or byproducts can be slowed down in a number of ways including but not limited to addition of packing beads in the wells, and/or the use of polymers such as polyethylene glycol in the wells (e.g., PEG attached to the capture beads and/or to packing beads). The packing beads also tend to increase the concentration of reagents and/or byproducts at the chemFET surface, thereby increasing the potential for signal. The presence of packing beads generally allows a greater time to sample (e.g., by 2- or 4-fold).
(451) Data capture rates can vary and be for example anywhere from 10-100 frames per second and the choice of which rate to use will be dictated at least in part by the well size and the presence of packing beads or other diffusion limiting techniques. Smaller well sizes generally require faster data capture rates.
(452) In some aspects of the invention that are flow-based and where the top face of the well is open and in communication with fluid over the entirety of the chip, it is important to detect the released PPi or other byproduct (e.g., H.sup.+) prior to its diffusion out of the well. Diffusion of reaction byproducts out of the well will lead to false negatives (because the byproduct is not detected in that well) and potential false positives in adjacent or downstream wells (where the byproduct may be detected), and thus should be avoided. Packing beads and/or polymers such as PEG may also help reduce the degree of diffusion and/or cross-talk between wells.
(453) Thus in some embodiments packing beads are used in addition to the nucleic acid-loaded beads. The packing beads may be magnetic (including superparamagnetic) but they are not so limited. In some embodiments the packing beads and the capture beads are made of the same material (e.g., both are magnetic, both are polystyrene, etc.), while in other embodiment they are made of different materials (e.g., the packing beads are polystyrene and the capture beads are magnetic). The packing beads are generally smaller than the capture beads. The difference in size may be vary and may be 5-fold, 10-fold, 15-fold, 20-fold or more. As an example, 0.35 μm diameter packing beads can be used with 5.91 μm capture beads. Such packing beads are commercially available from sources such as Bang Labs. The placement of the packing beads relative to the capture bead may vary. For example, the packing beads may surround the capture bead and thereby prevent the capture bead from contacting the chemFET surface. As another example, the packing beads may be loaded into the wells following the capture beads in which case the capture bead is in contact with the chemFET surface. The presence of packing beads between the capture bead and the chemFET surface slows the diffusion of the sequencing byproducts such as PPi, thereby facilitating data capture.
(454) The invention further contemplates the use of packing beads or modifications to the chemFET surface (as described herein) to prevent contact and thus interference of the chemFET surface with the template nucleic acids bound to the capture beads. A layer of packing beads that is 0.1-0.5 μm in depth or height would preclude this interaction.
(455) The sequencing reaction may be preceded by an analysis of the arrays to determine the location of beads. It has been found that in the absence of flow the background signal (i.e., noise) is less than or equal to about 0.25 mV, but that in the presence of DNA-loaded capture beads that signal increases to about 1.0 mV=/−0.5 mV. This increase is sufficient to allow one to determine wells with beads.
(456) The invention further contemplates kits comprising the various reagents necessary to perform a sequencing reaction and instructions of use according to the methods set forth herein.
(457) One preferred kit comprises one or more containers housing wash buffer, one or more containers each containing one of the following reagents: dATP buffer, dCTP buffer, dGTP buffer or dTTP buffer, dATP, dCTP, dGTP and dTTP stocks, apyrase, SSB, polymerase, packing beads, and optionally pyrophosphatase. Importantly the kits may comprise only naturally occurring dNTPs. The kits may also comprise one or more wash buffers comprising components as described in the Examples, but are not so limited. The kits may also comprise instructions for use including diagrams that demonstrate the methods of the invention.
(458) It is to be understood that interactions between receptors and ligands or between two members of a binding pair or between components of a molecular complex can also be detected using the chemFET arrays. Examples of such interactions include hybridization of nucleic acids to each other, protein-nucleic acid binding, protein-protein binding, enzyme-substrate binding, enzyme-inhibitor binding, antigen-antibody binding, and the like. Any binding or hybridization event that causes a change of the semiconductor charge density at the FET interface and thus changes the current that flows from the source to the drain of the sensors described herein can be detected according to the invention.
(459) In these embodiments, the passivation layer (or possibly an intermediate layer coated onto the passivation layer) is functionalized with nucleic acids (e.g., DNA, RNA, miRNA, cDNA, and the like), antigens (which can be of any nature), proteins (e.g., enzymes, cofactors, antibodies, antibody fragments, and the like), and the like. Conjugation of these entities to the passivation layer can be direct or indirect (e.g., using bifunctional linkers that bind to both the passivation layer reactive group and the entity to be bound).
(460) As an example, reaction groups such as amine or thiol groups may be added to a nucleic acid at any nucleotide during synthesis to provide a point of attachment for a bifunctional linker. As another example, the nucleic acid may be synthesized by incorporating conjugation-competent reagents such as Uni-Link AminoModifier, 3′-DMT-C6-Amine-ON CPG, AminoModifier II, N-TFA-C6-AminoModifier, C6-ThiolModifier, C6-Disulfide Phosphoramidite and C6-Disulfide CPG (Clontech, Palo Alto, Calif.). Other methods for attaching nucleic acids are discussed below.
(461) In one aspect of the invention, the chemFET arrays are provided in combination with nucleic acid arrays. Nucleic acids in the form of short nucleic acids (e.g., oligonucleotides) or longer nucleic acids (e.g., full length cDNAs) can be provided on chemFET surfaces of the arrays described herein. Nucleic acid arrays generally comprise a plurality of physically defined regions on a planar surface (e.g., “spots”) each of which has conjugated to it one and more preferably more nucleic acids. The nucleic acids are usually single stranded. The nucleic acids conjugated to a given spot are usually identical. In the context of an oligonucleotide array, these nucleic acids may be on the order of less 100 nucleotides in length (including about 10, 20, 25, 30, 40, 50, 60, 70, 80, 90 or 100 nucleotides in length). If the arrays are used to detect certain genes (including mutations in such genes or expression levels of such genes), then the array may include a number of spots each of which contains oligonucleotides that span a defined and potentially different sequence of the gene. These spots are then located across the planar surface in order to exclude position related effects in the hybridization and readout means of the array.
(462) The arrays are contacted with a sample being tested. The sample may be a genomic DNA sample, a cDNA sample from a cell, a tissue or a mass (e.g., a tumor), a population of cells that are grown on the array, potentially in a two dimensional array that corresponds to the underlying sensor array, and the like. Such arrays are therefore useful for determining presence and/or level of a particular gene or of its expression, detecting mutations within particular genes (such as but not limited to deletions, additions, substitutions, including single nucleotide polymorphisms), and the like.
(463) The binding or hybridization of the sample nucleic acids and the immobilized nucleic acids is generally performed under stringent hybridization conditions as that term is understood in the art. (See for example Sambrook et al. “Maniatis”.) Examples of relevant conditions include (in order of increasing stringency): incubation temperatures of 25° C., 37° C., 50° C. and 68° C.; buffer concentrations of 10×SSC, 6×SSC, 4×SSC, 1×SSC, 0.1×SSC (where SSC is 0.15 M NaCl and 15 mM citrate buffer) and their equivalents using other buffer systems; formamide concentrations of 0%, 25%, 50%, and 75%; incubation times from 5 minutes to 24 hours; 1, 2, or more washing steps; wash incubation times of 1, 2, or 15 minutes; and wash solutions of 6×SSC, 1×SSC, 0.1×SSC, or deionized water. By way of example hybridization may be performed at 50% formamide and 4×SSC followed by washes of 2×SSC/formamide at 50° C. and with 1×SSC.
(464) Nucleic acid arrays include those in which already formed nucleic acids such as cDNAs are deposited (or “spotted”) on the array in a specific location. Nucleic acids can be spotted onto a surface by piezoelectrically deposition, UV crosslinking of nucleic acids to polymer layers such as but not limited to poly-L-lysine or polypyrrole, direct conjugation to silicon coated SiO.sub.2 as described in published US patent application 2003/0186262, direct conjugation to a silanised chemFET surface (e.g., a surface treated with 3-aminopropyltriethoxysilane (APTES) as described by Uslu et al. Biosensors and Bioelectronics 2004, 19:1723-1731, for example.
(465) Nucleic acid arrays also include those in which nucleic acids (such as oligonucleotides of known sequence) are synthesized directly on the array. Nucleic acids can be synthesized on arrays using art-recognized techniques such as but not limited to printing with fine-pointed pins onto glass slides, photolithography using pre-made masks, photolithography using dynamic micromirror devices (such as DLP mirrors), ink-jet printing, or electrochemistry on microelectrode arrays. Reference can also be made to Nuwaysir et al. 2002 “Gene expression analysis using oligonucleotide arrays produced by maskless photolithography.”. Genome Res 12: 1749-1755. Commercial sources of this latter type of array include Agilent, Affymetrix, and NimbleGen.
(466) Thus the chemFET passivation layer may be coated with an intermediate layer of reactive molecules (and therefore reactive groups) to which the nucleic acids are bound and/or from which they are synthesized.
(467) The invention contemplates combining such nucleic acid arrays with the chemFET arrays and particularly the “large scale” chemFET arrays described herein. The chemFET/nucleic acid array can be used in a variety of applications, some of which will not require the wells (or microwells or reaction chambers, as they are interchangeably referred to herein). Since analyses may still be carried out in flow, including in a “closed” system (i.e., where the flow of reagents and wash solutions and the like is automated), there will be one or more flow chambers situated above and in contact with the array. The use of multiple flow chambers allows multiple samples (or nucleic acid libraries) to be analyzed simultaneously. There may be 2, 3, 4, 5, 6, 7, 8, 9, 10 or more flow chambers. This configuration applies equally to other biological arrays including those discussed herein such as protein arrays, antibody arrays, enzyme arrays, chemical arrays, and the like.
(468) Since the binding event between binding partners or between components of a complex is detected electronically via the underlying FET, such assays may be carried out without the need to manipulate (e.g., extrinsically label) the sample being assayed. This is advantageous since such manipulation invariably results in loss of sample and generally requires increased time and work up. In addition, the present method allows binding interactions to be studied in real time.
(469) Protein arrays used in combination with the chemFET arrays of the invention are also contemplated. Protein arrays comprise proteins or peptides or other amino acid comprising biological moiety bound to a planar surface in an organized and predetermined manner. Such proteins include but are not limited to enzymes, antibodies and antibody fragments or antibody mimics (e.g., single chain antibodies).
(470) In one embodiment, a protein array may comprise a plurality of different proteins (or other amino acid containing biological moieties). Each protein, and preferably a plurality of proteins, is present in a predetermined region or “cell” of the array. The regions (or cells) are aligned with the sensors in the sensor array such that there is one sensor for each region (or cell). The plurality of proteins in a single region (or cell) may vary depending on the size of the protein and the size of the region (or cell) and may be but is not limited to at least 10, 50, 100, 500, 10.sup.3, 10.sup.4 or more. The array itself may have any number of cells, including but not limited to at least 10, 10.sup.2, 10.sup.3, 10.sup.4, 10.sup.5, 10.sup.6, 10.sup.7, or more. In one application, the array is exposed to a sample that is known to contain or is suspected of containing an analyte that binds to the protein. The analyte may be a substrate or an inhibitor if the protein is an enzyme. The analyte may be any molecule that binds to the protein including another protein, a nucleic acid, a chemical species (whether synthetic or naturally occurring), and the like.
(471) It is to be understood that, like the nucleic acid arrays contemplated herein, the readout from the protein arrays will be a change in current through the chemFET and thus no additional step of labeling and/or label detection is required in these array methods.
(472) In another embodiment, the protein array may comprise a plurality of identical proteins (or other amino acid containing biological moieties). The identical proteins may be uniformly distributed on a planar surface or they may be organized into discrete regions (or cells) on that surface. In these latter embodiments, the regions (or cells) are aligned with the sensors in the sensor array such that there is one sensor for each region (or cell).
(473) The proteins may be synthesised off-chip, then purified and attached to the array. Alternatively they can be synthesised on-chip, similarly to the nucleic acids discussed above. Synthesis of proteins using cell-free DNA expression or chemical synthesis is amenable to on-chip synthesis. Using cell-free DNA expression, proteins are attached to the solid support once synthesized. Alternatively, proteins may be chemically synthesized on the solid support using solid phase peptide synthesis. Selective deprotection is carried out through lithographic methods or by SPOT-synthesis. Reference can be made to at least MacBeath and Schreiber, Science, 2000, 289:1760-1763, or Jones et al. Nature, 2006, 439:168-174. Reference can also be made to U.S. Pat. No. 6,919,211 to Fodor et al.
(474) Chemical compound microarrays in combination with chemFET arrays are also envisioned. Chemical compound microarrays can be made by covalently immobilizing the compounds (e.g., organic compounds) on the solid surface with diverse linking techniques (may be referred to in the literature as “small molecule microarray”), by spotting and drying compounds (e.g., organic compounds) on the solid surface without immobilization (may be referred to in the literature as “micro arrayed compound screening (μARCS)”), or by spotting organic compounds in a homogenous solution without immobilization and drying effect (commercialized as DiscoveryDot™ technology by Reaction Biology Corporation).
(475) Tissue microarrays in combination with chemFET arrays are further contemplated by the invention. Tissue microarrays are discussed in greater detail in Battifora Lab Invest 1986, 55:244-248; Battifora and Mehta Lab Invest 1990, 63:722-724; and Kononen et al. Nat Med 1998, 4:844-847.
(476) The configurations of the chemFET arrays and the biological or chemical arrays are similar in each instance and the discussion of one combination array will apply to others described herein or otherwise known in the art.
(477) In yet another aspect, the invention contemplates analysis of cell cultures (e.g., two-dimensional cells cultures) (see for example Baumann et al. Sensors and Actuators B 55 1999 77:89), and tissue sections placed in contact with the chemFET array. As an example, a brain section may be placed in contact with the chemFET array of the invention and changes in the section may be detected either in the presence or absence of stimulation such as but not limited to neurotoxins and the like. Transduction of neural processes and/or stimulation can thereby be analyzed. In these embodiments, the chemFETs may operate by detecting calcium and/or potassium fluxes via the passivation layer itself or via receptors for these ions that are coated onto the passivation layer.
(478) In yet another aspect, the invention contemplates the use of chemFET arrays, functionalized as described herein or in another manner, for use in vivo. Such an array may be introduced into a subject (e.g., in the brain or other region that is subject to ion flux) and then analyzed for changes based on the status of the subject.
(479) Methods for attaching nucleic acids, proteins, molecules, and the like to solid supports, particularly in the context of an array, have been described in the art. See for example Lipshutz et al. Nat. Genet. (supplement) 1999 21:20-24; Li et al. Proc. Natl. Acad. Sci., 2001, 98:31-36; Lockhart et al. Nat. Biotechnol. 1996 14:1675-1680; Wodicka et al. Nat. Biotechnol. 1997 15: 1359-1367; Chen et al. Journal of Biomedical Optics 1997 2:364; Duggan et al. Nat Genet 1991 21(1 Suppl):10-4; Marton et al. Nat Med. 1998 4(11):1293-301; Kononen et al. Nat Med 1998 4(7):844-847; MacBeath et al., Science 2000 289(5485):1760-1763; Haab et al. Genome Biology 2001 2(2); Pollack et al. Nat Genet 1999 23(1):41-6; Wang D G et al. Science 1998 280(5366):1077-82; Fodor et al. Science 1991 251:767-773; Fodor et al. Nature 1993 364:555-556; Pease et al. Proc. Natl Acad. Sci. USA 1994 91:5022-5026; Fodor Science 1997 277:393-395; Southern et al. Genomics 1992 13: 1008-1017; Schena et al. Science 1995 270(5235):467-70; Shalon et al. Genome Res 1996 6(7):639-45; Jongsma Proteomics 2006, 6:2650-2655; Sakata, Biosensors and Bioelectronics 2007, 22: 1311-1316.
(480) The systems described herein can be used for sequencing unlabeled biological polymers without optical detection.
(481) In some embodiments, the invention encompasses a sequencing apparatus adapted for sequencing unlabeled biological polymers without optical detection and comprising an array of at least 100 reaction chambers.
(482) Typically, each reaction chamber is capacitively coupled to a chemFET.
(483) Preferably, each reaction chamber is no greater than about 0.39 pL in volume and about 49 μm.sup.2 surface aperture, and more preferably has an aperture no greater than about 16 μm.sup.2 and volume no greater than about 0.064 pL. Preferably, the array has at least 1,000, 10,000, 100,000, or 1,000,000 reaction chambers.
(484) Typically, the reaction chambers comprise microfluidic wells.
(485) Preferably, the apparatus is adapted to sequence at least 10.sup.6 base pairs per hour, more preferably at least 10.sup.7 base pairs per hour, and most preferably at least 10.sup.8 base pairs per hour.
(486) In another embodiment, the invention encompasses a method for sequencing a biological polymer with the above-described apparatus comprising measuring time of incorporation of individual monomers into an elongating polymer.
(487) Typically, the biological polymer is a nucleic acid template and the monomer is a nucleotide. Preferably, the nucleic acid template has 200-700 base pairs. Preferably, the nucleic acid template is amplified prior to determining the sequence.
(488) Typically, the measuring is performed under diffusion limited conditions. The measuring may comprise detecting an electrical change, detecting an ion pulse, or detecting the release of inorganic pyrophosphate (“PPi”). Preferably, the incorporation takes place at an ionic strength no greater than 400 μM, wherein the concentration of Mg.sup.2+ or Mn.sup.2+ or other divalent cation is no greater than 100 μM. Preferably, at least 10.sup.6 base pairs are sequenced per hour, more preferably at least 10.sup.7 base pairs are sequenced per hour, and most preferably at least 10.sup.8 base pairs are sequenced per hour using the above-described method. Preferably the detecting is accomplished by an array of chemFETs as above described.
(489) In some embodiments, the invention encompasses a method for determining a nucleic acid sequence comprising: performing at least 100 sequencing reactions simultaneously; and determining the sequence without the use of labels and without the use of optical detection.
(490) The nucleic acid template used in this and other methods of the invention may be derived from a variety of sources by a variety of methods, all known to those of ordinary skill in the art. Templates may be derived from, but are not limited to, entire genomes of varying complexity, cDNA, mRNA or siRNA samples, or may represent entire populations, as in the various environmental and metabiome sequencing projects. Template nucleic acids may also be generated from specific subsets of nucleic acid populations including but not limited to PCR products, specific exons or regions of interest, or 16S or other diagnostic or identifying genomic regions.
(491) Typically, the required starting material for sequencing is less than 3 μg of nucleic acids. Preferably, the template nucleic acid is amplified prior to sequencing. The template nucleic acid may be optionally bound to one or more beads prior to sequencing.
(492) Typically, the determining step comprises measuring the amount of time it takes for one or more of a first monomer to incorporate into an elongating sequence. Typically, the incorporation of the monomer is measured under diffusion limited conditions. The measuring may comprise detecting an electrical change, detecting an ion pulse, or detecting the release of inorganic pyrophosphate (“PPi”).
(493) Preferably, at least 10.sup.6 base pairs are sequenced per hour, more preferably at least 10.sup.7 base pairs are sequenced per hour, and most preferably at least 10.sup.8 base pairs are sequenced per hour using the above-described method. Thus, the method may be used to sequence an entire human genome within about 24 hours, more preferably within about 20 hours, even more preferably within about 15 hours, even more preferably within about 10 hours, even more preferably within about 5 hours, and most preferably within about 1 hour. These rates may be achieved using multiple ISFET arrays as shown herein, and processing their outputs in parallel.
(494) The following examples are included for purposes of illustration and are not intended to limit the scope of the invention.
Example 1. Bead Preparation
(495) Binding of Single-Stranded Oligonucleotides to Streptavidin-Coated Magnetic Beads.
(496) Single-stranded DNA oligonucleotide templates with a 5′ Dual Biotin tag (HPLC purified), and a 20-base universal primer were ordered from IDT (Integrated DNA Technologies, Coralville, Ind.). Templates were 60 bases in length, and were designed to include 20 bases at the 3′ end that were complementary to the 20-base primer (Table 3, italics). The lyophilized and biotinylated templates and primer were re-suspended in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8) as 40 μM stock solutions and as a 400 μM stock solution, respectively, and stored at −20° C. until use.
(497) For each template, 60 μl of magnetic 5.91 μm (Bangs Laboratories, Inc. Fishers, Ind.) streptavidin-coated beads, stored as an aqueous, buffered suspension (8.57×10.sup.4 beads/μL), at 4° C., were prepared by washing with 120 μl bead wash buffer three times and then incubating with templates 1, 2, 3 and 4 (T1, T2, T3, T4 : Table 3) with biotin on the 5′ end, respectively.
(498) Due to the strong covalent binding affinity of streptavidin for biotin (Kd.sup.˜10-15), these magnetic beads are used to immobilize the templates on a solid support, as described below. The reported binding capacity of these beads for free biotin is 0.650 pmol/μL of bead stock solution. For a small (<100 bases) biotinylated ssDNA template, it was conservatively calculated that 9.1×10.sup.5 templates could be bound per bead. The beads are easily concentrated using simple magnets, as with the Dynal Magnetic Particle Concentrator or MPC-s (Invitrogen, Carlsbad, Calif.). The MPC-s was used in the described experiments.
(499) An MPC-s was used to concentrate the beads for 1 minute between each wash, buffer was then added and the beads were resuspended. Following the third wash the beads were resuspended in 120 μL bead wash buffer plus 1 μl of each template (40 μM). Beads were incubated for 30 minutes with rotation (Labquake Tube Rotator, Barnstead, Dubuque, Iowa). Following the incubation, beads were then washed thee times in 120 μL Annealing Buffer (20 mM Tris-HCl, 5 mM magnesium acetate, pH 7.5), and re-suspended in 60 μL of the same buffer.
(500) TABLE-US-00004 TABLE 3 Sequences for Templates 1, 2, 3, and 4 T1: 5′/52Bio/GCA AGT GCC CTT AGG CTT CAG TTC AAA AGT CCT AAC TGG GCA AGG CAC ACA GGG GAT AGG-3′ (SEQ ID NO: 1) T2: 5′/52Bio/CCA TGT CCC CTT AAG CCC CCC CCA TTC CCC CCT GAA CCC CCA AGG CAC ACA GGG GAT AGG-3′ (SEQ ID NO: 2) T3: 5′/52Bio/AAG CTC AAA AAC GGT AAA AAA AAG CCA AAA AAC TGG AAA ACA AGG CAC ACA GGG GAT AGG-3′ (SEQ ID NO: 3) T4: 5′/52Bio/TTC GAG TTT TTG CCA TTT TTT TTC GGT TTT TTG ACC TTT TCA AGG CAC ACA GGG GAT AGG-3′ (SEQ ID NO: 4)
(501) Annealing of Sequencing Primer. The immobilized templates, bound at the 5′ end to 5.91 μm magnetic beads, are then annealed to a 20-base primer complementary to the 3′ end of the templates (Table 3). A 1.0 μL aliquot of the 400 μM primer stock solution, representing a 20-fold excess of primer to immobilized template, is then added and then the beads plus template are incubated with primer for 15 minutes at 95° C. and the temperature was then slowly lowered to room temperature. The beads were then washed 3 times in 120 μL of 25 mM Tricine buffer (25 mM Tricine, 0.4 mg/ml PVP, 0.1% Tween 20, 8.8 mM Magnesium Acetate; ph 7.8) as described above using the MPC-s. Beads were resuspended in 25 mM Tricine buffer.
(502) Incubation of Hybridized Templates/Primer with DNA Polymerase.
(503) Template and primer hybrids are incubated with polymerase essentially as described by Margulies et al. Nature 2005 437(15):376-380 and accompanying supplemental materials.
(504) Loading of Prepared Test Samples onto the ISFET Sensor Array. The dimensions and density of the ISFET array and the microfluidics positioned thereon may vary depending on the application. A non-limiting example is a 512×512 array. Each grid of such an array (of which there would be 262144) has a single ISFET. Each grid also has a well (or as they may be interchangeably referred to herein as a “microwell”) positioned above it. The well (or microwell) may have any shape including columnar, conical, square, rectangular, and the like. In one exemplary conformation, the wells are square wells having dimensions of 7×7×10 μm. The center-to-center distance between wells is referred to herein as the “pitch”. The pitch may be any distance although it is preferably to have shorter pitches in order to accommodate as large of an array as possible. The pitch may be less than 50 μm, less than 40 μm, less than 30 μm, less than 20 μm, or less than 10 μm. In one embodiment, the pitch is about 9 μm. The entire chamber above the array (within which the wells are situated) may have a volume of equal to or less than about 30 μL, equal to or less than about 20 μL, equal to or less than about 15 μL, or equal to or less than 10 μL. These volumes therefore correspond to the volume of solution within the chamber as well.
(505) Loading of Beads in an ‘Open’ System. Beads with templates 1-4 were loaded on the chip (10 μL of each template). Briefly, an aliquot of each template was added onto the chip using an Eppendorf pipette. A magnet was then used to pull the beads into the wells.
(506) Loading of Beads in a ‘Closed’ System. Both the capture beads the packing beads are loaded using flow. Microliter precision of bead solution volume, as well as positioning of the bead solution through the fluidics connections, is achieved as shown in
(507) The chip comprising the ISFET array and flow cell is seated in the ZIF (zero insertion force) socket of the loading fixture, then attaching a stainless steel capillary to one port of the flow cell and flexible nylon tubing on the other port. Both materials are microfluidic-type fluid paths (e.g., on the order of <0.01″ inner diameter). The bead loading fitting, consisting of the major and minor reservoirs, it attached to the end of the capillary. A common plastic syringe is filled with buffer solution, then connected to the free end of the nylon tubing. The electrical leads protruding from the bottom of the chip are inserted into a socket on the top of a fixture unit (not shown).
(508) The chip comprising the ISFET array and flow cell is seated in a socket such as a ZIF (zero insertion force) socket of the loading fixture, then a stainless steel capillary may be attached to one port of the flow cell and flexible nylon tubing on the other port. Both materials are microfluidic-type fluid paths (e.g., on the order of <0.01″ inner diameter). The bead loading fitting, consisting of the major and minor reservoirs, it attached to the end of the capillary. A common plastic syringe is filled with buffer solution, then connected to the free end of the nylon tubing. The electrical leads protruding from the bottom of the chip are inserted into a socket on the top of a fixture unit (not shown).
(509) It will be appreciated that there will be other ways of drawing the beads into the wells of the flow chamber, including centrifugation or gravity. The invention is not limited in this respect.
(510) DNA Sequencing Using the ISFET Sensor Array in an Open System. A illustrative sequencing reaction can be performed in an ‘open’ system (i.e., the ISFET chip is placed on the platform of the ISFET apparatus and then each nucleotide (5 μL resulting in 6.5 μM each) is manually added in the following order: dATP, dCTP, dGTP and dTT (100 mM stock solutions, Pierce, Milwaukee, Wis.), by pipetting the given nucleotide into the liquid already on the surface of the chip and collecting data from the chip at a rate of 2.5 mHz. This can result in data collection over 7.5 seconds at approximately 18 frames/second. Data may then analyzed using LabView.
(511) Given the sequences of the templates, it is expected that addition of dATP will result in a 4 base extension for template 4. Addition of dCTP will result in a 4 base extension in template 1. Addition of dGTP will cause template 1, 2 and 4 to extend as indicated in Table 4 and addition of dTTP will result in a run-off (extension of all templates as indicated).
(512) Preferably when the method is performed in a non-automated manner (i.e., in the absence of automated flow and reagent introduction), each well contains apyrase in order to degrade the unincorporated dNTPs, or alternatively apyrase is added into each well following the addition and incorporation of each dNTP (e.g., dATP) and prior to the addition of another dNTP (e.g., dTTP). It is to be understood that apyrase can be substituted, in this embodiment or in any other embodiment discussed herein, with another compound (or enzyme) capable of degrading dNTPs.
(513) TABLE-US-00005 TABLE 4 Set-up of experiment and order of nucleotide addition. dATP dCTP dGTP dTTP T1 0 (3:C; 1:A)4 1 Run-off (25) T2 0 0 4 Run-off (26) T3 0 0 0 Run-off (30) T4 4 0 2 Run-off (24)
(514) DNA Sequencing Using Microfluidics on Sensor Chip. Sequencing in the flow regime is an extension of open application of nucleotide reagents for incorporation into DNA. Rather than add the reagents into a bulk solution on the ISFET chip, the reagents are flowed in a sequential manner across the chip surface, extending a single DNA base(s) at a time. The dNTPs are flowed sequentially, beginning with dTTP, then dATP, dCTP, and dGTP. Due to the laminar flow nature of the fluid movement over the chip, diffusion of the nucleotide into the microwells and finally around the nucleic acid loaded bead is the main mechanism for delivery. The flow regime also ensures that the vast majority of nucleotide solution is washed away between applications. This involves rinsing the chip with buffer solution and apyrase solution following every nucleotide flow. The nucleotides and wash solutions are stored in chemical bottles in the system, and are flowed over the chip using a system of fluidic tubing and automated valves. The ISFET chip is activated for sensing chemical products of the DNA extension during nucleotide flow.
Example 2. Low Ionic Strength Effects on Polymerase Activity
(515) Experiments were conducted to test various polymerases for certain functionalities associated with the sequencing aspects and embodiments discussed herein, including in particular low ionic strength nucleotide incorporation reactions. These functionalities include the ability of a polymerase to extend a sequencing primer in low ionic strength, affinity of a polymerase for primer/template nucleic acid hybrids in low ionic strength, rate of primer extension, and ability of polymerase to extend primers in the context of a well situated above a chemFET sensor. Each of these functionalities will be discussed below.
(516) Polymerase Activity in Low Ionic Strength. Polymerase activity can be measured by determining the extent to which a primer is extended along a template. A fluorescent reporter assay was devised to measure polymerase activity (
(517) In particular, the data shown in
(518) Using this assay polymerase activity has been measured in both standard and low ionic strength reaction conditions. A number of polymerases are relatively active in low ionic strength conditions when compared to standard conditions (
(519) Polymerase Affinity in Low Ionic Strength. The affinity of a polymerase for a template/primer hybrid can be measured indirectly with our standard assay (
(520) Rates of Polymerases. The rate at which a polymerase can extend a primer along a template can be measured with a standard assay by measuring the fraction of template extended per unit time. T4 exo-polymerase (3′ to 5′ exonuclease deficient) extends primer in both standard reaction conditions (STD; 10 mM Tris-HCl, 50 mM NaCl, 10 mM MgCl.sub.2, 1 mM dithiothreitol, 20 μM each dNTP, pH 7.9 at 25° C.) and in low ionic strength conditions (LS; 0.5 mM Tris, 80 μM MgCl.sub.2, 20 μM each dNTP, pH 9.0 at 25° C.), but at different rates. Previous studies of T4 demonstrate that T4 can extend a primer along a template at a maximal steady state rate of 4 nucleotides per second in a standard reaction solution (i.e., at standard ionic strength). At this rate, extension in the standard reaction solution is expected to be complete in 8.75 seconds, and indeed extension is complete at about 10 seconds experimentally (
(521) Therminator and Bst polymerases can extend primer in both standard reaction conditions (STD; 20 mM Tris-HCl, 10 mM (NH2).sub.2SO.sub.4, 10 mM KCl, 2 mM MgSO.sub.4, 0.1% Triton X-100, pH 8.8 at 25° C.) and in low ionic strength conditions (LS; 0.5 mM Tris, 100 μM MgCl.sub.2, 20 μM each dNTP, pH 9.0 at 25° C.). The steady state rate in low ionic strength conditions is slower than extension by these polymerases in standard ionic strength conditions, and it is slower than the steady state rate of T4 primer extension in low ionic strength extension solutions.
(522) Primer Extension in Wells of Sensor. Using an assay similar to the standard activity assay (
Example 3. Increasing Speed and Throughput
(523) One advantage of the array architecture used on the ISFET/chemFET chips is that all microwells receive and process chemical materials in parallel and the outputs of the FETs are analyzed in parallel. While a single chip is limited to N wells and the number N will increase over time as semiconductor processes shrink devices ever smaller, no matter how many FETs and wells (also known as reaction chambers) are on a single chip, a system may readily be build with two or more chips in the same machine, operated in parallel; or multiple machines may be operated in parallel and as their FETs provide output, those outputs can be processed in parallel to reduce the time required to sequence a genome or perform another analysis. Parallel computing is well known and there are multiple known ways for computer scientists to configure the processing of the outputs of several chips to expedite the computation of the desired output. Hence, with multiple chips, each having 256 k or more microwells, it is readily achievable to extract from a sample a human genome in a span of one or just a few hours. As chip speeds increase and device sizes shrink, these times may be achievable with a single chip.
(524) Having explained the theory and objective of the ion pulse detection approach, above, it may be helpful to some to provide basic examples of conditions that appear useful, based on simulation, in some circumstances, for achieving detection of DNA sequencing reactions.
(525) The Ions in Solution.
(526) The ionic strength is a factor that is involved in determining the ion pulse response. It is helpful to limit the ions in the reaction buffer. Assuming 0.5 mM Tris and use of to balance the charge, the following may be determined:
(527) TABLE-US-00006 Species Concentration (mM) Tris 0.452 Tris-H.sup.+ 0.048 Cl.sup.− 0.048
(528) Consider the case of 0.5 mM Tris, as other concentrations of Tris (e.g. 5-20 mM) are given by scaling the above concentrations. Assume the pH is constant at 8.8. The other species in the solution include the reactants and products:
(529) TABLE-US-00007 Species Concentration dynamics (μM) dNTP.sup.−−−− 0 .fwdarw. 6.5 H.sup.+ 10.sup.−2.8 .fwdarw. transient .fwdarw. 10.sup.−2.8 PPi.sup.−−− 0 .fwdarw. transient .fwdarw. 0 Pi.sup.−− 0 (if no PPi .fwdarw. Pi enzyme present)
(530) Since the concentration changes are small, it is possible to carry out an expansion of the ionic strength in concentration changes:
I={119.227+6.67[PPi]+2.78[Pi]+0.79[H+released]−9.52(6.5μM−[dNTP])}×10.sup.−6 (5) The pH is 8.80.
(531) For BST reaction buffer at 0.5 mM Tris, which contains ions and sulfuric acid:
I={1452.46+6.5[PPi]+2.87[Pi]+0.98[H+released]−9.98(6.5μM−[dNTP])}×10.sup.−6 (5′)
(532) For a 5 mM Tris buffer, the ionic strength is given by
I={555.803+7.43[PPi]+2.92[Pi]+0.98[H+released]−9.90(6.5μM−[dNTP])}×10.sup.−6 (5″) The pH is 9.03.
(533) For a 0.5 mM Tris plus 1 mM MgCl.sub.2, the ionic strength is given by
I={3119.86+6.94[PPi]+2.78[Pi]+0.91[H+released]−9.49(6.5μM−[dNTP])}×10.sup.−6 (5′″) The pH is 8.81.
(534) For a 5 mM Tris plus 1 mM MgCl.sub.2, the ionic strength is given by
I={3556.67+7.71[PPi]+2.93[Pi]+0.98[H+released]−9.89(6.5μM−[dNTP])}×10.sup.−6 (5″″) The pH is 9.04.
(535) For 0.5 mM Tris, the ionic strength is given by
I={250.699+6.44[PPi]+2.83[Pi]+0.96[H+released]−9.67(20μM−[dNTP])}×10.sup.−6 (5′″″)
(536) The pH is 8.53 with the dNTP present and 8.97 without.
(537) For 0.5 mM Tris plus 50 μM MgCl.sub.2, the ionic strength is given by
I={400.759+6.46[PPi]+2.84[Pi]+0.96[H+released]−9.66(20μM−[dNTP])}×10.sup.−6 (5″″″)
(538) The pH is 8.53 with the dNTP present and 8.97 without.
(539) Finally, consider a canonical buffer with pH=9 and ionic strength=300 mM in the absence of dNTP. The equation would be
I={493.20+6.46[PPi]+2.84[Pi]+0.96[H+released]−9.66(20μM−[dNTP])}×10.sup.−6 (5′″″″)
(540) Note that the signal is proportional to I.sup.−1/2, so it is desirable to keep the ionic strength of the reaction buffer solution as low as possible. Note that (I.sub.0+dI).sup.−1/2=I.sub.0.sup.−1/2−½(dI/I.sub.0)I.sub.0.sup.−1/2. Thus, the relative change to the signal is −½ (dI/I.sub.0). By comparing Eqs. 5-5′″″″, one sees that dI is roughly independent of the buffer ionic strength. Thus, if the initial ionic strength, I.sub.0, is lower, the relative signal change is larger.
(541) Reaction/Diffusion Calculation
(542) The reactions that occur are
d[DNA]/dt=−k.sub.bst[dNTP][DNA]
d[PPi]/dt=k.sub.bst[dNTP][DNA]
d[dNTP]/dt=k.sub.bst[dNTP][DNA]
(543) If the enzyme for the reaction PPi+H.sub.2O.fwdarw.2 Pi is present, then this reaction is assumed to occur instantaneously. All of these species diffuse with individual diffusion coefficients, except the DNA. The H.sup.+ released by adding one more base of DNA, however, does diffuse.
(544) A finite element method may be used to calculate the reaction and diffusion of these species. The method preferably is encoded in cylindrical coordinates. Reaction in each element follows the above scheme, with the note below about the effect of Mg.sup.++ on the reaction rate. In addition, the diffusion of Mg.sup.++ is also tracked. Recall that in the current calculations, however, Mg.sup.++ is assumed to be constant.
(545) The numerical code was a check against the analytical solution for diffusion into a cylinder of diameter D and depth H:
C(r,z,t)=c0−c0(2/π)Σ.sub.n=0.sup.∞[e.sup.−D(n+1/2)2π{circumflex over ( )}2t/Ĥ2/(n+½)]sin[(n+½)πz/H]
(546) The results agreed to 0.1% accuracy if a grid of (D/2)/100 points in the r direction and H/100 points in the z direction are used.
(547) To determine the parameter k.sub.bst, the results from the Rothberg et al Nature paper were used. That paper [resented a simplified, mass-transfer calculation of the relevant diffusion, from which the relevant reactions are
d[DNA]/dt=−k.sub.bst[dNTP][DNA]
d[PPi]/dt=k.sub.bst[dNTP][DNA]−k.sub.c([PPi])
d[dNTP]/dt=k.sub.bst[dNTP][DNA]+k.sub.c([dNTP].sub.0−[dNTP])
(548) Summing the last two reactions, one finds
d([PPi]+[dNTP])/dt=−k.sub.bst[dNTP][DNA]−K.sub.c([dNTP].sub.0−[PPi]−[dNTP])
(549) These equations may be solved numerically to determine the values of k.sub.bst and k.sub.c, yielding k.sub.bst=2.38 (μM s).sup.−1 and k.sub.c=0.2 s.sup.−1 at 30° C.
(550) The PPi.fwdarw.Pi Enzyme
(551) The reaction PPi+H.sub.2O.fwdarw.2 Pi occurs at a finite rate. The rate constant in the presence of the enzyme and Mg.sup.++ is about 10.sup.3 s.sup.−1 (Biochemistry 41 (2002) 12025). So, this reaction is over in 1 ms. Thus, one can say this conversion is instantaneous. However, in the absence of any enzyme or ionic catalyst, the rate constant is about 3.sup.−1 years.sup.−1 (J. Chem. Soc. Farady Trans, 93 (1007) 4295). In this case, the reaction does not proceed at all. Many ions, however, do catalyze this reaction to some extent. Thus, if this reaction is desired not to proceed, care must be taken that the solution is chosen so that no ions catalyzing the reaction are present.
(552) Diffusion Coefficients
(553) The signal due to the pH change from the sequencing reaction depends on the diffusion coefficients of the species involved. Reference (Biophys. J 78 (2000) 1657) gives values for Pi, PPi, and dNTP:
(554) TABLE-US-00008 Species 10.sup.6 D (cm.sup.2/s) at 37° C. Pi 9.5 PPi 7.48 ATP 5.3
(555) Reference (BBA 1291 (1996) 115) gives values for ATP over a temperature range. At 37° C., the value is 4.5×10.sup.−6 cm.sup.2/s, which is used herein. The values for Pi and PPi are calculated at temperature T by scaling the value above by the factor D.sub.ATP(T)/D.sub.ATP(37° C.). The diffusion coefficients of the ions at 25° C. are (J. Chem. Phys. 119 (2003) 11342.)
(556) TABLE-US-00009 Species 10.sup.5 D (cm.sup.2/s) at 25° C. Mg++ 1.25 (J. Amer. Chem. Soc. 77 (1955) 265) Na+ 1.33 (J. Chem. Phys. 119 (2003) 11342.)
(557) The diffusion coefficient of H.sup.+ is assumed equal to that of Nat, because by charge neutrality the H.sup.+ must either diffusion with a negative ion as a pair, or the H.sup.+ must exchange with another positive ion, of which the N.sup.+ is the most rapidly diffusing.
(558) One preferred implementation includes the use of small spheres (e.g. 0.1 μm) to decrease the diffusion coefficients of all the species. This shifts the pH curves to larger time by roughly the ratio that the diffusion coefficient is reduced. The textbook formula for the reduction of the diffusion coefficient is
D.sub.eff/D.sub.0=2(1−φ)/(2+φ), (6)
(559) where φ is the volume fraction occupied by the spheres. The radius of the small spheres does not directly enter into this expression. For no spheres φ=0. For randomly packed spheres, φ≈0.6, and for close packed spheres, φ=0.74.
(560) Another approach is to use other methods to reduce the diffusion coefficient of the H.sup.+, Pi, PPi, and dNTP. For example, the viscosity of the solution may be increased.
(561) The Sequencing Reaction
(562) Now consider the effect of the Mg.sup.++ diffusing in. It is not known exactly how the enzyme reaction rate depends on [Mg.sup.++]. However, a divalent metal ion is required for enzymatic activity. It is assumed that the rate is linearly proportional up to a maximum rate at 1 mM of Mg.sup.++. Thus, it is assumed the reaction rate constant is equal to
k.sub.bst min(1,[Mg.sup.++]/1mM) (7)
(563) One may use a diffusion coefficient for Mg.sup.++ of 1.25×10.sup.−5 cm.sup.2/2. For Model 2, at t=0, the Mg.sup.++ is only at z=0, at concentration [Mg.sup.++].sub.0=1 mM. For t>0 it diffuses in and affects the reaction rate by Eq. (8). For Model 1, the concentration of Mg.sup.++ is 1 mM at all times in the well. Alternatively, we may assume that the enzyme is presoaked with Mg.sup.++, so that the concentration of free Mg.sup.++ in solution is zero.
(564) In the range of 7-9 pH, the reactions release H.sup.+ roughly as follows
(565) TABLE-US-00010 DNA.sub.n n DNA.sub.n+1 n + 1 dNTP 4 PPi 3 2 Pi 4
dNTP+DNA.sub.n.fwdarw.PPi+DNA.sub.n+1 (8)
(566) is to not change the pH.
(567) The net effect of the reaction
PPi+H.sub.2O.fwdarw.2Pi (9)
(568) is to change the pH by release of one net H.sup.+.
(569) The H.sup.+ released by the DNA on the bead surface will diffuse out to beyond the DNA, and since the pH is much higher than the pKa of the DNA, these H.sup.+ ions will diffuse out of the well during the equilibration time t<0. (All these H.sup.+ diffuse away.) That is, if the DNA on the bead is allowed to equilibrate before the dNTP is added, the solution should just be the usual buffer solution (with a bit of extra Na.sup.+ near the bead to balance the negative charge on the DNA, within the Debye length). When the extra base of DNA is added, a pulse of H.sup.+ will come off—with a diffusion coefficient of something like Nat
(570) Well Sizes
(571) Simulation Runs were Done for Small Wells:
(572) TABLE-US-00011 D(μm) H(μm) 2 R.sub.bead (μm) V.sub.well (pl) [DNA].sub.well (μm) 1.5 2.25 1.05 0.005063 4.62 4 6 2.8 0.096 1.73 6 8 5.91 0.29 2.55
(573) For D=1.5 μm, h=0.2 μm (1 layer of beads). For D=4 μm, h=1.1 μm. Otherwise h=1.5 μm. In addition, one run was done for the 4 μm×6 μm wells with a larger bead of radius 3.6 μm.
(574) dNTP Front Thickness in Device Flow
(575) The dNTP is assumed in the calculations to instantaneously change from 0 to 6.5 μM above the wells at t=0. This condition is approximately true if the dNTP flows quickly enough past the wells. Define
(576) V=velocity of flow of the fluid above the well
(577) L=length of device
(578) D=diffusivity of dNTP=3.91×10-6 cm2/s
(579) d=size of well
(580) t=time to sweep flow across device
(581) t0=cutoff time for “instantaneous change”
(582) Wd=width of front due to diffusion at end of device=(2 D t)½
(583) Wm=width of front due to mixing length of flush volume region
(584) U=flush volume=15 μl (currently 52 μl)
(585) u=flow rate=50 μl/s
(586) Flow Cell Geometries are as Follows for Two Models, Designated 313 and 314:
(587) TABLE-US-00012 axis Width Height Volume Cell Flow Cell (mm) (mm) (mm) (μl) Volumes/s 314 7.249 7.249 0.336 17.657 2.775 313 4.999 4.999 0.293 7.321 6.693
(588) Take L=5 mm, d=4 μm and assume t=1/6.69s. Q=50 μl/s. Then Wd=10.8 μm. Also, V=L/t=3.3 cm/s. D.sub.eff=D.sub.0 2 (1−φ)/(2+φ)=1.20×10.sup.−6 cm.sup.2/s for φ=0.6. The value of to depends on the well size. It should be much shorter than the time at which dNTP diffuses to the bottom of the well for N.sub.poly=0 (i.e. much shorter than the time at which this ion pulse response occurs), say 10% of that value. For the 4×6 μm well the time for the dNTP to diffuse to the bottom of the well, assuming no obstruction to the diffusion is H.sup.2/(2D), so that to =10%×H.sup.2/(2D)=0.015 s. Simulation reveals the sphere with the DNA on it slows down the diffusion, so that perhaps t.sub.0=10%*0.2 s=0.02 s. Define W=max(W.sub.d, W.sub.m). The time the leading edge of the front takes to cross the device is d/V. The time it takes the whole front (leading plus trailing edge) to cross a well is (W+d)/V. The condition of “instantaneous” rise of the dNTP concentration is (W+d)/V<<t.sub.0. (W+d)/V≈U/u. Since U/u=0.3 s, and t.sub.0=0.02 s, this condition is not satisfied by the above assumptions. Either U must be made smaller or u must be made larger, or both, by a factor of 10. The condition that W.sub.m is given by the length of the flush region may be a bit too strict: the width of the mixing region may be substantially smaller than the length of the entire flush region. This is a conservative condition. However, W.sub.d V<<t.sub.0, so diffusive spreading of the front can be ignored.
(589) Now, a different calculation. It is desirable for the concentration of dNTP at the wall of the channel to rise from 0 to 6.5 μM in less than to. The concentration C(x,y,z) at t=0 may be assumed to be c.sub.0=6.5 μM at t=0. Take 0<x<L, −B<y<B, 0<z<W. It is desired to know the concentration at y=B for any value of 0<x<L. Strictly speaking, the concentration at the wall stays 0, since the velocity there is zero, unless diffusion is considered. So, it is necessary to consider convection in the x direction and diffusion in the z direction. Calculating the convection
v=[P.sub.0−P.sub.L)B.sup.2/(2μL)][1−(y/B).sup.2] (10)
(590) The flow rate is Q=(⅔) (p.sub.0−p.sub.L) B.sup.3 W/(μL). Calculate the time it takes for the average flow at the end of the device to be 0.95c.sub.0. First, note that without diffusion, the concentration is given by
c(x,y,z)=0,t<t*
c.sub.0,t>t*
t*=x/v(y) (11)
The flow becomes 95% of the maximum concentration at t′=L/v(0.95 B). The time to clear one volume is t″=volume/Q, where volume=2B W L. It can be found that t″=2/3/0.0975 t″=6.87 t″. Thus it takes 7 device volumes to bring the average concentration up to 95% of the final value (not 3). Of concern is the value of the concentration at the device wall. The concentration obeys the equation
∂c/∂t+v(y)∂c/∂x=D∂.sup.2c/∂y.sup.2
c(x,y,0)=0
c(0,y,t)=c.sub.0
∂c(x,±B,t)/∂y=0 (12)
(591) ignoring diffusion in the y direction as dominated by convection. Now the concentration reaches the wall by diffusion in a time δt, but diffusion from the concentration profile that existed at a smaller x−δx and y=y−δy, since convection occurs during this time δt. Note (δy).sup.2=2 D t. From the above, it is known that the concentration can diffuse on the order of δy=10 μm by diffusion during a time t″=1/6.69 s. It is desired that it diffuse during a time t.sub.0=0.02 s, which implies a diffusion distance of 3.95 μm. So, as a conservative calculation, one could calculate from Eq. (11) t*=L/v(B−3.95 μm). By the time t*+t.sub.0, the concentration at the wall will be approximately co through a combination of diffusion and convection. Thus
v(¾)(1/B W)Q[1−(y/B){circumflex over ( )}2]
(592) Using this, t*=L/v(B−3.95 μm)=1.84 s. At a value to later, the concentration at the wall will rise approximately to the maximum. Unfortunately, this calculation ignores that during 1.84 s, material from velocity profiles from smaller values of y will also have had the chance to diffuse to the wall. It is also ignoring that the diffusion which starts for smaller values of x proceeds to the right with lower velocities as y gets larger.
(593) Thus, one can ask when diffusion over the boundary layer is faster than convection. By diffusion δ/t=2 D/δ. By convection δ/t=v(B−δ). These two are equal for δ=0.001 mm=1.0 μm. So, when the front has reached this value of δ=L/v(B−δ)=6.78 s, the progression of the material to the wall is entirely dominated by diffusion. And the concentration will rise to approximately co during the time δ.sup.2/(2 D)=0.01 s. So, for δ<1.0 convection is more important than diffusion in the x direction.
(594) A numerical solution of Eq. (12) is necessary to check the above approximate calculations. A finite element method can be used to numerically solve Eq. (12). We find that for the above device, the concentration at the wall rises in about 0.15 s at x=0.5 mm and about 0.5 s at x=5.0 mm. This is too slow. The Reynolds number is Re=1 g/cm3*0.5 cm/(⅙.69 s)*0.0293 cm/0.01 cp=9.8. If one sets the height to be 0.1 mm, one finds that the concentration at the wall rises in about 0.04 s at x=0.5 mm and about 0.15 s at x=5.0 mm. This is probably too slow. Re=9.8. If set alternatively set Q=240 μl/s, the concentration at the wall rises in about 0.05 s at x=0.5 mm and about 0.2 s at x=5.0 mm. This is probably too slow. Re=47. If one sets both the height to be 0.1 mm and Q=240 μl/s, the concentration at the wall rises in about 0.013 s at x=0.5 mm and about 0.07 s at x=5.0 mm. This is probably fast enough. Re=47. The flow in all cases should be laminar.
(595) The above may be a bit too strict, in that the biggest signal is coming from the breakthrough of the dNTP, and this breakthrough is quite sharp. It is possible that the dNTP rise at the top of the well need not be complete in 0.02 s. Perhaps the rise just has to be complete significantly before the dNTP breakthrough at the bottom of the well occurs. And this is on the order of seconds. So, one needs to solve the reaction/diffusion equations with a dNTP concentration that is rising at the top of the well. Let it rise from 0 to 6.5 μM linearly between 0 and t, where t is the value specified in the previous paragraph for the rise time at the end of the device.
(596) From these numerical simulations, it may be found that a finite rise time of the dNTP at the surface of the well reduces the performance of the ion pulse device only slightly. For a 6 μm×8 μm device, a finite rise time of 0.2 s affects the observed sensor response curves only slightly. A finite rise time of 0.5 s affects the observed sensor response curves by shifting to longer times the first peak as well as reducing the height of the peaks for all values of N.sub.poly. A surface with a slower response time would be less affected by this finite dNTP rise time. The effects are similar, if a bit more dramatic, for the 4 μm×6 μm, since the characteristic diffusion times for this device are smaller and therefore the finite rise time is a larger relative perturbation to the idealized case of instantaneous dNTP rise at the surface of the well.
(597) Parametric Sensitivity of Results
(598) To determine the sensitivity of the ion pulse response to experimental variables, the 4 μm×6 μm system was taken as a reference. Since this detection method relies on diffusion, the well geometry has a significant impact on the resolution of any N.sub.poly curves that are plotted.
(599) The dNTP concentration was halved to 3.25 μM and doubled to 13 μM. The DNA on the bead surface was increased by a factor of 3. A run was done with a larger 3.6 μm bead. Both the DNA on the bead surface and the dNTP were increased by a factor of 3. The reaction rate, k, of the enzyme was halved and doubled. Sensor response curves for the 5 μM Tris, no free Mg.sup.++, were produced for all of the above.
(600) Double the enzyme reaction rate constant makes the peaks very slightly higher.
(601) One possibility is that the DNA concentration on the beads is so high that the reaction is effectively diffusion limited, and so the enzyme rate constant is not a major factor affecting the concentration curves. Similarly, reducing the enzyme reaction rate by ½ has little effect on the sensor curves. It spreads them out in time only very slightly.
(602) The large bead works very well. Since there is more mass of DNA on the bead, consuming it takes longer. In addition, the bead is taking up more space in the volume, so that diffusion to the well bottom takes a bit longer. The peaks are spread out to larger times.
(603) A larger concentration of DNA on the bead spreads out the peaks in time. This is because the larger concentration of DNA means a greater mass of DNA on the beads. This greater mass of DNA takes a longer time to reaction. Thus the final breakthrough of the dNTP takes longer.
(604) Lowering the solution dNTP concentration to 3.25 μM spreads out the peaks in time, but also reduces the peak height. The two effects probably approximately cancel in terms of distinguishing the different N.sub.poly curves. Conversely, increasing the dNTP concentration to 13 μM shifts the peaks to smaller time. While the peaks are higher, they are also substantially closer in time. The curves are at such a shorter time that they are running into the response time of the sensor. It appears the 13 μM curves may be less distinguishable than the 6.5 μM case.
(605) Increasing the concentration of DNA on the beads three times and increasing the dNTP concentration in the solution to 13 μM makes the peaks about 3 times higher but does not spread them out in time. The time scale is relatively unchanged because while the amount of DNA on the beads is greater, the reaction rate is also greater due to the higher dNTP concentration. The peaks are higher because there is more mass released in the same time.
(606) Packing Beads—Φ (Phi)
(607) The desirability for packing beads (to slow down diffusion of all species) is seen by comparing the Ion Pulse Response between similar well dimensions and between the 4×6 micron (width×height) and 6×8 micron wells. The N.sub.poly curves resolve with increased depth of the well and with packing beads (Φ=0.6).
(608) The formula to see the effect of Φ is the effective diffusion coefficient
D.sub.eff=D.sub.0*2(1−Φ))/(2+Φ))
(609) For Φ=0.6, D.sub.0/D.sub.eff=3.25 e.g. 3.25× slowdown due to beads
(610) For Φ=0.9, D.sub.0/D.sub.eff=14.5 e.g. 14.5× slowdown due to this The other formula needed is the characteristic time in the well
H.sup.2=2D.sub.efft
(611) So, for a 6 um well at Φ=0.6
t=6.sup.2/(2/3.25*D.sub.0)=58.5/D.sub.0
(612) A 2.25 well at Φ=0.9 yields
t=2.25.sup.2/(2/14.5*D.sub.0)=36.7/D.sub.0
(613) So, the 2.25 um well at Φ=0.9 would have all times about 36.7/58.5=0.63 smaller than the 6 um well at Φ=0.6. The smaller well will still be a little bit faster. If one had Φ=0.936, the smaller well would operate just about exactly as the larger well (assuming the well diameter and the DNA bead are scaled the same way as the well height).
(614) Signal to Noise
(615) The signal-to-noise ratio is proportional to the difference in mV between a N.sub.poly and N.sub.poly+1 signal and also to the square root of the length of time over which the mV curves are different. Define a as the sensor noise and f as the number of frames per second. Define the average difference of the absolute value between signal 1 and signal 2 during 0<t<t.sub.max as Δ, i.e.
Δ=∫.sub.0.sup.tmax|v.sub.1−v.sub.2|/t.sub.max.
(616) Take the absolute value because the difference between the curves can have both signs. Thus the condition is
Δ>σ[2/(ft.sub.max)].sup.1/2
(617) Consider, as an example, for the D=4 μm, H=6 μm wells, φ=0, dNTP=20 μM, 10×DNA case. Take f=20 Hz, and t.sub.max=2 s. The value of A is 0.6861. Thus
0.6861>0.22 σ
(618) So, using the value of σ=0.25 mV, this device will work (1−P=10.sup.−35). Consider as another example the D=4 μm, H=6 μm wells, y=0, dNTP=6.5 μM, 10×DNA case. Take f=20 Hz, and t.sub.max=4 s. The value of Δ is 0.1862. Thus
0.1862>0.1581σ
(619) So, using the value of σ=0.25 mV, this device will probably work (1−P=2.4×10.sup.−6).
(620) Note there is an optimal value of t.sub.max: too small a value, and one does not have enough signal; too large a value, and one is seeing mostly noise. The value of t.sub.max is only roughly here.
EQUIVALENTS
(621) While several inventive embodiments have been described and illustrated herein, those of ordinary skill in the art will readily envision a variety of other means and/or structures for performing the function and/or obtaining the results and/or one or more of the advantages described herein, and each of such variations and/or modifications is deemed to be within the scope of the inventive embodiments described herein. More generally, those skilled in the art will readily appreciate that all parameters, dimensions, materials, and configurations described herein are meant to be exemplary and that the actual parameters, dimensions, materials, and/or configurations will depend upon the specific application or applications for which the inventive teachings is/are used. Those skilled in the art will recognize, or be able to ascertain using no more than routine experimentation, many equivalents to the specific inventive embodiments described herein. It is, therefore, to be understood that the foregoing embodiments are presented by way of example only and that, within the scope of the appended claims and equivalents thereto, inventive embodiments may be practiced otherwise than as specifically described and claimed. Inventive embodiments of the present disclosure are directed to each individual feature, system, article, material, kit, and/or method described herein. In addition, any combination of two or more such features, systems, articles, materials, kits, and/or methods, if such features, systems, articles, materials, kits, and/or methods are not mutually inconsistent, is included within the inventive scope of the present disclosure.
(622) All definitions, as defined and used herein, should be understood to control over dictionary definitions, definitions in documents incorporated by reference, and/or ordinary meanings of the defined terms.
(623) All references, patents and patent applications disclosed herein are incorporated by reference with respect to the subject matter for which each is cited, which in some cases may encompass the entirety of the document.
(624) The indefinite articles “a” and “an,” as used herein in the specification and in the claims, unless clearly indicated to the contrary, should be understood to mean “at least one.”
(625) The phrase “and/or,” as used herein in the specification and in the claims, should be understood to mean “either or both” of the elements so conjoined, i.e., elements that are conjunctively present in some cases and disjunctively present in other cases. Multiple elements listed with “and/or” should be construed in the same fashion, i.e., “one or more” of the elements so conjoined. Other elements may optionally be present other than the elements specifically identified by the “and/or” clause, whether related or unrelated to those elements specifically identified. Thus, as a non-limiting example, a reference to “A and/or B”, when used in conjunction with open-ended language such as “comprising” can refer, in one embodiment, to A only (optionally including elements other than B); in another embodiment, to B only (optionally including elements other than A); in yet another embodiment, to both A and B (optionally including other elements); etc.
(626) As used herein in the specification and in the claims, “or” should be understood to have the same meaning as “and/or” as defined above. For example, when separating items in a list, “or” or “and/or” shall be interpreted as being inclusive, i.e., the inclusion of at least one, but also including more than one, of a number or list of elements, and, optionally, additional unlisted items. Only terms clearly indicated to the contrary, such as “only one of” or “exactly one of,” or, when used in the claims, “consisting of,” will refer to the inclusion of exactly one element of a number or list of elements. In general, the term “or” as used herein shall only be interpreted as indicating exclusive alternatives (i.e. “one or the other but not both”) when preceded by terms of exclusivity, such as “either,” “one of,” “only one of,” or “exactly one of.” “Consisting essentially of,” when used in the claims, shall have its ordinary meaning as used in the field of patent law.
(627) As used herein in the specification and in the claims, the phrase “at least one,” in reference to a list of one or more elements, should be understood to mean at least one element selected from any one or more of the elements in the list of elements, but not necessarily including at least one of each and every element specifically listed within the list of elements and not excluding any combinations of elements in the list of elements. This definition also allows that elements may optionally be present other than the elements specifically identified within the list of elements to which the phrase “at least one” refers, whether related or unrelated to those elements specifically identified. Thus, as a non-limiting example, “at least one of A and B” (or, equivalently, “at least one of A or B,” or, equivalently “at least one of A and/or B”) can refer, in one embodiment, to at least one, optionally including more than one, A, with no B present (and optionally including elements other than B); in another embodiment, to at least one, optionally including more than one, B, with no A present (and optionally including elements other than A); in yet another embodiment, to at least one, optionally including more than one, A, and at least one, optionally including more than one, B (and optionally including other elements); etc.
(628) It should also be understood that, unless clearly indicated to the contrary, in any methods claimed herein that include more than one step or act, the order of the steps or acts of the method is not necessarily limited to the order in which the steps or acts of the method are recited.
(629) In the claims, as well as in the specification above, all transitional phrases such as “comprising,” “including,” “carrying,” “having,” “containing,” “involving,” “holding,” “composed of,” and the like are to be understood to be open-ended, i.e., to mean including but not limited to. Only the transitional phrases “consisting of” and “consisting essentially of” shall be closed or semi-closed transitional phrases, respectively, as set forth in the United States Patent Office Manual of Patent Examining Procedures, Section 2111.03.