METHODS AND APPARATUS FOR MEASURING ANALYTES USING LARGE SCALE FET ARRAYS
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.
This application is a divisional of U.S. application Ser. No. 13/966,184 filed Aug. 13, 2013 (now U.S. Pat. No. 9,458,502), which is a continuation of U.S. application Ser. No. 13/193,128 filed Jul. 28, 2011 (now U.S. Pat. No. 8,524,057), which is a continuation of U.S. application Ser. No. 13/001,182 filed May 23, 2011 (now U.S. Pat. No. 8,470,164), which is a national stage filing under 35 U.S.C. §371 of PCT International application PCT/US2009/003766 filed Jun. 25, 2009, PCT Article 21(2) in English, which claims priority to patent applications GB 0811657.6 and GB 0811656.8, both filed on Jun. 25, 2008, the entire contents of all of which are incorporated by reference.
FIELD OF THE DISCLOSUREThe present disclosure is directed generally to inventive methods and apparatus relating to detection and measurement of one or more analytes via electronic sensors.
BACKGROUNDElectronic devices and components have found numerous applications in chemistry and biology (more generally, “life sciences”), especially for detection and measurement of various chemical and biological reactions and identification, detection and measurement of various compounds. One such electronic device is referred to as an ion-sensitive field effect transistor, often denoted in the relevant literature as ISFET (or pHFET). ISFETs conventionally have been explored, primarily in the academic and research community, to facilitate measurement of the hydrogen ion concentration of a solution (commonly denoted as “pH”).
More specifically, an ISFET is an impedance transformation device that operates in a manner similar to that of a MOSFET (Metal Oxide Semiconductor Field Effect Transistor), and is particularly configured to selectively measure ion activity in a solution (e.g., hydrogen ions in the solution are the “analytes”). A detailed theory of operation of an ISFET is given in “Thirty years of ISFETOLOGY: what happened in the past 30 years and what may happen in the next 30 years,” P. Bergveld, Sens. Actuators, 88 (2003), pp. 1-20, which publication is hereby incorporated herein by reference (hereinafter referred to as “Bergveld”).
Like a MOSFET, the operation of an ISFET is based on the modulation of charge concentration caused by a MOS (Metal-Oxide-Semiconductor) capacitance constituted by the polysilicon gate 64, the gate oxide 65 and the region 60 of the n-type well 54 between the source and the drain. When a negative voltage is applied across the gate and source regions (VGS<0 Volts), a “p-channel” 63 is created at the interface of the region 60 and the gate oxide 65 by depleting this area of electrons. This p-channel 63 extends between the source and the drain, and electric current is conducted through the p-channel when the gate-source potential VGS is negative enough to attract holes from the source into the channel. The gate-source potential at which the channel 63 begins to conduct current is referred to as the transistor's threshold voltage VTH (the transistor conducts when VGS has an absolute value greater than the threshold voltage VTH). The source is so named because it is the source of the charge carriers (holes for a p-channel) that flow through the channel 63; similarly, the drain is where the charge carriers leave the channel 63.
In the ISFET 50 of
As also shown in
With respect to ion sensitivity, an electric potential difference, commonly referred to as a “surface potential,” arises at the solid/liquid interface of the passivation layer 72 and the analyte solution 74 as a function of the ion concentration in the sensitive area 78 due to a chemical reaction (e.g., usually involving the dissociation of oxide surface groups by the ions in the analyte solution 74 in proximity to the sensitive area 78). This surface potential in turn affects the threshold voltage VTH of the ISFET; thus, it is the threshold voltage VTH of the ISFET that varies with changes in ion concentration in the analyte solution 74 in proximity to the sensitive area 78.
ID=β(VGS−VTH−½VDS)VDS (1)
where VDS is the voltage between the drain and the source, and β is a transconductance parameter (in units of Amps/Volts2) given by:
where μ represents the carrier mobility, Cox is the gate oxide capacitance per unit area, and the ratio W/L is the width to length ratio of the channel 63. If the reference electrode 76 provides an electrical reference or ground (VG=0 Volts), and the drain current ID and the drain-to-source voltage VDS are kept constant, variations of the source voltage VS of the ISFET directly track variations of the threshold voltage VTH, according to Eq. (1); this may be observed by rearranging Eq. (1) as:
Since the threshold voltage VTH of the ISFET is sensitive to ion concentration as discussed above, according to Eq. (3) the source voltage VS provides a signal that is directly related to the ion concentration in the analyte solution 74 in proximity to the sensitive area 78 of the ISFET. More specifically, the threshold voltage VTH is given by:
where VFB is the flatband voltage, QB is the depletion charge in the silicon and φF is the Fermi-potential. The flatband voltage in turn is related to material properties such as workfunctions and charge accumulation. In the case of an ISFET, with reference to
where Eref is the reference electrode potential relative to vacuum, Ψ0 is the surface potential that results from chemical reactions at the analyte solution/passivation layer interface (e.g., dissociation of surface groups in the passivation layer), and χsol is the surface dipole potential of the analyte solution 74. The fourth term in Eq. (5) relates to the silicon workfunction (q is the electron charge), and the last term relates to charge densities at the silicon surface and in the gate oxide. The only term in Eq. (5) sensitive to ion concentration in the analyte solution 74 is Ψ0, as the ion concentration in the analyte solution 74 controls the chemical reactions (dissociation of surface groups) at the analyte solution/passivation layer interface. Thus, substituting Eq. (5) into Eq. (4), it may be readily observed that it is the surface potential Ψ0 that renders the threshold voltage VTH sensitive to ion concentration in the analyte solution 74.
Regarding the chemical reactions at the analyte solution/passivation layer interface, the surface of a given material employed for the passivation layer 72 may include chemical groups that may donate protons to or accept protons from the analyte solution 74, leaving at any given time negatively charged, positively charged, and neutral sites on the surface of the passivation layer 72 at the interface with the analyte solution 74. A model for this proton donation/acceptance process at the analyte solution/passivation layer interface is referred to in the relevant literature as the “Site-Dissociation Model” or the “Site-Binding Model,” and the concepts underlying such a process may be applied generally to characterize surface activity of passivation layers comprising various materials (e.g., metal oxides, metal nitrides, metal oxynitrides).
Using the example of a metal oxide for purposes of illustration, the surface of any metal oxide contains hydroxyl groups that may donate a proton to or accept a proton from the analyte to leave negatively or positively charged sites, respectively, on the surface. The equilibrium reactions at these sites may be described by:
AOHAO−+HS+ (6)
AOH2+AOH+HS+ (7)
where A denotes an exemplary metal, Hs+ represents a proton in the analyte solution 74, Eq. (6) describes proton donation by a surface group, and Eq. (7) describes proton acceptance by a surface group. It should be appreciated that the reactions given in Eqs. (6) and (7) also are present and need to be considered in the analysis of a passivation layer comprising metal nitrides, together with the equilibrium reaction:
ANH3+ANH2+H+, (7b)
wherein Eq. (7b) describes another proton acceptance equilibrium reaction. For purposes of the present discussion however, again only the proton donation and acceptance reactions given in Eqs. (6) and (7) are initially considered to illustrate the relevant concepts.
Based on the respective forward and backward reaction rate constants for each equilibrium reaction, intrinsic dissociation constants Ka (for the reaction of Eq. (6)) and Kb (for the reaction of Eq. (7)) may be calculated that describe the equilibrium reactions. These intrinsic dissociation constants in turn may be used to determine a surface charge density σ0 (in units of Coulombs/unit area) of the passivation layer 72 according to:
σ0=−qB, (8)
where the term B denotes the number of negatively charged surface groups minus the number of positively charged surface groups per unit area, which in turn depends on the total number of proton donor/acceptor sites per unit area NS on the passivation layer surface, multiplied by a factor relating to the intrinsic dissociation constants Ka and Kb of the respective proton donation and acceptance equilibrium reactions and the surface proton activity (or pHS). The effect of a small change in surface proton activity (pHS) on the surface charge density is given by:
where βint is referred to as the “intrinsic buffering capacity” of the surface. It should be appreciated that since the values of NS, Ka and Kb are material dependent, the intrinsic buffering capacity βint of the surface similarly is material dependent.
The fact that ionic species in the analyte solution 74 have a finite size and cannot approach the passivation layer surface any closer than the ionic radius results in a phenomenon referred to as a “double layer capacitance” proximate to the analyte solution/passivation layer interface. In the Gouy-Chapman-Stern model for the double layer capacitance as described in Bergveld, the surface charge density σ0 is balanced by an equal but opposite charge density in the analyte solution 74 at some position from the surface of the passivation layer 72. These two parallel opposite charges form a so-called “double layer capacitance” Cdl (per unit area), and the potential difference across the capacitance Cdl is defined as the surface potential Ψ0, according to:
σ0=CdlΨ0=−σdl (10)
where σdl is the charge density on the analyte solution side of the double layer capacitance. This charge density σdl in turn is a function of the concentration of all ion species or other analyte species (i.e., not just protons) in the bulk analyte solution 74; in particular, the surface charge density can be balanced not only by hydrogen ions but other ion species (e.g., Na+, K+) in the bulk analyte solution.
In the regime of relatively lower ionic strengths (e.g., <1 mole/liter), the Debye theory may be used to describe the double layer capacitance Cdl according to:
where k is the dielectric constant ∈/∈0 (for relatively lower ionic strengths, the dielectric constant of water may be used), and λ is the Debye screening length (i.e., the distance over which significant charge separation can occur). The Debye length λ is in turn inversely proportional to the square root of the strength of the ionic species in the analyte solution, and in water at room temperature is given by:
The ionic strength I of the bulk analyte is a function of the concentration of all ionic species present, and is given by:
I=½Σsz2cs, (13)
where zs is the charge number of ionic species s and cs is the molar concentration of ionic species s. Accordingly, from Eqs. (10) through (13), it may be observed that the surface potential is larger for larger Debye screening lengths (i.e., smaller ionic strengths).
The relation between pH values present at the analyte solution/passivation layer interface and in the bulk solution is expressed in the relevant literature by Boltzman statistics with the surface potential Ψ0 as a parameter:
From Eqs. (9), (10) and (14), the sensitivity of the surface potential Ψ0 particularly to changes in the bulk pH of the analyte solution (i.e., “pH sensitivity”) is given by:
where the parameter α is a dimensionless sensitivity factor that varies between zero and one and depends on the double layer capacitance Cdl and the intrinsic buffering capacity of the surface βint as discussed above in connection with Eq. (9). In general, passivation layer materials with a high intrinsic buffering capacity βint render the surface potential Ψ0 less sensitive to concentration in the analyte solution 74 of ionic species other than protons (e.g., a is maximized by a large βint). From Eq. (15), at a temperature T of 298 degrees Kelvin, it may be appreciated that a theoretical maximum pH sensitivity of 59.2 mV/pH may be achieved at α=1. From Eqs. (4) and (5), as noted above, changes in the ISFET threshold voltage VTH directly track changes in the surface potential Ψ0; accordingly, the pH sensitivity of an ISFET given by Eq. (15) also may be denoted and referred to herein as ΔVTH for convenience. In exemplary conventional ISFETs employing a silicon nitride or silicon oxynitride passivation layer 72 for pH-sensitivity, pH sensitivities ΔVTH (i.e., a change in threshold voltage with change in pH of the analyte solution 74) over a range of approximately 30 mV/pH to 60 mV/pH have been observed experimentally.
Another noteworthy metric in connection with ISFET pH sensitivity relates to the bulk pH of the analyte solution 74 at which there is no net surface charge density σ0 and, accordingly, a surface potential Ψ0 of zero volts. This pH is referred to as the “point of zero charge” and denoted as pHpzc. With reference again to Eqs. (8) and (9), like the intrinsic buffering capacity βint, pHpzc is a material dependent parameter. From the foregoing, it may be appreciated that the surface potential at any given bulk pHB of the analyte solution 74 may be calculated according to:
Table 1 below lists various metal oxides and metal nitrides and their corresponding points of zero charge (pHpzc), pH sensitivities (ΔVTH), and theoretical maximum surface potential at a pH of 9:
Prior research efforts to fabricate ISFETs for pH measurements based on conventional CMOS processing techniques typically have aimed to achieve high signal linearity over a pH range from 1-14. Using an exemplary threshold sensitivity of approximately 50 mV/pH, and considering Eq. (3) above, this requires a linear operating range of approximately 700 mV for the source voltage VS. As discussed above in connection with
While the foregoing discussion relates primarily to a steady state analysis of ISFET response based on the equilibrium reactions given in Eqs. (6) and (7), the transient or dynamic response of a conventional ISFET to an essentially instantaneous change in ionic strength of the analyte solution 74 (e.g., a stepwise change in proton or other ionic species concentration) has been explored in some research efforts. One exemplary treatment of ISFET transient or dynamic response is found in “ISFET responses on a stepwise change in electrolyte concentration at constant pH,” J. C. van Kerkof, J. C. T. Eijkel and P. Bergveld, Sensors and Actuators B, 18-19 (1994), pp. 56-59, which is incorporated herein by reference.
For ISFET transient response, a stepwise change in the concentration of one or more ionic species in the analyte solution in turn essentially instantaneously changes the charge density σdl on the analyte solution side of the double layer capacitance Cdl. Because the instantaneous change in charge density σdl is faster than the reaction kinetics at the surface of the passivation layer 72, the surface charge density σ0 initially remains constant, and the change in ion concentration effectively results in a sudden change in the double layer capacitance Cdl. From Eq. (10), it may be appreciated that such a sudden change in the capacitance Cdl at a constant surface charge density σ0 results in a corresponding sudden change in the surface potential Ψ0.
As indicated in the bottom graph of
where Ψ1 is an equilibrium surface potential at an initial ion concentration in the analyte solution, Cdl,1 is the double layer capacitance per unit area at the initial ion concentration, Ψ2 is the surface potential corresponding to the ion-step stimulus, and Cdl,2 is the double layer capacitance per unit area based on the ion-step stimulus. The time decay profile 81 associated with the response 79 is determined at least in part by the kinetics of the equilibrium reactions at the analyte solution/passivation layer interface (e.g., as given by Eqs. (6) and (7) for metal oxides, and also Eq. (7b) for metal nitrides). One instructive treatment in this regard is provided by “Modeling the short-time response of ISFET sensors,” P. Woias et al., Sensors and Actuators B, 24-25 (1995) 211-217 (hereinafter referred to as “Woias”), which publication is incorporated herein by reference.
In the Woias publication, an exemplary ISFET having a silicon nitride passivation layer is considered. A system of coupled non-linear differential equations based on the equilibrium reactions given by Eqs. (6), (7), and (7a) is formulated to describe the dynamic response of the ISFET to a step (essentially instantaneous) change in pH; more specifically, these equations describe the change in concentration over time of the various surface species involved in the equilibrium reactions, based on the forward and backward rate constants for the involved proton acceptance and proton donation reactions and how changes in analyte pH affect one or more of the reaction rate constants. Exemplary solutions, some of which include multiple exponential functions and associated time constants, are provided for the concentration of each of the surface ion species as a function of time. In one example provided by Woias, it is assumed that the proton donation reaction given by Eq. (6) dominates the transient response of the silicon nitride passivation layer surface for relatively small step changes in pH, thereby facilitating a mono-exponential approximation for the time decay profile 81 of the response 79 according to:
Ψ0(t)=ΔΨ0e−1/τ, (18)
where the exponential function essentially represents the change in surface charge density as a function of time. In Eq. (16), the time constant τ is both a function of the bulk pH and material parameters of the passivation layer, according to:
τ=τ0×10pH/2 (19)
where τ0 denotes a theoretical minimum response time that only depends on material parameters. For silicon nitride, Woias provides exemplary values for τ0 on the order of 60 microseconds to 200 microseconds. For purposes of providing an illustrative example, using τ0=60 microseconds and a bulk pH of 9, the time constant r given by Eq. (19) is 1.9 seconds. Exemplary values for other types of passivation materials may be found in the relevant literature and/or determined empirically.
Previous efforts to fabricate two-dimensional arrays of ISFETs based on the ISFET design of
As shown in
As also shown in
In the column design of Milgrew et al. shown in
It should also be appreciated that in the column design of Milgrew et al. shown in
Furthermore, in the design of Milgrew et al., the p-channel MOSFET required to implement the transmission gate S1 in each pixel (e.g., see S11P in
The array design of Milgrew et al. was implemented using a 0.35 micrometer (μm) conventional CMOS fabrication process. In this process, various design rules dictate minimum separation distances between features. For example, according to the 0.35 μm CMOS design rules, with reference to
In sum, the ISFET pixel design of Milgrew et al. is aimed at ensuring accurate hydrogen ion concentration measurements over a pH range of 1-14. To ensure measurement linearity, the source and body of each pixel's ISFET are electrically coupled together. To ensure a full range of pH measurements, a transmission gate S1 is employed in each pixel to transmit the source voltage of an enabled pixel. Thus, each pixel of Milgrew's array requires four transistors (p-channel ISFET, p-channel MOSFET, and two n-channel MOSFETs) and two separate n-wells (
As noted earlier, individual ISFETs and arrays of ISFETs similar to those discussed above have been employed as sensing devices in a variety of chemical and biological applications. In particular, ISFETs have been employed as pH sensors in the monitoring of various processes involving nucleic acids such as DNA. Some examples of employing ISFETs in various life-science related applications are given in the following publications, each of which is incorporated herein by reference: Massimo Barbaro, Annalisa Bonfiglio, Luigi Raffo, Andrea Alessandrini, Paolo Facci and Imrich Barák, “Fully electronic DNA hybridization detection by a standard CMOS biochip,” Sensors and Actuators B: Chemical, Volume 118, Issues 1-2, 2006, pp. 41-46; Toshinari Sakurai and Yuzuru Husimi, “Real-time monitoring of DNA polymerase reactions by a micro ISFET pH sensor,” Anal. Chem., 64(17), 1992, pp 1996-1997; S. Purushothaman, C. Toumazou, J. Georgiou, “Towards fast solid state DNA sequencing,” Circuits and Systems, vol. 4, 2002, pp. IV-169 to IV-172; S. Purushothaman, C. Toumazou, C. P. Ou, “Protons and single nucleotide polymorphism detection: A simple use for the Ion Sensitive Field Effect Transistor,” Sensors and Actuators B: Chemical, Vol. 114, no. 2, 2006, pp. 964-968; A. L. Simonian, A. W. Flounders, J. R. Wild, “FET-Based Biosensors for The Direct Detection of Organophosphate Neurotoxins,” Electroanalysis, Vol. 16, No. 22, 2004, pp. 1896-1906; C. Toumazou, S. Purushothaman, “Sensing Apparatus and Method,” United States Patent Application 2004-0134798, published Jul. 15, 2004; and T. W. Koo, S. Chan, X. Su, Z. Jingwu, M. Yamakawa, V. M. Dubin, “Sensor Arrays and Nucleic Acid Sequencing Applications,” United States Patent Application 2006-0199193, published Sep. 7, 2006.
In general, the development of rapid and sensitive nucleic acid sequencing methods utilizing automated DNA sequencers has significantly advanced the understanding of biology. The term “sequencing” refers to the determination of a primary structure (or primary sequence) of an unbranched biopolymer, which results in a symbolic linear depiction known as a “sequence” that succinctly summarizes much of the atomic-level structure of the sequenced molecule. “DNA sequencing” particularly refers to the process of determining the nucleotide order of a given DNA fragment. Analysis of entire genomes of viruses, bacteria, fungi, animals and plants is now possible, but such analysis generally is limited due to the cost and time required to sequence such large genomes. Moreover, present conventional sequencing methods are limited in terms of their accuracy, the length of individual templates that can be sequenced, and the rate of sequence determination.
Despite improvements in sample preparation and sequencing technologies, none of the present conventional sequencing strategies, including those to date that may involve ISFETs, has provided the cost reductions required to increase throughput to levels required for analysis of large numbers of individual human genomes. The ability to sequence many human genomes facilitates an analysis of the genetic basis underlying disease (e.g., such as cancer) and aging, for example. Some recent efforts have made significant gains in both the ability to prepare genomes for sequencing and to sequence large numbers of templates simultaneously. However, these and other efforts are still limited by the relatively large size of the reaction volumes needed to prepare templates that are detectable by these systems, as well as the need for special nucleotide analogues, and complex enzymatic or fluorescent methods to “read out” nucleotide sequence.
SUMMARY OF INVENTIONWe have recognized and appreciated that large arrays of ISFETs may be particularly configured and employed to facilitate DNA sequencing techniques based on monitoring changes in chemical processes, including DNA synthesis. More generally, Applicants have recognized and appreciated that large arrays of chemically-sensitive FETs (i.e., chemFETs) may be employed to detect and measure static and/or dynamic concentrations/levels of a variety of analytes (e.g., hydrogen ions, other ions, non-ionic molecules or compounds, etc.) in a host of chemical and/or biological processes (e.g., biological or chemical reactions, cell or tissue cultures or monitoring, neural activity, nucleic acid sequencing, etc.) in which valuable information may be obtained based on such analyte measurements.
Accordingly, various embodiments of the present disclosure are directed generally to inventive methods and apparatus relating to large scale FET arrays for measuring one or more analytes. In the various embodiments disclosed herein, FET arrays include multiple chemFETs, that act as chemical sensors. An ISFET, as discussed above, is a particular type of chemFET that is configured for ion detection, and ISFETs may be employed in various embodiments disclosed herein. Other types of chemFETs contemplated by the present disclosure include enzyme FETs (EnFETs) which employ enzymes to detect analytes. It should be appreciated, however, that the present disclosure is not limited to ISFETs and EnFETs, but more generally relates to any FET that is configured for some type of chemical sensitivity. As used herein, chemical sensitivity broadly encompasses sensitivity to any molecule of interest, including without limitation organic, inorganic, naturally occurring, non-naturally occurring, and synthetic chemical and biological compounds, such as ions, small molecules, polymers such as nucleic acids, proteins, peptides, polysaccharides, and the like.
According to yet other embodiments, the present disclosure is directed generally to inventive methods and apparati relating to the use of the above-described large scale chemFET arrays in the analysis of chemical or biological samples. These samples are typically liquid (or are dissolved in a liquid) and of small volume, to facilitate high-speed, high-density determination of analyte (e.g., ion or other constituent) presence, concentration or other measurements on the analyte.
For example, some embodiments are directed to a “very large scale” two-dimensional chemFET sensor array (e.g., greater than 256 sensors), in which one or more chemFET-containing elements or “pixels” constituting the sensors of such an array are configured to monitor one or more independent biological or chemical reactions or events occurring in proximity to the pixels of the array. In some exemplary implementations, the array may be coupled to one or more microfluidics structures that form one or more reaction chambers, or “wells” or “microwells,” over individual sensors or groups of sensors of the array, and apparatus which delivers analyte samples (i.e., analyte solutions) to the wells and removes them from the wells between measurements. Even when microwells are not employed, the sensor array may be coupled to one or more microfluidics structures for the delivery of one or more analytes to the pixels and for removal of analyte(s) between measurements. Accordingly, inventive aspects of this disclosure, which are desired to be protected, include the various microfluidic structures which may be employed to flow analytes and where appropriate other agents useful in for example the detection and measurement of analytes to and from the wells or pixels, the methods of manufacture of the array of wells, methods and structures for coupling the arrayed wells with arrayed pixels, and methods and apparatus for loading the wells with sample to be analyzed, including for example loading the wells with DNA-bearing beads when the apparatus is used for DNA sequencing or related analysis, as will be discussed in greater detail below.
In association with the microfluidics, unique reference electrodes and their coupling to the flow cell are also shown.
In various aspects of the invention, an analyte of particular interest is a byproduct of nucleic acid synthesis. Such a byproduct can be monitored as the readout of a sequencing-by-synthesis method. One particularly important byproduct is inorganic pyrophosphate (PPi) which is released upon the addition (or incorporation) of a deoxynucleotide triphosphate (also referred to herein as dNTP) to the 3′ end of a nucleic acid (such as a sequencing primer). PPi may be hydrolyzed to orthophosphate (Pi) and free hydrogen ion (H+) in the presence of water (and optionally and far more rapidly in the presence of pyrophosphatase). As a result, nucleotide incorporation, and thus a sequencing-by-synthesis reaction, can be monitored by detecting PPi, Pi and/or H+. Conventionally, PPi has not been detected or measured by chemFETs. Optically based sequencing-by-synthesis methods have detected PPi via its sulfurylase-mediated conversion to adenosine triphosphate (ATP), and then luciferase-mediated conversion of luciferin to oxyluciferin in the presence of the previously generated ATP, with concomitant release of light. Such detection is referred to herein as “enzymatic” detection of PPi. The invention provides methods for detecting PPi using non-enzymatic methods. As used herein, non-enzymatic detection of PPi is detection of PPi that does not require an enzyme other than any enzyme required to produce or release the PPi in the first instance (e.g., a polymerase). An example of non-enzymatic detection of PPi is a detection method that does not require conversion of PPi to ATP.
H+ has been detected by measuring pH changes using standard pH meters or in some instances ISFETs. Importantly, many if not all prior attempts to detect nucleotide incorporation using ISFETs have focused solely on pH changes and not on detection or measurement of released PPi.
The instant invention contemplates and thus provides methods for monitoring nucleic acid sequencing reactions and thus determining the nucleotide sequence of nucleic acids by detecting H+ (or changes in pH), PPi (or Pi) in the absence of presence of PPi (or Pi) specific receptors, unincorporated dNTP, or some combination thereof. Some aspects therefore are aimed at monitoring pH while others are aimed at monitoring (and detecting) ion pulses at the FET surface resulting from changes in ionic species in the solution in contact with such surface.
Thus, various aspects of the invention provide methods and apparati for directly detecting released PPi as an indicator of nucleotide incorporation into a nucleic acid. The invention does so through the use of the chemFET arrays described herein. In some embodiments, a nucleic acid synthesis reaction is performed in a solution that is in contact with the chemFET and the released PPi is detected (or sensed) by the chemFET surface. Importantly, given the ability of the chemFETs of the invention to detect PPi directly, such synthesis and/or sequencing reactions may be performed, and in some instances are preferably performed, in an environment that is substantially pH insensitive (i.e., an environment in which changes in pH are not detected due to for example the strong buffering capacity of the environment).
In other embodiments, the analyte of interest is hydrogen ion, and large scale ISFET arrays according to the present disclosure are specifically configured to measure changes in H+ concentration (i.e., changes in pH).
In other embodiments, other biological or chemical reactions may be monitored, and the chemFET arrays may be specifically configured to measure hydrogen ions and/or one or more other analytes that provide relevant information relating to the occurrence and/or progress of a particular biological or chemical process of interest.
In various aspects, the chemFET arrays may be fabricated using conventional CMOS (or biCMOS or other suitable) processing technologies, and are particularly configured to facilitate the rapid acquisition of data from the entire array (scanning all of the pixels to obtain corresponding pixel output signals).
With respect to analyte detection and measurement, it should be appreciated that in various embodiments discussed in greater detail below, one or more analytes measured by a chemFET array according to the present disclosure may include any of a variety of biological or chemical substances that provide relevant information regarding a biological or chemical process (e.g., binding events such as hybridization of nucleic acids to each other, antigen-antibody binding, receptor-ligand binding, enzyme-inhibitor binding, enzyme-substrate binding, and the like). In some aspects, the ability to measure absolute or relative as well as static and/or dynamic levels and/or concentrations of one or more analytes, in addition to merely determining the presence or absence of an analyte, provides valuable information in connection with biological and chemical processes. In other aspects, mere determination of the presence or absence of an analyte or analytes of interest may provide valuable information may be sufficient.
A chemFET array according to various inventive 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 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 embodiment, the first and second analytes may be related to each other. As an example, the first and second analytes may be byproducts of the same biological or chemical reaction/process and therefore they may be detected concurrently to confirm the occurrence of a reaction (or lack thereof). Such redundancy is preferably in some analyte detection methods. 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 optionally to monitor biological or chemical processes such as binding events. 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 thereby consist of chemFETs of substantially similar or identical type that detect and/or measure the same 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. In another embodiment, the sensors in an array may be configured to detect and/or measure a single type (or class) of analyte even though the species of that type (or class) detected and/or measured may be different between sensors. As an example, all the sensors in an array may be configured to detect and/or measure nucleic acids, but each sensor detects and/or measures a different nucleic acid.
In yet other aspects, Applicants have specifically improved upon the ISFET array design of Milgrew et al. discussed above in connection with
With respect to chemFET array fabrication, Applicants have further recognized and appreciated that various techniques employed in a conventional CMOS fabrication process, as well as various post-fabrication processing steps (wafer handling, cleaning, dicing, packaging, etc.), may in some instances adversely affect performance of the resulting chemFET array. For example, with reference again to
Accordingly, one embodiment of the present invention is directed to an apparatus, comprising an array of CMOS-fabricated sensors, each sensor comprising one chemically-sensitive field effect transistor (chemFET) and occupying an area on a surface of the array of 10 μm2 or less.
Another embodiment is directed to a sensor array, comprising a two-dimensional array of electronic sensors including at least 512 rows and at least 512 columns of the electronic sensors, each sensor comprising one chemically-sensitive field effect transistor (chemFET) configured to provide at least one output signal representing a presence and/or concentration of an analyte proximate to a surface of the two-dimensional array.
Another embodiment is directed to an apparatus, comprising an array of CMOS-fabricated sensors, each sensor comprising one chemically-sensitive field effect transistor (chemFET). The array of CMOS-fabricated sensors includes more than 256 sensors, and a collection of chemFET output signals from all chemFETs of the array constitutes a frame of data. The apparatus further comprises control circuitry coupled to the array and configured to generate at least one array output signal to provide multiple frames of data from the array at a frame rate of at least 1 frame per second. In one aspect, the frame rate may be at least 10 frames per second. In another aspect, the frame rate may be at least 20 frames per second. In yet other aspects, the frame rate may be at least 30, 40, 50, 70 or up to 100 frames per second.
Another embodiment is directed to an apparatus, comprising an array of CMOS-fabricated sensors, each sensor comprising a chemically-sensitive field effect transistor (chemFET). The chemFET comprises a floating gate structure, and a source and a drain having a first semiconductor type and fabricated in a region having a second semiconductor type, wherein there is no electrical conductor that electrically connects the region having the second semiconductor type to either the source or the drain.
Another embodiment is directed to an apparatus, comprising an array of electronic sensors, each sensor consisting of three field effect transistors (FETs) including one chemically-sensitive field effect transistor (chemFET).
Another embodiment is directed to an apparatus, comprising an array of electronic sensors, each sensor comprising three or fewer field effect transistors (FETs), wherein the three or fewer FETs includes one chemically-sensitive field effect transistor (chemFET).
Another embodiment is directed to an apparatus, comprising an array of electronic sensors, each sensor comprising a plurality of field effect transistors (FETs) including one chemically-sensitive field effect transistor (chemFET), and a plurality of electrical conductors electrically connected to the plurality of FETs, wherein the plurality of FETs are arranged such that the plurality of electrical conductors includes no more than four conductors traversing an area occupied by each sensor and interconnecting multiple sensors of the array.
Another embodiment is directed to an apparatus, comprising an array of CMOS-fabricated sensors, each sensor comprising a plurality of field effect transistors (FETs) including one chemically-sensitive field effect transistor (chemFET), wherein all of the FETs in each sensor are of a same channel type and are implemented in a single semiconductor region of an array substrate.
Another embodiment is directed to a sensor array, comprising a plurality of electronic sensors arranged in a plurality of rows and a plurality of columns. Each sensor comprises one chemically-sensitive field effect transistor (chemFET) configured to provide at least one and in some instances at least two output signals representing a presence and/or a concentration of an analyte proximate to a surface of the array. For each column of the plurality of columns, the array further comprises column circuitry configured to provide a constant drain current and a constant drain-to-source voltage to respective chemFETs in the column, the column circuitry including two operational amplifiers and a diode-connected FET arranged in a Kelvin bridge configuration with the respective chemFETs to provide the constant drain-to-source voltage.
Another embodiment is directed to a sensor array, comprising a plurality of electronic sensors arranged in a plurality of rows and a plurality of columns. Each sensor comprises one chemically-sensitive field effect transistor (chemFET) configured to provide at least one output signal and in some instances at least two output signals representing a concentration of ions in a solution proximate to a surface of the array. The array further comprises at least one row select shift register to enable respective rows of the plurality of rows, and at least one column select shift register to acquire chemFET output signals from respective columns of the plurality of columns.
Another embodiment is directed to an apparatus, comprising an array of CMOS-fabricated sensors, each sensor comprising a chemically-sensitive field effect transistor (chemFET). The chemFET comprises a floating gate structure, and a source and a drain having a first semiconductor type and fabricated in a region having a second semiconductor type, wherein there is no electrical conductor that electrically connects the region having the second semiconductor type to either the source or the drain. The array includes a two-dimensional array of at least 512 rows and at least 512 columns of the CMOS-fabricated sensors. Each sensor consists of three field effect transistors (FETs) including the chemFET, and each sensor includes a plurality of electrical conductors electrically connected to the three FETs. The three FETs are arranged such that the plurality of electrical conductors includes no more than four conductors traversing an area occupied by each sensor and interconnecting multiple sensors of the array. All of the FETs in each sensor are of a same channel type and implemented in a single semiconductor region of an array substrate. A collection of chemFET output signals from all chemFETs of the array constitutes a frame of data. The apparatus further comprises control circuitry coupled to the array and configured to generate at least one array output signal to provide multiple frames of data from the array at a frame rate of at least 20 frames per second.
Another embodiment is directed to a method for processing an array of CMOS-fabricated sensors, each sensor comprising a chemically-sensitive field effect transistor (chemFET). The method comprises: A) dicing a semiconductor wafer including the array to form at least one diced portion including the array; and B) performing a forming gas anneal on the at least one diced portion.
Another embodiment is directed to a method for manufacturing an array of chemFETs. The method comprises: fabricating an array of chemFETs; depositing on the array a dielectric material; applying a forming gas anneal to the array before a dicing step; dicing the array; and applying a forming gas anneal after the dicing step. The method may further comprise testing the semiconductor wafer between one or more deposition steps.
Another embodiment is directed to a method for processing an array of CMOS-fabricated sensors. Each sensor comprises a chemically-sensitive field effect transistor (chemFET) having a chemically-sensitive passivation layer of silicon nitride and/or silicon oxynitride deposited via plasma enhanced chemical vapor deposition (PECVD). The method comprises depositing at least one additional passivation material on the chemically-sensitive passivation layer so as to reduce a porosity and/or increase a density of the passivation layer.
Other aspects of the invention relate to methods for monitoring nucleic acid synthesis reactions, including but not limited to those integral to sequencing-by-synthesis methods. Thus, various aspects of the invention provide methods for monitoring nucleic acid synthesis reactions, methods for determining or monitoring nucleotide incorporation into a nucleic acid, methods for determining the presence or absence of nucleotide incorporation, methods for determining the number of incorporated nucleotides, and the like.
In one such aspect, the invention provides a method for sequencing a nucleic acid comprising disposing (e.g., placing or positioning) a plurality of template nucleic acids into a plurality of reaction chambers, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer (thereby forming a template/primer hybrid) and is bound to a polymerase, synthesizing a new nucleic acid strand (or extending the sequencing primer) by incorporating one or more known nucleotide triphosphates (also referred to generically herein as dNTP) sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by a change in voltage and/or current at the at least one chemFET. The array (and thus the plurality) of chemFETs may be comprised of at least two, at least three, at least four, at least five, at least six, at least seven, at least eight, at least nine, at least ten, at least 100, at least 200, at least 300, at least 400, at least 500, at least 1000, at least 104, at least 105, at least 106, at least 107, or more chemFET (or chemFET sensors, or sensors, as the terms are used interchangeably herein). Similarly, the plurality of reaction chambers may be at least two, at least three, at least four, at least five, at least six, at least seven, at least eight, at least nine, at least ten, at least 100, at least 200, at least 300, at least 400, at least 500, at least 1000, at least 104, at least 105, at least 106, at least 107, or more reaction chambers.
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, wherein the plurality of reaction chambers is in contact with a chemFET array, wherein at least one chemFET is in contact with each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer (thereby forming a template/primer hybrid) and is bound to a polymerase, synthesizing a new nucleic acid strand (or extending the sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by a change in voltage and/or current at the at least one chemFET, wherein the chemFET array is any of the foregoing arrays.
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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer (thereby forming a template/primer hybrid) and is bound to a polymerase, synthesizing a new nucleic acid strand (by extending the sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by a change in voltage and/or current at the at least one chemFET within the array, wherein a center-to-center distance between adjacent reaction chambers (or “pitch”) is 1-10 μm.
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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer (thereby forming a template/primer hybrid) and is bound to a polymerase, synthesizing a new nucleic acid strand (by extending the sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by the generation of sequencing reaction byproduct, wherein (a) the chemFET array comprises more than 256 sensors, or (b) a center-to-center distance between adjacent reaction chambers is 1-10 μm. In an important embodiment, the sequencing reaction byproduct is PPi.
In yet another aspect, the invention provides a method for sequencing a nucleic acid comprising fragmenting (or isolating, for example in the context of enriched exon isolation) a target nucleic acid to generate a plurality of fragmented (or isolated) nucleic acids, attaching each of the plurality of fragmented (or isolated) nucleic acids to individual beads to generate a plurality of beads each attached to a single fragmented (or isolated) nucleic acid, amplifying the number of fragmented (or isolated) nucleic acids on each bead, delivering the plurality of beads attached to amplified fragmented (or isolated) nucleic acids to a chemFET array having a separate reaction chamber for each sensor in the array, wherein only one bead is situated in each reaction chamber, and performing simultaneous sequencing reactions in the plurality of reaction chambers.
It is to be understood, as described in greater detail herein, that the nucleic acids to be sequenced may be derived from longer nucleic acids that are then subsequently fragmented (i.e., converted into shorter nucleic acids) or they may be isolated at a length that is suitable to the reactions contemplated herein and thus would not have to be shortened (or fragmented). It should therefore be understood that for every aspect and limitation discussed herein that refers to the process of fragmenting a target nucleic acid, the method can also be carried out by isolating a target nucleic acid in the absence of fragmenting.
In still another aspect, the invention provides a method for sequencing nucleic acids comprising fragmenting a target nucleic acid to generate a plurality of fragmented nucleic acids, amplifying each fragmented nucleic acid separately in the presence of a bead and binding amplified copies of the fragmented nucleic acid to the bead, thereby producing a plurality of beads each attached to a plurality of identical fragmented nucleic acids, delivering the plurality of beads each attached to a plurality of identical fragmented nucleic acids to a chemFET array having a separate reaction chamber for each chemFET sensor in the array, wherein only one bead is situated in each reaction chamber, and performing simultaneous sequencing reactions in the plurality of reaction chambers.
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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer (thereby forming a template/primer hybrid) and is bound to a polymerase, synthesizing a new nucleic acid strand (by extending the sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting a change in the level of a sequencing byproduct as an indicator of incorporation of the one or more known nucleotide triphosphates.
The change in the level may an increase or a decrease in a level relative to a level prior to incorporation of the one or more known nucleotide triphosphates. The change in the level may be read as a change in voltage and/or current at a chemFET sensor or a change in pH, but it is not so limited.
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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, synthesizing a new nucleic acid strand by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, directly detecting release of PPi as an indicator of incorporation of the one or more known nucleotide triphosphates.
In another aspect, the invention provides a method for sequencing nucleic acids comprising fragmenting a template nucleic acid to generate a plurality of fragmented nucleic acids, attaching one strand from each of the plurality of fragmented nucleic acids individually to beads to generate a plurality of beads each having a single stranded fragmented nucleic acid attached thereto, delivering the plurality of beads having a single stranded fragmented nucleic acid attached thereto to a chemFET array having a separate reaction chamber for each sensor in the area, and wherein only one bead is situated in each reaction chamber, and performing simultaneous sequencing reactions in the plurality of chambers.
In one aspect, the invention provides a method for sequencing a nucleic acid comprising sequencing a plurality of identical template nucleic acids in a reaction chamber in contact with a chemFET, in an array which comprises at least 3 (and up to millions) of such assemblies of reaction chambers and chemFETs.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising detecting incorporation of one or more known nucleotide triphosphates to a 3′ end of a sequencing primer hybridized to a template nucleic acid in a reaction chamber in contact with a chemFET of a chemFET array that comprises at least three chemFET.
In one aspect, the invention provides a method for sequencing a nucleic acid comprising fragmenting a target nucleic acid to generate a plurality of fragmented nucleic acids, individually amplifying one or more of the plurality of fragmented nucleic acids, and sequencing the individually amplified fragmented nucleic acids using a chemFET array. In one embodiment, the chemFET array comprises at least three chemFETs. In some embodiments, the chemFET array comprises at least 500 chemFETs, or at least 100,000 chemFETs. In some embodiments, the plurality of fragmented nucleic acids is individually amplified using a water in oil emulsion amplification method.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising fragmenting a target nucleic acid to generate a plurality of fragmented nucleic acids, attaching each of the plurality of fragmented nucleic acids to individual beads to generate a plurality of beads each attached to a single fragmented nucleic acid, amplifying the fragmented nucleic acids on each bead resulting in a plurality of identical fragmented nucleic acids on each bead, delivering a plurality of beads attached to fragmented nucleic acids to a chemFET array having a separate reaction chamber for each sensor in the array, wherein only one bead is situated in each reaction chamber, and performing sequencing reactions simultaneously in the plurality of reaction chambers.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising fragmenting a target nucleic acid to generate a plurality of fragmented nucleic acids, individually amplifying one or more of the plurality of fragmented nucleic acids, and sequencing the individually amplified fragmented nucleic acids in a plurality of reaction chambers having a center-to-center distance of about 1-10 μm using a chemFET array. In various embodiments, the center-to-center distance is about 9 μm, about 5 μm, or about 2 μm.
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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, 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 by a change in voltage at the at least one chemFET within the array in the absence of a detectable pH change.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of template nucleic acids into a reaction chamber, wherein the plurality of template nucleic acids is attached to a single bead, each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, and the reaction chamber is in contact with a chemFET, 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 by detection of a first and a second voltage pulse at the chemFET.
In still 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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, 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 by non-enzymatically detecting released inorganic pyrophosphate. In one embodiment, the detecting step occurs in the absence of a detectable pH change.
In still 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, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, 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 by detecting released inorganic pyrophosphate and unincorporated nucleotide triphosphates.
In one embodiment, the released inorganic pyrophosphate is detected at t0 and unincorporated nucleotide triphosphates are detected at time t1. In a related embodiment, a time difference of t1−t0 indicates a number of known nucleotide triphosphates incorporated.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising contacting a template nucleic acid with a sequencing primer and a polymerase for times and conditions sufficient to allow the template nucleic acid to bind to the sequencing primer for form a template/primer hybrid and to allow the polymerase to bind to the template/primer hybrid, and synthesizing a new nucleic acid strand by incorporating nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting incorporation of the one or more known nucleotide triphosphates by detecting released inorganic pyrophosphate and unincorporated nucleotide triphosphates.
In one embodiment, the released inorganic pyrophosphate is detected at to and unincorporated nucleotide triphosphates are detected at time t1. In a related embodiment, the time difference between t1 and t0 (i.e., t1−t0) indicates the number of known nucleotide triphosphates incorporated.
In still another aspect, the invention provides a method for sequencing a nucleic acid comprising contacting a template nucleic acid with a sequencing primer and a polymerase for times and conditions sufficient to allow the template nucleic acid to bind to the sequencing primer for form a template/primer hybrid and to allow the polymerase to bind to the template/primer hybrid, and synthesizing a new nucleic acid strand (or extending sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer in a low ionic strength environment, and detecting incorporation of the one or more known nucleotide triphosphates by detection of one or more voltage pulses at the chemFET.
Depending on the embodiment, the low ionic strength environment may comprise less than 1 mM MgCl2, less than 0.5 mM MgCl2, less than 100 μM MgCl2, or less than 50 μM MgCl2. In other embodiments, the low ionic strength environment may comprise less than 1 mM MnCl2, less than 0.5 mM MnCl2, less than 100 μM MnCl2, or less than 50 μM MnCl2. In still other embodiments, the salt may be another Mg2+ or Mn2+ containing salt, or it may be a Ca2+ or Co2+ containing salt, or it may be a combination of one or more such salts.
In another aspect, the invention provides a method for determining incorporation of a nucleotide triphosphate into a newly synthesized nucleic acid (or onto a primer such as a sequencing primer) comprising combining a nucleotide triphosphate, a template/primer hybrid, and a polymerase, in a solution in contact with a chemFET, and detecting voltage pulses at the chemFET, wherein detection of a first and a second voltage pulse indicates incorporation of a nucleotide triphosphate and wherein detection of a first but not a second voltage pulse indicates lack of incorporation of a nucleotide triphosphate.
In one embodiment, the voltage pulses are detected at the chemFET independent of a binding event at the passivation layer of the chemFET.
In another aspect, the invention provides a method for determining incorporation of a nucleotide triphosphate into a newly synthesized nucleic acid (or onto a primer such as a sequencing primer) comprising combining a nucleotide triphosphate, a template/primer hybrid, and a polymerase, in a solution in contact with a chemFET, and detecting voltage pulses at the chemFET independent of a binding event at the passivation layer of the chemFET, wherein detection of a first and a second voltage pulse indicates incorporation of at least one nucleotide triphosphate.
In various embodiments of the some of the foregoing aspect, the first voltage pulse occurs at time t0 and the second voltage pulse occurs at time t1, and t1−t0 indicates the number of nucleotide triphosphates incorporated. In various embodiments, the incorporated nucleotide triphosphate is known. In various embodiments, the nucleotide triphosphate is a plurality of identical nucleotide triphosphates, the template/primer hybrid is a plurality of identical template/primer hybrids, and the polymerase is a plurality of identical polymerases. Alternatively, the polymerase may be a plurality of polymerases that are not identical and that rather may be comprised of 2, 3, or more types of polymerase. In some instances, a mixture of two polymerases may be used with one having suitable processivity and the other having suitable rate of incorporation. The ratio of the different polymerases can vary and the invention is not to be limited in this regard. Similarly the template/primer hybrid may be a plurality of template/primer hybrids which may not be identical to each other, provided that any template/primer hybrids in a single reaction chamber or attached to a single bead are identical to each other.
In still another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing (e.g., placing or positioning) a plurality of template nucleic acids into a plurality of reaction chambers, wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, synthesizing a new nucleic acid strand (or extending the sequencing primer) 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 by detecting at the at least one chemFET a first voltage pulse at time t0 having height h0 and a second pulse at time t1 having height h1, wherein h0 and h1 are each at least about 5 mV over baseline, and t1−t0 is at least 1 millisecond. In some embodiments, t1−t0 is at least 5 millisecond, at least 10 milliseconds, at least 20 milliseconds, at least 30 milliseconds, at least 40 milliseconds, or at least 50 milliseconds.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of beads into a plurality of reaction chambers, wherein each reaction chamber comprises a single bead, each bead is attached to a plurality of identical template nucleic acids, each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, and wherein the plurality of reaction chambers is in contact with a chemFET array comprising at least one chemFET for each reaction chamber, synthesizing a new nucleic acid strand (or extending the sequencing primer) by incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by detecting a first and a second voltage pulse at at least one chemFET within the array, wherein the chemFET array comprises at least three chemFET.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of beads into a plurality of reaction chambers, wherein each reaction chamber comprises a single bead, each bead is attached to a plurality of identical template nucleic acids, each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, and wherein the plurality of reaction chambers is in contact with a chemical-sensitive field effect transistor (chemFET) array comprising at least one chemFET for each reaction chamber, initiating synthesis of a new nucleic acid strand (or extension of the sequencing primer) by introducing a plurality of known identical nucleotide triphosphates into each reaction chamber, 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 the incorporation of the one or more known nucleotide triphosphates by detecting a first and a second voltage pulse at at least one chemFET within the array.
In one embodiment, each reaction chamber comprises adenosine triphosphate prior to introducing a plurality of known triphosphates into each reaction chamber.
In another aspect, the invention provides a method for sequencing a nucleic acid comprising (a) disposing a plurality of beads into a plurality of reaction chambers, each reaction chamber comprising a single bead, each bead attached to a plurality of identical template nucleic acids, each of the template nucleic acids hybridized to a sequencing primer and bound to a polymerase, and each reaction chamber in contact with at least one chemFET, (b) introducing a plurality of known identical nucleotide triphosphates into each reaction chamber, (c) detecting sequential incorporation at the 3′ end of the sequencing primer of one or more nucleotide triphosphates if complementary to corresponding nucleotides in the template nucleic acid, (d) washing unincorporated nucleotide triphosphates from the reaction chambers, and (e) repeating steps (b) through (d) in the same reaction chamber using a different plurality of known nucleotide triphosphates.
In one embodiment, sequential incorporation of one or more nucleotide triphosphates is detected by a first and a second voltage pulse at the chemFET. In some embodiments, step (e) comprises repeating steps (b) through (d) by separately introducing each different plurality of known nucleotide triphosphates into each reaction chamber
In another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of beads into a plurality of reaction chambers, wherein each reaction chamber comprises a single bead, each bead is attached to a plurality of identical template nucleic acids, each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, and wherein the each of the plurality of reaction chambers is in contact with at least one chemFET in a chemFET array, initiating synthesis of a new nucleic acid strand by introducing a divalent cation into each reaction chamber, 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 the incorporation of the one or more known nucleotide triphosphates by detecting a first and a second voltage pulse at at least one chemFET within the array.
In one embodiment, the divalent cation is Mg2+. In another embodiment, the divalent cation is Mn2+. In still a further embodiment, the divalent cation is a mixture of Mg2+ and Mn2+. Depending on the embodiment, the divalent cation is at a concentration of less than 1 mM, or less than 100 μM, or about 50 μM.
In yet another aspect, the invention provides a method for sequencing a nucleic acid comprising disposing a plurality of identical template nucleic acids into a reaction chamber, wherein the reaction chamber is in contact with a chemFET, and wherein each of the template nucleic acids is hybridized to a sequencing primer and is bound to a polymerase, incorporating one or more known nucleotide triphosphates sequentially at the 3′ end of the sequencing primer, and detecting the incorporation of the one or more known nucleotide triphosphates by a voltage pulse at the chemFET, wherein the incorporation of 10-1000 nucleotide triphosphates is detected.
Various embodiments may be embraced in the various foregoing aspects of the invention and these are recited below once for convenience and brevity.
It is to be understood that although various of the foregoing aspects and embodiments of the invention recite hybridization (or binding) of a sequencing primer to a template, the invention also contemplates the use of sequencing primers and/or template nucleic acids that hybridize to themselves (i.e., intramolecularly) thereby giving rise to free 3′ ends onto which nucleotide triphosphates may be incorporated. Such templates, referred to herein as self-priming templates, may be used in any of the foregoing methods.
Similarly, the invention equally contemplates the use of double stranded templates that are engineered to have particular sequences at their free ends that can be acted upon by nicking enzymes such as DNA nickase. In this way, the polymerase incorporates nucleotide triphosphates at the nicked site. In these instances, there is no requirement for a separate sequencing primer.
In some embodiments, the incorporation of 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, or at least 100 nucleotide triphosphates is detected. In other embodiments, the incorporation of 100-1000 nucleotide triphosphates is detected. In still other some embodiments, the incorporation of 250-750 nucleotide triphosphates is detected.
In some embodiments, the reaction chamber comprises a plurality of packing beads. In some embodiments, the detecting step occurs in the presence of a plurality of packing beads.
In some embodiments, the reaction chamber comprises a soluble non-nucleic acid polymer. In some embodiments, the detecting step occurs in the presence of a soluble non-nucleic acid polymer. In some embodiments, the soluble non-nucleic acid polymer is polyethylene glycol, or PEA, or a dextran, or an acrylamide, or a cellulose (e.g., methyl cellulose). In some embodiments, the non-nucleic acid polymer such as polyethylene glycol is attached to the single bead. In some embodiments, the non-nucleic acid polymer is attached to one or more (or all) sides of a reaction chamber, except in some instances the bottom of the reaction chamber which is the FET surface In some embodiments, the non-nucleic acid polymer is biotinylated such as but not limited to biotinylated polyethylene glycol.
In some embodiments, the method is carried out at a pH of about 7-9, or at about 8.5 to 9.5, or at about 9.
In some embodiments, the synthesizing and/or detecting step is carried out in a weak buffer. In some embodiments, the weak buffer comprises Tris-HCl, boric acid or borate buffer, acetate, morpholine, citric acid, carbonic acid, or phosphoric acid as a buffering agent.
In some embodiments, the synthesizing and/or detecting step is carried out in about 1 mM Tris-HCl. In some embodiments, the synthesizing and/or detecting step is carried out in less than 1 mM Tris-HCl. In some embodiments, the synthesizing and/or detecting step is carried out in about 0.9 mM Tris-HCl, about 0.8 mM Tris-HCl, about 0.7 mM Tris-HCl, about 0.6 mM Tris-HCl, about 0.5 mM Tris-HCl, about 0.4 mM Tris-HCl, about 0.3 mM Tris-HCl, or about 0.2 mM Tris-HCl.
In some embodiments, the synthesizing and/or detecting step is carried out in about 1 mM borate buffer. In some embodiments, the synthesizing and/or detecting step is carried out in less than 1 mM borate buffer. In some embodiments, the synthesizing and/or detecting step is carried out in about 0.9 mM borate buffer, about 0.8 mM borate buffer, about 0.7 mM borate buffer, about 0.6 mM borate buffer, about 0.5 mM borate buffer, about 0.4 mM borate buffer, about 0.3 mM borate buffer, or about 0.2 mM borate buffer.
In some embodiments, the detection step occurs in the absence of a detectable pH change. In some embodiments, the detection step is carried out in an environment of constant pH. In some embodiment, the chemFET is relatively pH insensitive. In some embodiments, the synthesizing and/or detecting step is carried out in about 0.5 mM Tris-HCl.
In some embodiments, the synthesizing and/or detecting step is carried out in a low ionic strength environment. In some embodiments, synthesizing and/or detecting step is carried out in a low ionic strength environment that comprises less than 1 mM MgCl2 or MnCl2, less than 0.5 mM MgCl2 or MnCl2, or less than 100 μM MgCl2 or MnCl2. In some embodiments, the synthesizing and/or detecting step is carried out in a low ionic strength environment that comprises about 100 μM MgCl2 or MnCl2, about 75 μM MgCl2 or MnCl2, about 50 μM MgCl2 or MnCl2, about 40 μM MgCl2 or MnCl2, about 30 μM MgCl2 or MnCl2, about 20 μM MgCl2 or MnCl2, or about 10 μM MgCl2 or MnCl2. It is to be understood that the invention contemplates similarly amounts of other divalent cations and it is therefore not limited to MgCl2 or MnCl2 only. Similarly it should be understood that the invention contemplates the use of two or more divalent cations (and/or their corresponding salts) provided the total ion concentration is at the levels stated above.
In some embodiments, the synthesizing and/or detecting step is carried out in about 0.5 mM TRIS and about 50 μM MgCl2 or MnCl2.
In various embodiments, the nucleotide triphosphates are unblocked. As used herein, an unblocked nucleotide triphosphate is a nucleotide triphosphate with an unmodified end that can be incorporated into a nucleic acid (at its 3′ end) and once it is incorporated can be attached to the following nucleotide triphosphate being incorporated. Blocked dNTP in contrast either cannot be added to a nucleic acid or their incorporation into a nucleic acid prevents any further nucleotide incorporation. In various embodiments, the nucleotide triphosphates are deoxynucleotide triphosphates (dNTPs).
In various embodiments, the chemFET comprises a silicon nitride passivation layer. In some embodiments, the chemFET comprises a passivation layer attached to inorganic pyrophosphate (PPi) receptors. In some embodiments, the chemFET comprises a passivation layer that is not bound to a nucleic acid.
In some embodiments, each reaction chamber is in contact with a single chemFET.
In some embodiments, the reaction chamber has a volume of equal to or less than about 1 picoliter (pL), including less than 0.5 pL, less than 0.1 pL, less than 0.05 pL, less than 0.01 pL, less than 0.005 pL. In some embodiment, the reaction chambers are separated by a center-to-center spacing of 1-10 μm. In some embodiments, the reaction chambers are separated by a center-to-center spacing of about 9 μm, about 5 microns, or about 2 microns. The reaction chambers may have a square cross section, for example at their base or bottom. Examples include an 8 μm by 8 μm cross section, a 4 μm by 4 μm cross section, or a 1.5 μm by 1.5 μm cross section. Alternatively, they may have a rectangular cross section, for example at their base or bottom. Examples include an 8 μm by 12 μm cross section, a 4 μm by 6 μm cross section, or a 1.5 μm by 2.25 μm cross section.
In some embodiments, the nucleotide triphosphates are pre-soaked in Mg2+ (e.g., in the presence of MgCl2) or Mn2+ (for example in the presence of MnCl2). In some embodiments, the polymerase is pre-soaked in Mg2+ (e.g., in the presence of MgCl2) or Mn2+ (for example in the presence of MnCl2).
In some embodiments, the method is carried out in a reaction chamber comprising a single capture bead, wherein a ratio of reaction chamber width to single capture bead diameter is at least 0.7, at least 0.8, or at least 0.9.
In some embodiments, the polymerase is free in solution. In some embodiments, the polymerase is immobilized to a bead. In some embodiments, the polymerase is immobilized to a capture bead. In some embodiments, the template nucleic acids are attached to capture beads.
In some embodiments, inorganic pyrophosphate (PPi) is not significantly hydrolyzed prior to detecting nucleotide incorporation. In some embodiments, inorganic pyrophosphate (PPi) is not significantly hydrolyzed prior to detecting a first voltage pulse.
In one aspect, the invention provides an apparatus comprising a chemFET array comprising a plurality of chemFET sensors each having a PPi receptor disposed on its surface, wherein the array comprises at least three chemFET sensors.
In one aspect, the invention provides an apparatus comprising a chemFET array comprising a plurality of chemFET sensors each having a PPi receptor disposed on its surface, wherein adjacent chemFET sensors in the array are separated by a center-to-center distance of about less than about 10 μm (e.g., using 0.18 μm CMOS fabrication, the pitch may be about 2.8 μm).
In another aspect, the invention provides an apparatus comprising a chemical-sensitive field effect transistor (chemFET) having disposed on its surface a PPi selective receptor. The PPi receptor may be immobilized to the chemFET surface.
In still another aspect, the invention provides an apparatus comprising a chemical-sensitive field effect transistor (chemFET) array comprising a plurality of chemFET sensors each having a PPi receptor disposed on its surface, wherein the plurality is a subset of the chemFET sensors in the array. In a related embodiment, the array comprises at least three chemFET sensors. The subset may represent 10%, 25%, 33%, 50%, 75% or more of the sensors in the array. The PPi selective receptor may be immobilized to the surface of each chemFET.
In some embodiments, the center-to-center distance between adjacent reaction chambers is about 1-9 μm, or about 2-9 μm, about 1 μm, about 2 μm, about 3 μm, about 4 μm, about 5 μm, about 6 μm, about 7 μm, about 8 μm, or about 9 μm. In some embodiments, the chemFET array comprises more than 256 sensors (and optionally more than 256 corresponding reaction chambers (or wells)), more than 300 sensors (and optionally more than 300 corresponding reaction chambers), more than 400 sensors (and optionally more than 400 corresponding reaction chambers), more than 500 sensors (and optionally more than 500 corresponding reaction chambers), more than 600 sensors (and optionally more than 600 corresponding reaction chambers), more than 700 sensors (and optionally more than 700 corresponding reaction chambers), more than 800 sensors (and optionally more than 800 corresponding reaction chambers), more than 900 sensors (and optionally more than 900 corresponding reaction chambers), more than 103 sensors (and optionally more than 103 corresponding reaction chambers), more than 104 sensors (and optionally more than 104 corresponding reaction chambers), more than 105 sensors (and optionally more than 105 corresponding reaction chambers), or more than 106 sensors (and optionally more than 106 corresponding reaction chambers). In some embodiments, the chemFET array comprises at least 512 rows and at least 512 columns of sensors.
In some embodiments, the sequencing byproduct is PPi. In some embodiments, PPi is measured directly. In related embodiments, PPi is detected by its binding to a PPi selective receptor immobilized on the surface of a chemFET sensor. In another embodiment, PPi is detected in the absence of a PPi selective receptor. In other embodiments, PPi is detected in a pH independent or insensitive environment.
In some embodiments, the sequencing reaction byproduct is hydrogen ion. In some embodiments, the sequencing reaction byproduct is Pi. In still other embodiments, the chemFET detects changes in any combination of the byproducts, optionally in combination with other parameters, as described herein.
In some embodiments, the PPi selective receptor is Compound 1, Compound 2, Compound 3, Compound 4, Compound 5, Compound 6, Compound 7, Compound 8, Compound 9 or Compound 10 as shown in
In another aspect, the invention provides an apparatus comprising a chemFET array having disposed on its surface a biological array or a chemical array.
The biological array may be a nucleic acid array, a protein array including but not limited to an enzyme array, an antibody array and an antibody fragment array, a cell array, and the like. The chemical array may be an organic small molecule array, or an inorganic molecule array, but it is not so limited. The chemFET array may comprise at least 5, at least 10, at least 102, at least 103, at least 104, at least 105, at least 106, or more sensors. The biological or chemical array may be arranged into a plurality of “cells” or spatially defined regions, and each of these regions is situated over a different sensor in the chemFET array, in some embodiments.
In yet another aspect, the invention provides a method for detecting a nucleic acid comprising contacting a nucleic acid array disposed on a chemFET array with a sample, and detecting binding of a nucleic acid from the sample to one or more regions on the nucleic acid array.
In another aspect, the invention provides a method for detecting a protein comprising contacting a protein array disposed on a chemFET array with a sample, and detecting binding of a protein from the sample to one or more regions on the protein array.
In yet another aspect, the invention provides a method for detecting a nucleic acid comprising contacting a protein array disposed on a chemFET array with a sample, and detecting binding of a nucleic acid from the sample to one or more regions on the protein array.
In another aspect, the invention provides a method for detecting an antigen comprising contacting an antibody array disposed on a chemFET array with a sample, and detecting binding of an antigen from the sample to one or more regions on the antibody array.
In another aspect, the invention provides a method for detecting an enzyme substrate or inhibitor comprising contacting an enzyme array disposed on a chemFET array with a sample, and detecting binding of an entity from the sample to one or more regions on the enzyme array.
It should be appreciated that all combinations of the foregoing concepts and additional concepts discussed in greater detail below (provided such concepts are not mutually inconsistent) are contemplated as being part of the inventive subject matter disclosed herein. In particular, all combinations of claimed subject matter appearing at the end of this disclosure are contemplated as being part of the inventive subject matter disclosed herein. It should also be appreciated that terminology explicitly employed herein that also may appear in any disclosure incorporated by reference should be accorded a meaning most consistent with the particular concepts disclosed herein.
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.
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.
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.
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.
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.
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.
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 μm3 (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 22, 32, 42, 52, 62, 72, 82, 92 or 102 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.
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).
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.
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.
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.
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:
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.
Preferably, at least 106 base pairs are sequenced per hour, more preferably at least 107 base pairs are sequenced per hour, even more preferably at least 108 base pairs are sequenced per hour, even more preferably at least 109 base pairs are sequenced per hour, and most preferably at least 1010 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.
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.
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.
where cB is the molar concentration of ion B (mol dm−3), zB is the charge number of that ion, and the sum is taken over all ions in the solution. The Mg+2 or Mn2+ concentration may be 1000 μM, 500 μM, 200 μM, 100 μM, 50 μM, 10 μM, 5 μM, or even 1 μM.
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
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.
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.
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.
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.
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
As illustrated in
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.
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.
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.
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.
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.
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.
In sum, the flow cell 200 in the system of
In the system 1000 of
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
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/
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.
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.:
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 λ(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:
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:
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:
from which is obtained:
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.
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
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.
In
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
More generally, Applicants have recognized and appreciated that for any of the ISFET output signals in
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 Npoly=0 in
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.
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.
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 ΔΨ0 of the ion pulse).
First, it may be appreciated from Eq. (17) that a higher initial equilibrium surface potential Ψ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 ΔΨ0 for an ion pulse. In one aspect, with reference again to Eq. (16), Applicants have appreciated that a higher Ψ1 is facilitated by the selection of a wash buffer having a pH that is significantly different than the pHpzc 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.
From Eq. (15), Applicants also have appreciated that a relatively high pH sensitivity (ΔΨ0/ΔpH) of the ISFET also may contribute to a higher Ψ1; in this manner, as discussed above in connection with Eq. (15), a passivation layer material with a large intrinsic buffering capacity βint may be desirable, as the dimensionless sensitivity factor a given in Eq. (15) is increased by a large β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 ΔΨ0 for an ion pulse and hence an increased SNR of the ISFET output signal.
Additionally, Applicants have appreciated from Eq. (17) that reducing or minimizing the initial double layer capacitance Cdl,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 Cdl,2 arising from the concentration change stimulus (i.e., PPi generation or dNTP increase), also results in a higher amplitude ΔΨ0 for an ion pulse. Accordingly, in various implementations, low ionic strength buffers may be appropriately selected to facilitate an increased amplitude ΔΨ0 for an ion pulse and hence an increased SNR of the ISFET output signal.
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.
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.
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.
With respect to the ISFET array 100 of the system 1000 shown in
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
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
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.
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.
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.
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.
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 micrometers2) 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.
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 micrometers2), 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.
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).
In other aspects of the system shown in
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.
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.
As noted above, various inventive embodiments disclosed herein specifically improve upon the ISFET array design of Milgrew et al. discussed above in connection with
To this end,
In one aspect of the embodiment shown in
As illustrated in
By employing the diode-connected MOSFET Q6 in the bias/readout circuitry 110j of
In
In another aspect of the embodiment shown in
In yet another aspect of the embodiment shown in
By not tying the body connection of each ISFET to its source, the possibility of some non-zero source-to-body voltage VSB may give rise to the “body effect,” as discussed above in connection with
In the top view of
With reference now to the cross-sectional view of
In the composite cross-sectional view of
Above the substrate, gate oxide, and polysilicon layers shown in
As indicated above,
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
embodiment, pixels are designed such that there is a greater vertical distance between the signal lines 1121, 1141, 1161 and 1181, and the topmost metal layer 304 constituting the floating gate structure 170.
In the embodiment described immediately above, with reference again to
To this end, in another embodiment some via and metal layers are reconfigured such that the signal lines 1121, 1141, 1161 and 1181 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.
In
1161 and 1181) required to operate the pixel. One noteworthy difference between
With reference now to the cross-sectional view of
More specifically, as in the embodiment of
which the correspondence between the lettered top views of respective layers and the cross-sectional view of
Accordingly, by consolidating the signal lines 1121, 1141, 1161 and 1181 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
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
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.
Thus, the pixel chip layout designs respectively shown in
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 Cox of 4.5 fF/μm2. 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 cm2/V·s (i.e., 1.9E10 μm2/V·s). From Eq. (2) above, this results in an ISFET transconductance parameter β on the order of approximately 170 to 300 μA/V2. 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 IDj 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
With respect to the analyte-sensitive passivation layer 172 shown in
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.
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 VTH, making accurate measurements challenging.
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.
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
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 (Al2O3), tantalum oxide (Ta3O5), tin oxide (SnO2) and silicon dioxide (SiO2). 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.
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.
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
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.
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).
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.
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.
In some embodiments, the passivation layer and/or the layers and/or molecules coated thereon dictate the analyte specificity of the array readout.
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 (Si3N4), silicon oxynitride (Si2N2O), silicon oxide (SiO2), aluminum oxide (Al2O3), tantalum pentoxide (Ta2O5), tin oxide or stannic oxide (SnO2), and the like.
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.
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.
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.
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.
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).
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.
As another example, receptors for neurotoxins are described in Simonian Electroanalysis 2004, 16: 1896-1906.
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
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).
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
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.
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.
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.
Alternatively, the receptor may be covalently conjugated or crosslinked to the polymer (e.g., it may be “grafted” onto a functionalized polymer).
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.
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
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.
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.
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
In another embodiment, the topmost metal layer 304 of the ISFETs floating gate structure 170 shown in
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).
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).
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.
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.
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
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
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
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 μm2 to about 50 μm2 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.
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”).
In one aspect of the array 100 shown in
In yet another implementation of an array similar to that shown in
In
Regarding the column select shift registers 1941 and 1942, 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,
With reference again for the moment to
In the embodiment of
In one exemplary implementation, the switches of both the even and odd output drivers 1981 and 1982 (e.g., the switches 1912, 1914, . . . 191512 shown in
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 1981 and 1982 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.
For example, in one embodiment, with reference to
While the example above describes asymmetric switches 191 for the output drivers 1981 and 1982 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
In yet another embodiment directed to facilitating rapid settling of the bus 175 shown in
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.
In one aspect of the array design discussed above in connection with
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
As illustrated in
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 VBODY 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.
As also shown in
Regarding data acquisition from the array 100, in one embodiment the array controller 250 of
The array controller 250 of
In the embodiment of
In the example of
For example, with respect to the method for detecting nucleotide incorporation, discussed above in connection with
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
In
As discussed above in connection with
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
Regarding pixel settling time, with reference again to
ΔVPIX(t)=A(1−e−1/k), (PP)
where A is the difference (VCOLj−VCOLj-1) between two pixel voltage values and k is a time constant associated with a capacitance of the bus 175.
For purposes of the present discussion, pixel “settling time” tsettle is defined as the time t at which ΔVPIX(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 np, then the voltage at the settling time tsettle is given by ΔVPIX (tsettle)=A[1−(np/A)]. Substituting in Eq. (PP) and solving for tsettle yields
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.
Subsequently, in block 506 of
In block 510 of
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
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.
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.
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.
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.
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.
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.
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.
Having discussed several aspects of an exemplary ISFET array and an array controller according to the present disclosure,
In one aspect of the embodiment shown in
In particular,
For each of the first and second groups of rows, the array 100A of
In one exemplary implementation of the array 100A of
Like the array 100 of
As noted in
While the exemplary arrays discussed above in connection with
The array 100D of
The array 100E of
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 micrometers2), providing for significantly dense ISFET arrays.
In the embodiments of ISFET arrays discussed above, array pixels employ a p-channel ISFET, as discussed above in connection with
For example,
One of the primary differences between the n-channel ISFET pixel design of
In addition to the pixel designs shown in
In
In
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,
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
One example of a method for forming a microwell array is now discussed, starting with reference to
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.
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.
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 (
Other photolithographic approaches may be used for formation of the microwell array, of course, the foregoing being only one example.
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.
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.
Naturally, in a production environment, modifications will be appropriate.
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.
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.
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 SiO2 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.
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.”
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.
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
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
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
Using a 6 um (micron) thick layer of tetra-methyl-ortho-silicate (TEOS) as a SiO2-like layer for microwell formation,
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.
The tapered slopes to the sides of the microwells also may be used to advantage. Referring to
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.
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.
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 Ta2O5 may be preferred as an etch stop over Al2O3 because the PZC of Al2O3 is right at the pH being used (i.e., about 8.8) and, hence, right at the point of zero charge. In addition Ta2O5 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.
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.
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.
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., SiO2), 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 102 microns or less.
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., SiO2), 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 SiO2. 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.
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.
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.
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:
have connections suitable for interconnecting with a fluidics delivery system—e.g., via appropriately-sized tubing;
have appropriate head space above wells;
minimize dead volumes encountered by fluids;
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);
be configured to achieve uniform transit time of the flow over the array;
generate or propagate minimal bubbles in the flow over the wells;
be adaptable to placement of a removable reference electrode inside or as close to the flow chamber as possible;
facilitate easy loading of beads;
be manufacturable at acceptable cost; and
be easily assembled and attached to the chip package.
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.
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).
A first example of a suitable experiment apparatus 3410 incorporating such a fluidic interface is shown in
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).
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.
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.
The chip 3412 will in turn be mounted to a carrier 3430, for packaging and connection to connector pins 3432.
For ease of description, to discuss fabrication starting with
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
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
The other alignment mark or set of marks produced by pattern 4022 is used for alignment with a subsequent layer to be discussed.
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.
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
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.
Though not illustrated in
Another way to introduce the reference electrode is shown in
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.
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
Another configuration, shown in
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.
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.
By applying a metallization to bottom 42F280 of member 42F200, the reference electrode may be formed.
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
A variation on this idea is depicted in
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.
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.
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
Yet another alternative for a fluidics assembly, as shown in
Some of the foregoing alternative embodiments also may be implemented in a hybrid plastic/PDMS configuration. For example, as shown in
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.
An example is shown in
Nanoports may be secured over the nanoport fluidic holes to facilitate connection of input and output tubing.
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.
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
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.
By contrast, an externally inserted electrode requires extra fluid tubing to the inlet port, which requires additional fluid flow between cycles.
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.
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.
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.
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.
The resulting fluidic system connections to the ISFET device thus incorporate shortened input and output fluidic lines, which is desirable.
Still another example embodiment for a fluidic assembly is shown in
Still further examples of flow cell structures are shown in
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
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.
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.
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 100K, 500K, or 1M 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.
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.
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.
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.
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
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.
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.
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.
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+.
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
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.
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 t1 peaks are not distinguishable). This limit, for the case Npoly=0 vs Npoly=1 is given by the following: Consider a first peak Npoly=0 which has a width=W0 and standard deviation σW
In certain embodiments the chemFET (or ISFET) readouts are independent of (or insensitive to) changes in the concentration of Pi and/or H+. In these embodiments, the rate of PPi hydrolysis to Pi with the concomitant release of H+ 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−1 years−1, as compared to Biochemistry, 2002, 41:12025 which reports the rate in the presence of enzyme and Mg2+ to be on the order of about 103 s−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.
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.
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.
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.
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.
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.
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.
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.
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), co-polymers 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.
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.
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.
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.
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.
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.
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%.
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.
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.
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.
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.
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 MgCl2 or MnCl2 is 100 μM or less. In some embodiments the concentration of Tris is 0.5 mM or less.
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 MgCl2. The MgCl2 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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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 104, at least 105, at least 106, 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.
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.
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.
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 Mg2+ (for example in the form of MgCl2) or Mn2+ (for example in the form of MnCl2). Other divalent cations can also be used including but not limited to Ca2+, Co2+. 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.
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.
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., t0) and the second (unincorporated nucleotides) ion pulse (i.e., t1). 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.
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.
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.
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.
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.
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).
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.
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+) 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.
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.
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.
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.
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.
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.
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.
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).
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.
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.
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.
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.
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 SiO2 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.
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.
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.
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.
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.
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).
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, 103, 104 or more. The array itself may have any number of cells, including but not limited to at least 10, 102, 103, 104, 105, 106, 107, 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.
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.
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).
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.
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).
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.
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.
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.
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.
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.
The systems described herein can be used for sequencing unlabeled biological polymers without optical detection.
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.
Typically, each reaction chamber is capacitively coupled to a chemFET.
Preferably, each reaction chamber is no greater than about 0.39 pL in volume and about 49 μm2 surface aperture, and more preferably has an aperture no greater than about 16 μm2 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.
Typically, the reaction chambers comprise microfluidic wells.
Preferably, the apparatus is adapted to sequence at least 106 base pairs per hour, more preferably at least 107 base pairs per hour, and most preferably at least 108 base pairs per hour.
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.
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.
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 Mg2+ or Mn2+ or other divalent cation is no greater than 100 μM. Preferably, at least 106 base pairs are sequenced per hour, more preferably at least 107 base pairs are sequenced per hour, and most preferably at least 108 base pairs are sequenced per hour using the above-described method. Preferably the detecting is accomplished by an array of chemFETs as above described.
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.
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.
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.
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”).
Preferably, at least 106 base pairs are sequenced per hour, more preferably at least 107 base pairs are sequenced per hour, and most preferably at least 108 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.
The following examples are included for purposes of illustration and are not intended to limit the scope of the invention.
Examples Example 1. Bead PreparationBinding of Single-Stranded Oligonucleotides to Streptavidin-Coated Magnetic Beads.
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.
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×104 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.
Due to the strong covalent binding affinity of streptavidin for biotin (Kd˜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×105 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.
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.
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.
Incubation of Hybridized Templates/Primer with DNA Polymerase.
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.
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.
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.
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
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).
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).
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.
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.
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).
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.
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 ActivityExperiments 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.
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 (
In particular, the data shown in
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 (
Polymerase Affinity in Low Ionic Strength.
The affinity of a polymerase for a template/primer hybrid can be measured indirectly with our standard assay (
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 MgCl2, 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 MgCl2, 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 (
Therminator and Bst polymerases can extend primer in both standard reaction conditions (STD; 20 mM Tris-HCl, 10 mM (NH2)2SO4, 10 mM KCl, 2 mM MgSO4, 0.1% Triton X-100, pH 8.8 at 25° C.) and in low ionic strength conditions (LS; 0.5 mM Tris, 100 μM MgCl2, 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.
Primer Extension in Wells of Sensor.
Using an assay similar to the standard activity assay (
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.
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.
The Ions in Solution.
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:
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:
This is the situation for the reaction buffer, scaled to 0.5 mM Tris-HCl. It is assumed that all the Mg++ is complexed to the enzyme by presoaking. From the solution of the reaction and diffusion equation for the sequencing reaction, the ionic strength is calculated, I=½Σs qs2Cs, and the sensor dynamics is given by relationships set out above.
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−6 (5)
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−6 (5′)
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−6 (5″)
For a 0.5 mM Tris plus 1 mM MgCl2, 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−6 (5′″)
For a 5 mM Tris plus 1 mM MgCl2, 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−6 (5″″)
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−6 (5′″″)
The pH is 8.53 with the dNTP present and 8.97 without.
For 0.5 mM Tris plus 50 μM MgCl2, 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−6 (5″″″)
The pH is 8.53 with the dNTP present and 8.97 without.
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−6 (5′″″″)
Note that the signal is proportional to I−1/2, so it is desirable to keep the ionic strength of the reaction buffer solution as low as possible. Note that (I0+dI)−1/2=I0−1/2−½(dI/I0) I0−1/2. Thus, the relative change to the signal is −½ (dI/I0). By comparing Eqs. 5-5′″″″, one sees that dI is roughly independent of the buffer ionic strength. Thus, if the initial ionic strength, I0, is lower, the relative signal change is larger.
Reaction/Diffusion Calculation
The reactions that occur are
d[DNA]/dt=−kbst[dNTP][DNA]
d[PPi]/dt=kbst[dNTP][DNA]
d[dNTP]/dt=kbst[dNTP][DNA]
If the enzyme for the reaction PPi+H20→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+ released by adding one more base of DNA, however, does diffuse.
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++ on the reaction rate. In addition, the diffusion of Mg++ is also tracked. Recall that in the current calculations, however, Mg++ is assumed to be constant.
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/π)Σn=0∞[e−D(n+1/2)̂2π̂2t/Ĥ2/(n+1/2)] sin [(n+1/2)πz/H]
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.
To determine the parameter kbst, 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=−kbst[dNTP][DNA]
d[PPi]/dt=kbst[dNTP][DNA]−kc([PPi])
d[dNTP]/dt=kbst[dNTP][DNA]+kc([dNTP]0−[dNTP])
Summing the last two reactions, one finds
d([PPi]+[dNTP])/dt=−kbst[dNTP][DNA]−kc([dNTP]0−[PPi]−[dNTP])
These equations may be solved numerically to determine the values of kbst and kc, yielding kbst=2.38 (μM s)−1 and kc=0.2 s−1 at 30° C.
The PPi→Pi Enzyme
The reaction PPi+H20→2 Pi occurs at a finite rate. The rate constant in the presence of the enzyme and Mg++ is about 103 s−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−1 years−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.
Diffusion Coefficients
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:
Reference (BBA 1291 (1996) 115) gives values for ATP over a temperature range. At 37° C., the value is 4.5×10−6 cm2/s, which is used herein. The values for Pi and PPi are calculated at temperature T by scaling the value above by the factor DATP(T)/DATP(37° C.). The diffusion coefficients of the ions at 25° C. are (J. Chem. Phys. 119 (2003) 11342.)
The diffusion coefficient of H+ is assumed equal to that of Na+, because by charge neutrality the H+ must either diffusion with a negative ion as a pair, or the H+ must exchange with another positive ion, of which the N+ is the most rapidly diffusing.
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
Deff/D0=2(1−φ)/(2+φ), (6)
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.
Another approach is to use other methods to reduce the diffusion coefficient of the H+, Pi, PPi, and dNTP. For example, the viscosity of the solution may be increased.
The Sequencing Reaction
Now consider the effect of the Mg++ diffusing in. It is not known exactly how the enzyme reaction rate depends on [Mg++]. However, a divalent metal ion is required for eyzymatic activity. It is assumed that the rate is linearly proportional up to a maximum rate at 1 mM of Mg++. Thus, it is assumed the reaction rate constant is equal to
kbstmin(1,[Mg++]/1 mM) (7)
One may use a diffusion coefficient for Mg++ of 1.25×10−5 cm2/2. For Model 2, at t=0, the Mg++ is only at z=0, at concentration [Mg++]0=1 mM. For t>0 it diffuses in and affects the reaction rate by Eq. (8). For Model 1, the concentration of Mg++ is 1 mM at all times in the well. Alternatively, we may assume that the enzyme is presoaked with Mg++, so that the concentration of free Mg++ in solution is zero.
In the range of 7-9 pH, the reactions release H+ roughly as follows
This is because all the H+ with pKa values lower than the pH of the solution will come off, to a good approximation. All of these species lower the pH if the solution is pH 7 to 9. So the net effect of the reaction
dNTP+DNAn→PPi+DNAn+1 (8)
is to not change the pH.
The net effect of the reaction
PPi+H2O→2Pi (9)
is to change the pH by release of one net H+.
The H+ 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+ ions will diffuse out of the well during the equilibration time t<0. (All these H+ 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+ 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+ will come off - - - with a diffusion coefficient of something like Na+.
Well Sizes
Simulation runs were done for small wells:
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.
dNTP Front Thickness in Device Flow
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
-
- V=velocity of flow of the fluid above the well
- L=length of device
- D=diffusivity of dNTP=3.91×10−6 cm2/s
- d=size of well
- t=time to sweep flow across device
- t0=cutoff time for “instantaneous change”
- Wd=width of front due to diffusion at end of device=(2 D t)1/2
- Wm=width of front due to mixing length of flush volume region
- U=flush volume=15 μl (currently 52 μl)
- u=flow rate=50 μl/s
Flow cell geometries are as follows for two models, designated 313 and 314:
Take L=5 mm, d=4 μm and assume t=1/6.69 s. Q=50 μl/s. Then Wd=10.8 μm. Also, V=L/t=3.3 cm/s. Deff=D0 2 (1−φ)/(2+φ)=1.20×10−6 cm2/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 Npoly=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 H2/(2D), so that to =10%×H2/(2D)=0.015 s. Simulation reveals the sphere with the DNA on it slows down the diffusion, so that perhaps t0=10%*0.2 s=0.02 s. Define W=max(Wd, Wm). 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<<t0. (W+d)/V≈U/u. Since U/u=0.3 s, and t0=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 Wm 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, Wd V<<t0, so diffusive spreading of the front can be ignored.
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 t0. The concentration C(x,y,z) at t=0 may be assumed to be c0=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 [6]
v=[(P0−PL)B2/(2 μL)][1−(y/B)2] (10)
The flow rate is Q=(2/3) (p0−pL) B3 W/(μL). Calculate the time it takes for the average flow at the end of the device to be 0.95c0. First, note that without diffusion, the concentration is given by
c(x,y,z)=0, t<t*
c0, 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∂2c/∂y2 (12)
c(x,y,0)=0
c(0,y,t)=c0
∂c(x,±B,t)/∂y=0
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)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 t0=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*+t0, the concentration at the wall will be approximately c0 through a combination of diffusion and convection. Thus
v=(3/4)(1/BW)Q[1−(y/B)̂2]
Using this, t*=L/v(B−3.95 μm)=1.84 s. At a value t0 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.
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 c0 during the time δ2/(2 D)=0.01 s. So, for δ<1.0 μm, convection is more important than diffusion in the x direction.
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/(1/6.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.
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.
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 Npoly. 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.
Parametric Sensitivity of Results
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 Npoly curves that are plotted.
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++, were produced for all of the above.
Double the enzyme reaction rate constant makes the peaks very slightly higher.
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.
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.
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.
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 Npoly 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.
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.
Packing Beads—Φ (Phi)
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 Npoly curves resolve with increased depth of the well and with packing beads (Φ=0.6).
The formula to see the effect of Φ is the effective diffusion coefficient
Deff=D0*2(1−Φ)/(2+Φ)
For Φ=0.6, D0/Deff=3.25 e.g. 3.25× slowdown due to beads
For Φ=0.9, D0/Deff=14.5 e.g. 14.5× slowdown due to this
The other formula needed is the characteristic time in the well
H2=2Defft
So, for a 6 um well at Φ=0.6
t=62/(2/3.25*D0)=58.5/D0
A 2.25 well at Φ=0.9 yields
t=2.252/(2/14.5*D0)=36.7/D0
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).
Signal to Noise
The signa-to-noise ratio is proportional to the difference in mV between a Npoly and Npoly+1 signal and also to the square root of the length of time over which the mV curves are different. Define σ 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<tmax as Δ, i.e.
Δ=∫0tmax|v1−v2|/tmax.
Take the absolute value because the difference between the curves can have both signs. Thus the condition is
Δ>σ[2/(ftmax)]1/2
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 tmax=2 s. The value of Δ is 0.6861. Thus
0.6861>0.22 σ
So, using the value of σ=0.25 mV, this device will work (1−P=10−35). Consider as another example the D=4 μm, H=6 μm wells, φ=0, dNTP=6.5 μM, 10×DNA case. Take f=20 Hz, and tmax=4 s. The value of Δ is 0.1862. Thus
0.1862>0.1581σ
So, using the value of σ=0.25 mV, this device will probably work (1−P=2.4×10−6).
Note there is an optimal value of tmax: too small a value, and one does not have enough signal; too large a value, and one is seeing mostly noise. The value of tmax is only roughly here.
EQUIVALENTSWhile 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.
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.
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.
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.”
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.
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.
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.
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.
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.
Claims
1.-23. (canceled)
24. An apparatus comprising:
- a substrate including an array of sensors;
- an array of wells disposed over the substrate, a well of the array of wells corresponding to a sensor of the array of sensors; and
- a cover disposed over the array of wells and defining a flow volume between the array of wells and the cover, the cover defining a fluid inlet and a fluid outlet in fluid communication with the flow volume, the flow volume at a diffuser including a non-flat wall defining a curved boundary of the flow volume, the non-flat wall being the top or the bottom wall, the fluid flow at the diffuser being more restricted in the center of the flow volume than at the edges of the flow volume when viewed in cross-section perpendicular to the fluid flow.
25. The apparatus of claim 24, wherein the non-flat wall defines a bottom boundary of the flow volume.
26. The apparatus of claim 24, wherein the non-flat wall defines a top boundary of the flow volume.
27. The apparatus of claim 24, wherein the sensor is a field effect transistor.
28. The apparatus of claim 27, wherein the field effect transistor is an ion sensitive field effect transistor.
29. The apparatus of claim 24, wherein the fluid inlet corresponds with a first corner of the array of wells and the fluid outlet corresponds with a second corner of the array of wells.
30. The apparatus of claim 24, wherein the cover further comprises flow disruptors extending into the flow volume.
31. The apparatus of claim 24, wherein a width of the flow volume is smaller at a first position in proximity to the fluid inlet than at a second position further from the fluid inlet than the first position.
32. The apparatus of claim 31, wherein a width of the flow volume is smaller at a third position in proximity to the fluid outlet than at the second position.
33. A system comprising:
- a fluidic system including a fluid source;
- a computational system; and
- an apparatus comprising: a substrate including an array of sensors; an array of wells disposed over the substrate, a well of the array of wells corresponding to a sensor of the array of sensors; and a cover disposed over the array of wells and defining a flow volume between the array of wells and the cover, the cover defining a fluid inlet and a fluid outlet in fluid communication with the flow volume, the flow volume at a diffuser including a non-flat wall defining a curved boundary of the flow volume, the non-flat wall being the top or the bottom wall, the fluid flow at the diffuser being more restricted in the center of the flow volume than at the edges of the flow volume when viewed in cross-section perpendicular to the fluid flow.
34. The system of claim 33, wherein the non-flat wall defines a top boundary of the flow volume.
35. The system of claim 33, wherein the sensor is a field effect transistor.
36. The system of claim 35, wherein the field effect transistor is an ion sensitive field effect transistor.
37. The system of claim 33, wherein the fluid inlet corresponds with a first corner of the array of wells and the fluid outlet corresponds with a second corner of the array of wells.
38. The system of claim 33, wherein the cover further comprises flow disruptors extending into the flow volume.
39. The system of claim 33, wherein a width of the flow volume is smaller at a first position in proximity to the fluid inlet than at a second position further from the fluid inlet than the first position.
40. The system of claim 39, wherein a width of the flow volume is smaller at a third position in proximity to the fluid outlet than at the second position.
41. A method of sequencing a nucleic acid, the method comprising:
- flowing a reagent solution into a fluid inlet of an apparatus, the apparatus comprising: a substrate including an array of sensors; an array of wells disposed over the substrate, a well of the array of wells corresponding to a sensor of the array of sensors; and a cover disposed over the array of wells and defining a flow volume between the array of wells and the cover, the cover defining a fluid inlet and a fluid outlet in fluid communication with the flow volume, the flow volume at a diffuser including a non-flat wall defining a curved boundary of the flow volume, the non-flat wall being the top or the bottom wall, the fluid flow at the diffuser being more restricted in the center of the flow volume than at the edges of the flow volume when viewed in cross-section perpendicular to the fluid flow; and
- measuring a signal from the sensor of the array of sensors.
42. The method of claim 41, wherein a nucleic acid is disposed within the well and the reagent solution includes a nucleotide, wherein measuring the signal includes measuring with the sensor a byproduct released in response to incorporation of the nucleotide.
43. The method of claim 41, wherein a width of the flow volume is smaller at a first position in proximity to the fluid inlet than at a second position further from the fluid inlet than the first position.
Type: Application
Filed: Sep 26, 2016
Publication Date: May 25, 2017
Inventors: Jonathan M. ROTHBERG (Guilford, CT), Wolfgang HINZ (Killingworth, CT), Kim L. JOHNSON (Carlsbad, CA), James BUSTILLO (Castro Valley, CA), John LEAMON (Stonington, CT), Jonathan SCHULTZ (Guilford, CT)
Application Number: 15/276,709