# Patent

## Title
Multiplexed digital assays with combinatorial use of signals

## Abstract

System, including methods, apparatus, and compositions, for performing a multiplexed digital assay on a greater number of targets through combinatorial use of signals.

## Background
CROSS-REFERENCE TO PRIORITY APPLICATION
This application is based upon and claims the benefit under 35 U.S.C. §119(e) of U.S. Provisional Patent Application Ser. No. 61/454,373, filed Mar. 18, 2011, which is incorporated herein by reference in its entirety for all purposes.


CROSS-REFERENCES
This application incorporates by reference in their entireties for all purposes the following materials: U.S. Pat. No. 7,041,481, issued May 9, 2006; U.S. Patent Application Publication No. 2010/0173394 A1, published Jul. 8, 2010; PCT Patent Application Publication No. WO 2011/120006 A1, published Sep. 29, 2011; PCT Patent Application Publication No. WO 2011/120024 A1, published Sep. 29, 2011; U.S. patent application Ser. No. 13/287,120, filed Nov., 1, 2011; U.S. Provisional Patent Application Ser. No. 61/507,082, filed Jul. 12, 2011; U.S. Provisional Patent Application Ser. No. 61/510,013, filed Jul. 20, 2011; and Joseph R. Lakowicz, PRINCIPLES OF FLUORESCENCE SPECTROSCOPY (2nd Ed. 1999).
INTRODUCTION
Digital assays generally rely on the ability to detect the presence or activity of individual copies of an analyte in a sample. In an exemplary digital assay, a sample is separated into a set of partitions, generally of equal volume, with each containing, on average, less than about one copy of the analyte. If the copies of the analyte are distributed randomly among the partitions, some partitions should contain no copies, others only one copy, and, if the number of partitions is large enough, still others should contain two copies, three copies, and even higher numbers of copies. The probability of finding exactly 0, 1, 2, 3, or more copies in a partition, based on a given average concentration of analyte in the partitions, is described by a Poisson distribution. Conversely, the average concentration of analyte in the partitions may be estimated from the probability of finding a given number of copies in a partition.
Estimates of the probability of finding no copies and of finding one or more copies may be measured in the digital assay. Each partition can be tested to determine whether the partition is a positive partition that contains at least one copy of the analyte, or is a negative partition that contains no copies of the analyte. The probability of finding no copies in a partition can be approximated by the fraction of partitions tested that are negative (the “negative fraction”), and the probability of finding at least one copy by the fraction of partitions tested that are positive (the “positive fraction”). The positive fraction or the negative fraction then may be utilized in a Poisson equation to determine the concentration of the analyte in the partitions.
Digital assays frequently rely on amplification of a nucleic acid target in partitions to enable detection of a single copy of an analyte. Amplification may be conducted via the polymerase chain reaction (PCR), to achieve a digital PCR assay. The target amplified may be the analyte itself or a surrogate for the analyte generated before or after formation of the partitions. Amplification of the target can be detected optically with a fluorescent probe included in the reaction. In particular, the probe can include a dye that provides a fluorescence signal indicating whether or not the target has been amplified.
A digital PCR assay can be multiplexed to permit detection of two or more different targets within each partition. Amplification of the targets can be distinguished by utilizing target-specific probes labeled with different dyes, which produce fluorescence detected in different detection channels, namely, at different wavelengths or wavelength regions (“colors”) of emission (and/or excitation). If a detector for a digital PCR assay can distinguishably measure the fluorescence emitted by R different dyes, then the assay is effectively capable of measuring R different targets. However, instruments with more detection channels, to detect more colors, are more expensive than those with fewer detection channels. Also, increasing the number of distinguishable dyes is expensive and becomes impractical beyond a certain number. On the other hand, many applications, especially where sample is limited, could benefit greatly from higher degrees of multiplexing.
A new approach is needed to increase the multiplex levels of digital assays.
SUMMARY
The present disclosure provides a system, including methods, apparatus, and compositions, for performing a multiplexed digital assay on a greater number of targets through combinatorial use of signals.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a flowchart of an exemplary method of performing a digital assay, in accordance with aspects of the present disclosure.
FIG. 2 is a schematic view of an exemplary system for performing the digital assay of FIG. 1, in accordance with aspects of the present disclosure.
FIG. 3 is a schematic view of a pair of targets and corresponding probes capable of reporting the presence or absence of target amplification via emitted light that may be detected to create a dedicated signal for each target in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 4 is a pair of exemplary graphs of respective dedicated signals that may be created by detecting light emitted from the probes of FIG. 3 in a digital amplification assay performed in droplets, with each signal created from light detected over the same time interval from a fluid stream containing the droplets, in accordance with aspects of the present disclosure.
FIG. 5 is a schematic representation of how copies of the pair of targets of FIG. 3 are distributed among the droplets from which light is detected in FIG. 4, based on the intensity of the respective dedicated signals of FIG. 4, in accordance with aspects of the present disclosure.
FIG. 6 is a schematic view of three targets and corresponding exemplary probes capable of reporting the presence or absence of target amplification via emitted light that may be detected to create a pair of composite signals in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 7 is a pair of exemplary graphs of a pair of composite signals that may be created by detecting fluorescence emission from the three probes of FIG. 6 in a digital amplification assay performed in droplets, with emitted light detected in two different wavelength regimes over the same time interval from a fluid stream containing the droplets, in accordance with aspects of the present disclosure.
FIG. 8 is a schematic representation of how copies of the three targets of FIG. 6 are distributed among the droplets from which light is detected in FIG. 7, based on the intensity of the respective composite signals of FIG. 7, in accordance with aspects of the present disclosure.
FIG. 9A is a schematic view of the third target of FIG. 6 and another exemplary probe configuration capable of reporting the presence or absence of third target amplification via emitted light, which may be used in conjunction with the first and second target probes of FIG. 6 to create only a pair of composite signals representing amplification of the three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 9B is a schematic view of the third target of FIG. 6 and yet another exemplary probe configuration specific for the third target, which may be used in conjunction with the first and second target probes of FIG. 6 to create only a pair of composite signals representing the three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 10 is a schematic view of the third target of FIG. 6 and still another exemplary probe configuration capable of reporting the presence or absence of third target amplification via emitted light, which may be used in conjunction with the first and second target probes of FIG. 6 to create only a pair of composite signals representing amplification of the three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 11 is a schematic view of three targets and corresponding exemplary primers that enable use of only two probes to report amplification of the three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 12 is a schematic view of the third target of FIG. 6 and another exemplary probe and primer configuration that enables use of only two probes to report amplification of three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 13 is a schematic view of the three targets of FIG. 6 with yet another exemplary configuration of only two probes that enables use of the two probes to report amplification of three targets in a digital amplification assay, in accordance with aspects of the present disclosure.
FIG. 14 is a schematic view of a population of fragments containing a pair of unlinked targets, T1 and T2, in accordance with aspects of the present disclosure.
FIG. 15 is a schematic view of a population of fragments taken as in FIG. 14, but with the pair of targets always linked to each other on the same individual fragments, in accordance with aspects of the present disclosure.
FIG. 16 is a schematic view of a population of fragments taken as in FIG. 14, but with the pair of targets only partially linked to each other within the population, in accordance with aspects of the present disclosure.
FIG. 17 is a schematic representation of a set of exemplary multi-labeled probes for use in digital amplification assays, in accordance with aspects of the present disclosure.
FIG. 18 is a schematic illustration of a template molecule being copied by DNA polymerase during target amplification in the presence of a multi-labeled probe molecule and depicting probe degradation by the polymerase to separate a quencher from fluorophores of the probe molecule, in accordance with aspects of the present disclosure.
FIG. 19 is an exemplary two-dimensional histogram of droplet intensities, showing clusters that may be obtained in a multiplexed digital amplification assay for three targets performed with a combination of single-labeled and dual-labeled probes each labeled with FAM, VIC, or both FAM and VIC, in accordance with aspects of the present disclosure.
FIG. 20 is another exemplary two-dimensional histogram of droplet intensities, showing clusters that may be obtained in the assay of FIG. 19, with partial resolution of the cluster for target-1+2-positive droplets from the cluster for target-3-positive droplets, in accordance with aspects of the present disclosure.
FIG. 21 is an exemplary two-dimensional intensity histogram of droplet intensities, showing clusters that may be obtained in a digital amplification assay performed with only a multi-labeled FAM, VIC probe, in accordance with aspects of the present disclosure.
FIG. 22 is another exemplary two-dimensional intensity histogram of droplet intensities, showing clusters that may be obtained in a digital amplification assay performed as in FIG. 21, but with the assay supplemented with pair of single-labeled FAM or VIC probes in addition to the multi-labeled FAM, VIC probe, in accordance with aspects of the present disclosure.



DETAILED DESCRIPTION
The present disclosure provides a system, including methods, apparatus, and compositions, for performing a multiplexed digital assay on a greater number of targets through combinatorial use of signals. The method may be described as a color-based approach to multiplexing.
A method of performing a multiplexed digital amplification assay, such as a PCR assay, is provided. In the method, more than R targets may be amplified in partitions. R signals may be created. The signals may be representative of light detected in R different wavelength regimes from the partitions, where R≧2. An average level of each target in the partitions may be calculated based on the R signals, with the level calculated accounting for a coincidence, if any, of different targets in the same individual partitions.
Another method of performing a multiplexed digital amplification assay is provided. In the method, more than R targets may be amplified in droplets. R signals may be created, with the signals representative of light detected in R different wavelength regimes from the droplets, where R≧2. An average level of each of the more than R targets may be calculated by finding solutions to a set of simultaneous equations.
Yet another method of performing a multiplexed digital amplification assay is provided. In the method, R targets may be amplified in droplets. R signals may be created, where R≧2, with the signals representative of light detected in R different wavelength regimes from the droplets. Each of the signals may report amplification of a different combination of at least two of the targets. An average level of each target in the droplets may be calculated based on the R signals, without determining which of the at least two targets for each signal amplified in individual amplification-positive droplets for such signal.
A composition is provided. The composition may comprise a droplet containing a probe. The probe may include an oligonucleotide, a first fluorophore, a second fluorophore, and an energy transfer moiety. The energy transfer moiety may be a quencher and/or an energy transfer partner for one or both of the first and second fluorophores.
Further aspects of the present disclosure are presented in the following sections: (I) system overview, and (II) examples.
I. SYSTEM OVERVIEW
FIG. 1 shows a flowchart of an exemplary method 40 of performing a digital assay. The steps presented for method 40 may be performed in any suitable order and in any suitable combination. Furthermore, the steps may be combined with and/or modified by any other suitable steps, aspects, and/features of the present disclosure.
A sample may be prepared for the assay, indicated at 42. Preparation of the sample may include any suitable manipulation of the sample, such as collection, dilution, concentration, purification, lyophilization, freezing, extraction, combination with one or more assay reagents, performance of at least one preliminary reaction to prepare the sample for one or more reactions in the assay, or any combination thereof, among others. Preparation of the sample may include rendering the sample competent for subsequent performance of one or more reactions, such as one or more enzyme catalyzed reactions and/or binding reactions.
In some embodiments, preparation of the sample may include combining the sample with reagents for amplification and for reporting whether or not amplification occurred. Such reagents may include any combination of primers for the targets (e.g., a forward primer and a reverse primer for each target), reporters (e.g., probes) for detecting amplification of the targets, dNTPs and/or NTPs, at least one enzyme (e.g., a polymerase, a ligase, a reverse transcriptase, or a combination thereof, each of which may or may not be heat-stable), or the like. Accordingly, preparation of the sample may render the sample (or partitions thereof) capable of amplification of each of one or more targets, if present, in the sample (or a partition thereof).
The sample may be separated into partitions, indicated at 44. Separation of the sample may involve distributing any suitable portion or all of the sample to the partitions. Each partition may be and/or include a fluid volume that is isolated from the fluid volumes of other partitions. The partitions may be isolated from one another by a carrier fluid, such as a continuous phase of an emulsion, by a solid phase, such as at least one wall of a container, or a combination thereof, among others. In some embodiments, the partitions may be droplets disposed in a continuous phase, such that the droplets and the continuous phase collectively form an emulsion.
The partitions may be formed by any suitable procedure, in any suitable manner, and with any suitable properties. For example, the partitions may be formed with a fluid dispenser, such as a pipette, with a droplet generator, by agitation of the sample (e.g., shaking, stirring, sonication, etc.), or the like. Accordingly, the partitions may be formed serially, in parallel, or in batch. The partitions may have any suitable volume or volumes. The partitions may be of substantially uniform volume or may have different volumes. Exemplary partitions having substantially the same volume are monodisperse droplets. Exemplary volumes for the partitions include an average volume of less than about 100, or 1 μL, less than about 100, 10, or 1 nL, or less than about 100, 10, or 1 pL, among others.
The partitions, when formed, may be competent for performance of one or more reactions in the partitions. Alternatively, one or more reagents may be added to the partitions after they are formed to render them competent for reaction. The reagents may be added by any suitable mechanism, such as a fluid dispenser, fusion of droplets, or the like. Any of the reagents may be combined with the partitions (or a bulk phase sample) in a macrofluidic or microfluidic environment.
One or more reactions may be performed in the partitions, indicated at 46. Each reaction performed may occur selectively (and/or substantially) in only a subset of the partitions, such as less than about one-half, one-fourth, or one-tenth of the partitions, among others. The reaction may involve a target, which may, for example, be a template and/or a reactant (e.g., a substrate), and/or a binding partner, in the reaction. The reaction may occur selectively (or selectively may not occur) in partitions containing at least one copy of the target.
The reaction may or may not be an enzyme-catalyzed reaction. In some examples, the reaction may be an amplification reaction, such as a polymerase chain reaction and/or ligase chain reaction. Accordingly, a plurality of amplification reactions for a plurality of targets may be performed simultaneously in the partitions.
Performing a reaction may include subjecting the partitions to one or more conditions that promote occurrence of the reaction. The conditions may include heating the partitions and/or incubating the partitions at a temperature above room temperature. In some examples, the conditions may include thermally cycling the partitions to promote a polymerase chain reaction and/or ligase chain reaction.
R signals may be created that are representative of light detected from the partitions, indicated at 48. The R signals may be 2, 3, 4, or more signals. In some examples, light corresponding to each signal may be detected with a distinct sensor, and/or light corresponding to at least two signals may be detected at different times with the same sensor. The R signals may correspond to light detected in respective wavelength regimes that are different from one another. Each wavelength regime may be characterized by a wavelength(s) or and/or wavelength range(s) at which the partitions are illuminated (e.g., with excitation light) and/or a wavelength(s) or and/or wavelength range(s) at which light from the partitions is detected (e.g., emitted light). The light detected may be light emitted from one or more fluorophores.
Each of the R signals may be created in a distinct detection channel. Accordingly, R signals may be created in R detection channels.
Each signal may be a composite signal that represents two, three, four, or more reactions/assays and thus two, three, four, or more targets of the reactions/assays. The composite signal may include two or more integral signal portions that each represent a different reaction/assay and target. Analysis of one of the composite signals by itself, without the benefit of data from the other composite signals, may (or may not) permit estimation of a collective concentration, but not individual concentrations, for two or more targets represented by the composite signal. Instead, as described further below, analysis of the composite signals together permits calculation of the concentration of each target. (The terms “estimation” and “calculation” are used interchangeably.)
The R composite signals (and/or R detection channels) may represent more than R reactions and/or targets, with the number of reactions/assays and targets depending on configuration. For example, the R signals may be or include two composite signals (arbitrarily termed α and β) collectively representing three reactions/assays and/or three targets (arbitrarily termed 1, 2, and 3), with each composite signal representing a different combination of two reactions/assays/targets (e.g., targets 1 and 2 for α and targets 1 and 3 for β). Alternatively, the R signals may be three composite signals (arbitrarily termed α, β, and γ) collectively representing up to seven reactions/assays/targets (1 to 7), if each composite signal represents a different combination of up to four reactions/assays/targets each (e.g., targets 1, 2, 3, and 4 for α; targets 2, 4, 5, and 6 for β; and targets 3, 4, 6, and 7 for γ). The R signals may be four composite signals representing up to fifteen reactions/assays/targets, if each composite signal represents a different combination of up to eight reactions/assays/targets each.
More generally, 2R−1 targets can be assayed with R composite signals (or wavelength regimes). To assay 2R−1 targets in R wavelength regimes, each target may be represented by a different wavelength regime or combination of wavelength regimes than every other target. A set of 2R−1 targets can be represented and assayed when all of the wavelength regimes have been utilized individually and in all possible combinations.
Each composite signal may be created based on detected light emitted from one or more probes in the partitions. The one or more probes may report whether at least one of two or more particular reactions represented by the signal has occurred in a partition and thus whether at least one copy of at least one of two or more particular targets corresponding to the two or more particular reactions is present in the partition. The intensity of a composite signal corresponding to the probes may be analyzed to determine whether or not at least one of the particular reactions has occurred and at least one copy of one of the particular targets is present. The intensity may vary among the partitions according to whether at least one of the particular reactions occurred or did not occur (e.g., above a threshold extent) and at least one of the particular targets is present in or absent from each individual partition.
The probes represented by a composite signal may include different fluorophores. In other words, light emitted from different fluorophores may be detected to create at two different integral portions of the composite signal for a particular wavelength regime. Alternatively, or in addition, the same fluorophore may be included in one probe or two or more probes for at least two targets represented by the composite signal. In some cases, the same fluorophore may be included in a probe for each target represented by the composite signal.
Each probe may include a nucleic acid (e.g., an oligonucleotide) and at least one fluorophore. Different probes with different oligonucleotide sequences and/or different fluorophores (or fluorophore combinations) may be used to create at least two different integral portions of the signal. Alternatively, or in addition, the same probe may be used as a reporter for at least two different targets represented by the composite signal (e.g., see Examples 3-5). In some cases, the same probe may be used as a reporter for each target represented by the composite signal.
The composite signal detected from each partition, and the partition itself, may be classified as being positive or negative for the reactions/assays/targets represented by the signal or corresponding wavelength regime. Classification may be based on the strength of the signal. If the signal/partition is classified as positive, at least one of the reactions/assays represented by the signal is deemed to have occurred and at least one copy of at least one of the targets represented by the signal is deemed to be present in the partition. In contrast, if the signal/partition is classified as negative, none of the reactions/assays represented by the signal is deemed to have occurred and no copy of any of the targets represented by the signal is deemed to be present in the partition.
The composite signals collectively permit estimation of target concentrations by representation of a different combination of targets in each detection channel. Accordingly, each target, when present without any of the other targets in a partition, may produce a unique target signature among the wavelength regimes. For example, some of the targets, if present alone in a partition, may selectively change the signal strength for only one wavelength regime. Others of the targets, if present alone in a partition, may selectively change the signal strength for a unique combination of two of the wavelength regimes, still other targets may selectively change the signal strength for a unique combination of three of the wavelength regimes, and so on, optionally up to the number of wavelength regimes/detection channels.
A fraction of the partitions may have a coincidence of different targets, where each of these partitions contains a copy of each of two or more targets in the same individual partitions. Moreover, each of these partitions may contain a copy of each of two or more distinct targets, which, for a particular partition, collectively may produce a signature that is the same as that of a target not present in the partition. However, the fraction of partitions containing two or more distinct targets may (or may not) be kept relatively low, by working in a dilute regime, such as with less than about one-half, one-fifth, or one-tenth, among others, of the partitions containing more than one target molecule when the partitions are formed. In any event, a suitable estimation of concentration, as described below, may take into account the occurrence of two or more target molecules, representing the same target or different targets, in the same individual partitions. Alternatively, if working in a sufficiently dilute regime, the occurrence of two or more target molecules per partition may be sufficiently rare to ignore for a desired accuracy of concentration.
A number of partitions that are positive may be determined for each signal, indicated at 50. For example, a number of partitions that are positive only for each particular composite signal or corresponding wavelength regime/detection channel may be determined individually (e.g., counted) for each signal or channel (i.e., a number for each channel). Also, a number of partitions that are positive only for each particular combination (or at least one combination or each of two or more combinations) of composite signals or corresponding wavelength regimes may be determined individually (e.g., counted) for each combination of signals or channels (i.e., a number for each combination, and particularly each combination corresponding to a particular target).
A distinct fraction of the partitions positive for each signal alone and for each signal combination may be determined. The fraction for each signal or signal combination may be determined by dividing the number of partitions for the signal/combination, determined at 50, by the total number of partitions from which signals are detected. The total number of partitions may be counted or estimated.
A level of each target may be calculated, indicated at 52. The level may be an average level, such as an average concentration of molecules of the target per partition. Generally, if R signals are detected from the partitions in R wavelength regimes, the average level of each of more than R targets (e.g., up to 2R−1 targets) may be calculated. The level of each target may be calculated based on the respective numbers of partitions positive for each signal alone and signal combination. The calculation may be based on copies of each target having a Poisson distribution among the partitions. The concentrations may, for example, be calculated by finding solutions to a series of simultaneous equations (interchangeably termed a set of simultaneous equations), each having the same variables. The simultaneous equations may be linear equations. Alternatively, or in addition, each equation may contain at least 2R−1 variables. The solutions may be found by numerical analysis, also termed numerical approximation. Further aspects of calculating average target levels are described elsewhere in the present disclosure, such as in Examples 6 and 7.
FIG. 2 shows an exemplary system 60 for performing the digital assay of FIG. 1. System 60 may include a partitioning assembly, such as a droplet generator 62 (“DG”), a thermal incubation assembly, such as a thermocycler 64 (“TC”), a detection assembly (a detector) 66 (“DET”), and a data processing assembly (a processor) 68 (“PROC”), or any combination thereof, among others. The data processing assembly may be, or may be included in, a controller that communicates with and controls operation of any suitable combination of the assemblies. The arrows between the assemblies indicate optional movement or transfer, such as movement or transfer of fluid (e.g., a continuous phase of an emulsion) and/or partitions (e.g., droplets) or signals/data. Any suitable combination of the assemblies may be operatively connected to one another, and/or one or more of the assemblies may be unconnected to the other assemblies, such that, for example, material/data is transferred manually.
System 60 may operate as follows. Droplet generator 62 may form droplets disposed in a carrier fluid, such as a continuous phase. The droplets may be cycled thermally with thermocycler 64 to promote amplification of targets in the droplets. Composite signals may be detected from the droplets with detector 66. The signals may be processed by processor 68 to calculate levels of the targets.
Further aspects of sample preparation, droplet generation, assays, reagents, reactions, thermal cycling, detection, and data processing, among others, that may be suitable for the methods and systems disclosed herein, are described below and in the documents listed above under Cross-References, which are incorporated herein by reference, particularly U.S. Patent Application Publication No. 2010/0173394 A1, published Jul. 8, 2010; PCT Patent Application Publication No. WO 2011/120006 A1, published Sep. 29, 2011; PCT Patent Application Publication No. WO 2011/120024 A1, published Sep. 29, 2011; and U.S. patent application Ser. No. 13/287,120, filed Nov., 1, 2011.
II. EXAMPLES
This section presents selected aspects and embodiments of the present disclosure related to methods and compositions for performing multiplexed digital assays. These aspects and embodiments are included for illustration and are not intended to limit or define the entire scope of the present disclosure.
Example 1
Digital Amplification Assays with Dedicated Signals and Composite Signals
This example compares and contrasts exemplary digital amplification assays utilizing (i) a pair of dedicated signals for two targets, see FIGS. 3-5, and (ii) a pair of composite signals for three targets, see FIGS. 6-8. The principles explained here may be extended to R signals for 2R−1 targets.
FIG. 3 shows a pair of nucleic acid targets 80, 82 (“Target 1” and “Target 2”) and corresponding probes 84, 86 (“Probe 1” and “Probe 2”) that may be used to create a dedicated signal for amplification of each target in a digital amplification assay. Each probe may include an oligonucleotide 88, 90, a fluorophore 92, 94, and a quencher 96. The fluorophore(s) and quencher are associated with and/or attached to the oligonucleotide, such as attached covalently. The probe also or alternatively may include a binding moiety (a minor groove binder) for the minor groove of a DNA duplex, which may be conjugated to the oligonucleotide and may function to permit a shorter oligonucleotide to be used in the probe. The probe may be a TaqMan probe, a molecular beacon probe, a scorpion probe, a locked nucleic acid probe, or the like.
Each oligonucleotide may provide target specificity by hybridization predominantly or at least substantially exclusively to an amplicon produced by amplification of only one of the two targets. Hybridization of the oligonucleotide to its corresponding target/amplicon is illustrated schematically at 98.
Fluorophores 92, 94 may be optically distinguishable from each other, as illustrated schematically by a distinct hatch pattern for each fluorophore. For example, the fluorophores may have distinct absorption spectra and/or maxima, and/or distinct emission spectra and/or emission maxima. Proper selection of the wavelength regime used for detection allows the fluorophores to be distinguished. In other words, the wavelength or wavelength band of excitation light used for each wavelength regime and/or the wavelength or wavelength band of emitted light received and sensed by the sensor for the wavelength regime provides selective detection of light from only one of the fluorophores in the detection channel. Exemplary fluorophores that may be suitable include FAM, VIC, HEX, ROX, TAMRA, JOE, etc.
Quencher 96 is configured to quench the signal produced by fluorophore 92 or 94 in a proximity-dependent fashion. Light detected from the fluorophore may increase when the associated oligonucleotide 88 or 90 binds to the amplified target, to increase the separation distance between the fluorophore and the quencher, or when the probe is cleaved during target amplification, among others. In some cases, the quencher may be replaced by, or supplemented with, a fluorophore that is capable of energy transfer with fluorophore 92 or 94.
FIG. 4 shows a pair of exemplary graphs 102, 104 of data collected in an exemplary digital amplification assay for Target 1 and Target 2 performed in droplets. Each graph plots a dedicated signal 106 (“Signal 1”) or signal 108 (“Signal 2”) that represents light detected from respective probes 84, 86 (and/or one or more modified (e.g., cleavage) products thereof) (see FIG. 3). Each dedicated signal is created from light detected over the same time period from a fluid stream containing the droplets and flowing through an examination region of a channel. Signal 1 reports whether or not Target 1 is present in each droplet, and Signal 2 reports whether or not Target 2 is present in each target. In particular, if the strength of Signal 1 (or Signal 2) increases above a threshold 110, then Target 1 (or Target 2) is deemed to be present (and amplified) in a corresponding droplet. In the present illustration, each droplet, whether positive or negative for each target, produces an increase in signal strength above the baseline signal that forms an identifiable peak 112. Accordingly, each signal may vary in strength with the presence or absence of a droplet and with the presence or absence of a corresponding target.
Each target is present here at an average level (or frequency) of about 0.2 in the droplets. In other words, each target is amplified and detected on average about once every five droplets. Accordingly, the expected frequency of droplets containing both targets is the product of the two droplet frequencies, or about 0.04 (1 out of every 25 droplets). Consistent with this frequency, a droplet that is positive for both targets is present only once on the twenty droplets analyzed here, and is indicated by a dashed box at 114 extending around the signal peaks for the droplet in graphs 102, 104.
FIG. 5 schematically represents the distribution of Targets 1 and 2 in a set of droplets 116 corresponding to and in the same order as the droplet signal peaks of FIG. 4. Droplets positive for Signal 1, such as the droplet indicated at 118, are hatched according to fluorophore 92, and droplets positive for Signal 2, such as the droplet indicated at 120, are hatched according to fluorophore 94 (see FIG. 3). A double-positive droplet 122 containing both Target 1 and Target 2 is double-hatched and indicated by dashed box 114.
FIG. 6 shows three targets 80, 82, and 140 and corresponding exemplary probes 84, 86, and 142, respectively, that may be used to create a pair of composite signals for the three targets in a digital amplification assay. Two of the targets and probes, namely, targets 80 and 82 (Target 1 and Target 2) and probes 84 and 86 (Probe 1 and Probe 2) are the same targets and probes shown and utilized in FIGS. 3-5. Target 140 (Target 3) and its corresponding probe 142 (Probe 3) may be introduced into the multiplexed assay of FIGS. 3-5 to increase the level of multiplexing and the amount of target information that can be extracted from the assay without increasing the number of detection channels.
Amplification of Target 3 is reported by Probe 3. The probe includes an oligonucleotide 144 that hybridizes specifically to Target 3 (and/or an amplicon thereof), relative to Targets 1 and 2. The probe may be double-labeled with the same fluorophores (92, 94) present individually in Probe 1 and Probe 2 for reporting respective Target 1 and Target 2 amplification. The probes for the three targets may be selected to permit detection of target amplification in only two detection channels, rather than the three detection channels that would be necessary with the use of a dedicated detection channel for each target. Examples 2-5 describe other exemplary probe configurations that may be suitable to increase the level of multiplexing.
FIG. 7 shows a pair of exemplary graphs 152, 154 of a pair of composite signals 156, 158 that may be detected in a pair of wavelength regimes/detection channels. The composite signals, arbitrarily designated α and β, are representative of light detected from the three probes of FIG. 6 in a digital amplification assay performed in droplets. Each composite signal is created from light detected over the same time period from a fluid stream containing the droplets. To simplify the presentation, Target 1 and Target 2 are present at the same frequency and in the same droplets as in FIGS. 4 and 5.
Each composite signal, α or β (156 or 158), represents a pair of targets. Signal a (graph 152) has a strength for each droplet that indicates whether the droplet is positive or negative for at least one member of a first pair of targets, namely, Target 1 and Target 3. Signal β (graph 154) has a strength for each droplet that indicates whether the droplet is positive or negative for at least one member of a different second pair of targets, namely, Target 2 and Target 3. Accordingly, each composite signal analyzed by itself may provide no information about how frequently a given member of the pair of targets is present in droplets.
The composite signals analyzed in combination provide additional information about target frequency that cannot be deduced from the composite signals in isolation from one another. Each target, when present without other targets in a droplet, produces a signal signature that is distinct from the signatures of each other target individually. The signature for Target 1 in a droplet is indicated at 160: positive for Signal α and negative for Signal β. The signature for Target 2 in a droplet is indicated at 162: negative for Signal α and positive for Signal β. Furthermore, the signature for Target 3 in a droplet is outlined by dashed boxes at 164: positive for both Signal α and Signal β. Finally, the signature for none of Targets 1, 2, and 3 in a droplet is indicated at 166: negative for both Signal α and Signal β.
FIG. 8 schematically represents the distribution of Targets 1, 2, and 3 in a set of droplets 168 corresponding to and in the same order as the droplet signal peaks of FIG. 7. Single-positive droplets that are positive for Signal α only, such as the droplet indicated at 118, are hatched according to fluorophore 92, and droplets positive for Signal β, such as the droplet indicated at 120, are hatched according to fluorophore 94 (see FIG. 6). Each double-positive droplet 170 is double-hatched and indicated by dashed box 164.
The single-positive signatures indicated at 160 and 162 unambiguously identify corresponding droplets 118, 120 as containing at least one copy of Target 1 or Target 2, respectively, and no copy of Target 3. Accordingly, the number of each type of single-positive droplet may be used, in a ratio with the total number of droplets, to calculate an average level of Target 1 and of Target 2. However, this estimate may not be accurate enough if droplets contain multiple target molecules, because the estimate ignores any droplets containing Target 1 and/or Target 2, but having the signature of Target 3. These droplets can produce the same Target 3 signature while containing Targets 1+2, Targets 1+3, or Targets 2+3. If the concentration of each target is low enough, the frequency of droplets containing at least two different targets may be negligible and/or ignored. Here, the concentrations of Targets 1, 2, and 3 are high enough to produce, on average, only about one droplet with both of Targets 1 and 2 per twenty droplets (see FIGS. 4 and 5). Accordingly, one of double-positive droplets 170 is expected to contain both of Targets 1 and 2, and the other two double-positive droplets are expected to contain Target 3.
It is not necessary to know the target composition of each double-positive droplet. Instead, it is sufficient to know the frequency of droplets with each target signature. Poisson statistics then may be utilized to calculate the average level of each target. Calculation may be performed, in some cases by finding solutions to a series of simultaneous equations, such as numerically, to obtain a best-fit, or by a closed-form approach, among others.
Example 2
Exemplary Target-Specific Probes for Composite Signals
This example describes additional, exemplary target-specific probes that may be utilized in any suitable multiplexed digital assay; see FIGS. 9A, 9B, and 10. The principles explained here may be extended to any number of signals and/or targets.
FIGS. 9A, 9B, and 10 shows third target 140 of FIG. 6 and other exemplary probe configurations of probes 180, 182 (Probes 3A and 3B) each specific for Target 3 (and/or an amplicon thereof). Probes 3A and 3B may be used together, in place of Probe 3 of FIG. 6, and in conjunction with Probe 1 and Probe 2 of FIG. 6, to create only a pair of composite signals for Targets 1 to 3 in a digital amplification assay.
One of the probes (e.g., Probe 3A) may include fluorophore 92 and the other probe (e.g., Probe 3B) may include fluorophore 94. Accordingly, light emitted by Probe 3A can be detected in the same detection channel as light from Probe 1 of FIG. 6, and light emitted by Probe 3B in the same channel as Probe 2.
FIGS. 9A and 9B show Probes 3A and 3B binding specifically to distinct sites of Target 3. The probes may bind to non-overlapping (or only partially overlapping) sites on the same strand of the target (FIG. 9A). Alternatively, the probes may bind to opposite strands of Target 3 (FIG. 9B).
FIG. 10 shows Probes 3A and 3B that are capable of binding to the same site of Target 3 (and/or an amplicon thereof). Each of the probes may contain the same oligonucleotide 144 and thus differ only by the particular fluorophore (92 or 94) attached to oligonucleotide 144. The probes of this example may be blended in any suitable ratio or set of different ratios for a multiplexed assay.
Example 3
Exemplary Tailed Primers and Shared Probes for Composite Signals
This example describes exemplary shared probes that enable creation of composite signals, and exemplary tailed primers forming binding sites for the shared probes; see FIG. 11. The principles explained here may be extended to any suitable number of composite signals and/or targets.
FIG. 11 shows three targets 80, 82, and 140 (i.e., Targets 1 to 3 of FIG. 6) and corresponding primers that enable assay of the three targets with only two probes, namely, probe 190 (Probe A) and probe 192 (Probe B), in a digital amplification assay. Probe A and Probe B may include respective fluorophores and a quencher (e.g., see FIG. 3).
Target 1 may be amplified with a pair of forward and reverse primers 194, 196. Forward primer 194 may be a tailed primer with a 3′ binding portion 198 that is complementary to a region of Target 1 and a 5′ tail portion 200 that is not. The tail portion may introduce a binding site for Probe A into the resulting amplicon, indicated by a dashed line at 202, such that Probe A (like Probe 1 of FIG. 6) can report amplification of Target 1.
Target 2 may be amplified with a pair of forward and reverse primers 204, 206. Reverse primer 206 may be structured analogously to forward primer 194 for Target 1, with a 3′ binding portion complementary to a region of Target 2 and a 5′ tail portion 207 that is not. The tail portion may introduce a binding site for Probe B into the resulting amplicon, such that Probe B (like Probe 2 of FIG. 6) can report amplification of Target 2.
Target 3 may be amplified with a pair of forward and reverse primers 208, 210. Both of the primers may be structured analogously to forward primer 194, with a 3′ binding portion complementary to a region of Target 3 and a 5′ tail portion that is not. Tail portions 200, 207 of primers 208, 210 may introduce respective binding sites for Probes A and B into the amplified product, such that a combination of both Probe A and Probe B (like Probe 3 of FIG. 6) reports amplification of Target 3.
Example 4
Exemplary Ligation Strategy to Enable Shared Probes
This example describes an exemplary ligation strategy to enable the use of shared probes; see FIG. 12. The principles explained here may be extended to any number of composite signals and/or targets.
FIG. 12 shows Target 3 (at 140) of FIG. 6 and another exemplary probe and primer configuration that enables assay of Targets 1 to 3 of FIG. 6 in a digital amplification assay with only Probes 1 and 2 (at 84 and 86) of FIG. 6. The ligation, extension, and digestion steps presented below may be performed in any suitable order and before or after a sample providing Target 3 is separated into partitions.
A template 220 (and/or a complementary strand and/or amplicon thereof) may be designed to bind each of Probes 1 and 2.
Template 220 also may be designed to bind to adjacent regions of Target 3 via opposing end regions 222, 224 of the template. (The template may be described as a molecular inversion “probe,” but is generally not attached to a fluorophore.) The 5′ and 3′ termini of the template may be joinable directly to one another by ligation when bound to Target 3 or may form a gap 226 of one or more nucleotides between the aligned 5′ and 3′ termini of the template. The gap may be closed by extending the 3′ terminus of the template, while bound to Target 3, before ligation of the template to form a closed loop. After ligation, and optional extension, copies of the template that fail to ligate (and thus have not found a copy of Target 3 for binding), may be degraded by an exonuclease. Ligated copies of the template may be resistant to this degradation, such that the number of ligated template molecules corresponds to the number of Target 3 molecules.
Ligated template 220 (and/or a complementary strand thereof) may provide one or more sites for binding of at least one primer 228 (or 230). The primer may amplify the ligated template by rolling circle amplification. Alternatively, or in addition, a pair of primers 228, 230 (forward and reverse) may be included to produce a cascade of amplification. In some embodiments, ligated template 220 may be linearized by cleavage at a predefined site 232 before amplification with primers 228, 230. In any event, the presence of Target 3 in partitions is reported by a combination of both Probes 1 and 2 in this embodiment.
Example 5
Exemplary Shared Probes with Multi-target Specificity
This example describes exemplary shared probes each capable of binding to sequence regions of two different targets; see FIG. 13. The principles explained here may be extended to any number of composite signals and/or targets.
FIG. 13 shows targets 80, 82, 140 (i.e., Targets 1 to 3) bound by probes with multi-target specificity. In particular, a probe 240 (Probe 1/3) includes an oligonucleotide 242 capable of binding to a sequence region present in Target 1 and another sequence region present in Target 3. Also, another probe 244 (Probe 2/3) includes an oligonucleotide 246 capable of binding to a sequence region present in Target 2 and another sequence region present in Target 3.
Example 6
Increasing Multiplex Levels in Digital Amplification without Additional Detection Channels
This example describes an exemplary approach for increasing the multiplex level of a multiplexed digital amplification assay.
A. Introduction
The ability to measure multiple targets simultaneously (multiplexing) within every partition of a digital amplification system is often limited by the detection approach. Commonly one measures fluorescence to classify partitions as positive (if the measured fluorescence is high) or negative (if the measured fluorescence is low). Some chemistries, such as TaqMan, allow measurements of several targets simultaneously by utilizing target-specific probes labeled with different dyes. If the detector can measure the fluorescence emitted by R different dyes, then the digital amplification system is effectively capable of measuring R different targets. Typically, instruments that are capable of detecting more colors are more expensive than those with fewer colors. Increasing the number of detectable dyes is expensive and is impractical beyond a certain number. On the other hand, many applications especially where sample is limited could benefit greatly from higher degrees of multiplexing.
This example presents an approach that, given an instrument capable of detecting multiple colors, can dramatically increase the number of simultaneously measured targets without requiring any changes to the detection optics of the instrument. The standard approach of designing assays is that a given target is assessed based on fluorescence produced from a single probe with a single dye. Thus, if the instrument is capable of detecting two colors, such as the light emitted from the dyes FAM and VIC, one measures the concentration of one target by counting the number of partitions with positive FAM signals and another target by counting the number of partitions with positive VIC signals.
One can design assays that produce fluorescence on multiple channels simultaneously. If processed on a digital amplification platform with a large number of partitions these assays can be multiplexed together with single channel assays and can be measured by counting the number of partitions with fluorescence on both channels.
B. Example with Two-color FAM-VIC Detection
Assuming we are looking at two unlinked loci (target 1 and target 2), and given some number of FAM-only positive droplets as well as some number of VIC-only positive droplets, we can estimate how many FAM-VIC double-positive droplets we expect. If we are operating at low concentrations this number should be small and can be worked out in a straightforward fashion.
If we set up a third assay (target 3) such that it has two additional probes—one labeled with FAM and one labeled with VIC, we can estimate the concentration of this third target locus by how many excess FAM-VIC double-positive droplets we have compared to the expectation. This would reduce the overall precision, but not much, and basically not at all if we are operating in a dilute regime (i.e., the total number of droplets is much larger than the number of positive droplets). Below is an example of an algorithm that can be used to determine the concentration of the excess FAM-VIC double positive droplets (or other partitions).
The use of multiple probes labeled with the same dye will increase the fluorescence of the negative droplets, which can present a challenge in extreme cases if fluorescence of the negative droplets starts approaching that of the positive droplets. This challenge can be addressed effectively by using sufficiently robust assays. One can also use common probes (e.g., see Examples 3-5) and avoid the elevation of negative fluorescence altogether. For the above example, we can consider using a common FAM probe for target 1 and target 3 and a common VIC probe for target 2 and target 3 by utilizing tailed primers or locked nucleic acid probes.
C. Additional Considerations
One gains an ever-larger advantage from this approach when one uses four or more colors. There are six combinations of two colors if one has four to choose from. Together with single colors, this would give a total of ten reporters. If we go further and use triplets of colors we would end up with 13 reporters.
The advantage of using this multi-color scheme becomes more pronounced with higher numbers of partitions. For that reason, this approach is of particular utility when combined with more recent implementations of digital amplification such as digital PCR in droplets where thousands or millions of partitions can be produced in an easy and cost effective manner.
Several assay schemes can be employed to assess a target with multiple colors simultaneously. One could design a multi-labeled probe—e.g., a single probe can be labeled with both FAM and VIC on the same molecule (e.g., see FIG. 6). As another example, the same oligonucleotide may be labeled separately with FAM and VIC, to produce a FAM-labeled version and a VIC-labeled version of the same probe, and then the two versions mixed (e.g., see FIG. 9B). In other cases, such as for a TaqMan assay, two probes can be designed to bind to different regions of the same amplicon strand (e.g., see FIG. 9A). Alternatively, the probes can bind to opposite strands of the amplicon (e.g., see FIG. 10), which may position the dyes away from the quenchers and facilitate the fluorescence increase from the bound probes.
This approach is general and can be used with a range of chemistries including ligation chain reaction, molecular beacons, scorpion probes, molecular inversion probes, or the like.
D. Mathematical Approach for Estimating Excess FAM-VIC Droplets
The following is an example of an algorithm that can be used to estimate concentrations of a joint FAM-VIC species (e.g., target 3 of FIGS. 6-8).

    
    
        1. Get 2×2 table of FAM versus VIC counts.
        2. Compute concentration of distinct FAM and joint FAM-VIC as if there are 1 species.
        3. Compute concentration of distinct VIC and joint FAM-VIC as if there are 1 species.
        4. Try out different concentrations of joint FAM-VIC (from which the concentration of distinct FAM and distinct VIC can be found), and find the best fit of the probability table (Table 1) with the observed counts.
    
    












TABLE 1







FAM−
FAM+




















VIC+
(1 − f) v (1 − c)
1 − sum of others



VIC−
(1 − f) (1 − v) (1 − c)
f (1 − v) (1 − c)










E. MATLAB Implementation of the Algorithm
Below is an example of a MATLAB implementation of the algorithm. Note that the algorithm can be expanded in a straightforward fashion to high order multiplexes.

    % Consider three types of DNA fragments: Fam-Vic together,
    % Fam fragment, Vic fragment. We observe some probabilities (counts in
    % FAM-VIC cross plot), and the goal is to infer the concentrations.
    % First let us do forward. Given concentrations, compute counts. Then to do
    % inverse, we simply try out different values of concentrations and select
    % one which gives actual counts.
    N=20000;
    A=10000;
    B=20000;
    AB=10000; % Joined together
    cA=A/N;
    cB=B/N;
    cAB=AB/N;
    fprintf(1, ‘% f % f % fn’, cAB, cA, cB);
    pA=1−exp(−cA);
    pB=1−exp(−cB);
    pAB=1−exp(−cAB);
    % A is X and B is Y in cross plot
    p(2,1)=(1−pA)*(1−pB)*(1−pAB); % Bottom left
    p(2,2)=pA*(1−pB)*(1−pAB); % Bottom right
    p(1,1)=(1−pA)*pB*(1−pAB); % Top Left
    p(1,2)=1−p(2,1)−p(2,2)−p(1,1); % Top Right
    disp(round(p*N));
    % Also compute marginals directly
    cAorAB=(A+AB)/N; %=c_A+c_AB;
    cBorAB=(B+AB)/N; %=c_B+c_AB;
    pAorAB=1−exp(−cAorAB); % Can be computed from p too
    pBorAB=1−exp(−cBorAB);
    % Inverse
    H=p*N; % We are given some hits
    % H=[08000; 20000];
    % Compute prob
    estN=sum(H(:));
    i_p=H/estN;
    i_pAorAB=i_p(1,2)+i_p(2,2);
    i_pBorAB=i_p(1,1)+i_p(1,2);
    i_cAorAB=−log(1−i_pAorAB);
    i_cBorAB=−log(1−i_pBorAB);
    maxVal=min(i_cAorAB, i_cBorAB);
    delta=maxVal/1000;
    errArr=[ ];
    gcABArr=[ ];
    forgcAB=0:delta:maxVal
    gcA=i_cAorAB−gcAB;
    gcB=i_cBorAB−gcAB;
    gpA=1−exp(−gcA);
    gpB=1−exp(−gcB);
    gpAB=1−exp(−gcAB);
    gp(2,1)=(1−gpA)*(1−gpB)*(1−gpAB); % Bottom left
    gp(2,2)=gpA*(1−gpB)*(1−gpAB); % Bottom right
    gp(1,1)=(1−gpA)*gpB*(1−gpAB); % Top Left
    gp(1,2)=1−gp(2,1)−gp(2,2)−gp(1,1); % Top Right
    gH=gp*estN;
    err=sqrt(sum((H(:)−gH(:)).^2));
    errArr=[errArr; err];
    gcABArr=[gcABArr; gcAB];
    end
    figure, plot(gcABArr, errArr);
    minidx=find(errArr==min(errArr(:)));
    minidx=minidx(1);
    estAB=gcABArr(minidx);
    estA=i_cAorAB−estAB;
    estB=i_cBorAB−estAB;
    fprintf(1, ‘% f % f % f\n’, estAB, estA, estB);
    gpA=1−exp(−estA);
    gpB=1−exp(−estB);
    gpAB=1−exp(−estAB);
    gp(2,1)=(1−gpA)*(1−gpB)*(1−gpAB); % Bottom left
    gp(2,2)=gpA*(1−gpB)*(1−gpAB); % Bottom right
    gp(1,1)=(1−gpA)*gpB*(1−gpAB); % Top Left
    gp(1,2)=1−gp(2,1)−gp(2,2)−gp(1,1); % Top Right
    gH=gp*estN;
    disp(round(gH));
    % Confirm the results using simulation
    numMolA=round(estA*estN);
    numMolB=round(estB*estN);
    numMolAB=round(estAB*estN);
    A=unique(randsample(estN, numMolA, 1));
    B=unique(randsample(estN, numMolB, 1));
    AB=unique(randsample(estN, numMolAB, 1));
    U=1:estN;
    notA=setdiff(U, A);
    notB=setdiff(U, B);
    notAB=setdiff(U, AB);
    AorBorAB=union(A, union(B, AB));
    none=setdiff(U, AorBorAB);
    simcount(2,1)=length(none);
    simcount(2,2)=length(intersect(A, intersect(notB, notAB)));
    simcount(1,1)=length(intersect(B, intersect(notA, notAB)));
    simcount(1,2)=length(AorBorAB)−simcount(2,2)−simcount(1,1);
    disp(simcount);


Example 7
Algorithm for Computation of DNA Fragmentation or for Digital Amplification Multiplexing
This example describes an exemplary algorithm to compute a level of DNA fragmentation and/or levels of target in a multiplexed amplification assay.
A. Introduction
1. Totally Fragmented Targets
Consider two DNA targets T1 and T2 corresponding to two dyes FAM and VIC, respectively. Let T1 and T2 be always on separate DNA fragments, as illustrated schematically in FIG. 14. Let the number of DNA fragments with T1 and T2 targets be M1 and M2, respectively.
Let the counts of FAM and VIC positive partitions be N1 and N2, respectively. Note that N1 and N2 will be smaller than M1 and M2, respectively, as there can be multiple DNA fragments in a partition. Let the total number of partitions be N. We will refer to digital amplification partitions simply as partitions. In this case, we can expect to see the counts of partitions listed in Table 2.











TABLE 2







VIC Negative
VIC positive
Total



















FAM Positive
N1 · (N − N2)/N
N1 · N2/N
N1


FAM Negative
(N − N1) · (N − N2)/N
(N − N1) · N2/N
N − N1


Total
N − N2
N2
N









If we denote the probability of seeing a partition to be FAM positive as p1=N1/N, and of seeing a partition to be VIC positive as p2=N2/N, then the probability table is given by Table 3.











TABLE 3







VIC Negative
VIC positive
Probability



















FAM Positive
p1 · (1 − p2)
p1 · p2
p1


FAM Negative
(1 − p1) · (1 − p2)
(1 − p1) · p2
1 − p1


Probability
1 − p2
p2
1









In this case, we can say that 100% fragmentation exists.

We can compute the number of T1 and T2 molecules, M1 and M2, respectively as follows, wherein (where log=loge):

M1=−N log(1−p1)

M2=−N log(1−p2)

(Given N digital partitions in which P are positive, the number of molecules is M=−N log(1−P/N).)
2. No Fragmentation
Now consider the other extreme. Both targets T1 and T2 are always found together on the same DNA fragments (e.g., see FIG. 15). They are linked, perhaps because their loci are quite close to each other on the same part of a chromosome, and fragmentation of the chromosome during DNA isolation did not separate the loci. Therefore, N1=N2. The expected counts and probabilities are listed in Tables 4 and 5, respectively.











TABLE 4







VIC Negative
VIC positive
Total





















FAM Positive
0
N1
N1



FAM Negative
N − N1
0
N − N1



Total
N − N1
N1
N





















TABLE 5







VIC Negative
VIC positive
Probability



















FAM Positive
0
p1
p1


FAM Negative
(1 − p1)
0
1 − p1


Probability
1 − p1
p1
1









In this case, we can say that 0% fragmentation exists.

We can compute the number of T1 and T2 molecules as follows, where p1=N1/N:

M1=−N log(1−p1)

M2=−N log(1−p1)

3. Partial Fragmentation
In the intermediate situation, where the targets are together on some fragments, but also happen to be on separate fragments, then we have partial fragmentation (e.g., see FIG. 16).
Suppose we have M3 molecules of linked T1 and T2 fragments, M1 molecules of separate T1 fragments, and M2 molecules of separate T2 fragments.
We can make a table of counts of partitions (Table 6).











TABLE 6







VIC Negative
VIC positive
Total





















FAM Positive
N01
N11
N1



FAM Negative
N00
N10
N − N1



Total
N − N2
N2
N










4. Problem Statement
How can we find the number of molecules M1, M2 and M3, and thereby getting extent of fragmentation? For example if M1=M2=M3, then we would say that there is 50% fragmentation, as 50% of linked molecules got fragmented into separate fragments and 50% remained intact.
B. Multiplexing of PCR to Many Targets Using Few Colors
We have an algorithm that provides a solution to the above problem, and there is another interesting application for the algorithm. With this algorithm and by using a FAM probe for target T1, a VIC probe for target T2, and both FAM and VIC probes placed close to each other for target T3, we can achieve multiplexing of 3 targets by using 2 colors. Basically, we get a third color for “free” using the algorithm.
Now consider the case if there are 3 dyes. Thus, we will have a 2×2×2 table of 8 observed counts.
There are seven different kinds of targets: T1, T2, T3, T12, T23, T13, and T123. Here, for example, T12 means a target that is amplified to produce an amplicon bound by Dye 1 and Dye 2 probes, and T123 means a target that is amplified to produce an amplicon bound by one or more probes containing all three dyes, and likewise for the others.
If we have 4 dyes, then we have 24=16 counts, and we can now multiplex quantitation of 24−1=15 target genes. In general, with R colors, we have 2R observed counts and we can have 2R−1 targets.
C. Solution for 2 Dyes and 3 Targets
Suppose we are given the counts listed in Table 7.











TABLE 7







VIC Negative
VIC positive
Total





















FAM Positive
N01
N11
N1



FAM Negative
N00
N10
N − N1



Total
N − N2
N2
N










Consider the following three cases:

    
    
        1. We turn off VIC, as if VIC cannot be seen at all.
        2. We turn off FAM, as if FAM cannot be seen at all.
        3. We consider both FAM and VIC as if they are only one color.
    
    


We will have the following three observations:

    
    
        1. Turning off VIC: We will be able to see T1 or T12, together and indistinguishably, as if there were one target. It gives the total number of molecules of T1 and T12, as if there were one target species.
        2. Turning off FAM: We will be able to count T2 or T12, indistinguishably, as if there were one target. It gives the total number of molecules of T2 and T12, as if there were one target species.
        3. Considering both FAM and VIC indistinguishably: We will be able to count T1, T2, or T12, indistinguishably. It gives the total number of molecules of T1, T2, and T12, as if there were one target species.
    
    


This allows us to step 3 equations in 3 unknowns. We can show these three cases in the form of a table.













TABLE 8





Visible
Invisible
Distinct
Indistinct
Target
Positive



Dyes
Dyes
Dyes
Dyes
Detected
Count
Molecules







FAM
VIC
—
—
{T1, T12}
N1
−log(1 − N1/N)


VIC
FAM
—
—
{T2, T12}
N2
−log(1 − N2/N)


FAM,
—
—
{FAM,
{T1, T2, T12}
N01 + N10 +
−log(1 − (N01 + N10 +


VIC


VIC}

N11
N11)/N)









The three linear equations in three unknowns M1, M2 and M12 are:




  
    
      [
      
        
          
            101
          
        
        
          
            011
          
        
        
          
            111
          
        
      
      ]
    
    ⁡
    
      [
      
        
          
            
              M
              ⁢
              
                
              
              ⁢
              1
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              2
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              12
            
          
        
      
      ]
    
  
  =
  
    
      [
      
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
        
          
            
              M
              
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
      
      ]
    
    =
    
      [
      
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        N
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        N
                        ⁢
                        
                          
                        
                        ⁢
                        2
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          01
                        
                        +
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          10
                        
                        +
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          11
                        
                      
                      N
                    
                  
                  )
                
              
            
          
        
      
      ]
    
  




We solve the above equations to get values of M1, M2, and M12.
We can then compute the extent of fragmentation in % as follows:

M=(M1+M2)/2

F=M/(M+M12)*100

D. Algorithm Steps
Now we write down the steps of the algorithm clearly using the input of Table 9.











TABLE 9







VIC Negative
VIC positive
Total





















FAM Positive
N11
N12
N1



FAM Negative
N21
N22
N − N1



Total
N − N2
N2
N










Algorithm:

Step 1. Compute the three entities:

M1or12=−N·log(1−N1/N)

M2or12=−N·log(1−N2/N)

M1or2or12=−N·log(1−(N01+N10+N11)/N)

Step 2. Solve the following linear equations:




  
    
      [
      
        
          
            101
          
        
        
          
            011
          
        
        
          
            111
          
        
      
      ]
    
    ⁡
    
      [
      
        
          
            
              M
              ⁢
              
                
              
              ⁢
              1
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              2
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              12
            
          
        
      
      ]
    
  
  =
  
    
      [
      
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
        
          
            
              M
              
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
      
      ]
    
    =
    
      [
      
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        N
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        N
                        ⁢
                        
                          
                        
                        ⁢
                        2
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          01
                        
                        +
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          10
                        
                        +
                        
                          N
                          ⁢
                          
                            
                          
                          ⁢
                          11
                        
                      
                      N
                    
                  
                  )
                
              
            
          
        
      
      ]
    
  




Step 3. Compute extent of fragmentation.

M=(M1+M2)/2

F=M/(M+M12)*100.

Step 4. Compute confidence intervals based on the concentration and expected number of positive counts. To compute confidence interval of F, we note that it is ratio of two random variables. Then we apply techniques to estimate confidence interval of ratio of two random variables for F.
The output is M1, M2, and F and their confidence intervals.
E. Outline of an Alternative Solution
Now we present an alternative solution to the problem of partial fragmentation phrased in terms of optimization of an objective criterion. First let us make the following table. Depending upon what type of molecules we have in a partition we will have corresponding FAM and VIC colors of the partition.
Table 10 maps molecules into partition colors.








TABLE 10







Presence of Molecules in a Partition
Partition Color











T1
T2
T12
FAM (T1)
VIC (T2)





0
0
0
Neg
Neg


0
1
0
Neg
Pos


1
0
0
Pos
Neg


0
0
1
Pos
Pos


0
1
1


1
0
1


1
1
0


1
1
1









Suppose we make a “guess” of M1, M2, and M12. We can compute the probabilities of a partition having a copy of the T1, T2, or T12 target using inverse equations:

p12=1−exp(−M12/N)

p1=1−exp(−M1/N)

p2=1−exp(−M2/N).

From these probabilities we can compute the predicted counts as shown in Table 11.










TABLE 11







VIC Negative
VIC positive


















FAM Positive
p1(1 − p2)(1 − p12)N
(1 − (sum of other 3




cells in this table))N


FAM
(1 − p1)(1 − p2)(1 − p12)N
(1 − p1)p2(1 − p12)N


Negative









The probabilities above can be filled using Table 10, which maps presence of molecules into partition colors.
If our “guess” is really correct, then the predicted counts will “match” closely with our expected counts. Thus an optimization algorithm under the above objective criterion that needs to be minimized could also be used to solve the problem.

    
    
        M1, M2, M12=best guess=least deviation of predicted counts and actual counts
    
    


We can then compute the extent of fragmentation in % as follows:

M=(M1+M2)/2

F=M/(M+M12)*100.

F. Solution for 3 Dyes and 7 Targets
To solve the problem for greater number of dyes, we need more equations as there are more unknowns. For 3 dyes and 7 targets, we need 7 equations to find the concentration of these 7 targets.
We have Table 12 with 23=8 rows of counts of partitions depending upon which dyes are positive.
Denote the three dyes by D1, D2 and D3.










TABLE 12









Presence of color in a partition













D1
D2
D3
Count of Partitions







Neg
Neg
Neg
N000



Neg
Pos
Neg
N010



Pos
Neg
Neg
N100



Neg
Neg
Pos
N001



Neg
Pos
Pos
N011



Pos
Neg
Pos
N101



Pos
Pos
Neg
N110



Pos
Pos
Pos
N111








Total
N









Consider Table 13.













TABLE 13





Visible
Invisible
Distinct
Indistinct
Target



Dyes
Dyes
Dyes
Dyes
Detected
Positive Count







D1, D2,
—
—
{D1, D2,
{T1, T2,
C1 =


D3


D3}
T3, T12,
N − N000






T13, T23,






T123}


D1
D2, D3
—
—
{T1, T12,
C2 =






T13, T123}
N100 +







N101 +







N110 + N111


D2
D1, D3
—
—
{T2, T12,
C3 =






T23, T123}
N010 +







N011 +







N110 + N111


D3
D1, D2
—
—
{T3, T13,
C4 =






T23, T123}
N001 +







N011 +







N101 + N111


D1, D2
D3
D1, D2
—
{T1, T13}
Compute







concentration







M1or13 by







solving 2 dye







problem*


D1, D3
D2
D1, D3
—
{T1, T12}
Compute







concentration







M1or12 by







solving 2 dye







problem*


D2, D3
D1
D2, D3
—
{T2, T12}
Compute







concentration







M2or12 by







solving 2 dye







problem*





*Note that for row 5, we could have also detected {T2, T23} or {T12, T123}, and for row 6, {T3, T23}, and for row 7, {T3, T13} and {T23, T123}. This choice does not matter and all lead to the same results.






We make 3 recursive calls, in rows 5, 6 and 7.
For example, in row 5, we are really solving 2 dyes, D1 and D2, case with 3 targets: {T1, T13} as one target, {T2, T23} as second target, and {T12, T123} as third linked target.
Similarly, in rows 6 and 7, we make recursive calls to solve the simpler problem of 2 dyes and 3 targets.
The system of equation looks like as follows:




  
    
      [
      
        
          
            1111111
          
        
        
          
            1001101
          
        
        
          
            0101011
          
        
        
          
            0010111
          
        
        
          
            1000100
          
        
        
          
            1001000
          
        
        
          
            0101000
          
        
      
      ]
    
    ⁡
    
      [
      
        
          
            
              M
              ⁢
              
                
              
              ⁢
              1
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              2
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              3
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              12
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              13
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              23
            
          
        
        
          
            
              M
              ⁢
              
                
              
              ⁢
              123
            
          
        
      
      ]
    
  
  =
  
    
      [
      
        
          
            
              M
              all
            
          
        
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                13
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                123
              
            
          
        
        
          
            
              M
              
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
                ⁢
                
                  
                
                ⁢
                o
                ⁢
                
                  
                
                ⁢
                23
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                123
              
            
          
        
        
          
            
              M
              
                3
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                13
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                23
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                123
              
            
          
        
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                13
              
            
          
        
        
          
            
              M
              
                1
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
        
          
            
              M
              
                2
                ⁢
                
                  
                
                ⁢
                or
                ⁢
                
                  
                
                ⁢
                12
              
            
          
        
      
      ]
    
    =
    
      [
      
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        1
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        2
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        3
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              
                -
                N
              
              ⁢
              
                
              
              ⁢
              
                log
                ⁡
                
                  (
                  
                    1
                    -
                    
                      
                        C
                        ⁢
                        
                          
                        
                        ⁢
                        4
                      
                      N
                    
                  
                  )
                
              
            
          
        
        
          
            
              Recursive
              ⁢
              
                
              
              ⁢
              call
            
          
        
        
          
            
              Recursive
              ⁢
              
                
              
              ⁢
              call
            
          
        
        
          
            
              Recursive
              ⁢
              
                
              
              ⁢
              call
            
          
        
      
      ]
    
  




The above system of linear equations can be solved to find concentrations of 7 targets.
G. R Dyes and 2R−1 Targets
Now we give the steps to generalize the algorithm to an arbitrary number of dyes.
This is a recursive algorithm. Since we have a solution for the case when R=2, we can solve for any number of R through recursion.
Input:
For R dyes, we are given a table of counts of partitions which has 2R rows, for all possible combinations of presence or absence of colors in a partition.
Algorithm:
Step 1. Set up a system of 2R−1 linear equations. To obtain these equations consider different ways of making some colors invisible. Also consider the case when all colors are indistinguishable, which gives the first row. It can be shown that we can obtain 2R−1 equations in 2R−1 target concentrations (unknowns). When we make certain colors invisible, then we reduce the problem to a case with fewer colors, which can be solved recursively, all the way down to case when there are 2 dyes and 3 targets.
Step 2. Solve the equations.
Step 3. Compute confidence intervals based on the concentration and expected number of positive counts.
Output:
Concentrations of 2R−1 targets along with confidence intervals.
Example 8
Multi-Labeled Probes
This example describes exemplary multi-labeled probes and methods of using the probes in a multiplexed digital amplification assay; see FIGS. 17-22. The multi-labeled probes may be utilized for the color-based multiplexed assays described in Section I and in Examples 6 and 7. Alternatively, or in addition, the multi-labeled probes may be utilized for intensity-based multiplexing, as described in U.S. Provisional Patent Application Ser. No. 61/507,082, filed Jul. 12, 2011; and U.S. Provisional Patent Application Ser. No. 61/510,013, filed Jul. 20, 2011, which are incorporated herein by reference.
Typical 5′ nuclease assay (TaqMan) probes have a single fluorophore and a quencher, or two fluors where the second fluor acts to quench the first fluorophore through FRET (e.g., a TAMRA “quencher”). The fluorophore and quencher are typically attached to the 5′- and 3′-most bases/nucleotides of the oligonucleotide probe.
Oligonucleotide synthesis chemistry allows fluorophores to be added to internal nucleotides/bases of a probe. Attaching multiple fluorophores to one oligonucleotide probe allows creation of a wider range of probes, which can be used to enhance multiplexing capabilities. The ability to put multiple fluors on one probe allows the resulting emission fluorescence to be “tuned” to achieve more fluorescence signatures than are possible through single fluorophores. Probe spectra will generally be a composite result of the multiple fluorophores included. For fluorophores that can resonate or otherwise interact, the proximity of the fluorophores and location (internal or end label) will allow additional tuning possibilities. It is not necessary that the multiple fluorophores (or quenchers) be different. Addition of multiple molecules of the same fluor also may allow different fluorescence output on a per-probe basis. For some applications, it may be beneficial to put the quencher on the 5′-end of the oligonucleotide so that probe degradation removes the quencher, but leaves the other fluorophores attached on the remaining portion of the oligonucleotide. That way, when the quencher is cleaved off, the same signal is obtained from the probe, but each probe could have a different signature (e.g., probe 1—FAM; probe 2—FAM+Cy3; probe 3—FAM+HEX; probe 4—FAM+FAM; etc.).
The approach is also applicable to multiplexing for other target molecules. Although detection of target nucleic acids is clearly an area of great interest, it is also possible to apply this concept to other types of probes (for example, antibody probes for protein targets).
FIG. 17 show a schematic representation of an unlabeled oligonucleotide (A) and a set of exemplary multi-labeled probes (B-G) containing the oligonucleotide. Each solid circle and open circle represents a quencher (also termed a quencher moiety) and a fluorophore, respectively, attached to the oligonucleotide. F1 to F4 identify structurally distinct fluorophores. The oligonucleotide may have any suitable length, such as 5 to 500, 10 to 200, or 15 to 100 nucleotides, among others.
The multi-labeled probes of FIG. 17 are as follows: B) oligonucleotide with a 5′ quencher and fluorophores (F1 and F2) at respective internal and 3′-end positions; C) oligonucleotide with a 5′ quencher moiety and fluorophores (F1 to F4) at 3′ and multiple internal positions; D) oligonucleotide with a 5′ quencher moiety and fluorophores (3×F1, F2) at 3′ and multiple internal positions; E) oligonucleotide with a 5′ quencher moiety and an internal quencher moiety plus fluorophores at 3′ and internal positions; F) oligonucleotide with a 5′ quencher moiety and fluorophores at 3′ and internal positions, with greater separation between fluorophores; and G) oligonucleotide with fluorophores at the 5′-end and an internal position and a quencher at the 3′-end position.
FIG. 18 depicts separation of the quencher and fluorophores of a probe molecule during target amplification. A template molecule is shown being copied (complementarily) by DNA polymerase in the presence of an exemplary multi-labeled probe molecule (FIG. 17B). A) Polymerase is extending a nucleic acid strand upstream of a probe binding site. B) 5′ nuclease activity of the polymerase removes the 5′-most nucleotides (including the quencher in this example) from the probe, allowing the quencher and fluorophores to be spatially separated. C) Polymerase continues to remove 5′ nucleotides from the probe until binding of the probe becomes unstable and the remaining probe fragment dissociates from the template strand. Multiple fluorophores still may be connected to each other in the probe fragment, but are separated from the quencher.
FIG. 19 shows an exemplary two-dimensional histogram of droplet intensities, showing droplet clusters that may be obtained in a multiplexed digital amplification assay for three targets performed with a combination of single-labeled and dual-labeled probes each labeled with FAM, VIC, or both FAM and VIC. Here, the droplet cluster produced by target-3 positives is not resolved from droplets double-positive for both target 1 and 2.
FIG. 20 shows another exemplary two-dimensional histogram of droplet intensities that may be obtained in the assay of FIG. 19, with partial resolution of a cluster for target-1+2-positive droplets from a cluster for target-3-positive droplets.
Detector optics, particularly excitation sources and optical filters, may be selected to optimize the separation of clusters. For example, two clusters might be substantially overlapping (hard to separate) at a first wavelength condition but substantially non-overlapping at a second wavelength condition. Detectors with spectrophotometer gratings or exchangeable filter sets could provide greater flexibility in wavelength selection.
FIG. 21 shows an exemplary two-dimensional histogram of droplet intensities, illustrating two clusters of data points that may be obtained in a digital amplification assay performed with only a multi-labeled probe.
FIG. 22 shows another exemplary two-dimensional histogram of droplet intensities, illustrating various clusters of data points that may be obtained in a digital amplification assay performed as in FIG. 21, but with the assay supplemented with a pair of single-labeled probes (for measuring target 1 and target 2 amplification), in addition to the multi-labeled probe (for measuring target 3 amplification). Putting multiple fluors on the same oligonucleotide can be more advantageous when combined with different levels of probe and/or with different fluors that have different spectral overlap in two (or more) detection channels. The example of FIG. 22 shows the dual fluor probe having slightly higher fluorescence than the FAM-only and VIC-only probes, but it could be either higher or lower (or the same). If about the same (or not), the color-based multiplexing approach disclosed herein may be utilized advantageously to determine average target levels.
Several chemistries can yield the multi-labeled probes disclosed herein, including labeled deoxynucleotides, click chemistry, and various other linker chemistries. The probes could be produced through custom synthesis by an oligonucleotide supplier. References describing exemplary synthetic routes for internal fluorophore incorporation are listed below and are incorporated herein by reference.

    Haralambidis J., Chai M. and Tregear G. W. (1987) Preparation of base-modified nucleosides suitable for non-radioactive label attachment and their incorporation into synthetic oligodeoxyribonucleotides. Nucleic Acids Res. 15, 4857-4876.
    Randolph J. B., and Waggoner A. S. (1997) Stability, specificity and fluorescence brightness of multiply-labelled fluorescence DNA probes. Nucleic Acids Research; 25:2923-2929.
    Brumbaugh J. A., Middendorf L. R., Grone D. L., and Ruth J. L. Proc. Natl. Acad. Sci. USA 1988; 85:5610-5614.
    Singh D., Vijayanti K., Ganesh K. N. (1990) Oligonucleotides, part 5+: synthesis and fluorescence studies of DNA oligomers d(AT) 5 containing adenines covalently linked at C-8 with dansyl fluorophore. Nucleic Acids Res.; 18:3339-3345.
    Tae Seok Seo, Zengmin Li, Hameer Ruparel, and Jingyue Ju (2003) Click Chemistry to Construct Fluorescent Oligonucleotides for DNA Sequencing. J. Org. Chem.; 68: 609-612.


Example 9
Selected Embodiments
This example describes selected aspects and embodiments related to digital assays with combinatorial use of signals, presented without limitation as a series of numbered paragraphs. Each of these paragraphs can be combined with one or more other paragraphs, and/or with disclosure from elsewhere in the present disclosure, in any suitable manner. Some of the paragraphs below expressly refer to and further limit other paragraphs, providing without limitation examples of some of the suitable combinations.
A1. A method of performing a digital assay, comprising: (a) creating R signals representative of light detected from each of a plurality of partitions of a sample; and (b) estimating a concentration of more than R different targets in the partitions based on how the R signals vary relative to one another among the partitions.
B1. A method of performing a digital assay for more than R targets, comprising: (a) separating a sample into partitions; (b) creating R signals representative of light detected from the partitions; and (c) estimating a concentration of more than R targets in the partitions based on the R signals.
B2. The method of paragraph B1, wherein at least one of the R signals reports the presence or absence of a target in each partition independently of every other one of the R signals.
B3. The method of paragraph B1 or B2, wherein each of the R signals reports the presence or absence of a different target in each partition independently of every other one of the R signals.
B4. The method of any of paragraphs B1 to B3, wherein a combination of two or more of the R signals collectively reports the presence or absence an R+1 target in each partition.
B5. The method of any of paragraphs B1 to B4, further comprising a step of determining a number of the partitions that are positive for each of the R signals alone and a number that are positive for at least one combination of two or more of the R signals.
B6. The method of paragraph B5, wherein the step of estimating a concentration includes a step of estimating concentrations of each of the more than R targets that collectively correspond to the determined numbers of positives, if each target has a Poisson distribution among the partitions.
B7. The method of any of paragraphs B1 to B6, wherein the step of estimating includes a step of finding solutions to a set of simultaneous equations, and wherein the equations each have the same variables.
B8. The method of paragraph B7, wherein the simultaneous equations are linear equations.
B9. The method of paragraph B7 or B8, wherein the solutions are obtained by numerical analysis.
B10. The method of any of paragraphs B7 to B9, wherein the step of finding solutions includes a step of finding solutions to at least 2R−1 equations.
B11. The method of any of paragraphs B7 to B10, wherein each equation is based on copies of each target having a Poisson distribution among the partitions.
B12. The method of any of paragraphs B1 to B11, wherein each of the R signals is a composite signal that includes two or more integral signal portions corresponding to the presence or absence of different targets in individual partitions.
B13. The method of any of paragraphs B1 to B12, wherein the step of separating a sample forms the partitions with an average concentration per partition of less than about one copy of each of the more than R targets.
B14. The method of any of paragraphs B1 to B13, wherein the step of separating a sample forms one or more partitions containing no copies of a target for each of the more than R targets.
B15. The method of any of paragraphs B1 to B14, wherein the partitions are droplets.
B16. The method of any of paragraphs B1 to B15, wherein each of the R signals is representative of fluorescence emission that is detected.
B17. The method of any of paragraphs B1 to B16, further comprising a step of performing an amplification reaction in one or more of the partitions before the step of creating R signals.
B18. The method of paragraph B17, wherein the step of performing an amplification reaction includes a step of thermally cycling partitions.
B19. The method of any of paragraphs B1 to B18, wherein each target is a nucleic acid.
B20. The method of any of paragraphs B1 to B19, wherein the partitions include a first probe reporting a presence or absence of a first target molecule in individual partitions, a second probe reporting a presence or absence of a second target molecule in individual partitions, and at least one third probe reporting a presence or absence of a third target molecule in individual partitions.
B21. The method of claim B20, wherein the at least one third probe is a single third probe.
B22. The method of any of paragraphs B1 to B19, wherein the partitions include a first probe reporting a presence or absence of a first target molecule in individual partitions, a second probe reporting a presence or absence of a second target molecule in individual partitions, and wherein the first and second probes collectively report a presence or absence of a third target molecule in individual partitions.
B23. The method of any of paragraphs B1 to B22, wherein each of the R signals is representative at least predominantly of light detected from a different fluorophore.
B24. The method of any of paragraphs B1 to B23, wherein the partitions all have substantially same volume.
B25. The method of any of paragraphs B1 to B24, wherein at least one of the more than R targets is a linked version of at least two of the other R targets.
C1. A method of performing a digital assay, comprising: (a) separating a sample into partitions, with each partition capable of amplifying more than R targets, if present in the partition; (b) creating R signals representative of light detected from the partitions, wherein, for at least one of the more than R targets, amplification in a partition of the at least one target selectively changes only one of the R signals, and wherein, for at least one other of the more than R targets, amplification in a partition of at least one other target coordinately changes two or more of the signals; and (c) estimating a concentration of each of the more than R targets based on the R signals created.
C2. The method of paragraph C1, wherein the step of amplifying includes a step of thermally cycling the partitions before the step of creating.
C3. The method of paragraph C1 or C2, wherein the R signals are representative of light detected from fluid carrying droplets through an examination region.
C4. The method of any of paragraphs C1 to C3, wherein the step of creating includes a step of creating two or more signals that represent different wavelengths and/or wavelength ranges of detected light.
C5. The method of any of paragraphs C1 to C4, wherein the step of creating includes a step of creating two or more signals that are representative of a same wavelength range of detected light produced by illumination with a different wavelength or wavelength range for each of the two or more signals.
D1. A method of performing a multiplexed digital amplification assay, the method comprising: (a) amplifying more than R targets in partitions; (b) creating R signals representative of light detected in R different wavelength regimes from the partitions, where R≧2; and (c) calculating an average level of each target in the partitions based on the R signals, wherein the level calculated accounts for a coincidence, if any, of different targets in the same individual partitions. For example, if T denotes the number of targets, then T>R.
D2. The method of paragraph D1, wherein amplification of each target is reported by a different signal or combination of the signals than any of the other targets individually.
D3. The method of paragraph D1 or D2, wherein each of the signals reports amplification of a different combination of at least two of the targets.
D4. The method of paragraph D3, wherein the partitions are droplets, further comprising a step of determining a number of droplets exhibiting amplification of any of the at least two targets for each signal, and wherein the step of calculating is based on the number determined for each of the R signals.
D5. The method of any of paragraphs D1 to D4, wherein the step of calculating includes a step of finding solutions to a set of simultaneous equations.
D6. The method of paragraph D5, wherein there are T targets, wherein the step of finding solutions includes a step of finding solutions to T simultaneous equations, and wherein R<T≦2R−1.
D7. The method of paragraph D5 or D6, wherein the solutions are obtained by numerical analysis.
D8. The method of any of paragraphs D1 to D7, wherein there are three targets, and wherein the signals representative of light are detected from two fluorophores associated with probes that bind to amplicons of respective targets during amplification.
D9. The method of paragraph D8, wherein the two fluorophores are VIC and FAM.
D10. The method of paragraph D8 or D9, wherein the partitions contain a first probe for a first of the three targets, a second probe for a second of the three targets, and a third probe for a third of the three targets, and wherein the first probe is labeled exclusively with VIC, the second probe is labeled exclusively with FAM, and the third probe is labeled with both VIC and FAM.
D11. The method of paragraph D10, wherein the third probe includes a FAM-labeled probe that is not labeled with VIC and a VIC-labeled probe that is not labeled with FAM.
D12. The method of any of paragraphs D1 to D11, wherein the average level is a concentration.
D13. The method of any of paragraphs D1 to D12, wherein the step of calculating an average level includes a step of determining a total number of amplification-positive partitions for each type of target and a step of determining a total number of partitions.
D14. The method of any of paragraphs D1 to D13, further comprising a step of distributing copies of the more than R targets among the partitions such that some partitions contain more than one copy of a given target.
D15. The method of any of paragraphs D1 to D13, wherein there are at least four targets.
E1. A method of performing a multiplexed digital amplification assay, the method comprising: (a) amplifying more than R targets in droplets; (b) creating R signals representative of light detected in R different wavelength regimes from the droplets, where R≧2; and (c) calculating an average level of each of the more than R targets by finding solutions to a series of simultaneous equations. For example, if T denotes the number of targets, then T>R.
E2. The method of paragraph E1, wherein amplification of each target is reported by a different signal or combination of the signals than any of the other targets individually.
E3. The method of paragraph E1 or E2, wherein each of the signals reports amplification of a different combination of at least two of the targets.
E4. The method of paragraph E3, further comprising a step of determining a total number of droplets that are amplification-positive for any of the at least two targets reported on by each signal, and wherein the step of calculating is based on the total number determined for each of the R signals.
E5. The method of any of paragraphs E1 to E4, wherein there are T targets, wherein the step of calculating includes a step of finding solutions to T simultaneous equations, and wherein R<T≦2R−1.
E6. The method of any of paragraphs E1 to E5, wherein the solutions are obtained by numerical analysis.
E7. The method of any of paragraphs E1 to E6, wherein the level accounts for any coincidence of different targets in the same individual droplets.
F1. A method of performing a multiplexed digital amplification assay, the method comprising: (a) amplifying more than R targets in droplets; (b) creating R signals representative of light detected in R different wavelength regimes from the droplets, wherein R≧2 and each of the signals reports amplification of a different combination of at least two of the targets; and (c) calculating an average level of each target in the droplets based on the R signals and without determining which of the at least two targets for each signal amplified in individual amplification-positive droplets for such signal. For example, if T denotes the number of targets, then T>R.
F2. The method of paragraph F1, further comprising a step of determining a total number of droplets that are amplification-positive for any of the at least two targets reported on by each signal, and wherein the step of calculating is based on the total number determined for each of the R signals.
F3. The method of paragraph F1 or F2, wherein the step of calculating includes a step of finding solutions to a set of simultaneous equations.
G1. A composition, comprising: a droplet containing a probe, the probe including an oligonucleotide, a first fluorophore, a second fluorophore, and an energy transfer moiety, wherein the energy transfer moiety is a quencher and/or an energy transfer partner for one or both of the first and second fluorophores.
G2. The composition of paragraph G1, wherein the first fluorophore, the second fluorophore, and the energy transfer moiety are each covalently attached to the oligonucleotide.
G3. The composition of paragraph G1 or G2, further comprising a plurality of droplets containing the probe and disposed in a carrier fluid.
G4. The composition of any of paragraphs G1 to G3, wherein the droplet contains a template molecule and amplification reagents capable of amplifying at least a region of the template molecule, and wherein the probe is capable of binding to amplicons generated by amplification of the region of the template molecule.
G5. The composition of any of paragraphs G1 to G4, wherein the probe is a first probe, and wherein the droplet further comprises a second probe including the first fluorophore or the second fluorophore, but not both the first fluorophore and the second fluorophore.
G6. The composition of any of paragraphs G1 to G5, wherein the energy transfer moiety is a quencher attached to a nucleotide at the 5′-end of the oligonucleotide.
G7. The composition of any of paragraphs G1 to G6, wherein each of the first fluorophore, the second fluorophore, and the energy transfer moiety is attached to a different nucleotide of the oligonucleotide.
G8. The composition of paragraph G7, wherein a pair of the first fluorophore, the second fluorophore, and the energy transfer moiety are attached to the same nucleotide of the oligonucleotide.
G9. The composition of any of paragraphs G1 to G8, wherein at least one of the first fluorophore, the second fluorophore, and the energy transfer moiety is attached to a nucleotide of the oligonucleotide via another of the first fluorophore, the second fluorophore, and the energy transfer moiety.
G10. The composition of any of paragraphs G1 to G9, wherein a fluorophore or a quencher is attached to the 5′-end of the oligonucleotide, wherein a fluorophore or a quencher is attached to the 3′-end of the oligonucleotide, and wherein a fluorophore or a quencher is attached to a nucleotide intermediate the 5′-end and the 3′-end.
G11. The composition of any of paragraphs G1 to G10, wherein one or more of the first fluorophore, the second fluorophore, and the energy transfer moiety are attached to one or more internal nucleotides disposed intermediate the 5′-end and the 3′-end of the oligonucleotide.
The disclosure set forth above may encompass multiple distinct inventions with independent utility. Although each of these inventions has been disclosed in its preferred form(s), the specific embodiments thereof as disclosed and illustrated herein are not to be considered in a limiting sense, because numerous variations are possible. The subject matter of the inventions includes all novel and nonobvious combinations and subcombinations of the various elements, features, functions, and/or properties disclosed herein. The following claims particularly point out certain combinations and subcombinations regarded as novel and nonobvious. Inventions embodied in other combinations and subcombinations of features, functions, elements, and/or properties may be claimed in applications claiming priority from this or a related application. Such claims, whether directed to a different invention or to the same invention, and whether broader, narrower, equal, or different in scope to the original claims, also are regarded as included within the subject matter of the inventions of the present disclosure. Further, ordinal indicators, such as first, second, or third, for identified elements are used to distinguish between the elements, and do not indicate a particular position or order of such elements, unless otherwise specifically stated.

## Claims

1. A method of performing a multiplexed digital amplification assay, the method comprising:
amplifying more than R targets in partitions;
creating R signals representative of light detected in R different wavelength regimes from the partitions, where R≧2; and
calculating an average level of each target in the partitions based on the R signals, wherein the level calculated accounts for a coincidence of all possible combinations of the more than R targets in the same individual partitions;
wherein the more than R targets include three targets, and wherein the average levels of the three targets are calculated based on light detected from only two fluorophores associated with probes that bind to amplicons of the three targets during amplification.

2. The method of claim 1, wherein amplification of each target of the more than R targets is reported by a different signal or combination of the signals than any of the other more than R targets individually.

3. The method of claim 1, wherein each of the signals reports amplification of a different combination of at least two of the more than R targets.

4. The method of claim 3, wherein the partitions are droplets, further comprising a step of determining a number of droplets exhibiting amplification of any of the at least two targets for each signal, and wherein the step of calculating is based on the number determined for each of the R signals.

5. The method of claim 1, wherein the two fluorophores are VIC and FAM.

6. The method of claim 5, wherein the partitions contain a first probe for a first of the three targets, a second probe for a second of the three targets, and a third probe for a third of the three targets, and wherein the first probe is labeled exclusively with VIC, the second probe is labeled exclusively with FAM, and the third probe is labeled with both VIC and FAM.

7. The method of claim 6, wherein the third probe includes a FAM-labeled probe that is not labeled with VIC and a VIC-labeled probe that is not labeled with FAM.

8. The method of claim 1, wherein the average level is a concentration.

9. The method of claim 1, wherein the step of calculating an average level includes a step of determining a total number of amplification-positive partitions for each type of the more than R targets and a step of determining a total number of partitions.

10. The method of claim 1, further comprising a step of distributing copies of the more than R targets among the partitions such that some partitions contain more than one copy of a given target.

11. The method of claim 1, wherein there are at least four targets.

12. A method of performing a multiplexed digital amplification assay, the method comprising:
amplifying targets in partitions;
detecting light from the partitions; and
calculating levels of at least three of the targets based on light detected from only two fluorophores, wherein the levels calculated account for a coincidence of all possible combinations of the at least three targets in the same individual partitions.

13. The method of claim 12, wherein the partitions contain a first probe for a first target of the at least three targets, a second probe for a second target of the at least three targets, and a third probe for a third target of the at least three targets, and wherein the first probe is labeled exclusively with a first fluorophore, the second probe is labeled exclusively with a second fluorophore, and the third probe is labeled with the first fluorophore and the second fluorophore.

14. The method of claim 13, wherein the third probe includes a probe labeled with the first fluorophore and not the second fluorophore and another probe labeled with the second fluorophore and not the first fluorophore.

15. The method of claim 13, wherein the first fluorophore is VIC and the second fluorophore is FAM.