Method and system for identifying a fluorochrome panel - Patents.com

Abstract

A method is provided for identifying a fluorochrome panel suitable for use in a flow cytometry protocol. The method of interest includes receiving an instrument identifier and a request for N fluorochrome identifiers, selecting two or more subsets of the N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier, and identifying a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. Systems and non-transitory computer readable storage media for carrying out the invention are further provided.

Description

[Background technology]

[0001] Characterization of analytes in biological fluids has become an important part of biological research, medical diagnosis, and the evaluation of a patient's overall health and wellness. Detecting analytes in biological fluids, such as human blood or blood-derived products, can provide results that may be involved in determining treatment protocols for patients with various medical conditions.

[0002] Particle analysis (e.g., flow cytometry) is a technique used to characterize and often sort biological material, such as cells in a blood sample or particles of interest in another type of biological or chemical sample. A flow cytometer typically includes a sample reservoir for receiving a fluid sample, such as a blood sample, and a sheath reservoir containing a sheath fluid. The flow cytometer carries particles (including cells) in the fluid sample as a cell stream to a flow cell, while further directing the sheath fluid to the flow cell. To characterize the components in the flow stream, light is illuminated into the flow stream. Changes in the biological material in the flow stream, such as morphology or the presence of fluorescent labels, may change the light observed, and these changes allow for characterization and separation. To characterize the components in the flow stream, light must be directed into the flow stream and collected. The light source in the flow cytometer can be a variety of light sources and may include one or more broad spectrum lamps, light emitting diodes, and single wavelength lasers. The light source is aligned with the flow stream, and optical responses from the illuminated particles are collected and quantified.

[0003] Parameters measured using particle analyzers typically include light at the excitation wavelength that is scattered by the particles at narrow angles, primarily in the forward direction, referred to as forward scatter (FSC), excitation light that is scattered by the particles in a direction perpendicular to the excitation laser, referred to as side scatter (SSC), and light emitted from fluorescent molecules or dyes. Different cell types can be distinguished by their light scattering properties and by fluorescent emission by labeling different cellular proteins or other components with fluorochrome-labeled antibodies or other fluorescent probes. The forward scattered light, side scattered light, and fluorescent light are detected by photodetectors located within the particle analyzer. Summary of the Invention [Problem to be solved by the invention]

[0004] When a flow cytometry protocol involves the detection of fluorescence, the experimental design typically involves the identification of a fluorochrome panel, i.e., a collection of fluorochromes that will be used together in a given flow cytometry workflow. The process of fluorochrome panel design is necessary because biological resolution (i.e., the ability to distinguish different components of interest within or between particles of interest) is directly affected by both the measurement variance of the flow cytometry "raw" data and the mathematical process of spectral correction or unmixing. Both of these factors are highly dependent on the choice of fluorochromes in the panel. First, measurement variance (i.e., noise) in flow cytometry arises from a wide range of sources, including a constant baseline measurement noise in the flow cytometer electronics, optical shot noise that varies linearly with signal intensity, and multiplicative measurement noise that varies quadratically with signal intensity and arises from random fluctuations in the flow cytometer lasers and fluids. The measurement noise itself is determined by the choice of fluorochromes. For example, bright fluorochromes produce more shot noise than dimmer fluorochromes, and the signal magnitude of dimmer fluorochromes is smaller compared to the constant "noise floor" of the instrument's optics and electronics. The raw measurement noise in "detector space" (with a number of dimensions equal to the number of detectors in the instrument) is then propagated to the final biological data in "correction space" or "unmixing space" (with a number of dimensions equal to the number of fluorochromes in the sample) by the mathematical process of fluorescence correction (in the case of conventional cytometers) or spectral unmixing (in the case of full-spectrum cytometers). The variance in "unmixing space" is important because it is the space in which the final biological analysis of interest (e.g. gating, clustering, sorting, marker quantification, etc.) is performed. This mathematical mapping of noise into biological space is heavily driven by the spectral characteristics of the fluorochromes themselves.

[0005] Conventional flow cytometry, in which separate photodetectors are dedicated to dye-specific fluorescence emission bands, places severe limitations on the number of fluorochromes that may be used simultaneously in a flow experiment, which may not exceed the number of fluorescence detection channels of the instrument. In contrast, full-spectrum flow cytometers naturally use many more detectors than fluorochromes, with over 180 fluorescence channels available in commercially available full-spectrum flow cytometers.

[0006] The present inventors have recognized that there is a practical limit to the number of fluorochromes that may be used simultaneously in a flow cytometry experiment. Despite the commercial availability of approximately 100 distinct fluorochromes for flow cytometry, the size of the panel is still limited. This practical limit arises from the inevitable spectral overlap and similarity of the fluorochromes used. Thus, a method and system for selecting a suitable fluorochrome panel is desirable. Embodiments of the present invention fulfill this need. [Means for solving the problem]

[0007] An aspect of the invention includes a method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol. The method of interest includes receiving an instrument identifier and a request for N fluorochrome identifiers, where N is an integer of at least 2. The instrument identifier may be associated with a spectral matrix associated with the instrument on which the given flow cytometry protocol is to be performed, for example. The spectral matrix described herein has spectral features associated with each fluorochrome identifier in the collection of fluorochrome identifiers. In some cases, the spectral matrix further has autofluorescence spectral features (i.e., spectral features associated with unique fluorescent signals generated by particles such as cells when measured in a flow cytometer). The spectral features of the spectral matrix may include, for example, spillover values ​​and may be derived from multiple sources (e.g., experimental and simulation data). The number of panels (i.e., N) may have any convenient value and may be in the range of 2 to 100, for example, 2 to 75, in some cases. The disclosed method further comprises selecting, from the set of fluorochrome identifiers in the spectral matrix, two or more subsets of N fluorochrome identifiers and identifying a fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. In some embodiments, the method is performed without evaluating antigen data of the sample.

[0008] In some cases, the numerical stability indicator is a condition number. In some cases, the condition number may be calculated by determining the singular values ​​associated with a spectral matrix that includes the spectral features associated with the N fluorochrome identifiers of the subset, and determining the ratio of the maximum singular value to the minimum singular value. In some embodiments, the condition number may be calculated by determining the matrix norm (e.g., L 2For example, embodiments of the subject method may calculate the Frobenius norm for each of two or more subsets of N fluorochrome identifiers, determine the Frobenius norm of the pseudoinverse of each of the two or more subsets of N fluorochrome identifiers, and determine, for each of the two or more subsets of N fluorochrome identifiers, the product of the Frobenius norm of the subset of N fluorochrome identifiers and the Frobenius norm of the pseudoinverse of the subset of N fluorochrome identifiers.

[0009] In some embodiments, the fluorochrome panel has an optimized numerical stability index (e.g., minimum condition number). For example, in one variation, identifying the fluorochrome panel includes a constrained optimization method (e.g., a minimization algorithm). In an embodiment, the constrained optimization method is performed sequentially. For example, when selecting two or more subsets of the N fluorochrome identifiers, a first subset of the N fluorochrome identifiers and a second subset of the N fluorochrome identifiers may be selected. The method may further include calculating a numerical stability index for each of the first subset of the N fluorochrome identifiers and the second subset of the N fluorochrome identifiers, and comparing the numerical stability index of the first subset of the N fluorochrome identifiers to the numerical stability index of the second subset of the N fluorochrome identifiers. The first subset of the N fluorochrome identifiers may be obtained, possibly randomly, from the spectral matrix. In one variation, when generating the second subset of the N fluorochrome identifiers, fluorochrome identifiers in the first subset of the N fluorochrome identifiers are replaced with alternative fluorochrome identifiers in the spectral matrix. The disclosed method may further include iteratively generating subsets of fluorochrome identifiers, calculating a numerical stability index for each iteratively generated subset of N fluorochrome identifiers, and comparing the calculated numerical stability index for each iteratively generated subset of N fluorochrome identifiers with a numerical stability index of a previously generated subset of N fluorochrome identifiers to determine an optimized numerical stability index. In some cases, the optimized numerical stability index is a minimized condition number.

[0010] The method of interest may further include evaluating a maximum number of fluorochromes that may be used simultaneously in a given fluorochrome panel. In some embodiments, when evaluating the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel, an optimized numerical stability index is determined for each of a plurality of values ​​of N. For example, the plurality of values ​​of the number of panels may include i values, where i is, for example, the number of fluorochrome identifiers in the spectral matrix. In some cases, when evaluating the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel, a maximum value of the number of panels that does not result in the optimized numerical stability index exceeding a threshold is determined.

[0011] Aspects of the invention further include a system. The system of interest comprises a processor configured to perform the method of the invention, namely, receiving an instrument identifier and a request for N fluorochrome identifiers, selecting two or more subsets of the N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier, and identifying a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. In some embodiments, the system comprises a particle analyzer, such as a flow cytometer (e.g., a full spectrum flow cytometer). A non-transitory computer readable storage medium having instructions stored thereon for identifying a fluorochrome panel suitable for use in a flow cytometry protocol is further provided. [Brief description of the drawings]

[0012] The invention may be best understood from the following detailed description when read in conjunction with the accompanying drawing figures, in which:

[0013] [Figure 1A] 1 is a chart showing plots of flow cytometer data characterized by a fluorochrome panel with low condition numbers. [Figure 1B]1 is a chart showing plots of flow cytometer data characterized by a fluorochrome panel with relatively high condition numbers. [Diagram 2] 1 is a flow chart for implementing a method for identifying a fluorochrome panel according to an embodiment. [Figure 3-1] FIG. 1 illustrates the identification of fluorochrome panels using an exemplary spectral matrix. [Figure 3-2] FIG. 1 illustrates the identification of fluorochrome panels using an exemplary spectral matrix. [Figure 4] FIG. 1 is a functional block diagram illustrating a flow cytometry system according to an embodiment. [Diagram 5] FIG. 1 illustrates a sorting control system according to an embodiment. [Figure 6A] FIG. 1 is a schematic diagram illustrating a particle sorting system according to an embodiment. [Figure 6B] FIG. 1 is a schematic diagram illustrating a particle sorting system according to an embodiment. [Figure 7] FIG. 1 is a block diagram illustrating a computing system according to an embodiment. [Figure 8] 1 is a chart showing a comparison of numerical stability indices obtained for identified fluorochrome panels according to an embodiment of the present invention and for randomly determined fluorochrome panels; [Figure 9] 13 is a chart showing the improvement of numerical stability metrics by iteratively generating intermediate subsets of fluorochrome identifiers. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0014] A method is described for identifying a fluorochrome panel suitable for use in a flow cytometry protocol. The method of interest includes receiving an instrument identifier and a request for N fluorochrome identifiers, selecting two or more subsets of the N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier, and identifying a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. Systems and non-transitory computer readable storage media for implementing the method are also described.

[0015] Before the present invention is described in more detail, it should be understood that the invention is not limited to the particular embodiments described, as such may, of course, vary. The scope of the present invention will be limited only by the appended claims, and it should be further understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.

[0016] When a range of values ​​is given, it is understood that each intervening value between the upper and lower limits of that range, to one-tenth of the unit of the lower limit unless the context clearly indicates otherwise, and any other stated or intervening value in that stated range, is encompassed within the invention. The upper and lower limits of these smaller ranges may be independently included in the smaller ranges and are further encompassed within the invention, subject to any specifically excluded limit in the stated range. Where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also encompassed within the invention.

[0017] In this specification, a range is provided with the term "about" before the numerical values. The term "about" is used herein to literally support the exact number that it precedes, as well as a number that is close to or approximately the number that it precedes. When determining whether a number is close to or approximately a specifically stated number, the unstated number that is close or approximately may be a number that, in the context in which the specifically stated number is presented, provides a substantial equivalent to the specifically stated number.

[0018] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present invention, representative and illustrative methods and materials are described.

[0019] All publications and patents cited herein are incorporated by reference as if each individual publication or patent was specifically and individually indicated to be incorporated by reference, and are incorporated by reference herein to disclose and describe the methods and / or materials in connection with which the publications are cited. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that the present invention is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided may be different from the actual publication dates which may need to be independently confirmed.

[0020] It should be noted that, as used in this specification and the appended claims, the singular forms "a," "an," and "the" include plural referents unless the context clearly dictates otherwise. It should be further noted that the claims may be drafted to exclude any optional element. As such, this statement is intended to serve as a prelude to the use of exclusive terminology, such as "solely," "only," or the use of a "negative" limitation in relation to the recitation of claim elements.

[0021] As will be apparent to one of ordinary skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein comprises separate components and features which may be readily separated from, or combined with, any of the features of the other multiple embodiments without departing from the scope or spirit of the invention. Any recited method may be carried out in the order of events recited or in any other order which is logically possible.

[0022] Although the systems and methods have been described or will be described for grammatical fluidity with functional descriptions, the claims should not be construed as necessarily limited in any way by limitations of "means" or "step" construction unless expressly recited under 35 U.S.C. 112, but should be accorded the full scope of meaning and equivalents of the definitions given by the claims under the judicial theory of equivalents, with the express understanding that if a claim is expressly recited under 35 U.S.C. 112, then the full statutory equivalents under 35 U.S.C. 112 should be accorded.

[0023] How to Identify Fluorochrome Panels As discussed above, aspects of the invention include methods of identifying a fluorochrome panel. As described herein, a "fluorochrome panel" refers to a collection of different fluorescent molecular substances (i.e., dyes) that may be used to identify particles in a sample or specific moieties or components that bind to particles in a sample. The distinct fluorochromes in a fluorochrome panel may differ with respect to properties such as absorption spectrum, extinction coefficient, emission spectrum and quantum efficiency (i.e., the number of photons emitted per absorption of one photon), or a combination thereof. Thus, the distinct or distinct fluorochromes may differ from each other in terms of chemical composition and / or in terms of one or more properties of the dye. For example, fluorochromes of a given pair may be considered distinct when they differ from each other in terms of excitation maxima and / or emission maxima, the magnitude of such difference being, in some cases, 5 nm or more, such as 10 nm or more, such as 15 nm or more, and in some cases the magnitude of such difference being within the range of 5-400 nm, such as 10-200 nm, such as 15-100 nm, such as 25-50 nm.

[0024] The fluorochrome panels of interest are suitable for use in flow cytometry protocols. A fluorochrome panel may be described as "suitable for use" in a flow cytometry protocol if it generates easy-to-understand flow cytometer data that provides reliable insight into the characteristics of interest of the sample under investigation. In some embodiments, a fluorochrome panel is suitable for use in a flow cytometry protocol if it improves biological resolution. "Biological resolution" refers to the ability to distinguish different entities of interest (e.g., molecules, antigens, moieties, epitopes, etc.) in a biological specimen. In some cases, the fluorochrome panels identified herein provide maximum biological resolution despite measurement variance and variance in the flow cytometer data space (e.g., flow cytometer data that has undergone fluorescence compensation or spectral unmixing). "Maximum" biological resolution is, in some variations, evaluated relative to the biological resolution achieved using one or more other sets of fluorochromes that are different from the fluorochromes identified as described herein (i.e., that include one or more different fluorochromes relative to the fluorochromes identified as described herein).

[0025] In some cases, the method is performed without evaluating antigen data. By "antigen data" is meant information regarding the antigens to be evaluated in a given flow cytometry protocol for which a fluorochrome panel is desired. In other words, in contrast to methods that consider the properties of the antigens to which they bind and select the fluorochromes accordingly, embodiments of the subject method are instead directed to providing a collection of fluorochromes (i.e., a "color palette") suitable for use in the same flow cytometry protocol. In these cases, the antigenicity of the sample is not taken into account during the identification of the fluorochrome panel. In one variation, after a fluorochrome panel of interest has been identified, the method may optionally further assign fluorochromes in the identified panel to specific antigens (e.g., via antibodies targeted to the relevant antigens).

[0026] As described above, the methods of the present invention include receiving an instrument identifier. By "instrument identifier" is meant information or data that refers to a particular instrument (e.g., particle analyzer, flow cytometer). In some cases, the instrument identified by the instrument identifier is a flow cytometer. Any convenient flow cytometer configured to analyze fluorescent particle-modulated light may be used. In some cases, flow cytometers of interest include flow cytometers manufactured by BD Biosciences.Exemplary flow cytometers include BD Biosciences FACSCanto™ flow cytometer, BD Biosciences FACSCanto™ II flow cytometer, BD Accuri™ flow cytometer, BD Accuri™ C6 Plus flow cytometer, BD Biosciences FACSCelesta™ flow cytometer, BD Biosciences FACSLyric™ flow cytometer, BD Biosciences FACSVerse™ flow cytometer, BD Biosciences FACSymphony™ flow cytometer, BD Biosciences LSRFortessa™ flow cytometer, BD Biosciences LSRFortessa™ X-20 flow cytometer, BD Biosciences FACSPresto™ flow cytometer, BD Biosciences FACSVia™ flow cytometer, as well as BD Biosciences FACSCalibur™ cell sorter, BD Biosciences FACSCount™ cell sorter, BD Biosciences FACSLyric™ cell sorter, BD Biosciences These include the Via™ cell sorter, BD Biosciences Influx™ cell sorter, BD Biosciences Jazz™ cell sorter, BD Biosciences Aria™ cell sorter, BD Biosciences FACSAria™ II cell sorter, BD Biosciences FACSAria™ III cell sorter, BD Biosciences FACSAria™ Fusion cell sorter, and BD Biosciences FACSMelody™ cell sorter, BD Biosciences FACSymphony™ S6 cell sorter, etc.

[0027] An instrument identifier may be associated with, for example, a spectral matrix associated with the instrument referenced by the instrument identifier. As described herein, a "spectral matrix" refers to a matrix that includes information about the spectral properties of the fluorochromes of a collection. In an embodiment, a spectral matrix of interest includes a collection of spectral features associated with each fluorochrome identifier in a collection of fluorochrome identifiers. A "spectral feature" refers to a characteristic of the fluorescence spectrum of an individual fluorochrome, expressed as one or more numerical values. In other words, a spectral matrix associated with an instrument identifier (also referred to herein as an "input" or "original" spectral matrix) may represent a "palette" of fluorochromes or dyes from which a fluorochrome panel of interest may be selected. In some cases, a spectral matrix includes information related to the spectral properties of each fluorophore of an instrument of interest (e.g., a flow cytometer). Such information may vary from machine to machine. For example, differences in the number and arrangement of lasers and detection channels between distinct instruments may result in different spectral features associated with each instrument. In some embodiments, the spectral features are obtained from experimental data (i.e., the results of experiments performed on a particular instrument). In other cases, the spectral features are obtained from simulated data. The input spectral matrix may be associated with any suitable number of fluorochrome identifiers. In some embodiments, the number of fluorochrome identifiers in the input spectral matrix is ​​in the range of 2-150, such as 2-140, such as 2-130, such as 2-120, such as 2-110, such as 2-100, such as 2-90, such as 2-80, such as 2-70, such as 2-60, such as 2-50. The collection of fluorochrome identifiers used in a given embodiment of the invention may be referred to as a fluorochrome palette, which collectively refers to the fluorochrome identifiers that may be selected for a given candidate fluorochrome panel.

[0028] In some cases, the spectral features have one or more spillover values. By "spillover value" is meant the relative amount of signal that a given fluorochrome emits into each detector band. In some cases, the spillover value is normalized to the detector with the maximum signal for that fluorochrome (i.e., the "peak" detector). In some cases, particle modulated light indicative of a particular fluorochrome is obtained at one or more detectors in the particle analyzer (e.g., a flow cytometer) that are not the peak detector or detectors for that fluorochrome. Thus, the light may "spillover" and be detected by the off-peak detector. In other words, the particular fluorescent label used in an experiment and the associated fluorescent emission band of the fluorescent label may be selected to generally coincide with a particular detector. However, as more detectors are provided and more labels are utilized, perfect correspondence between a particular detector and a fluorescent emission spectrum may not be possible. Although the peak of the emission spectrum of a particular fluorescent molecule may lie within the window of one particular detector, it is generally true that a portion of the emission spectrum of that label also overlaps with the window of one or more other detectors. This may be referred to as spillover.

[0029] In some embodiments, the spectral matrices described herein may include one or more autofluorescence spectral features. Autofluorescence is an intrinsic fluorescent signal generated by particles such as cells when measured by a flow cytometer. Autofluorescence arises from endogenous molecules with fluorescent activity, such as metabolites within the cell. Different cells of the same type (e.g., lymphocytes) may have the same autofluorescence spectrum, but at different intensities, e.g., larger cells generally tend to have a larger autofluorescence signal. In some cases, different types of particles are associated with different autofluorescence spectra. For example, different types of cells (e.g., lymphocytes vs. monocytes) may have different autofluorescence spectra as well as various autofluorescence levels (e.g., the spectral features of lymphocyte autofluorescence may be different from the spectral features of monocyte autofluorescence). In some cases, e.g., in the case of spectral cytometry, the spectral features of autofluorescence are measured by looking at unstained cells, and when multiple autofluorescence spectra are included, they are included in the spectral unmixing process as an additional one or more "fluorochrome" parameters.

[0030] The method of interest further includes receiving a request for N fluorochrome identifiers. N, also referred to herein as the "panel number", is an integer of at least 2 and defines the size of the fluorochrome panel to be identified. The size of the fluorochrome panel to be identified may vary and is determined based on the requirements of a particular experimental design. For example, N may be determined by one of ordinary skill in the relevant art (e.g., molecular biology, flow cytometry, fluorescence microscopy, etc.) based on the number of analytes to be investigated. In some cases, the size of the fluorochrome panel (and the value of N) may be in the range of 2-150 fluorochromes, such as 2-140 fluorochromes, such as 2-130 fluorochromes, such as 2-120 fluorochromes, such as 2-110 fluorochromes, such as 2-100 fluorochromes, such as 2-110 fluorochromes, such as 2-12 ... the number of fluorescent dyes may be within the range of, for example, 2 to 95 fluorescent dyes, for example, 2 to 90 fluorescent dyes, for example, 2 to 85 fluorescent dyes, for example, 2 to 80 fluorescent dyes, for example, 2 to 75 fluorescent dyes, for example, 2 to 70 fluorescent dyes, for example, 2 to 65 fluorescent dyes, for example, 2 to 60 fluorescent dyes, for example, 2 to 55 fluorescent dyes, for example, 2 to 50 fluorescent dyes, for example, 2 to 45 fluorescent dyes, for example, 2 to 40 fluorescent dyes, for example, 2 to 40 fluorescent dyes, for example, 2 to 35 fluorescent dyes, for example, 2 to 30 fluorescent dyes, for example, 2 to 25 fluorescent dyes, for example, 2 to 20 fluorescent dyes, for example, 2 to 15 fluorescent dyes, for example, 2 to 10 fluorescent dyes, for example, 2 to 5 fluorescent dyes. In some cases, the fluorochrome panel has 2 or more fluorochromes, 3 or more fluorochromes, 4 or more fluorochromes, 5 or more fluorochromes, 6 or more fluorochromes, 7 or more fluorochromes, 8 or more fluorochromes, 9 or more fluorochromes, 10 or more fluorochromes, 11 or more fluorochromes, 12 or more fluorochromes, 13 or more fluorochromes, 14 or more fluorochromes, 15 or more fluorochromes, 16 or more fluorochromes, 17 or more fluorochromes, 18 or more fluorochromes, 19 or more fluorochromes, or 20 or more different fluorochromes.

[0031] The embodiment of the invention further comprises selecting two or more subsets of N fluorochrome identifiers from the set of fluorochrome identifiers in the spectral matrix associated with the instrument identifier, i.e., the fluorochrome palette. In other words, the two or more subsets each have N fluorochrome identifiers from the spectral matrix, where the value of N is the same in each case. Since the two or more subsets have N fluorochrome identifiers, each subset has the same number of fluorochrome identifiers as the number of panels (i.e., N). In an example where a request for five fluorochrome identifiers (i.e., N=5) is received, the two or more subsets each have five fluorochrome identifiers. In other words, the number of fluorochrome identifiers in each of the two or more subsets of fluorochrome identifiers corresponds to the number of fluorochrome identifiers requested and the number of fluorochromes in the fluorochrome panel being identified. Neither the two or more subsets of fluorochrome identifiers nor the fluorochrome panel includes a number of fluorochrome identifiers different from the number of panels (i.e., N).

[0032] Each of the fluorochrome identifiers in the two or more subsets of the N fluorochrome identifiers is associated with a spectral feature in the same way that the fluorochrome identifiers in the spectral matrix are associated with a spectral feature. Two entities may be described as "associated" with each other if they are linked and / or associated with each other in data space. Thus, in some cases, each of the two or more subsets of the N fluorochrome identifiers may be associated with a spectral matrix having spectral features that are associated with only the fluorochromes in the subset of the total set of possible fluorochromes in the original spectral matrix. In other words, in some embodiments, the two or more subsets of the N fluorochrome identifiers are associated with a submatrix of the original (i.e., input) spectral matrix associated with the instrument identifier. A "submatrix" is used herein in the conventional sense to describe a matrix obtained by deleting a certain combination of rows and / or columns of another matrix. In some cases, the two or more subsets of the N fluorochrome identifiers are associated with a submatrix of the original spectral matrix (i.e., the spectral matrix associated with the instrument identifier as described above) in which there are no rows and / or columns corresponding to non-selected fluorochromes.

[0033] The method further includes, after selecting the two or more subsets of the N fluorochrome identifiers, identifying a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on the numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers. As used herein, a "numerical stability index" refers to a quantitative measure of the numerical stability of a given set or subset of fluorochrome identifiers, i.e., a measure of the stability and accuracy of flow cytometer data obtained when fluorochromes corresponding to the fluorochrome identifiers of a given subset are analyzed with a particular instrument.

[0034] As described above, embodiments of the spectral matrix of interest include autofluorescence spectral features. If the autofluorescence characteristics of a particular sample to be analyzed via a fluorochrome panel of interest are known, the method of interest may identify a fluorochrome panel that takes said autofluorescence characteristics into account. In other words, such a method may identify a fluorochrome panel having fluorochromes characterized by excitation and / or emission spectra that overlap minimally with the autofluorescence spectrum. For example, the method may include spectral features related to autofluorescence (e.g., spillover values) when calculating a numerical stability index. The resulting numerical stability index may then be used to select a fluorochrome panel suitable for use in a flow cytometry protocol analyzing samples with a particular fluorescence profile (i.e., with improved biological resolution).

[0035] In some embodiments, the numerical stability index is a condition number. The "condition number" provides an index by which the similarity of the fluorescence spectra and the impact of this similarity on the biological resolution may be evaluated. In an embodiment, a completely stable set or subset of fluorochrome identifiers has a condition number of 1, and a completely unstable set or subset of fluorochrome identifiers has a condition number of infinity. Thus, if two subsets of N fluorochrome identifiers are associated with different condition numbers, the subset of N fluorochrome identifiers with the lower condition number will have less spectral overlap. Therefore, when the numerical stability index is a condition number, the method calculates the condition number for each of two or more subsets of N fluorochrome identifiers.

[0036] The condition number may be computed via any convenient protocol. For example, given a linear equation Ax=b, the condition number may be defined as the maximum ratio of the relative error in x to the relative error in b. Assuming A is a regular matrix and e is defined as the error in b, then the solution A -1 The error of b is A -1 e. The ratio of the relative error of the solution to the relative error of b is calculated as follows (from https: / / en(dot)Wikipedia(dot)org / wiki / Condition_number#Matrices):

[0037]

number

[0038] As a result, the maximum non-zero values ​​of e and b are the product of the two operator norms:

[0039]

number

[0040] In some cases, the numerical stability index takes the form of a condition number of the spectral matrix, calculated with respect to the matrix norm. Any suitable matrix norm may be used to calculate the condition number. For example, the above definition may be used for any consistent norm, i.e. a norm that satisfies:

[0041]

number

[0042] In embodiments, the calculation of the condition number includes singular value decomposition (SVD). In such embodiments, the condition number is calculated by determining the singular values ​​of the spectral features (i.e., the singular values ​​of the spectral matrix associated with the relevant subset of fluorochrome identifiers). In some embodiments, the numerical stability index is calculated by the derivative L (also referred to as the Euclidean norm) as follows: 2 It takes the form of the condition number of the spectral matrix computed with respect to the norm. κ(A)=(σ max (A)) / (σ min (A) Here, σ max (A) and σ min (A) are the maximum and minimum singular values ​​of A, respectively.

[0043] In further embodiments, the numerical stability indicator takes the form of a condition number calculated with respect to any other induced norm, including but not limited to the 1-norm, 2-norm, or infinity-norm. In some embodiments, the calculation of the condition number includes the calculation of the Frobenius norm. As understood in the field of linear algebra, the Frobenius norm is a matrix norm defined as the square root of the sum of the squares of the elements of a matrix. The Frobenius norm can be found, for example, in Golub, GH and Van Loan, CF (1996) Matrix Computations, 3, vol. 14, No. 1, pp. 1171-1175, 2002, which is incorporated herein by reference in its entirety. rded. Baltimore, MD: Johns Hopkins. In some cases, the Frobenius norm is calculated as follows (adapted from Golub and Van Loan):

[0044]

number

[0045] where A is an m×n matrix. In an embodiment, when calculating the condition number, the spectral features associated with the N fluorochrome identifiers of the two or more subsets are used to determine the Frobenius norm of each of the two or more subsets of the N fluorochrome identifiers. For example, if each of the two or more subsets of the N fluorochrome identifiers is included in a submatrix of the original spectral matrix (e.g., as described above), the method obtains the Frobenius norm of each submatrix.

[0046] In some cases, calculating the condition number further involves determining the Frobenius norm of a pseudoinverse for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. As is known in the art of linear algebra, a "pseudoinverse" is a generalization of the inverse of an m by n (i.e., non-square) matrix. In some cases, the pseudoinverse may be calculated as follows: A + =(A T A) -1 A T where A is an m × n matrix, and A T is the transpose of A, and A + is the pseudoinverse. After computing the pseudoinverse, one may obtain the Frobenius norm of said pseudoinverse (e.g., computed as shown above). In a further case, computing the condition number further comprises, for each of two or more subsets of N fluorochrome identifiers, determining the product of the Frobenius norm of the subset of N fluorochrome identifiers and the Frobenius norm of the pseudoinverse of the subset of N fluorochrome identifiers.

[0047] Methods for calculating condition numbers are further described in Spidlen et al. (2016) Data File Standard for Flow Cytometry. International Society for Advancement of Cytometry. Version 4.0, the entire disclosure of which is incorporated herein by reference. In some cases, a spectral matrix (e.g., a submatrix) may be decomposed into a product of three matrices (adapted from Spidlen et al.) as follows: WM=UwV T where U is an a×b rectangular matrix (a>b), M is an a×b rectangular matrix, and W is an a×a diagonal matrix (with diagonal components W i is the weight, e.g. 1 / the estimated standard error of the detector (W i =1 / σ i )), V T are U and V T is column-orthogonal, i.e., a b × b square matrix (source Spidlen et al.), U T U=VV T =1 w is the diagonal component w i The inverse of w is also a diagonal matrix with diagonal entries equal to the reciprocals of the singular values ​​(w i -1 =1 / w). The general solution may therefore be written as (adapted from Spidlen et al.): F = (Vw -1 U T W)D where W can be omitted or set to 1 for no weights, and the unmixing matrix (V -1 U T W) can be computed once and applied to many measurement vectors D. In an embodiment, when computing the condition number, the ratio of the largest singular value to the smallest singular value is determined.

[0048] In some cases, more than one of the above methods (e.g., L 2In some cases, the hybrid numerical stability index is calculated by combining multiple numerical stability indices calculated using a set of standardized norms (e.g., 1-norm, 2-norm, infinity-norm, Frobenius norm, etc.). The hybrid numerical stability index may be calculated via any convenient protocol. For example, in some cases, the hybrid numerical stability index is calculated by averaging multiple numerical stability indices.

[0049] 1A and 1B are diagrams that demonstrate how flow cytometer data generated using a fluorochrome panel with a low condition number are characterized by lower spread compared to fluorochrome panels with a relatively higher condition number. When a set of overly similar dyes is used simultaneously in a spectral flow cytometry experiment, the variance of the data can become "bulky". This manifests as severe broadening in the data regardless of dye co-expression, and is most pronounced in the negative population (where no dye is present). This "unmixing broadening" is a consequence of the numerical stability of the spectral matrix, which may be described by the condition number. As shown in FIG. 1A, flow cytometer data unmixed with a spectral matrix with a low condition number (i.e., 3.7) (i.e., a submatrix of the original input matrix associated with the instrument identifier) ​​has a low spread, especially in the negative population. On the other hand, FIG. 1B shows that flow cytometer data unmixed with a spectral matrix with a high condition number has a high spread in the negative population.

[0050] In some embodiments, the fluorochrome panel has an optimized numerical stability index. By "optimized" it is meant that the method identifies a subset of fluorochrome identifiers associated with fluorochromes that have a small (e.g., smallest) amount of spectral overlap given certain constraints (e.g., number of panels, runtime, etc.). In some embodiments, the method identifies a fluorochrome panel that is associated with an absolute minimum of possible numerical stability indices given a palette of fluorochromes present in the original spectral matrix and the associated spectral features of the fluorochromes. In other embodiments, the method identifies a fluorochrome panel that is associated with a numerical stability index that approximates the lowest possible numerical stability index, but is not necessarily the lowest possible numerical stability index. For example, in some embodiments, the method identifies multiple fluorochrome panels that are each suitable for use in a particular flow cytometry protocol. In some cases, the method identifies a number of fluorochrome panels in the range of 1-20, such as 1-10, such as 1-5, such as 1-3. When identifying multiple fluorochrome panels, not all fluorochrome panels may be associated with the lowest possible numerical stability index. Alternatively, one or more of the fluorochrome panels may approximate the lowest possible numerical stability index. In this case, a numerical stability index that approximates the lowest possible numerical stability index may still be considered to be "optimal".

[0051] In some embodiments, constrained optimization methods are used in identifying the fluorochrome panel. "Constrained optimization" is used herein in the conventional sense to describe the process of optimizing variables when those variables are constrained. Any suitable constrained optimization method may be used. In some cases, the constrained optimization method is a minimization algorithm. By "minimization algorithm" is meant a type of constrained optimization method that seeks to minimize a particular variable. When the constrained optimization method is a minimization algorithm, the variable that is minimized is a numerical stability index. For example, in some cases, the variable that is minimized is a condition number (e.g., calculated as described above). Thus, in some embodiments, the constrained optimization method seeks to obtain a subset of N fluorochrome identifiers from within the original (i.e., input) spectral matrix that has a minimized numerical stability index (e.g., condition number). Examples of constrained optimization techniques include local search, local repair, backtracking, and constraint propagation. These may be combined with minimization techniques such as simulated annealing and genetic (evolutionary) algorithms.

[0052] In some embodiments, the constrained optimization method is performed sequentially. In other words, the constrained optimization method selects one subset of N fluorochrome identifiers, and selects another subset of N fluorochrome identifiers based on the calculation results performed on the first subset of N fluorochrome identifiers, and then performs the same. For example, in an embodiment, when selecting two or more subsets of N fluorochrome identifiers, a first subset of N fluorochrome identifiers and a second subset of N fluorochrome identifiers are selected. The first subset of N fluorochrome identifiers is selected before the second subset of N fluorochrome identifiers and may be obtained via any convenient protocol. In some embodiments, the first subset of N fluorochrome identifiers is randomly obtained from the spectral matrix associated with the instrument identifier (i.e., the original input spectral matrix). The method of interest further includes, after identifying the first subset of N fluorochrome identifiers, calculating a numerical stability index associated with the selected first subset of N fluorochrome identifiers based on the spectral features associated with the N fluorochrome identifiers of the first subset. In other words, a numerical stability index (eg, condition number) of a spectral matrix having spectral features associated with the N fluorochrome identifiers of the first subset is calculated.

[0053] After selecting and analyzing the first subset of N fluorochrome identifiers, a second subset of N fluorochrome identifiers may then be generated. In an embodiment, the second subset of N fluorochrome identifiers is a variation of the first subset of N fluorochrome identifiers. For example, in some cases, when generating the second subset of N fluorochrome identifiers, fluorochrome identifiers in the first subset of N fluorochrome identifiers are replaced with alternative fluorochrome identifiers in the spectral matrix associated with the instrument identifier. In other words, the constrained optimization method may switch fluorochrome identifiers for different fluorochrome identifiers in the original (i.e., input) spectral matrix representing the fluorochromes or palette of dyes from which to select. In some cases, fluorochromes in the first subset of N fluorochrome identifiers are not replaced with fluorochrome identifiers already in the first subset of N fluorochrome identifiers. The method of interest further includes, after generating the second subset of N fluorochrome identifiers, calculating a numerical stability index based on the spectral features associated with the N fluorochrome identifiers of the second subset. In other words, calculate a numerical stability index (eg, condition number) of a spectral matrix having spectral features associated with the N fluorochrome identifiers of the second subset.

[0054] The method of interest further comprises comparing the numerical stability metric of the first subset of N fluorochrome identifiers to the numerical stability metric of the second subset of N fluorochrome identifiers. For example, if the numerical stability metric is a condition number, the method may determine whether the first subset of N fluorochrome identifiers or the second subset of N fluorochrome identifiers is associated with a lower condition number. If the condition number associated with the first subset of N fluorochrome identifiers is less than the condition number associated with the second subset of N fluorochrome identifiers, discard the second subset of N fluorochrome identifiers as a less optimal combination compared to the first subset of N fluorochrome identifiers. Alternatively, if the condition number associated with the second subset of N fluorochrome identifiers is less than the condition number associated with the first subset of N fluorochrome identifiers, discard the first subset of N fluorochrome identifiers as a less optimal combination of fluorochromes compared to the second subset of N fluorochrome identifiers.

[0055] In some cases, the above sequential process of computing and evaluating subsets of N fluorochrome identifiers is performed iteratively. In such cases, the constrained optimization method further generates subsets of fluorochrome identifiers iteratively. Each iteration of the constrained optimization may include generating any convenient number of different subsets of the N fluorochrome identifiers. In some cases, a given iteration of the method may include generating one subset of the N fluorochrome identifiers, generating two subsets of the N fluorochrome identifiers, generating three subsets of the N fluorochrome identifiers, generating four subsets of the N fluorochrome identifiers, or generating five or more subsets of the N fluorochrome identifiers. For example, in some embodiments, either the first or second (or third or fourth, etc., if applicable) of the N fluorochrome identifiers associated with a relatively low condition number (e.g., as determined above) may serve as a seed for the next portion of the iterative process. "Seed" refers to a subset of N fluorochrome identifiers that is determined in one iteration of the method to be associated with a more optimized numerical stability index compared to one or more subsets of N fluorochrome identifiers that have been slightly altered (e.g., swapped as described above). For example, if the condition number associated with a first subset of N fluorochrome identifiers is smaller than the condition number associated with a second subset of N fluorochrome identifiers, then the first subset of N fluorochrome identifiers may serve as a seed for the next part of the iterative process.

[0056] The method of interest further comprises calculating a numerical stability index for each iteratively generated subset of the N fluorochrome identifiers, and comparing the calculated numerical stability index for each iteratively generated subset of the N fluorochrome identifiers with the numerical stability index of the subset of the N fluorochrome identifiers already generated to determine an optimized numerical stability index. In other words, after generating each seed subset and one or more subsets of the N fluorochrome identifiers generated by replacing the fluorochrome identifiers of each seed subset with alternative fluorochrome identifiers in the spectral matrix, the method calculates a numerical stability index for each subset, and then compares the calculated numerical stability indexes. The subset with the more optimal numerical stability index serves as a seed for the next part (i.e., iteration) of the iterative process.

[0057] In some embodiments, the iterative process repeats itself until a condition is met. Any suitable condition may be used to terminate the iterative process. In some cases, the iterative process terminates after a certain runtime has elapsed. In other cases, the iterative process terminates when a certain numerical stability index threshold is exceeded. For example, if the numerical stability index is a condition number, the iterative process may terminate when the constrained optimization method generates a subset of N fluorochrome identifiers that is associated with a condition number smaller than a predetermined value. In still other cases, the iterative process terminates when a numerical stability index is determined and evaluated for each possible combination of N fluorochrome identifiers or for a predetermined number of combinations of N fluorochrome identifiers. In still other cases, the iterative process terminates when the numerical stability index calculated for the iteratively generated subset of N fluorochrome identifiers converges. In other words, the iterative process terminates when only small changes in the numerical stability index are observed between subsequent subsets of N fluorochrome identifiers.

[0058] FIG. 2 shows a sample flow chart for identifying a fluorochrome panel according to an embodiment of the present invention. The method, in step 201, receives an input spectral matrix (201a) and a panel number (201b). The input spectral matrix is ​​associated with an instrument identifier and thus relates to a particular machine (e.g., particle analyzer, flow cytometer) as described above. In step 202, a constrained optimization method is used to prepare two or more subsets of N fluorochrome identifiers and evaluate a numerical stability metric for each of the two or more subsets of N fluorochrome identifiers. For example, in step 202a, a random subset of N fluorochrome identifiers is selected from the input spectral matrix. Then, in step 202b, a numerical stability metric (i.e., condition number) associated with the subset of N fluorochrome identifiers selected in step 202a is calculated. In step 202c, a minimization algorithm is used to select a second subset of the N fluorochrome identifiers, a numerical stability index is calculated for the second subset of the N fluorochrome identifiers, and the numerical stability index calculated for the second subset of the N fluorochrome identifiers is compared to the numerical stability index calculated in step 202b. This process is repeated (i.e., iterated) in step 202d until the numerical stability index converges. The resulting fluorochrome panel is then output in step 203.

[0059] FIG. 3 illustrates the identification of a fluorochrome panel according to the flow chart shown in FIG. 2. In step 301, an input spectral matrix 301a and a number of panels 301b are received. In the example of FIG. 3, the number of panels (i.e., N) is 4. The spectral matrix 301a includes the fluorochrome identifiers (top row: FL1, FL2, FL3...) and the detectors (leftmost column: U10, U9, U8...) that may detect light emitted from the fluorochromes. The spectral features (i.e., spillover values) associated with the fluorochrome identifiers are included in the spectral matrix 301a. In step 302, a subset of four fluorochrome identifiers (in this example, [FL2, FL3, FL5, FL6]) is selected. This subset is associated with a spectral matrix 302a, which is a submatrix of the input spectral matrix 301a that includes only data related to the selected subset of four fluorochrome identifiers (in this example, [FL2, FL3, FL5, FL6]). A numerical stability index (i.e., condition number 302b) is calculated based on the spectral matrix 302a associated with the subset of four fluorochrome identifiers. In step 302c, a minimization algorithm is used to select a second subset of N fluorochrome identifiers, a numerical stability index is calculated for the second subset of N fluorochrome identifiers, and the numerical stability index calculated for the second subset of N fluorochrome identifiers is compared to the numerical stability index calculated in step 302b. This process is repeated (i.e., iterated) in step 302d until the numerical stability index converges. In the example of FIG. 3, this occurs for the subset of fluorochrome identifiers including [FL2, FL3, FL5, FL6]. This subset is then output as fluorochrome panel 303.

[0060] In embodiments, the method further evaluates the maximum number of fluorochromes that may be used simultaneously in the fluorochrome panel. Such evaluation may be desirable when one seeks to use a large number of fluorochromes (e.g., 50 or more, 75 or more, 100 or more) at once in a flow cytometry protocol and wants to understand the feasibility of this experiment. In these embodiments, the method may determine an optimized numerical stability index for each of a plurality of values ​​of N. In other words, the above process for identifying fluorochrome panels may be repeated multiple times. In each case of the process, the number of panels (i.e., the value of N) varies. In some cases, the plurality of values ​​of N includes i values. Any convenient value of i may be used. For example, in some embodiments, i may be in the range of 2-150, such as 2-140, such as 2-130, such as 2-120, such as 2-110, such as 2-100, such as 2-90, such as 2-80, such as 2-70, such as 2-60, such as 2-50. In some embodiments, i is the number of fluorochrome identifiers in the spectral matrix associated with the instrument identifier.

[0061] In embodiments, when evaluating the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel, a maximum value of N is determined that does not result in an optimized numerical stability index exceeding a threshold. For example, if the numerical stability index is a condition number, the method determines an optimized condition number for each of a plurality of values ​​of N, and compares each determined optimized condition number to a threshold. If the determined optimized condition number for a given value of N is greater than the threshold, that value of N exceeds the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel. Any convenient threshold may be used. In some cases, the threshold is entered by a user. In some embodiments, the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel is N-1, and the condition number of N is the number that first exceeds the threshold for the i values ​​of N in an ascending series (e.g., N=2, N=3...N=i). In some embodiments, the method outputs the maximum number of fluorochromes that may be used simultaneously in a fluorochrome panel, and the fluorochrome panel identified with that panel number.

[0062] A fluorochrome panel of interest may include any suitable collection of fluorochromes. Fluorescent dyes of interest according to certain embodiments have excitation maxima in the range of 100 nm to 800 nm, such as 150 nm to 750 nm, for example 200 nm to 700 nm, for example 250 nm to 650 nm, for example 300 nm to 600 nm, for example 400 nm to 500 nm. Fluorescent dyes of interest according to certain embodiments have emission maxima in the range of 400 nm to 1000 nm, for example 450 nm to 950 nm, for example 500 nm to 900 nm, for example 550 nm to 850 nm, for example 600 nm to 800 nm. In some cases, the fluorescent dye is a luminescent dye such as a fluorescent dye having a peak emission wavelength of 200 nm or more, for example, 250 nm or more, for example, 300 nm or more, for example, 350 nm or more, for example, 400 nm or more, for example, 450 nm or more, for example, 500 nm or more, for example, 550 nm or more, for example, 600 nm or more, for example, 650 nm or more, for example, 700 nm or more, for example, 750 nm or more, for example, 800 nm or more, for example, 850 nm or more, for example, 900 nm or more, for example, 950 nm or more, for example, 1000 nm or more, for example, 1050 nm or more. For example, the fluorescent dye may be a fluorescent dye having a peak emission wavelength within a range of 200 nm to 1200 nm, for example, 300 nm to 1100 nm, for example, 400 nm to 1000 nm, for example, 500 nm to 900 nm, and may include a fluorescent dye having a peak emission wavelength of 600 nm to 800 nm.

[0063] Fluorescent dyes of interest include, but are not limited to, bodipy dyes, coumarin dyes, rhodamine dyes, acridine dyes, anthraquinone dyes, arylmethane dyes, diarylmethane dyes, chlorophyll-containing dyes, triarylmethane dyes, azo dyes, diazonium dyes, nitro dyes, nitroso dyes, phthalocyanine dyes, cyanine dyes, asymmetric cyanine dyes, quinone-imine dyes, azine dyes, eurhodin dyes, safranine dyes, indamines, indophenol dyes, fluorine dyes, oxazine dyes, oxazone dyes, thiazine dyes, thiazole dyes, xanthene dyes, fluorene dyes, pyronine dyes, fluorine dyes, rhodamine dyes, phenanthridine dyes, squaraines, bodipy, squaraine roxitanes, naphthalenes, coumarins, oxadiazoles, anthracenes, pyrenes, acridines, arylmethines, or tetrapyrroles, and combinations thereof. In some embodiments, the conjugate may include two or more dyes, such as two or more dyes selected from bodipy dyes, coumarin dyes, rhodamine dyes, acridine dyes, anthraquinone dyes, arylmethane dyes, diarylmethane dyes, chlorophyll-containing dyes, triarylmethane dyes, azo dyes, diazonium dyes, nitro dyes, nitroso dyes, phthalocyanine dyes, cyanine dyes, asymmetric cyanine dyes, quinone-imine dyes, azine dyes, eurydine dyes, safranine dyes, indamines, indophenol dyes, fluorine dyes, oxazine dyes, oxazone dyes, thiazine dyes, thiazole dyes, xanthene dyes, fluorene dyes, pyronine dyes, fluorine dyes, rhodamine dyes, phenanthridine dyes, squaraines, bodipy, squaraine rotaxanes, naphthalenes, coumarins, oxadiazoles, anthracenes, pyrenes, acridines, arylmethines, or tetrapyrroles, and combinations thereof.

[0064] In certain embodiments, fluorescent dyes of interest include, but are not limited to, fluorescein isothiocyanate (FITC), phycoerythrin (PE) dyes, peridinin chlorophyll protein-cyanine dyes (e.g., PerCP-Cy5.5), phycoerythrin-cyanine (PE-Cy) dyes (PE-Cy7), allophycocyanin (APC) dyes (e.g., APC-R700), allophycocyanin-cyanine dyes (e.g., APC-Cy7), and coumarin dyes (e.g., V450 or V500). In some cases, the fluorescent dye may include one or more of 1,4-bis-(o-methylstyryl)-benzene (bis-MSB 1,4-bis[2-(2-methylphenyl)ethenyl]-benzene), C510 dye, C6 dye, Nile Red dye, T614 dye (e.g., N-[7-(methanesulfonamido)-4-oxo-6-phenoxycyclomen-3-yl]formamide), LDS 821 dye (2-(6-(p-dimethylaminophenyl)-2,4-neopentylene-1,3,5-hexatrienyl)-3-ethylbenzothiazolium perchlorate), mFluor dye (e.g., mFluor Red dye such as mFluor 780NS).

[0065] Fluorescent dyes of interest include fluorescein, hydroxycoumarin, aminocoumarin, methoxycoumarin, Cascade Blue, Pacific Blue, Pacific Orange, Lucifer Yellow, NBD, R-Phycoerythrin (PE), PE-Cy5 conjugate, PE-Cy7 conjugate, Red 613, PerCP, TruRed, FluorX, BODIPY-FL, TRITC, X-Rhodamine, Lissamine Rhodamine B, Texas Red, allophycocyanin (APC), APC-Cy7 conjugate, Cy2, Cy3, Cy3B, Cy3.5, Cy5, Cy5.5, Cy7, Hoechst 33342, DAPI, Hoechst 33258, SYTOX Blue, Chromomycin A3, Mithramycin, YOYO-1, ethidium bromide, acridine orange, SYTOX Green, TOTO-1, TO-PRO-1, Thiazole Orange, Propidium Iodide (PI), LDS 751, 7-AAD, SYTOX Orange, TOTO-3, TO-PRO-3, DRAQ5, Indo-1, Fluo-3, DCFH, DHR, SNARF, Y66H, Y66F, EBFP, EBFP2, Azurite, GFPuv, T-Sapphire, TagBFP, Cerulean, mCFP, ECFP, CyPet, Y66W, dKeima-Red, mKeima-Red, TagCFP, AmCyan1, mTFP1 (Teal), S65A, Midoriishi-Cyan, Wild Type GFP, S65C, TurboGFP, TagGFP, TagGFP2, AcGFP1, S65L, Emerald, S65T, EGFP, Azami-Green, ZsGreen1, Dronpa-Green, TagYFP, EYFP, Topaz, Venus, mCitrine, YPet, TurboYFP, PhiYFP, PhiYFP-m, ZsYellow1, mBanana, Kusabira-Orange, mOrange, mOrange2, mKO, TurboRFP, tdTomato, DsRed-Express2, TagRFP, DsRed monomer, DsRed2("RFP"), mStrawberry, TurboFP602, AsRed2, mRFP1, J-Red, mCherry, HcRed1, mKate2, Katushka (TurboFP635), mKate (TagFP635), TurboFP635, mPlum, mRaspberry, mNeptune, E2-Crimson, Monochlorobimane, Calcein, Alexa Fluor 350, Alexa Fluor 405, Alexa Fluor 430, Alexa Fluor 488, Alexa Fluor 500, Alexa Fluor 514, Alexa Fluor 532, Alexa Fluor 546, Alexa Fluor 555, Alexa Fluor 568, Alexa Fluor 594, Alexa Fluor 610, Alexa Fluor 633, Alexa Fluor 647, Alexa Fluor 660, Alexa Fluor 680, Alexa Fluor 700, Alexa Fluor 750, Alexa Fluor 790, and Hyper.

[0066] In some cases, the fluorescent dye panel includes one or more polymer dyes (e.g., fluorescent polymer dyes). The fluorescent polymer dyes used in the subject methods and systems can vary. In some cases of the methods, the polymer dye includes a conjugated polymer. A conjugated polymer (CP) is characterized by a delocalized electronic structure that includes a backbone of alternating unsaturated bonds (e.g., double and / or triple bonds) and saturated bonds (e.g., single bonds), where π electrons can move from one bond to the other. Thus, the conjugated backbone may impart an extended linear structure to the polymer dye, with limited bond angles between the repeat units of the polymer. For example, proteins and nucleic acids are also macromolecules, but in some cases, they do not form extended rod structures, but rather fold into higher order three-dimensional shapes. In addition, CPs may form "rigid rod" shaped polymer backbones, exhibiting limited twist (e.g., torsion) angles between monomer repeat units along the polymer backbone chain. In some cases, the polymer dye includes a CP with a rigid rod-like structure. The structural properties of the polymer dye may affect the fluorescent properties of the molecule.

[0067] Any convenient polymer dye may be utilized in the subject devices and methods. In some cases, the polymer dye is a multi-chromophore with a structure capable of harvesting light to amplify the fluorescent output of the fluorophore. In some cases, the polymer dye can harvest light and efficiently convert it to emitted light of a longer wavelength. In some cases, the polymer dye has a light-harvesting multi-chromophore system that can efficiently transfer energy to a nearby light-emitting species (e.g., a "signaling chromophore"). Energy transfer mechanisms include, for example, resonance energy transfer (e.g., Förster (or fluorescence) resonance energy transfer (FRET)), quantum charge exchange (Dexter energy transfer), and the like. In some cases, these energy transfer mechanisms are relatively short-range, i.e., the proximity of the light-harvesting multi-chromophore system to the signaling chromophore allows for efficient energy transfer. Under conditions for efficient energy transfer, when the number of individual chromophores in the light-harvesting multi-luminophore system is large, amplification of the emission from the signaling chromophore occurs, i.e., the emission from the signaling chromophore is stronger when the incident light (excitation light) is at a wavelength that is absorbed by the light-harvesting multi-luminophore system than when the signaling chromophore is directly excited by the pump light.

[0068] The multiluminophores may be conjugated polymers. Conjugated polymers (CPs) are characterized by a delocalized electronic structure and can be used as highly responsive optical reporters for chemical and biological targets. Because the effective conjugation length is substantially shorter than the length of the polymer chain, the backbone contains many conjugated segments in close proximity. Thus, conjugated polymers are efficient at harvesting light and allow light amplification by Förster energy transfer.

[0069] Polymeric dyes of interest include those disclosed in U.S. Pat. Nos. 7,270,956, 7,629,448, 8,158,444, 8,227,187, 8,455,613, 8,575,303, 8,802,450, 8,969,509, 9,139,869, 9,371,559, 9,547,008, 10,094,838, 10,302,648, 10,458,989, 10,641,775, and 10,962,546, the entire disclosures of which are incorporated herein by reference, and Gaylord et al., J. Am. Chem. Soc., 2001, 123, the entire disclosures of which are incorporated herein by reference. (26), pp 6417-6418; Feng et al., Chem. Soc. Rev., 2010,39, 2411-2419; and Traina et al., J. Am. Chem. Soc., 2011, 133 (32), pp 12600-12607. Particular polymeric dyes that may be used include, but are not limited to, BD Horizon Brilliant™ Dye, such as BD Horizon Brilliant™ Violet Dye (e.g., BV421, BV510, BV605, BV650, BV711, BV786); BD Horizon Brilliant™ Ultraviolet Dye (e.g., BUV395, BUV496, BUV737, BUV805); and BD Horizon Brilliant™ Blue Dye (e.g., BB515) (BD Biosciences, San Jose, Calif.). Any fluorescent dye known to one of skill in the art may be used in the subject methods, including, but not limited to, those fluorescent dyes listed above or fluorescent dyes yet to be discovered.

[0070] The fluorochromes in the fluorochrome panel of interest and / or the fluorochromes referenced in the spectral matrix may or may not be bound to a biological molecule, such as a biopolymer. A biopolymer may be a biopolymer. A "biopolymer" is a polymer of one or more types of repeating units. Biopolymers are commonly found in living systems and specifically include polysaccharides (such as carbohydrates), peptides (the term is used to include polypeptides and proteins whether or not bound to polysaccharides), and polynucleotides, and analogs thereof, such as compounds composed of or containing amino acid analogs or non-amino acid groups, or nucleotide analogs or non-nucleotide groups. This includes polynucleotides in which the conventional backbone is replaced by a non-naturally occurring or synthetic backbone, and nucleic acids (or synthetic or naturally occurring analogs) in which one or more of the conventional bases are replaced by groups (natural or synthetic) that can participate in Watson-Crick type hydrogen bonding interactions. Polynucleotides include single-stranded or multi-stranded configurations, where one or more of the strands may or may not be perfectly aligned with another strand. In particular, "biopolymer" includes DNA (including cDNA), RNA and oligonucleotides, regardless of source. Thus, biomolecules may include polysaccharides, nucleic acids and polypeptides. For example, a nucleic acid may be an oligonucleotide, truncated or full-length DNA or RNA. In embodiments, the oligonucleotide, truncated and full-length DNA or RNA is comprised of 10 or more nucleotide monomers, such as 15 or more, such as 25 or more, such as 50 or more, such as 100 or more, such as 250 or more, such as 500 or more nucleotide monomers. For example, the length of the oligonucleotide, truncated and full-length DNA or RNA of interest may range from 10 nucleotides to 10 nucleotides. 8 Nucleotides, e.g. 10 2 Nucleotide ~10 7 Nucleotides, e.g. 10 3 Nucleotide ~10 6The length of the nucleic acid may be within the range of nucleotides. In embodiments, the biopolymer is not a single nucleotide or a short oligonucleotide (e.g., less than 10 nucleotides). By "full length" is meant that the DNA or RNA is a nucleic acid polymer having 70% or more of the complete sequence (as found in nature), e.g., 75% or more, such as 80% or more, such as 85% or more, such as 90% or more, such as 95% or more, such as 97% or more, such as 99% or more, e.g., 100% of the full length sequence of the DNA or RNA (as found in nature).

[0071] The polypeptide may optionally be a truncated or full length protein, enzyme or antibody. In embodiments, the polypeptide, truncated and full length protein, enzyme or antibody is composed of 10 or more amino acid monomers, such as 15 or more, such as 25 or more, such as 50 or more, such as 100 or more, such as 250 or more, such as 500 or more amino acid monomers. For example, the length of the polypeptide, truncated and full length protein, enzyme or antibody of interest may be between 10 amino acids and 10 8 Amino acids, e.g. 10 2 Amino acids ~10 7 Amino acids, e.g. 10 3 Amino acids ~10 6 The length of the biopolymer may be in the range of amino acids. In embodiments, the biopolymer is not a single amino acid or a short polypeptide (e.g., less than 10 amino acids). By "full length" is meant a polypeptide polymer having 70% or more of the complete sequence of the protein, enzyme or antibody (as found in nature), such as a polypeptide polymer having 75% or more, such as 80% or more, such as 85% or more, such as 90% or more, such as 95% or more, such as 97% or more, such as 99% or more, such as 100% of the full length sequence of the protein, enzyme or antibody (as found in nature).

[0072] In some cases, the fluorescent dye is conjugated to a specific binding member. The specific binding member and the fluorescent dye may be conjugated (e.g., covalently bonded) to one another at any convenient location on the two molecules via any linker. As used herein, the term "specific binding member" refers to one member of a pair of molecules that have binding specificity for one another. One member of the pair of molecules may have a surface or cavity region that specifically binds to a surface region or cavity of the other member of the pair of molecules. Thus, the members of the pair have the property of specifically binding to one another to generate a binding complex. In some embodiments, the affinity between the specific binding members in the binding complex is greater than 10 -6 M or less, e.g. 10 -7 M or less, e.g. 10 -8 M or less, e.g. 10 -9 M or less, 10 -10 M or less, 10 -11 M or less, 10 -12 M or less, 10 -13 M or less, 10 -14 M or less, e.g. 10 -15 K below M d In some embodiments, a specific binding member specifically binds with high avidity. High avidity means that the binding member specifically binds with a dissociation constant greater than 10×10 -9 M or less, e.g. 1×10 -9 M or less, 3×10 -10 M or less, 1×10 -10 M or less, 3×10 -11 M or less, 1×10 -11 M or less, 3×10 -12 M or less, or 1 x 10 -12 Apparent K below M d This means that the antibody specifically binds with an apparent affinity characterized by

[0073] The specific binding member may be proteinaceous. As used herein, the term "proteinaceous" refers to a moiety consisting of amino acid residues. The proteinaceous moiety may be a polypeptide. In some cases, the proteinaceous specific binding member is an antibody. In some embodiments, the proteinaceous specific binding member is an antibody fragment, e.g., a binding fragment of an antibody that specifically binds to a polymeric dye. As used herein, the terms "antibody" and "antibody molecule" are used interchangeably and refer to a protein consisting of one or more polypeptides substantially encoded by all or part of a recognized immunoglobulin gene. For example, the recognized human immunoglobulin genes include the kappa (k), lambda (l) and heavy chain gene loci, which include a myriad of variable region genes, as well as the constant region genes μ(u), δ(d), γ(g), σ(e) and α(a), which encode the IgM, IgD, IgG, IgE and IgA isotypes, respectively. The variable region of an immunoglobulin light or heavy chain is composed of a "framework" region (FR) interrupted by three hypervariable regions, also called "complementarity determining regions" or "CDRs". The extent of the framework regions and CDRs has been precisely defined (see "Sequences of Proteins of Immunological Interest," E. Kabat et al., US Department of Health and Human Services, (1991)). The numbering of all antibody amino acid sequences described herein is according to the Kabat system. The sequences of the framework regions of different light or heavy chains are relatively conserved within a species. The framework region of an antibody, i.e. the combined framework regions of the constituent light and heavy chains, serves to position and align the CDRs. The CDRs are primarily responsible for binding to an epitope of an antigen. The term "antibody" is meant to include full length antibodies and may refer to natural antibodies from any organism, artificial antibodies, or antibodies that have been recombinantly produced for experimental, therapeutic or other purposes as further defined below.Antibody fragments of interest include, but are not limited to, Fab, Fab', F(ab')2, Fv, ScFv, or other antigen-binding subsequences of antibodies, either produced by modification of whole antibodies or synthesized de novo using recombinant DNA technology. The antibodies may be monoclonal or polyclonal and may have other specific activities on cells (e.g., antagonist, agonist, neutralizing, inhibitory or stimulatory). It is understood that the antibodies may have additional conservative amino acid substitutions that do not substantially affect antigen binding function or other antibody functions. In some embodiments, the specific binding member is a Fab fragment, a F(ab')2 fragment, a ScFv, a bispecific antibody, or a trispecific antibody. In some embodiments, the specific binding member is an antibody. Optionally, the specific binding member is a murine antibody or binding fragment thereof. Optionally, the specific binding member is a recombinant antibody or binding fragment thereof.

[0074] In an embodiment, a sample is analyzed using a panel of fluorescent dyes of interest. In some cases, the sample analyzed is a biological sample. The term "biological sample" is used in the conventional sense to refer to a whole organism, a plant, a fungus, or a subset of animal tissues, cells, or components, as may be found in blood, mucus, lymph, synovial fluid, cerebrospinal fluid, saliva, bronchoalveolar lavage fluid, amniotic fluid, amniotic cord blood, urine, vaginal fluid, and semen, or the like. Thus, a "biological sample" refers to both a naturally occurring organism or a subset of its tissues, and homogenates, lysates, or extracts prepared from an organism or a subset of its tissues, including, but not limited to, plasma, serum, spinal fluid, lymph, skin sections, respiratory tract, gastrointestinal tract, cardiovascular, genitourinary tract, tears, saliva, milk, blood cells, tumors, organs. A biological sample may be any type of organism's tissue, including both healthy and diseased tissue (e.g., cancerous, malignant, necrotic, etc.). In certain embodiments, the biological sample is a liquid sample, such as blood or a derivative thereof, such as plasma, tears, urine, semen, etc., and in some cases, the sample is a blood sample, including whole blood, such as blood obtained from a venipuncture or finger prick (which may or may not be combined with any reagents, such as preservatives, anticoagulants, etc., prior to assay).

[0075] In some embodiments, the sample source is a "mammal" or "mammalian," which terms are used broadly to describe organisms belonging to the class Mammalia, including the orders Carnivora (e.g., dogs and cats), Rodentia (e.g., mice, guinea pigs, and rats), and Primates (e.g., humans, chimpanzees, and monkeys). In some cases, the subject is a human. The method may be applied to samples obtained from human subjects of both genders and at any stage of development (i.e., neonate, infant, juvenile, adolescent, adult), and in some embodiments, the human subject is a juvenile, adolescent, or adult. While the present invention may be applied to samples from human subjects, it should be understood that the method may also be performed on samples from other animal subjects (i.e., "non-human subjects"), including, but not limited to, birds, mice, rats, dogs, cats, livestock, and horses.

[0076] The fluorochromes in the fluorochrome panel may be configured to target different types of cells (e.g., via antibodies targeted to the cells). A variety of cells may be characterized using the subject methods. Target cells of interest include, but are not limited to, stem cells, T cells, dendritic cells, B cells, granulocytes, leukemia cells, lymphoma cells, viral cells (e.g., HIV cells), NK cells, macrophages, monocytes, fibroblasts, epithelial cells, endothelial cells, and red blood cells. Target cells of interest include cells bearing convenient cell surface markers or antigens that may be internalized or labeled by convenient affinity agents or conjugates thereof. For example, the target cell may comprise cell surface antigens such as CD11b, CD123, CD14, CD15, CD16, CD19, CD193, CD2, CD25, CD27, CD3, CD335, CD36, CD4, CD43, CD45RO, CD56, CD61, CD7, CD8, CD34, CD1c, CD23, CD304, CD235a, T cell receptor alpha / beta, T cell receptor gamma / delta, CD253, CD95, CD20, CD105, CD117, CD120b, Notch4, Lgr5 (N-terminus), SSEA-3, TRA-1-60 antigen, disialoganglioside GD2, and CD71. In some embodiments, the target cells are selected from HIV-containing cells, Treg cells, antigen-specific T cell populations, tumor cells, or hematopoietic progenitor cells (CD34+) from whole blood, bone marrow, or umbilical cord blood.

[0077] In certain embodiments, the fluorochrome panels identified by the methods may be used in a flow cytometry protocol (e.g., to analyze a sample as described above). In carrying out such methods, a sample (e.g., in a flow stream of a flow cytometer) is illuminated with light from a light source. In some embodiments, the light source is a broadband light source, emitting light having a wide range of wavelengths, e.g., 50 nm or more, e.g., 100 nm or more, e.g., 150 nm or more, e.g., 200 nm or more, e.g., 250 nm or more, e.g., 300 nm or more, e.g., 350 nm or more, e.g., 400 nm or more, e.g., 500 nm or more. For example, one suitable broadband light source emits light having a wavelength between 200 nm and 1500 nm. Another example of a suitable broadband light source includes a light source that emits light having a wavelength between 400 nm and 1000 nm. When the method involves illumination with a broadband light source, broadband light source protocols of interest may include, but are not limited to, halogen lamps, deuterium arc lamps, xenon arc lamps, stabilized fiber-coupled broadband light sources, broadband LEDs with continuous spectra, superluminescent light emitting diodes, semiconductor light emitting diodes, broad spectrum LED white light sources, multi-LED integrated white light sources, or any combination thereof, among other broadband light sources.

[0078] In other embodiments, the methods of the present invention are illuminated using a narrowband light source emitting at a specific wavelength or narrow range of wavelengths, e.g., irradiated using a light source emitting light in a narrow range of wavelengths, such as 50 nm or less, e.g., 40 nm or less, e.g., 30 nm or less, e.g., 25 nm or less, e.g., 20 nm or less, e.g., 15 nm or less, e.g., 10 nm or less, e.g., 5 nm or less, e.g., 2 nm or less, e.g., irradiated using a light source emitting light of a specific wavelength (i.e., monochromatic light). When the methods are illuminated using a narrowband light source, narrowband light source protocols of interest include, but are not limited to, narrow wavelength LEDs, laser diodes, or broadband light sources coupled to one or more optical bandpass filters, diffraction gratings, monochromators, or any combination thereof.

[0079] In some embodiments, the method uses one or more lasers to irradiate the sample. As discussed above, the type and number of lasers will depend on the sample and the desired light to be collected, and may be a gas laser, such as a helium neon laser, an argon laser, a krypton laser, a xenon laser, a nitrogen laser, a CO2 laser, a CO2 laser, an argon fluorine (ArF) excimer laser, a krypton fluorine (KrF) excimer laser, a xenon chlorine (XeCl) excimer laser, or a xenon fluorine (XeF) excimer laser, or a combination thereof. In other cases, the method uses a dye laser, such as a stilbene laser, a coumarin laser, or a rhodamine laser, to irradiate the flow stream. In still other cases, the method irradiates the flowstream using a metal vapor laser, such as a helium-cadmium (HeCd), helium-mercury (HeHg), helium-selenium (HeSe), helium-silver (HeAg), strontium, neon-copper (NeCu), copper, or gold laser, and combinations thereof. In still other cases, the method irradiates the flowstream using a solid state laser, such as a ruby, Nd:YAG, NdCrYAG, Er:YAG, Nd:YLF, Nd:YVO4, Nd:YCa4O(BO3)3, Nd:YCOB, titanium sapphire, thulium YAG, ytterbium YAG, Yb2O3, or cerium doped laser, and combinations thereof.

[0080] The sample may be illuminated using one or more of the light sources described above, such as 2 or more, such as 3 or more, such as 4 or more, such as 5 or more, such as 10 or more. The light source may include any combination of light source types. For example, in some embodiments, the method illuminates the sample in the flow stream using an array of lasers, such as an array having one or more gas lasers, one or more dye lasers, and one or more solid state lasers.

[0081] The sample may be irradiated with a wavelength in the range of 200 nm to 1500 nm, for example, 250 nm to 1250 nm, for example, 300 nm to 1000 nm, for example, 350 nm to 900 nm, for example, 400 nm to 800 nm. For example, if the light source is a broadband light source, the sample may be irradiated with a wavelength in the range of 200 nm to 900 nm. In other cases, if the light source has multiple narrowband light sources, the sample may be irradiated with a specific wavelength in the range of 200 nm to 900 nm. For example, the light source may be multiple narrowband LEDs (1 nm to 25 nm) each independently emitting light having a wavelength in the range of 200 nm to 900 nm. In other embodiments, the narrowband light source includes one or more lasers (for example, a laser array), and the sample is irradiated with a specific wavelength in the range of 200 nm to 700 nm using a laser array having, for example, gas lasers, excimer lasers, dye lasers, metal vapor lasers, and solid-state lasers as described above.

[0082] When two or more light sources are used, the light sources may be used to irradiate the sample simultaneously or sequentially or a combination thereof. For example, each of the light sources may be used to irradiate the sample simultaneously. In other embodiments, each of the light sources is used to irradiate the flow stream sequentially. When two or more light sources are used to irradiate the sample sequentially, the time for which each light source irradiates the sample may independently be 0.001 microseconds or more, such as 0.01 microseconds or more, such as 0.1 microseconds or more, such as 1 microseconds or more, such as 5 microseconds or more, such as 10 microseconds or more, such as 30 microseconds or more, such as 60 microseconds or more. For example, the method may irradiate the sample using a light source (e.g., a laser) for a duration in the range of 0.001 microseconds to 100 microseconds, such as 0.01 microseconds to 75 microseconds, such as 0.1 microseconds to 50 microseconds, such as 1 microsecond to 25 microseconds, such as 5 microseconds to 10 microseconds. In embodiments in which two or more light sources are used to sequentially illuminate the sample, the duration for which the sample is illuminated by each light source may be the same or different.

[0083] The time between illumination by each light source may further vary as desired and may be independently separated by a delay of 0.001 microseconds or more, such as 0.01 microseconds or more, such as 0.1 microseconds or more, such as 1 microseconds or more, such as 5 microseconds or more, such as 10 microseconds or more, such as 15 microseconds or more, such as 30 microseconds or more, such as 60 microseconds or more. For example, the time between illumination by each light source may be in the range of 0.001 microseconds to 60 microseconds, such as 0.01 microseconds to 50 microseconds, such as 0.1 microseconds to 35 microseconds, such as 1 microseconds to 25 microseconds, such as 5 microseconds to 10 microseconds. In an embodiment, the time between illumination by each light source is 10 microseconds. In an embodiment in which the sample is illuminated sequentially by more than two (i.e., three or more) light sources, the delay between illumination by each light source may be the same or different.

[0084] The sample may be illuminated continuously or at discrete intervals. In some cases, the method uses a light source to illuminate the sample continuously. In other cases, the light source illuminates the sample at discrete intervals, such as every 0.001 milliseconds, 0.01 milliseconds, 0.1 milliseconds, 1 millisecond, 10 milliseconds, 100 milliseconds, 1000 milliseconds, or other intervals.

[0085] Depending on the light source, the sample may be illuminated from various distances, for example 0.01 mm or more, such as 0.05 mm or more, for example 0.1 mm or more, such as 0.5 mm or more, for example 1 mm or more, such as 2.5 mm or more, for example 5 mm or more, such as 10 mm or more, for example 15 mm or more, such as 25 mm or more, for example 50 mm or more. Furthermore, the angle or illumination may vary within the range of 10°-90°, for example 15°-85°, for example 20°-80°, for example 25°-75°, for example 30°-60°, for example an angle of 90°.

[0086] In embodiments, light from the illuminated sample is conveyed to a light detection system and measured by one or more photodetectors. In practicing the subject methods, light from the sample is conveyed to three or more wavelength separators, each configured to pass light having a predetermined spectral region. The spectral regions of light from each of the wavelength separators are conveyed to one or more photodetection modules having optics configured to convey light having the predetermined sub-spectral regions to a photodetector.

[0087] The light may be measured continuously or at discrete intervals using a light detection system. In some cases, the method measures the light continuously. In other cases, the light is measured at discrete intervals, such as every 0.001 milliseconds, every 0.01 milliseconds, every 0.1 milliseconds, every 1 millisecond, every 10 milliseconds, every 100 milliseconds, every 1000 milliseconds, or other intervals.

[0088] Measurements of collected light may be made one or more times during the subject methods, such as two or more times, such as three or more times, such as five or more times, such as ten or more times, in some embodiments, light propagation is measured two or more times, and the data is optionally averaged.

[0089] In some embodiments, the method conditions the light before detecting it with the optical detection system of interest. For example, light from the sample source may pass through one or more lenses, mirrors, pinholes, slits, gratings, optical refractors, and any combination thereof. In some cases, the collected light passes through one or more focusing lenses to reduce the profile of the light that is directed to the optical detection system or optical collection system, for example as described above. In other cases, the light emitted from the sample passes through one or more collimators to reduce the divergence of the light beam transmitted to the optical detection system.

[0090] system Aspects of the invention further include a system configured to perform the above method. The system of interest includes a processor configured to receive an instrument identifier and a request for N fluorochrome identifiers, where N is an integer of at least two. The processor of interest is further configured to select two or more subsets of the N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier, and identify a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. In some embodiments, the processor is configured to identify the fluorochrome panel without evaluating antigen data. In some cases, the processor is further configured to evaluate a maximum number of fluorochromes that may be used simultaneously in the fluorochrome panel. The system may optionally be a particle analyzer (e.g., a flow cytometer). In some embodiments, the system is a full spectrum flow cytometer.

[0091] In an embodiment, the subject processor operates with programmable logic, which may be implemented in hardware, software, firmware, or any combination thereof, to identify the fluorochrome panel. For example, when the programmable logic is implemented in software, the identification of the fluorochrome panel may be at least partially achieved by a computer-readable data storage medium having program code including instructions configured, when executed, to receive an instrument identifier and a request for N fluorochrome identifiers, select two or more subsets of the N fluorochrome identifiers from the set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier, and identify the fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using the spectral features associated with the N fluorochrome identifiers of the two or more subsets.

[0092] The programmable logic may use any convenient mechanism for calculating the numerical stability index. In some cases, the numerical stability index is a condition number. In such cases, the condition number may be calculated by determining singular values ​​for two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets, and determining a ratio of the maximum singular value to the minimum singular value. In other embodiments, the programmable logic is configured to calculate the condition number by determining a Frobenius norm for each of the two or more subsets of N fluorochrome identifiers, determining a Frobenius norm of a pseudoinverse for each of the two or more subsets of N fluorochrome identifiers, and determining, for each of the two or more subsets of N fluorochrome identifiers, a product of the Frobenius norm of the subset of N fluorochrome identifiers and the Frobenius norm of the pseudoinverse for the subset of N fluorochrome identifiers. In some cases, the programmable logic is configured to identify the fluorochrome panel by a constrained optimization method (e.g., a minimization algorithm). In some cases, the constrained optimization method is performed iteratively (e.g., by generating a first subset, generating a second subset by modifying the first subset, and comparing the numerical stability index of both the first and second subsets). In some cases, the constrained optimization method further includes iteratively generating subsets of fluorochrome identifiers, calculating a numerical stability index (e.g., a condition number) for each iteratively generated subset of N fluorochrome identifiers, and comparing the calculated numerical stability index for each iteratively generated subset of N fluorochrome identifiers with the numerical stability index of the already generated subset of N fluorochrome identifiers to determine an optimized numerical stability index.

[0093] In a further embodiment, the programmable logic is configured to evaluate a maximum number of fluorochromes that may be simultaneously used in the fluorochrome panel, in such an embodiment, the multiple values ​​of N include i values, where i is, for example, the number of fluorochrome identifiers in the spectral matrix, and determine a maximum value for N that does not cause the optimized numerical stability index to exceed a threshold value when evaluating the maximum number of fluorochromes that may be simultaneously used in the fluorochrome panel.

[0094] The programmable logic of interest may be implemented in any of a variety of devices, such as a specially programmed event processing computer, a wireless communication device, an integrated circuit device, etc. In some embodiments, the programmable logic may be executed by a specially programmed processor, which may include one or more processors, such as one or more digital signal processors (DSPs), configurable microprocessors, application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), or other equivalent integrated or discrete logic circuits. Combinations of computing devices, such as a combination of a DSP and a microprocessor, multiple microprocessors, one or more microprocessors in combination with a DSP core, or any other such configuration in at least partial data connectivity, may implement one or more of the features described herein.

[0095] In some cases, the system is a particle analyzer. Particle analyzers of interest may include a flow cell for carrying particles in a flow stream, a light source for illuminating particles in the flow stream at an interrogation point, and a particle modulated light detector for detecting the particle modulated light. In some embodiments, the particle analyzer is a flow cytometer. When the particle analyzer is a flow cytometer, the flow cytometer is a full spectrum flow cytometer.

[0096] As used herein, "flow cell" is used in the conventional sense to refer to a component, such as a cuvette, that includes a flow path having a liquid flow stream for carrying particles within a sheath fluid. A cuvette of interest includes a container with a passageway therethrough. The flow stream may include a liquid sample injected from a sample tube. A flow cell of interest includes a flow path accessible to light. In some cases, the flow cell includes a transparent material (e.g., quartz) that allows the passage of light. In some embodiments, the flow cell is a stream-in-air flow cell in which optical interrogation of particles occurs outside of the flow cell (i.e., in free space).

[0097] In some cases, the flow stream is configured to illuminate the light from the light source at the investigation point. The flow stream in which the flow path is configured may include a liquid sample injected from a sample tube. In some embodiments, the flow stream may include a narrow, rapidly flowing liquid stream, such stream arranged such that linearly separated particles carried within such stream are in line and separated from each other. As used herein, an "investigation point" refers to an area in the flow cell where particles are illuminated by light from a light source, for example for analysis. The size of the investigation point may vary as desired. For example, if 0 μm represents the axis of the light emitted by the light source, the investigation point may be in the range of -100 μm to 100 μm, such as -50 μm to 50 μm, such as -25 μm to 40 μm, such as -15 μm to 30 μm.

[0098] Particle modulated light may be observed after illuminating the particles in the flow cell. "Particle modulated light" refers to light received from particles in a flow stream after illuminating the particles with light from a light source. In some cases, the particle modulated light is side scattered light. As described herein, side scattered light refers to light that is refracted and reflected off the surface and internal structures of a particle. In further embodiments, the particle modulated light includes forward scattered light (i.e., light that travels through or around the particle primarily in a forward direction). In still other cases, the particle modulated light includes fluorescent light (i.e., light emitted from a fluorescent dye after illumination with excitation wavelength light).

[0099] As mentioned above, aspects of the invention further include a light source configured to illuminate particles passing through the flow cell at the interrogation point. Any convenient light source may be used as the light source described herein. In some embodiments, the light source is a laser. In embodiments, the laser may be any convenient laser, such as a continuous wave laser. For example, the laser may be a diode laser, such as an ultraviolet diode laser, a visible diode laser, and a near infrared diode laser. In other embodiments, the laser may be a helium neon (HeNe) laser. In some cases, the laser is a gas laser, such as a helium neon laser, an argon laser, a krypton laser, a xenon laser, a nitrogen laser, a CO2 laser, a CO2 laser, an argon fluorine (ArF) excimer laser, a krypton fluorine (KrF) excimer laser, a xenon chlorine (XeCl) excimer laser, or a xenon fluorine (XeF) excimer laser, or a combination thereof. In other cases, the flow cytometer of interest has a dye laser, such as a stilbene laser, a coumarin laser, or a rhodamine laser. In still other cases, the lasers of interest include metal vapor lasers, such as helium-cadmium (HeCd), helium-mercury (HeHg), helium-selenium (HeSe), helium-silver (HeAg), strontium, neon-copper (NeCu), copper, or gold lasers, and combinations thereof. In still other cases, flow cytometers of interest include solid-state lasers, such as ruby, Nd:YAG, NdCrYAG, Er:YAG, Nd:YLF, Nd:YVO4, Nd:YCa4O(BO3)3, Nd:YCOB, titanium sapphire, thulium YAG, ytterbium YAG, Yb2O3, or cerium-doped lasers, and combinations thereof.

[0100] The laser light source according to an embodiment may further include one or more optical adjustment components. In an embodiment, the optical adjustment components are disposed between the light source and the flow cell and may include any device capable of changing the spatial width of the illumination or other characteristics of the illumination from the light source, such as the illumination direction, wavelength, beam width, beam intensity, and focus. The optical adjustment protocol may include any convenient device that adjusts one or more characteristics of the light source, including, but not limited to, lenses, mirrors, filters, optical fibers, wavelength separators, pinholes, slits, collimation protocols, and combinations thereof. In an embodiment, the flow cytometer of interest includes one or more focusing lenses. The focusing lens may be, in one example, a reduction lens. In yet another embodiment, the flow cytometer of interest includes optical fibers.

[0101] If the optical adjustment component is configured to move, it may be configured to move continuously or at discrete intervals, for example in increments of 0.01 μm or more, such as 0.05 μm or more, for example 0.1 μm or more, such as 0.5 μm or more, for example 1 μm or more, such as 10 μm or more, for example 100 μm or more, such as 500 μm or more, for example 1 mm or more, such as 5 mm or more, for example 10 mm or more, such as 25 mm or more.

[0102] Any movement protocol may be used to move the optical adjustment component structure, for example coupled to a movable support stage or directly coupled to a motor-driven translation stage, a lead screw translation assembly, a geared translation device, for example using stepper motors, servo motors, brushless electric motors, brushed DC motors, microstep drive motors, and high resolution stepper motors, among other types of motors.

[0103] The light source may be positioned at any suitable distance from the flow cell, for example the light source and the flow cell are separated by 0.005 mm or more, such as 0.01 mm or more, for example 0.05 mm or more, such as 0.1 mm or more, for example 0.5 mm or more, such as 1 mm or more, for example 5 mm or more, such as 10 mm or more, for example 25 mm or more, for example 100 mm or more. In addition, the light source may be positioned at any suitable angle relative to the flow cell, for example at an angle between 10 degrees and 90 degrees, for example 15 degrees and 85 degrees, for example 20 degrees and 80 degrees, for example 25 degrees and 75 degrees, for example 30 degrees and 60 degrees, for example at an angle of 90 degrees.

[0104] In some embodiments, the light source of interest may include multiple lasers, e.g., 2 or more lasers, e.g., 3 or more lasers, e.g., 4 or more lasers, e.g., 5 or more lasers, e.g., 10 or more lasers, e.g., 15 or more lasers, configured to provide laser light for separate illumination of the flow stream. Depending on the desired wavelength of light for illumination of the flow stream, each laser may have a different specific wavelength within the range of 200 nm to 1500 nm, e.g., 250 nm to 1250 nm, e.g., 300 nm to 1000 nm, e.g., 350 nm to 900 nm, e.g., 400 nm to 800 nm. In some embodiments, the lasers of interest may include one or more of a 405 nm laser, a 488 nm laser, a 561 nm laser, and a 635 nm laser.

[0105] As discussed above, particle analyzers of interest may further include one or more particle modulation light detectors for detecting particle modulation light intensity data. In some embodiments, the one or more particle modulation light detectors include one or more forward scattered light detectors configured to detect forward scattered light. For example, particle analyzers of interest may include one forward scattered light detector or multiple forward scattered light detectors, such as two or more, such as three or more, such as four or more, such as five or more forward scattered light detectors. In some embodiments, the particle analyzer includes one forward scattered light detector. In other embodiments, the particle analyzer includes two forward scattered light detectors.

[0106] Any convenient detector for detecting collected light may be used in the forward scattered light detectors described herein. Detectors of interest may include, but are not limited to, light sensors or photodetectors, such as active pixel sensors (APS), avalanche photodiodes, image sensors, charge-coupled devices (CCD), intensified charge-coupled devices (ICCD), light-emitting diodes, photon counters, bolometers, pyroelectric detectors, photoresistors, photocells, photodiodes, photomultiplier tubes (PMT), phototransistors, quantum dot photoconductors or photodiodes, and combinations thereof, among other detectors. In some embodiments, the collected light is measured with a charge-coupled device (CCD), a semiconductor charge-coupled device (CCD), an active pixel sensor (APS), a complementary metal-oxide semiconductor (CMOS) image sensor, or an N-type metal-oxide semiconductor (NMOS) image sensor. In some embodiments, the detector is a photomultiplier tube, e.g., a 0.01 cm 2 ~10cm 2 , e.g. 0.05 cm 2 ~9 cm 2 , e.g. 0.1 cm 2 ~8 cm 2 , e.g. 0.5 cm 2 ~7 cm 2 , e.g. 1 cm 2 ~5 cm 2The photomultiplier tube has an active detection surface area of ​​each region within the range of .

[0107] In embodiments, the forward scattered light detector is configured to measure light continuously or at discrete intervals. In some cases, the detector of interest is configured to measure the collected light continuously. In other cases, the detector of interest is configured to measure at discrete intervals, such as every 0.001 milliseconds, every 0.01 milliseconds, every 0.1 milliseconds, every 1 millisecond, every 10 milliseconds, every 100 milliseconds, such as every 1000 milliseconds, or other intervals.

[0108] In further embodiments, the one or more particle modulation light detectors may comprise one or more side scatter light detectors for detecting side scattered wavelengths of light (i.e., light that is refracted and reflected off the surface and internal structure of the particle). In some embodiments, the particle analyzer comprises one side scatter light detector. In other embodiments, the particle analyzer comprises multiple side scatter light detectors, for example 2 or more, for example 3 or more, for example 4 or more, for example 5 or more side scatter light detectors.

[0109] Any convenient detector for detecting collected light may be used for the side scatter light detectors described herein. Detectors of interest may include, but are not limited to, light sensors or photodetectors, such as active pixel sensors (APS), avalanche photodiodes, image sensors, charge-coupled devices (CCD), intensified charge-coupled devices (ICCD), light-emitting diodes, photon counters, bolometers, pyroelectric detectors, photoresistors, photocells, photodiodes, photomultiplier tubes (PMT), phototransistors, quantum dot photoconductors or quantum dot photodiodes, and combinations thereof, among other detectors. In some embodiments, the collected light is measured with a charge-coupled device (CCD), a semiconductor charge-coupled device (CCD), an active pixel sensor (APS), a complementary metal-oxide semiconductor (CMOS) image sensor, or an N-type metal-oxide semiconductor (NMOS) image sensor. In some embodiments, the detector is a photomultiplier tube, e.g., a 0.01 cm 2 ~10 cm 2 , e.g. 0.05m 2 ~9 cm 2 , e.g. 0.1 cm 2 ~8 cm 2 , e.g. 0.5 cm 2 ~7 cm 2 , e.g. 1 cm 2 ~5 cm 2 The photomultiplier tube has an active detection surface area of ​​each region within the range of .

[0110] In embodiments, a subject particle analyzer further comprises a fluorescence detector configured to detect light at one or more fluorescence wavelengths, hi other embodiments, the particle analyzer comprises a plurality of fluorescence detectors, e.g., 2 or more, e.g., 3 or more, e.g., 4 or more, 5 or more, 10 or more, 15 or more, e.g., 20 or more fluorescence detectors.

[0111] Any convenient detector for detecting collected light may be used in the fluorescence photodetectors described herein. Detectors of interest may include, but are not limited to, light sensors or photodetectors, such as active pixel sensors (APS), avalanche photodiodes, image sensors, charge-coupled devices (CCD), intensified charge-coupled devices (ICCD), light-emitting diodes, photon counters, bolometers, pyroelectric detectors, photoresistors, photocells, photodiodes, photomultiplier tubes (PMT), phototransistors, quantum dot photoconductors or quantum dot photodiodes, and combinations thereof, among other detectors. In some embodiments, the collected light is measured with a charge-coupled device (CCD), a semiconductor charge-coupled device (CCD), an active pixel sensor (APS), a complementary metal-oxide semiconductor (CMOS) image sensor, or an N-type metal-oxide semiconductor (NMOS) image sensor. In some embodiments, the detector is a photomultiplier tube, e.g., a 0.01 cm 2 ~10 cm 2 , e.g. 0.05 cm 2 ~9 cm 2 , e.g. 0.1 cm 2 ~8 cm 2 , e.g. 0.5 cm 2 ~7 cm 2 , e.g. 1 cm 2 ~5 cm 2 The photomultiplier tube has an active detection surface area of ​​each region within the range of .

[0112] When a particle analyzer of interest has multiple fluorescence detectors, each of the fluorescence detectors may be identical, or the collection of fluorescence detectors may be a combination of different types of detectors. For example, when a particle analyzer of interest has two fluorescence detectors, in some embodiments, the first fluorescence detector is a CCD type device and the second fluorescence detector (or image sensor) is a CMOS type device. In other embodiments, both the first and second fluorescence detectors are CCD type devices. In yet other embodiments, both the first and second fluorescence detectors are CMOS type devices. In yet other embodiments, the first fluorescence detector is a CCD type device and the second fluorescence detector is a photomultiplier tube (PMT). In yet other embodiments, the first fluorescence detector is a CMOS type device and the second fluorescence detector is a photomultiplier tube. In yet other embodiments, both the first and second fluorescence detectors are photomultiplier tubes.

[0113] In embodiments of the present disclosure, the fluorescent light detector of interest is configured to measure collected light at one or more wavelengths, e.g., two or more wavelengths, e.g., five or more different wavelengths, e.g., ten or more different wavelengths, e.g., twenty-five or more different wavelengths, e.g., fifty or more different wavelengths, e.g., one or more than one hundred different wavelengths, e.g., two or more than one hundred different wavelengths, e.g., two or more than one hundred different wavelengths, e.g., three or more than one hundred different wavelengths, e.g., four ...

[0114] In some embodiments, the fluorescent light detector of interest is configured to measure light collected over a range of wavelengths (e.g., 200 nm to 1000 nm). In certain embodiments, the detector of interest is configured to collect a spectrum of light over a range of wavelengths. For example, a particle analyzer may have one or more detectors configured to collect a spectrum of light over one or more of the wavelength ranges of 200 nm to 1000 nm. In yet other embodiments, the detector of interest is configured to measure light emitted from a sample in the flow stream at one or more particular wavelengths. For example, the particle analyzer may have one or more detectors configured to measure light at one or more of: 450 nm, 518 nm, 519 nm, 561 nm, 578 nm, 605 nm, 607 nm, 625 nm, 650 nm, 660 nm, 667 nm, 670 nm, 668 nm, 695 nm, 710 nm, 723 nm, 780 nm, 785 nm, 647 nm, 617 nm, and any combination thereof. In some embodiments, the one or more detectors may be configured to pair with a particular fluorophore, such as a fluorophore used with a sample in a fluorescence assay.

[0115] In some embodiments, the particle analyzer comprises one or more wavelength separators disposed between the flow cell and the one or more particle modulation light detectors. The term "wavelength separator" is used herein in its conventional sense to refer to an optical component configured to separate light collected from a sample into predetermined spectral regions. In some embodiments, the particle analyzer comprises one wavelength separator. In other embodiments, the particle analyzer comprises multiple wavelength separators, e.g., 2 or more wavelength separators, e.g., 3 or more, e.g., 4 or more, e.g., 5 or more, e.g., 6 or more, e.g., 7 or more, e.g., 8 or more, e.g., 9 or more, e.g., 10 or more, e.g., 15 or more, e.g., 25 or more, e.g., 50 or more, e.g., 75 or more, e.g., 100 or more wavelength separators. In some embodiments, the wavelength separator is configured to separate light collected from a sample into predetermined spectral regions by passing light having a predetermined spectral region and reflecting light in one or more remaining spectral regions. In other embodiments, the wavelength separator is configured to separate light collected from a sample into predetermined spectral regions by passing light having a predetermined spectral region and absorbing light in one or more remaining spectral regions. In yet other embodiments, the wavelength separator is configured to spatially diffract light collected from the sample into a predetermined spectral region. Each wavelength separator may be any convenient light separation protocol, such as one or more dichroic mirrors, bandpass filters, diffraction gratings, beam splitters, or prisms. In some embodiments, the wavelength separator is a prism. In other embodiments, the wavelength separator is a diffraction grating. In certain embodiments, the wavelength separator of the subject light detection system is a dichroic mirror.

[0116] Suitable flow cytometry systems include those described in Ormerod (ed.), Flow Cytometry: A Practical Approach, Oxford Univ. Press (1997);Jaroszeski et al. (eds.), Flow Cytometry Protocols, Methods in Molecular Biology No. 91, Humana Press (1997);Practical Flow Cytometry, 3rd ed., Wiley-Liss (1995);Virgo, et al. (2012) Ann Clin Biochem. Jan;49(pt 1):17-28;Linden, et. al., Semin Throm Hemost. 2004 Oct;30(5):502-11;Alison, et al. J Pathol, 2010 Dec; 222(4):335-344; and Herbig, et al. (2007) Crit Rev Ther Drug Carrier Syst. 24(3):203-255.In some cases, flow cytometry systems of interest include BD Biosciences FACSCanto™ flow cytometers, BD Biosciences FACSCanto™ II flow cytometers, BD Accuri™ flow cytometers, BD Accuri™ C6 Plus flow cytometers, BD Biosciences FACSCelesta™ flow cytometers, BD Biosciences FACSLyric™ flow cytometers, BD Biosciences FACSVerse™ flow cytometers, BD Biosciences FACSymphony™ flow cytometers, BD Biosciences LSRFortessa™ flow cytometers, BD Biosciences LSRFortessa™ X-20 flow cytometers, BD Biosciences FACSPresto™ flow cytometers, BD Biosciences FACSVia™ flow cytometers, and BD Biosciences FACSCalibur™ cell sorters, BD Biosciences FACSCount™ cell sorters, BD Biosciences FACSLyric™ cell sorters, BD These include BD Biosciences Via™ cell sorters, BD Biosciences Influx™ cell sorters, BD Biosciences Jazz™ cell sorters, BD Biosciences Aria™ cell sorters, BD Biosciences FACSAria™ II cell sorters, BD Biosciences FACSAria™ III cell sorters, BD Biosciences FACSAria™ Fusion cell sorters, BD Biosciences FACSMelody™ cell sorters, and BD Biosciences FACSymphony™ S6 cell sorters.

[0117] In some embodiments, the subject system is a flow cytometry system, e.g., as described in U.S. Pat. Nos. 10,663,476, 10,620,111, 10,613,017, 10,605,713, 10,585,031, 10,578,542, 10,578,469, 10,481,074, 10,302,545, 10,145,793, 10,113,967, 10,006,852, 9,952,076, 9,933,341, 9,726,527, 9,453,789, 9,200,334, 9,000,491, 9,100,620, 9,200,492, 9,300,820, 9,400,440, 9,500,460, 9,600,480, 9,700,480, 9,800,480, 9,900,480, 9,100,480, 9,200,480, 9,300,480, 9,400,480, 9,500,480, 9,600,480, 9,700,480, 9,800,480, 9,900,480, 9,100,480, 9,100,480, 9,200,480, 9,300,480, 9,400,480, 9,500,480, 9,1 97640, U.S. Pat. No. 9095494, U.S. Pat. No. 9092034, U.S. Pat. No. 8975595, U.S. Pat. No. 8753573, U.S. Pat. No. 8233146, U.S. Pat. No. 8140300, U.S. Pat. No. 7544326, U.S. Pat. No. 7201875, U.S. Pat. No. 7129505, U.S. Pat. No. 6821740, U.S. Pat. No. 6813017, U.S. Pat. No. 6809804, U.S. Pat. No. 6372506, U.S. Pat. No. 5700692, U.S. Pat. No. 5643796, U.S. Pat. No. 5627040, U.S. Pat. No. 5620842, U.S. Pat. No. 5602039, U.S. Pat. No. 4987086, and U.S. Pat. No. 4498766.

[0118] 4 shows a system 400 for flow cytometry according to an exemplary embodiment of the invention. System 400 includes a flow cytometer 410, a controller / processor 490, and a memory 495. Flow cytometer 410 includes one or more excitation lasers 415a-415c, a focusing lens 420, a flow chamber 425, a forward scatter detector 430, a side scatter detector 435, a fluorescence collection lens 440, one or more beam splitters 445a-445g, one or more bandpass filters 450a-450e, one or more longpass (LP) filters 455a-455b, and one or more fluorescence detectors 460a-460f.

[0119] The pump lasers 415a-415c emit light in the form of laser beams. The wavelengths of the laser beams emitted from the pump lasers 415a-415c are 488 nm, 633 nm, and 325 nm, respectively, in the exemplary system of FIG. 4. The laser beams are first directed through one or more of the beam splitters 445a and 445b. The beam splitter 445a transmits the 488 nm light and reflects the 633 nm light. The beam splitter 445b transmits the UV light (light having a wavelength in the range of 10-400 nm) and reflects the 488 nm and 633 nm lights.

[0120] The laser beam is then directed to a focusing lens 420, which focuses the laser beam onto the portion of the fluid stream in which the sample particles reside in a flow chamber 425. The flow chamber is the part of the fluid system that directs particles in the stream, usually one at a time, into the focused laser beam for investigation. The flow chamber can comprise a flow cell in a benchtop flow cytometer or a nozzle tip in a stream-in-air cytometer.

[0121] Light from the laser beam(s) interacts with particles in the sample by diffraction, refraction, reflection, scattering and absorption, and is re-emitted at a variety of different wavelengths depending on the particle's properties, such as the particle's size, internal structure, and the presence of one or more fluorescent molecules attached to or naturally present on or in the particle. The fluorescent radiation, as well as the diffracted, refracted, reflected and scattered light, may be sent through one or more of beam splitters 445c-445g, bandpass filters 450a-450e, longpass filters 455a-455b, and fluorescence collection lens 440 to one or more of forward scatter detector 430, side scatter detector 435, and one or more fluorescence detectors 460a-460f.

[0122] Fluorescence collection lens 440 collects light emitted by particle-laser beam interactions and directs the light toward one or more beam splitters and filters. Bandpass filters, such as bandpass filters 450a-450e, allow a narrow range of wavelengths to pass through the bandpass filter. For example, bandpass filter 450a is a 510 / 20 filter. The first number represents the center of the spectral band. The second number represents the range of the spectral band. Thus, a 510 / 20 filter extends 10 nm on either side of the center of the spectral band, from 500 nm to 520 nm. Shortpass filters transmit light with wavelengths equal to or less than a specified wavelength. Longpass filters, such as longpass filters 455a-455b, transmit light with wavelengths equal to or greater than a specified wavelength of light. For example, longpass filter 455b, a 670 nm longpass filter, transmits light above 670 nm. Filters are often selected to optimize the specificity of the detector for a particular fluorochrome. The filters may be configured so that the spectral band of light transmitted to the detector is close to the emission peak of the fluorescent dye.

[0123] The forward scatter detector 430 is positioned slightly off-axis from the direct beam through the flow cell and is configured to detect diffracted light, i.e., excitation light traveling primarily in a forward direction through or around the particle. The intensity of light detected by the forward scatter detector depends on the overall size of the particle. The forward scatter detector may include a photodiode. The side scatter detector 435 is configured to detect refracted and reflected light at the surface and internal structure of the particle, which tends to increase as the complexity of the particle's structure increases. Fluorescent emission from fluorescent molecules bound to the particle may be detected by one or more of the fluorescence detectors 460a-460f. The side scatter detector 435 and the fluorescence detector may include photomultiplier tubes. The signals detected by the forward scatter detector 430, the side scatter detector 435 and the fluorescence detector may be converted to electronic signals (voltage) by the detectors. This data may provide information about the sample.

[0124] Those skilled in the art will recognize that flow cytometers according to embodiments of the present invention are not limited to the flow cytometer shown in Figure 4, but may include any flow cytometer known in the art. For example, a flow cytometer may have any number of lasers, beam splitters, filters, and detectors of various wavelengths and in a variety of different configurations.

[0125] During operation, the operation of the flow cytometer is controlled by the controller / processor 490, and measurement data from the detectors may be stored in memory 495 and processed by the controller / processor 490. Although not explicitly shown, the controller / processor 490 is coupled to the detectors to receive output signals from the detectors, and may further be coupled to electrical and electromechanical components of the flow cytometer 410 to control lasers, fluid flow parameters, and the like. Input / output (I / O) functionality 497 may also be provided in the system. The memory 495, controller / processor 490, and I / O functionality 497 may be provided entirely as an integral part of the flow cytometer 410. In such an embodiment, a display may also form part of the I / O functionality 497 to present experimental data to a user of the flow cytometer 410. Alternatively, some or all of the memory 495 and the controller / processor 490 and I / O functionality may be part of one or more external devices, such as a general purpose computer. In some embodiments, memory 495 and some or all of controller / processor 490 may be in wireless or wired communication with flow cytometer 410. Controller / processor 490 in conjunction with memory 495 and I / O functionality 497 may be configured to perform a variety of functions associated with the preparation and analysis of flow cytometer experiments.

[0126] The system illustrated in FIG. 4 includes six different detectors that detect fluorescence in six different wavelength bands (which may be referred to herein as “filter windows” for a given detector) as determined by the configuration of filters and / or splitters in the beam path from the flow cell 425 to each detector. Different fluorescent molecules in a panel of fluorochromes used in a flow cytometer experiment emit light in their own unique wavelength bands. The particular fluorescent labels used in the experiment and the associated fluorescent emission bands of the fluorescent labels may be selected to approximately match the filter windows of the detectors. The I / O functionality 497 may be configured to receive data regarding a flow cytometer experiment having a panel of fluorescent labels and multiple cell populations having multiple markers, each cell population having a subset of multiple markers. The I / O functionality 497 may be further configured to receive biological data assigning one or more markers to one or more cell populations, marker concentration data, emission spectrum data, data assigning labels to one or more markers, and cytometer configuration data. Flow cytometer experiment data, such as label spectral characteristics and flow cytometer configuration data, may be further stored in memory 495. The controller / processor 490 may be configured to evaluate one or more assignments of labels to markers.

[0127] In some embodiments, the subject system is a particle sorting system configured to sort particles using an enclosed particle sorting module, such as described in U.S. Patent Application Publication No. 2017 / 0299493, filed March 28, 2017, the disclosure of which is incorporated herein by reference. In an embodiment, particles (e.g., cells) of a sample are sorted using a sorting determination module having a plurality of sorting determination units, such as described in U.S. Patent Application Publication No. 2020 / 0256781, filed December 23, 2019, the disclosure of which is incorporated herein by reference. In some embodiments, a system for sorting components of a sample includes a particle sorting module having deflection plates, such as described in U.S. Patent Application Publication No. 2017 / 0299493, filed March 28, 2017, the disclosure of which is incorporated herein by reference.

[0128] 5 is a functional block diagram of an example of a sorting control system for analyzing and displaying biological events, such as a processor 500. The processor 500 can be configured to perform various processes for controlling the graphical display of biological events.

[0129] A flow cytometer or sorting system 502 may be configured to acquire bio-event data. For example, a flow cytometer may generate flow cytometry event data (e.g., particle-modulated light data). The flow cytometer 502 may be configured to provide the bio-event data to the processor 500. A data communication channel may be provided between the flow cytometer 502 and the processor 500. The bio-event data may be provided to the processor 500 via the data communication channel.

[0130] The processor 500 may be configured to receive bio-event data from the flow cytometer 502. The bio-event data received from the flow cytometer 502 may include flow cytometry event data. The processor 500 may be configured to provide a graphical display including a first plot of the bio-event data to the display device 506. The processor 500 may be further configured to render a region of interest as a gate (e.g., a first gate) around a population of the bio-event data displayed by the display device 506, e.g., overlaid on the first plot. In some embodiments, the gate may be a logical combination of one or more illustrated regions of interest plotted on a histogram or bivariate plot of a parameter. In some embodiments, the display may be used to display particle parameters or saturation detector data.

[0131] The processor 500 may further be configured to display the bioevent data on the display device 506 within the gate differently than other events within the bioevent data outside the gate. For example, the processor 500 may be configured to render the color of the bioevent data included within the gate differently than the color of the bioevent data outside the gate. The display device 506 may be implemented as a monitor, tablet computer, smartphone, or other electronic device configured to provide a graphical interface.

[0132] The processor 500 may be configured to receive a gate selection signal identifying a gate from a first input device. For example, the first input device may be implemented as a mouse 510. The mouse 510 may transmit a gate selection signal identifying a gate to be displayed on or manipulated via the display device 506 (e.g., by clicking on or within the desired gate when the cursor is on the desired gate) to the processor 500. In some implementations, the first device may be implemented as a keyboard 508 or other means for providing input signals to the processor 500, such as a touch screen, a pen for input, a photodetector, or a voice recognition system. Some input devices may include multiple input functions. In such implementations, each input function may be considered an input device. For example, as shown in FIG. 5, the mouse 510 may include a right mouse button and a left mouse button, and the right mouse button and the left mouse button may each generate a trigger event.

[0133] A trigger event can cause the processor 500 to change how the data is displayed, what portions of the data are actually displayed on the display device 506, and / or provide input for further processing, such as selection of a population of interest for particle sorting.

[0134] In some embodiments, the processor 500 may be configured to detect when a gate selection is initiated by the mouse 510. The processor 500 may further be configured to automatically modify the visualization of the plot to facilitate gating. This modification may be based on a particular distribution of the bio-event data received by the processor 500. In some embodiments, the processor 500 extends the first gate such that a second gate is generated (e.g., as described above).

[0135] The processor 500 may be connected to a storage device 504. The storage device 504 may be configured to receive and store bio-event data from the processor 500. The storage device 504 may be further configured to receive and store flow cytometry event data from the processor 500. The storage device 504 may be further configured to enable retrieval of bio-event data, such as flow cytometry event data, by the processor 500.

[0136] The display device 506 may be configured to receive display data from the processor 500. The display data may include plots of the biological event data and gates outlining sections of the plots. The display device 506 may be further configured to modify the information displayed in response to input received from the processor 500 in conjunction with input from the flow cytometer 502, the storage device 504, the keyboard 508, and / or the mouse 510.

[0137] In some embodiments, the processor 500 can generate a user interface for receiving exemplary events for screening. For example, the user interface can include a mechanism for receiving exemplary events or exemplary images. The exemplary events or images, or exemplary gates, can be provided prior to collection of event data for a sample, or based on an initial set of events for a portion of a sample.

[0138] FIG. 6A is a schematic diagram illustrating a particle sorting system 600 (e.g., flow cytometer 502) according to one embodiment presented herein. In some embodiments, the particle sorting system 600 is a cell sorting system. As shown in FIG. 6A, a droplet forming transducer 602 (e.g., a piezoelectric oscillator) is coupled to a fluid conduit 601, which may be coupled to, may include, or may be a nozzle 603. Within the fluid conduit 601, a sheath fluid 604 hydrodynamically focuses a sample fluid 606 containing particles 609 into a moving fluid column 608 (e.g., a stream). Within the moving fluid column 608, the particles 609 (e.g., cells) are aligned in a single file across a monitoring region 611 (e.g., where a laser stream intersects) and are illuminated by an illumination source 612 (e.g., a laser). Vibration of the droplet forming transducer 602 causes the moving fluid column 608 to break up into multiple droplets 610 , some of which contain particles 609 .

[0139] In operation, a detection station 614 (e.g., an event detector) identifies when a particle of interest (or cell of interest) crosses the monitoring region 611. The detection station 614 feeds a timing circuit 628, which in turn feeds a flash charge circuit 630. At a droplet break-off point, signaled by a timed droplet delay (Δt), a flash charge can be applied to the moving fluid column 608 such that the droplet of interest carries a charge. The droplet of interest can contain one or more particles or cells to be sorted. The charged droplets can then be sorted by activating a deflection plate (not shown) to deflect the charged droplets into a receptacle, such as a collection tube or a multi-well or microwell sample plate, where a well or microwell can be specifically associated with the droplet of interest. As shown in FIG. 6A, the droplets can be collected in a drain receptacle 638.

[0140] The detection system 616 (e.g., a droplet boundary detector) serves to automatically determine the phase of the droplet drive signal when a particle of interest passes through the monitoring region 611. An exemplary droplet boundary detector is described in U.S. Pat. No. 7,679,039, the entire contents of which are incorporated herein by reference. The detection system 616 enables the instrument to accurately calculate the position of each detected particle in the droplet. The detection system 616 may provide an amplitude signal 620 and / or a phase signal 618, which are then provided (via an amplifier 622) to an amplitude control circuit 626 and / or a frequency control circuit 624. The amplitude control circuit 626 and / or the frequency control circuit 624 then control the droplet forming transducer 602. The amplitude control circuit 626 and / or the frequency control circuit 624 may be provided within a control system.

[0141] In some embodiments, the sorting electronics (e.g., detection system 616, detection station 614, and processor 640) may be coupled to a memory configured to store detected events and sorting decisions based on the detected events. The sorting decisions may be included in the event data for the particles. In some embodiments, the detection system 616 and detection station 614 may be implemented as a single detection unit or may be communicatively coupled such that event measurements may be collected by one of the detection system 616 or detection station 614 and provided to a non-collection element.

[0142] FIG. 6B is a schematic diagram showing a particle sorting system according to one embodiment presented herein. The particle sorting system 600 shown in FIG. 6B includes deflection plates 652 and 654. An electric charge can be applied via a stream charging wire in the barb. This generates a stream of droplets 610 containing particles 609 for analysis. The particles can be illuminated with one or more light sources (e.g., lasers) to generate light scattering and fluorescence information. The information about the particles is analyzed, such as by sorting electronics (not shown in FIG. 6B) or other detection systems. Deflection plates 652 and 654 can be independently controlled to attract or repel the charged droplets and direct them toward a desired collection vessel (e.g., one of 672, 674, 676, or 678). 6B, deflector plates 652 and 654 can be controlled to direct particles along a first path 662 toward a container 674 or along a second path 668 toward a container 678. If the particle is not a particle of interest (e.g., does not exhibit scattering or illumination information within a specified sorting range), the deflector plates may allow the particle to continue along flow path 664. Such uncharged droplets may be channeled into a waste container, such as via an aspirator 670.

[0143] Sorting electronics can be included to initiate measurement collection, receive fluorescent signals for the particles, and determine how to adjust the deflection plates to sort the particles. Exemplary implementations of the embodiment shown in Figure 6B include the BD FACSAria™ line of flow cytometers marketed by Becton, Dickinson and Company (Franklin Lakes, NJ).

[0144] Computer Control System Aspects of the present disclosure further include a computer control system, the computer control system comprising one or more computers for full or partial automation. In some embodiments, the system comprises a computer having a non-transitory computer readable storage medium having a computer program stored thereon, the computer program comprising instructions for receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2, when loaded into the computer. Additionally, the computer comprises instructions for selecting two or more subsets of the N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with the instrument identifier (each fluorochrome identifier in the spectral matrix is ​​associated with a spectral feature), and identifying a fluorochrome panel from the two or more subsets of the N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using the spectral features associated with the N fluorochrome identifiers of the two or more subsets.

[0145] The system may include a display and an operator input device. The operator input device may be, for example, a keyboard, a mouse, etc. The processing module has a processor that can access a memory in which instructions for performing the steps of the subject method are stored. The processing module may include an operating system, a graphical user interface (GUI) controller, a system memory, a memory storage device, an input / output controller, a cache memory, a data backup unit, and many other devices. The processor may be a commercially available processor or one of other processors that are available or become available. The processor executes an operating system, which interfaces with firmware and hardware in a well-known manner to facilitate the processor coordinating and executing the functions of various computer programs, which may be written in various programming languages, such as Java, Perl, C++, Python, other high-level or low-level languages, and combinations thereof, as is known in the art. The operating system generally cooperates with the processor to coordinate and execute the functions of the other elements of the computer. The operating system further provides scheduling, input / output control, file and data management, memory management, and communication control and related services, all according to known techniques. In some embodiments, the processor includes analog electronics that provide feedback control, such as negative feedback control.

[0146] The system memory may be any of a variety of known or future memory storage devices. Examples include any commonly available random access memory (RAM), magnetic media such as resident hard disks or tapes, optical media such as read-write compact disks, flash memory devices or other memory storage devices. The memory storage devices may be any of a variety of known or future devices including compact disk drives, tape drives, or diskette drives. Such types of memory storage devices typically read from and / or write to a program storage medium (not shown) such as a compact disk. Any of these program storage media, or other program storage media now in use or that may be developed in the future, may be considered a computer program product. As will be appreciated, these program storage media typically store computer software programs and / or data. Computer software programs, also referred to as computer control logic, are typically stored in the system memory and / or in program storage devices used in conjunction with the memory storage devices.

[0147] In some embodiments, a computer program product is described that includes a computer usable medium having stored thereon control logic (a computer software program including program code). The control logic, when executed by a processor, a computer, causes the processor to perform the functions described herein. In other embodiments, some functions are implemented primarily in hardware, for example using hardware state machines. The implementation of hardware state machines to perform the functions described herein will be apparent to one of ordinary skill in the relevant art.

[0148] The memory may be any suitable device from which the processor can store and retrieve data, such as a magnetic device, an optical device, or a solid-state storage device (including a magnetic or optical disk or tape or RAM or any other suitable device, fixed or portable). The processor may include a general-purpose digital microprocessor suitably programmed from a computer-readable medium storing the necessary program code. The program may be provided to the processor remotely via a communication channel, or may be pre-recorded in a computer program product such as a memory, or in other portable or fixed computer-readable storage media using either a device connected to the memory. For example, a magnetic or optical disk may store the program and may be read by a disk writer / reader. The system of the present invention further comprises a program, for example in the form of a computer program product, an algorithm for use in implementing the method as described above. The program of the present invention may be recorded on a computer-readable medium, such as any medium that can be directly read and accessed by a computer. Such media include, but are not limited to, magnetic storage media; optical storage media, such as CD-ROM; electrical storage media, such as RAM and ROM; portable flash drives; and hybrids of these categories, such as magnetic / optical storage media.

[0149] The processor may also have access to a communication channel to communicate with a user at a remote location, meaning that the user does not have direct contact with the system, but rather relays input information to the input manager from an external device, such as a computer connected to a wide area network ("WAN"), a telephone network, a satellite network, or any other suitable communication channel, including a cellular phone (i.e., a smartphone).

[0150] In some embodiments, a system according to the present disclosure may be configured to include a communication interface. In some embodiments, the communication interface includes a receiver and / or a transmitter for communicating with a network and / or another device. The communication interface may be configured for wired or wireless communication, including but not limited to radio frequency (RF) communication, such as radio frequency identification (RFID), ZigBee communication protocol, Wi-Fi, infrared, wireless universal serial bus (USB), ultra-wideband (UWB), Bluetooth® communication protocol, and cellular communication, such as code division multiple access (CDMA) or Global System for Mobile Communications (GSM).

[0151] In one embodiment, the communication interface is configured to include one or more communication ports, e.g., physical ports or interfaces, such as a USB port, a USB-C port, an RS-232 port, or any other suitable electrical connection port, to enable data communication between the subject system and other external devices, such as computer terminals (e.g., in a clinic or hospital environment) configured for similar complementary data communication.

[0152] In one embodiment, the communications interface is configured for infrared communications, Bluetooth® communications, or any other suitable wireless communications protocol, allowing the subject system to communicate with other devices, such as computer terminals and / or networks, communication-enabled mobile phones, personal digital assistants, or any other communications devices that a user may use in conjunction therewith.

[0153] In one embodiment, the communication interface is configured to provide a connection for data transfer utilizing Internet Protocol (IP) via a cellular network, Short Message Service (SMS), a wireless connection to a personal computer (PC) on a local area network (LAN) connected to the Internet, or a Wi-Fi connection to the Internet at a Wi-Fi hotspot.

[0154] In one embodiment, the subject system is configured to communicate wirelessly with the server device via a communication interface using common standards such as 802.11, Bluetooth RF protocol, or IrDA infrared protocol. The server device may be another portable device, such as a smart phone, a personal digital assistant (PDA), or a notebook computer; or a larger device, such as a desktop computer, an appliance, etc. In some embodiments, the server device has a display, such as a liquid crystal display (LCD), and an input device, such as a button, a keyboard, a mouse, or a touch screen.

[0155] In some embodiments, the communications interface is configured to automatically or semi-automatically communicate data stored in the subject system, e.g., any data storage unit, with a network or server device using one or more of the communications protocols and / or mechanisms described above.

[0156] The output controller may include a controller for any of a variety of known display devices for presenting information to a user, whether human or machine, local or remote. Where one of the display devices provides visual information, this information may generally be logically and / or physically organized as an array of pixels. A graphical user interface (GUI) controller may include any of a variety of known or future software programs for providing a graphical input / output interface between the system and the user, and for processing user input. The functional elements of the computer may communicate with each other via a system bus. Some of these communications may be accomplished using a network or other type of remote communication in alternative embodiments. The output manager may further provide information generated by the processing module to a user at a remote location, for example, via the Internet, telephone, or satellite network, according to known techniques. Presentation of data by the output manager may be performed according to a variety of known techniques. As some examples, the data may include SQL, HTML, or XML documents, emails, or other files, or other forms of data. The data may include Internet URL addresses so that the user may retrieve additional SQL, HTML, XML, or other documents or data from remote sources. The platform or platforms present in the subject system are typically a class of computers commonly referred to as servers, but may be any type of computer platform now known or later developed. However, the platforms may also be mainframe computers, workstations, or other computer types. The platforms may be connected via any type of cabling, now known or later, or other communication systems, including wireless systems, whether networked or not. The platforms may be co-located or physically separated.Various operating systems may be used on any of the computer platforms, depending in some cases on the type and / or configuration of the computer platform selected. Suitable operating systems include Windows NT, Windows XP, Windows 7, Windows 8, Windows 10, iOS, macOS, Linux, Ubuntu, Fedora, OS / 400, i5 / OS, IBM i, Android™, SGI IRIX, Oracle Solaris, and the like.

[0157] FIG. 7 illustrates a general configuration of an exemplary computing device 700 according to an embodiment. The general configuration of the computing device 700 illustrated in FIG. 7 includes an arrangement of computer hardware and software components. However, not all of these general conventional elements need to be illustrated to provide a useful disclosure. As illustrated, the computing device 700 includes a processing unit 710, a network interface 720, a computer-readable medium drive 730, an input / output device interface 740, a display 750, and an input device 760, all of which may communicate with each other via a communication bus. The network interface 720 may provide connectivity to one or more networks or computing systems. Thus, the processing unit 710 may receive information and instructions from other computing systems or services via a network. The processing unit 710 may further communicate with a memory 770 and may further provide output information for an optional display 750 via the input / output device interface 740. For example, analysis software (e.g., data analysis software or programs, such as FlowJo®), stored as executable instructions in non-transitory memory of the analysis system, can display flow cytometry event data to a user. Input / output device interface 740 may further accept input from any input device 760, such as a keyboard, mouse, digital pen, microphone, touch screen, gesture recognition system, voice recognition system, game pad, accelerometer, gyroscope, or other input device.

[0158] Memory 770 may include computer program instructions (grouped in some embodiments as modules or components) that processing unit 710 executes to implement one or more embodiments. Memory 770 typically includes RAM, ROM, and / or other persistent, secondary, or non-transitory computer-readable media. Memory 770 may store an operating system 772 that provides computer program instructions for use by processing unit 710 in the general management and operation of computing device 700. Data may be stored in data storage device 790. Memory 770 may further include computer program instructions and other information for implementing aspects of the present disclosure.

[0159] usefulness The methods, systems, and computer readable media of the present invention may be used when it is desirable to automatically determine a set of fluorochromes that can be used for particle analysis (e.g., flow cytometry). In some cases, the present invention is particularly useful for full-spectrum (i.e., "spectral") flow cytometry panel experimental design. In other words, the present invention serves as a first step in spectral panel design by identifying whether a set of dyes can be used simultaneously. In some cases, the methods, systems, and computer readable media described herein help determine which set or sets of fluorochromes are likely to provide the highest quality data (e.g., greatest biological resolution). In further cases, the present invention may be used to determine the maximum number of fluorochromes from a given set of fluorochromes (e.g., a given set of all commercially available flow cytometry fluorochromes) that may be measured simultaneously in a full-spectrum flow cytometer. The present invention accomplishes the above through an automated optimization algorithm, the use of the spectral characteristics of the fluorochromes as readily available and easy-to-measure inputs to the algorithm, and the use of a spectral matrix as a computationally efficient heuristic method for optimization.

[0160] The embodiments of the present invention are further used in applications where cells prepared from biological samples may be desired for use in research, clinical trials, or therapy. In some embodiments, the subject methods and devices may facilitate obtaining individual cells prepared from a target fluid or tissue biological sample. For example, the subject methods and systems may facilitate obtaining cells from a fluid or tissue sample used as a research or diagnostic specimen for a disease such as cancer. Similarly, the subject methods and systems may facilitate obtaining cells from a fluid or tissue sample used for therapy. The disclosed methods and devices allow for the separation and collection of cells from biological samples (e.g., organs, tissues, tissue slices, fluids) with high efficiency and low cost compared to conventional flow cytometry systems.

[0161] kit Aspects of the present disclosure further include kits, the kits comprising instructions and / or programmable logic for carrying out the claimed methods, for example, the kits comprising a program configured to identify a fluorochrome panel (e.g., as described above in the Method section), such as in the form of a computer readable medium (e.g., a flash drive, USB storage, compact disc, DVD, Blu-ray disc, etc.) or instructions for downloading the program from an Internet web protocol or cloud server.

[0162] The kits may further include instructions for practicing the subject methods. These instructions may be present in the subject kits in a variety of forms, one or more of which may be present in the kit. One form in which these instructions may be provided is information printed on a suitable medium or substrate, such as one or more pieces of paper on which the information is printed, kit packaging, a package insert, etc. Yet another form in which these instructions may be present is a computer readable medium on which the information is recorded, such as a diskette, a compact disc (CD), a portable flash drive, etc. Yet another form in which these instructions may be present is a website address that may be used via the internet to access the information at a remote location.

[0163] The following examples are offered by way of illustration and not by way of limitation.

[0164] experiment Experiments were performed to demonstrate the quality of fluorochrome panels identified by the claimed method versus randomized fluorochrome panels. Condition numbers were calculated for 100 sets of random fluorochrome panels and 10 / 100 sets of fluorochrome panels identified according to the method. These calculations were performed for each of a series of different panel numbers (i.e., N=29 to N=43). As shown in Figure 8, optimization results in an order of magnitude improvement over random fluorophore sets.

[0165] In addition, 100 different trials were performed to identify a fluorochrome panel with panel number 41. As shown in Figure 9, the condition numbers associated with the subset of N fluorochrome identifiers steadily decrease as the number of iterations of the minimization algorithm increases. The instrument identifiers for these trials were associated with a BD FACSymphony™ A5 SE cell analyzer.

[0166] Regardless of the scope of the appended claims, the present disclosure is further defined by the following caveats.

[0167] Appendix 1. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2; selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with an instrument identifier, each fluorochrome identifier in the spectral matrix being associated with a spectral feature; A method for identifying a fluorochrome panel from two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of the N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets.

[0168] Appendix 2. The method of appendix 1, in which the numerical stability index is the condition number.

[0169] Appendix 3. The method of Appendix 2, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0170] Appendix 4. The method of appendix 3, in which the ratio of the largest singular value to the smallest singular value is determined when calculating the condition number.

[0171] Appendix 5. The method of Appendix 2, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0172] Appendix 6. When calculating the condition number, for each of two or more subsets of N fluorescent dye identifiers, the Frobenius norm of the spectral matrix associated with a subset of N fluorochrome identifiers, and the Frobenius norm of the pseudoinverse of each spectral matrix associated with two or more subsets of N fluorochrome identifiers 6. The method of claim 5, further comprising determining

[0173] Appendix 7. The method of appendix 6, wherein when calculating the condition number, for each of two or more subsets of N fluorochrome identifiers, a product of the Frobenius norm of a spectral matrix associated with the subset of N fluorochrome identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of N fluorochrome identifiers is determined.

[0174] Addendum 8. The method of any one of Addendums 1-7, wherein the fluorescent dye panel comprises an optimized numerical stability index.

[0175] Appendix 9. The method of appendix 8, wherein a constrained optimization method is used to identify the fluorescent dye panel.

[0176] Addendum 10. The method of Addendum 9, wherein the constrained optimization method is a minimization algorithm.

[0177] Addendum 11. The method of Addendum 9 or 10, wherein the constrained optimization method is carried out sequentially.

[0178] Addendum 12. The method of Addendum 11, wherein when selecting two or more subsets of N fluorescent dye identifiers, a first subset of N fluorescent dye identifiers and a second subset of N fluorescent dye identifiers are selected.

[0179] Addendum 13. The method of Addendum 12, wherein a first subset of N fluorescent dye identifiers is randomly obtained from a spectral matrix associated with the instrument identifier.

[0180] Addendum 14. The method of Addendum 12 or 13, wherein the constrained optimization method further calculates a numerical stability index associated with a first subset of the selected N fluorescent dye identifiers.

[0181] Addendum 15. A method according to any one of Addendums 12 to 14, wherein when generating the second subset of N fluorescent dye identifiers, fluorescent dye identifiers in the first subset of N fluorescent dye identifiers are replaced with alternative fluorescent dye identifiers in a spectral matrix associated with the instrument identifier.

[0182] Addendum 16. The method of Addendum 14 or 15, wherein the constrained optimization method calculates a numerical stability index associated with a second subset of the N fluorescent dye identifiers.

[0183] Addendum 17. The method of Addendum 16, wherein the constrained optimization method compares a numerical stability metric associated with a first subset of the N fluorescent dye identifiers to a numerical stability metric associated with a second subset of the N fluorescent dye identifiers.

[0184] Addendum 18. A method according to any one of Addendums 11 to 17, wherein the constrained optimization method further comprises iteratively generating a subset of fluorescent dye identifiers.

[0185] Appendix 19. In the constrained optimization method, Calculate a numerical stability index for each subset of N fluorescent dye identifiers generated iteratively; The method described in Appendix 18, further comprising comparing a numerical stability index calculated for each subset of N fluorescent dye identifiers generated iteratively with a numerical stability index of a subset of N fluorescent dye identifiers already generated, and determining an optimized numerical stability index.

[0186] Addendum 20. The method of any one of Addendums 8 to 19, wherein the optimized numerical stability index is a minimized condition number.

[0187] Addendum 21. The method of any one of Addendums 1 to 20, wherein the spectral feature is a spillover value.

[0188] Addendum 22. The method of any one of Addendums 1 to 21, wherein the spectral features are obtained from experimental data.

[0189] Addendum 23. The method of any one of Addendums 1 to 21, wherein the spectral features are obtained from simulated data.

[0190] Appendix 24. The method of any one of appendices 1 to 23, performed without evaluation of antigen data.

[0191] Addendum 25. The method of any one of Addendums 1 to 24, wherein N is an integer in the range of 2 to 100.

[0192] Addendum 26. The method of Addendum 25, wherein N is an integer in the range of 2 to 75.

[0193] Addendum 27. The method of any one of Addendums 1 to 26, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 150.

[0194] Addendum 28. The method of Addendum 27, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 100.

[0195] Addendum 29. The method of any one of Addendums 1 to 28, comprising assessing the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0196] Addendum 30. The method of Addendum 29, wherein an optimized numerical stability index is determined for each of multiple values ​​of N when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0197] Addendum 31. The method of Addendum 30, wherein the multiple values ​​of N include i values, where i is the number of fluorescent dye identifiers in the spectral matrix associated with the instrument identifier.

[0198] Addendum 32. The method of Addendum 30 or 31, in which when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel, a maximum value of N is determined that does not result in an optimized numerical stability index exceeding a threshold.

[0199] Addendum 33. The method of any one of Addendums 1 to 32, wherein the spectral matrix associated with the device identifier includes autofluorescence spectral features.

[0200] Addendum 34. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: receiving a request for N fluorochrome identifiers, where N is an integer of at least two; A method comprising identifying a fluorochrome panel comprising a subset of the N fluorochrome identifiers, the method being performed without evaluating antigen data.

[0201] Appendix 35. Receiving a device identifier; 35. The method of claim 34, further comprising selecting two or more subsets of N fluorescent dye identifiers from a set of fluorescent dye identifiers in a spectral matrix associated with an instrument identifier, each fluorescent dye identifier in the spectral matrix being associated with a spectral feature.

[0202] Addendum 36. The method of Addendum 35, wherein the numerical stability index is a condition number.

[0203] Addendum 37. The method of Addendum 36, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0204] Addendum 38. The method of Addendum 37, wherein when calculating the condition number, the ratio of the largest singular value to the smallest singular value is determined.

[0205] Addendum 39. The method of Addendum 36, wherein when calculating the condition number, for each of two or more subsets of N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of N fluorescent dye identifiers is determined.

[0206] Addendum 40. When calculating the condition number, for each of two or more subsets of N fluorochrome identifiers, the Frobenius norm of the spectral matrix associated with a subset of N fluorochrome identifiers, and the Frobenius norm of the pseudoinverse of each spectral matrix associated with two or more subsets of N fluorochrome identifiers 40. The method of claim 39, further comprising determining

[0207] Addendum 41. The method of claim 40, wherein when calculating the condition number, for each of two or more subsets of N fluorochrome identifiers, a product of the Frobenius norm of a spectral matrix associated with the subset of N fluorochrome identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of N fluorochrome identifiers is determined.

[0208] Addendum 42. The method of any one of Addendums 35 to 41, wherein the fluorescent dye panel includes an optimized numerical stability index.

[0209] Addendum 43. The method of Addendum 42, wherein a constrained optimization method is used to identify the fluorescent dye panel.

[0210] Addendum 44. The method of Addendum 43, wherein the constrained optimization method is a minimization algorithm.

[0211] Addendum 45. The method of Addendum 43 or 44, wherein the constrained optimization method is carried out sequentially.

[0212] Addendum 46. The method of Addendum 45, wherein when selecting two or more subsets of N fluorescent dye identifiers, a first subset of N fluorescent dye identifiers and a second subset of N fluorescent dye identifiers are selected.

[0213] Addendum 47. The method of Addendum 46, wherein a first subset of N fluorescent dye identifiers is randomly obtained from a spectral matrix associated with the instrument identifier.

[0214] Addendum 48. A method as described in Addendum 46 or 47, wherein the constrained optimization method further calculates a numerical stability index associated with a first subset of the selected N fluorescent dye identifiers.

[0215] Addendum 49. A method according to any one of Addendums 46 to 48, wherein when generating the second subset of N fluorescent dye identifiers, fluorescent dye identifiers in the first subset of N fluorescent dye identifiers are replaced with alternative fluorescent dye identifiers in a spectral matrix associated with the instrument identifier.

[0216] Addendum 50. The method of Addendum 48 or 49, wherein the constrained optimization method calculates a numerical stability index associated with a second subset of the N fluorescent dye identifiers.

[0217] Addendum 51. The method of Addendum 50, wherein the constrained optimization method compares a numerical stability metric associated with a first subset of the N fluorescent dye identifiers to a numerical stability metric associated with a second subset of the N fluorescent dye identifiers.

[0218] Addendum 52. A method according to any one of Addendums 45 to 51, wherein the constrained optimization method further comprises iteratively generating a subset of fluorescent dye identifiers.

[0219] Note 53. In the constrained optimization method, Calculate a numerical stability index for each subset of N fluorescent dye identifiers generated iteratively; The method described in Appendix 52, further comprising comparing a numerical stability index calculated for each subset of N fluorescent dye identifiers generated iteratively with a numerical stability index of a subset of N fluorescent dye identifiers already generated, and determining an optimized numerical stability index.

[0220] Addendum 54. The method of any one of Addendums 42 to 53, wherein the optimized numerical stability index is a minimized condition number.

[0221] Addendum 55. The method of any one of Addendums 35 to 54, wherein the spectral feature is a spillover value.

[0222] Addendum 56. The method of any one of Addendums 35 to 55, wherein the spectral features are obtained from experimental data.

[0223] Addendum 57. The method of any one of Addendums 35 to 55, wherein the spectral features are obtained from simulated data.

[0224] Addendum 58. The method of any one of Addendums 34 to 57, wherein N is an integer in the range of 2 to 100.

[0225] Addendum 59. The method of Addendum 58, wherein N is an integer in the range of 2 to 75.

[0226] Addendum 60. The method of any one of Addendums 35 to 59, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 150.

[0227] Addendum 61. The method of Addendum 60, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 100.

[0228] Addendum 62. The method of any one of Addendums 34 to 61, comprising assessing the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0229] Addendum 63. The method of Addendum 62, wherein an optimized numerical stability index is determined for each of multiple values ​​of N when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0230] Addendum 64. The method of Addendum 63, wherein the multiple values ​​of N include i values, where i is the number of fluorescent dye identifiers in the spectral matrix associated with the instrument identifier.

[0231] Addendum 65. The method of Addendum 63 or 64, wherein when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel, a maximum value of N is determined that does not result in an optimized numerical stability index exceeding a threshold.

[0232] Addendum 66. The method of any one of Addendums 34 to 65, wherein the spectral matrix associated with the device identifier includes autofluorescence spectral features.

[0233] Addendum 67. Receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2; selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with an instrument identifier, each fluorochrome identifier in the spectral matrix being associated with a spectral feature; Identifying a fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. 16. A system comprising: a processor configured to:

[0234] Addendum 68. The system of Addendum 67, wherein the numerical stability index is a condition number.

[0235] Addendum 69. The system of Addendum 68, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0236] Addendum 70. The system of Addendum 69, wherein when calculating the condition number, a ratio of a maximum singular value to a minimum singular value is determined.

[0237] Addendum 71. The system of Addendum 68, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0238] Addendum 72. When calculating the condition number, for each of two or more subsets of N fluorochrome identifiers, the Frobenius norm of the spectral matrix associated with a subset of N fluorochrome identifiers, and the Frobenius norm of the pseudoinverse of each spectral matrix associated with two or more subsets of N fluorochrome identifiers 72. The system of claim 71, further comprising:

[0239] Addendum 73. The system of Addendum 72, further comprising, when calculating the condition number, determining, for each of two or more subsets of the N fluorescent dye identifiers, the product of the Frobenius norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of the N fluorescent dye identifiers.

[0240] Addendum 74. A system described in any one of Addendums 67 to 73, wherein the fluorescent dye panel includes an optimized numerical stability index.

[0241] Addendum 75. The system of Addendum 74, wherein a constrained optimization method is used to identify the fluorescent dye panel.

[0242] Addendum 76. The system of Addendum 75, wherein the constrained optimization method is a minimization algorithm.

[0243] Addendum 77. The system of any one of Addendum 75 to 76, which sequentially performs a constrained optimization method.

[0244] Addendum 78. The system of Addendum 77, wherein when selecting two or more subsets of N fluorescent dye identifiers, a first subset of N fluorescent dye identifiers and a second subset of N fluorescent dye identifiers are selected.

[0245] Addendum 79. The system of Addendum 78, wherein a first subset of N fluorescent dye identifiers is randomly obtained from a spectral matrix associated with the instrument identifier.

[0246] Addendum 80. The system of Addendum 78 or 79, wherein the constrained optimization method further calculates a numerical stability index associated with a first subset of the selected N fluorescent dye identifiers.

[0247] Addendum 81. A system described in any one of Addendums 78 to 80, wherein when generating the second subset of N fluorescent dye identifiers, fluorescent dye identifiers in the first subset of N fluorescent dye identifiers are replaced with alternative fluorescent dye identifiers in a spectral matrix associated with the instrument identifier.

[0248] Addendum 82. The system of Addendum 80 or 81, wherein the constrained optimization method calculates a numerical stability index associated with a second subset of the N fluorescent dye identifiers.

[0249] Addendum 83. The system of Addendum 82, wherein the constrained optimization method compares a numerical stability metric associated with a first subset of the N fluorescent dye identifiers to a numerical stability metric associated with a second subset of the N fluorescent dye identifiers.

[0250] Addendum 84. A system according to any one of Addendums 77 to 83, wherein the constrained optimization method further comprises iteratively generating a subset of fluorescent dye identifiers.

[0251] Note 85. In constrained optimization, Calculate a numerical stability index for each subset of N fluorescent dye identifiers generated iteratively; The system of claim 84, further comprising: comparing a numerical stability index calculated for each subset of N fluorescent dye identifiers iteratively generated with a numerical stability index of a subset of N fluorescent dye identifiers that have already been generated, and determining an optimized numerical stability index.

[0252] Addendum 86. The system of any one of Addendums 74 to 85, wherein the optimized numerical stability metric is a minimized condition number.

[0253] Addendum 87. The system of any one of Addendums 67 to 86, wherein the spectral feature is a spillover value.

[0254] Addendum 88. The system of any one of Addendums 67 to 87, wherein the spectral features are obtained from experimental data.

[0255] Addendum 89. The system of any one of Addendums 67 to 87, wherein the spectral features are obtained from simulated data.

[0256] Addendum 90. The system of any one of Addendums 67-89, wherein the processor is configured to identify the fluorochrome panel without evaluating antigen data.

[0257] Addendum 91. The system of any one of Addendums 67 to 90, wherein N is an integer in the range of 2 to 100.

[0258] Addendum 92. The system of Addendum 91, wherein N is an integer in the range of 2 to 75.

[0259] Addendum 93. A system described in any one of Addendums 67 to 92, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 150.

[0260] Addendum 94. The system of Addendum 93, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 100.

[0261] Addendum 95. A system described in any one of Addendums 67 to 94, wherein the processor is configured to estimate the maximum number of fluorescent dyes that may be used simultaneously in the fluorescent dye panel.

[0262] Addendum 96. The system of Addendum 95, which determines an optimized numerical stability index for each of multiple values ​​of N when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0263] The system of claim 96, wherein the plurality of values ​​of N includes i values, where i is the number of fluorescent dye identifiers in the spectral matrix associated with the instrument identifier.

[0264] Addendum 98. A system described in Addendum 96 or 97, which determines a maximum value of N that does not result in an optimized numerical stability index exceeding a threshold when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0265] Addendum 99. The system of any one of Addendums 67-98, comprising a particle analyzer.

[0266] Addendum 100. The system of Addendum 99, wherein the particle analyzer is a flow cytometer.

[0267] Addendum 101. The system of Addendum 100, wherein the flow cytometer is a full-spectrum flow cytometer.

[0268] Addendum 102. A system described in any one of Addendums 67 to 101, wherein the spectral matrix associated with the device identifier includes autofluorescence spectral features.

[0269] Clause 103. Receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2; selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with an instrument identifier, each fluorochrome identifier in the spectral matrix being associated with a spectral feature; Identifying a fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. A non-transitory computer readable storage medium having stored thereon instructions for identifying, by a method, a fluorochrome panel suitable for use in a flow cytometry protocol.

[0270] Clause 104. The non-transitory computer-readable storage medium of Clause 103, wherein the numerical stability index is a condition number.

[0271] Addendum 105. The non-transitory computer-readable storage medium of Addendum 104, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0272] Clause 106. The non-transitory computer-readable storage medium of clause 105, wherein when calculating the condition number, a ratio of a maximum singular value to a minimum singular value is determined.

[0273] Addendum 107. The non-transitory computer-readable storage medium of Addendum 106, wherein when calculating the condition number, for each of two or more subsets of N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of N fluorescent dye identifiers is determined.

[0274] Addendum 108. When calculating the condition number, for each of two or more subsets of N fluorescent dye identifiers, the Frobenius norm of the spectral matrix associated with a subset of N fluorochrome identifiers, and the Frobenius norm of the pseudoinverse of each spectral matrix associated with two or more subsets of N fluorochrome identifiers 108. The non-transitory computer-readable storage medium of claim 107, further comprising:

[0275] Addendum 109. The non-transitory computer-readable storage medium of Addendum 108, further comprising, when calculating the condition number, determining, for each of two or more subsets of the N fluorescent dye identifiers, a product of the Frobenius norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of the N fluorescent dye identifiers.

[0276] Addendum 110. A non-transitory computer-readable storage medium according to any one of Addendums 103 to 109, wherein the fluorescent dye panel includes an optimized numerical stability index.

[0277] Addendum 111. The non-transitory computer-readable storage medium of Addendum 110, wherein the identifying the fluorochrome panel uses a constrained optimization method.

[0278] Clause 112. The non-transitory computer-readable storage medium of Clause 111, wherein the constrained optimization method is a minimization algorithm.

[0279] Addendum 113. The non-transitory computer-readable storage medium of any one of Addendums 111 to 112, for sequentially executing a constrained optimization method.

[0280] Addendum 114. A non-transitory computer-readable storage medium as described in Addendum 113, wherein when selecting two or more subsets of N fluorescent dye identifiers, a first subset of N fluorescent dye identifiers and a second subset of N fluorescent dye identifiers are selected.

[0281] Addendum 115. The non-transitory computer-readable storage medium of Addendum 114, wherein a first subset of N fluorescent dye identifiers is randomly obtained from a spectral matrix associated with the instrument identifier.

[0282] Addendum 116. A non-transitory computer-readable storage medium as described in Addendum 114 or 115, wherein the constrained optimization method further calculates a numerical stability index associated with a first subset of the selected N fluorescent dye identifiers.

[0283] Addendum 117. A non-transitory computer-readable storage medium according to any one of Addendums 114 to 116, wherein when generating a second subset of N fluorescent dye identifiers, fluorescent dye identifiers in the first subset of N fluorescent dye identifiers are replaced with alternative fluorescent dye identifiers in a spectral matrix associated with the instrument identifier.

[0284] Addendum 118. The non-transitory computer-readable storage medium of Addendum 116 or 117, wherein the constrained optimization method calculates a numerical stability index associated with a second subset of the N fluorescent dye identifiers.

[0285] Addendum 119. The non-transitory computer-readable storage medium of Addendum 118, wherein the constrained optimization method compares a numerical stability metric associated with a first subset of the N fluorescent dye identifiers to a numerical stability metric associated with a second subset of the N fluorescent dye identifiers.

[0286] Addendum 120. A non-transitory computer-readable storage medium according to any one of Addendums 113 to 119, wherein the constrained optimization method further comprises iteratively generating a subset of fluorescent dye identifiers.

[0287] Appendix 121. In constrained optimization, Calculate a numerical stability index for each subset of N fluorescent dye identifiers generated iteratively; A non-transitory computer-readable storage medium as described in Appendix 120, which compares a numerical stability index calculated for each subset of N fluorescent dye identifiers generated iteratively with a numerical stability index of a subset of N fluorescent dye identifiers that have already been generated, and determines an optimized numerical stability index.

[0288] Addendum 122. The non-transitory computer-readable storage medium of any one of Addendums 110 to 121, wherein the optimized numerical stability metric is a minimized condition number.

[0289] Addendum 123. The non-transitory computer-readable storage medium of any one of Addendums 103 to 122, wherein the spectral feature is a spillover value.

[0290] Addendum 124. The non-transitory computer-readable storage medium of any one of Addendums 103 to 123, wherein the spectral features are obtained from experimental data.

[0291] Addendum 125. The non-transitory computer-readable storage medium of any one of Addendums 103 to 123, wherein the spectral features are obtained from simulation data.

[0292] Addendum 126. A non-transitory computer-readable storage medium according to any one of Addendums 103 to 125, wherein the method is performed without evaluating antigen data.

[0293] The non-transitory computer-readable storage medium of any one of claims 103 to 126, wherein N is an integer in the range of 2 to 100.

[0294] 128. The non-transitory computer-readable storage medium of claim 127, wherein N is an integer in the range of 2 to 75.

[0295] Addendum 129. A non-transitory computer-readable storage medium according to any one of Addendums 103 to 128, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 150.

[0296] Addendum 130. The non-transitory computer-readable storage medium of Addendum 129, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 100.

[0297] Addendum 131. A non-transitory computer-readable storage medium according to any one of Addendums 103 to 130, wherein the method further comprises estimating the maximum number of fluorescent dyes that may be used simultaneously in the fluorescent dye panel.

[0298] Addendum 132. A non-transitory computer-readable storage medium as described in Addendum 131, determining an optimized numerical stability index for each of multiple values ​​of N when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0299] The non-transitory computer-readable storage medium of claim 132, wherein the plurality of values ​​of N includes i values, where i is the number of fluorescent dye identifiers in the spectral matrix associated with the instrument identifier.

[0300] Addendum 134. A non-transitory computer-readable storage medium according to Addendum 132 or 133, which determines a maximum value of N that does not result in an optimized numerical stability index exceeding a threshold when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0301] Addendum 135. A non-transitory computer-readable storage medium according to any one of Addendums 103 to 134, wherein the spectral matrix associated with the device identifier includes autofluorescence spectral features.

[0302] Addendum 136. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: (a) inputting into a processor an instrument identifier and a request for N fluorochrome identifiers, where N is an integer of at least 2; The processor selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with an instrument identifier, each fluorochrome identifier in the spectral matrix being associated with a spectral feature; Identifying a fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets. It is structured as follows: (b) receiving the identified fluorochrome panel from the processor.

[0303] Addendum 137. The method of Addendum 136, wherein the numerical stability index is a condition number.

[0304] Addendum 138. The method of Addendum 137, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

[0305] Addendum 139. The method of Addendum 138, wherein when calculating the condition number, the ratio of the largest singular value to the smallest singular value is determined.

[0306] Addendum 140. The method of Addendum 137, wherein when calculating the condition number, for each of two or more subsets of N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of N fluorescent dye identifiers is determined.

[0307] Addendum 141. When calculating the condition number, for each of two or more subsets of N fluorochrome identifiers, the Frobenius norm of the spectral matrix associated with a subset of N fluorochrome identifiers, and the Frobenius norm of the pseudoinverse of each spectral matrix associated with two or more subsets of N fluorochrome identifiers 141. The method of claim 140, further comprising determining

[0308] Addendum 142. The method of claim 141, further comprising, when calculating the condition number, determining, for each of two or more subsets of N fluorescent dye identifiers, the product of the Frobenius norm of a spectral matrix associated with the subset of N fluorescent dye identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of N fluorescent dye identifiers.

[0309] Addendum 143. The method of any one of Addendums 136 to 142, wherein the fluorescent dye panel includes an optimized numerical stability index.

[0310] Addendum 144. The method of Addendum 143, wherein a constrained optimization method is used to identify the fluorescent dye panel.

[0311] Addendum 145. The method of Addendum 144, wherein the constrained optimization method is a minimization algorithm.

[0312] Addendum 146. The method of Addendum 144 or 145, wherein the constrained optimization method is carried out sequentially.

[0313] Addendum 147. The method of Addendum 146, wherein when selecting two or more subsets of N fluorescent dye identifiers, a first subset of N fluorescent dye identifiers and a second subset of N fluorescent dye identifiers are selected.

[0314] Addendum 148. The method of Addendum 147, wherein a first subset of N fluorescent dye identifiers is randomly obtained from a spectral matrix associated with the instrument identifier.

[0315] Addendum 149. A method as described in Addendum 147 or 148, wherein the constrained optimization method further calculates a numerical stability index associated with a first subset of the selected N fluorescent dye identifiers.

[0316] Addendum 150. A method according to any one of Addendums 147 to 149, wherein when generating the second subset of N fluorescent dye identifiers, fluorescent dye identifiers in the first subset of N fluorescent dye identifiers are replaced with alternative fluorescent dye identifiers in a spectral matrix associated with the instrument identifier.

[0317] Addendum 151. A method according to Addendum 149 or 150, wherein the constrained optimization method calculates a numerical stability index associated with a second subset of the N fluorescent dye identifiers.

[0318] Addendum 152. The method of Addendum 151, wherein the constrained optimization method compares a numerical stability metric associated with a first subset of the N fluorescent dye identifiers to a numerical stability metric associated with a second subset of the N fluorescent dye identifiers.

[0319] Addendum 153. A method according to any one of Addendums 146 to 152, wherein the constrained optimization method further comprises iteratively generating a subset of fluorescent dye identifiers.

[0320] Appendix 154. In constrained optimization, Calculate a numerical stability index for each subset of N fluorescent dye identifiers generated iteratively; The method described in Appendix 153, further comprising comparing a numerical stability index calculated for each subset of N fluorescent dye identifiers generated iteratively with a numerical stability index of a subset of N fluorescent dye identifiers already generated, and determining an optimized numerical stability index.

[0321] Addendum 155. The method of any one of Addendums 143 to 154, wherein the optimized numerical stability index is a minimized condition number.

[0322] Addendum 156. The method of any one of Addendums 136 to 155, wherein the spectral feature is a spillover value.

[0323] Addendum 157. The method of any one of Addendums 136 to 156, wherein the spectral features are obtained from experimental data.

[0324] Addendum 158. The method of any one of Addendums 136 to 156, wherein the spectral features are obtained from simulated data.

[0325] Addendum 159. The method of any one of Addendums 136-158, wherein the processor is configured to identify the fluorochrome panel without evaluating antigen data.

[0326] Addendum 160. The method of any one of Addendums 136-159, wherein N is an integer in the range of 2-100.

[0327] Addendum 161. The method of Addendum 160, wherein N is an integer in the range of 2 to 75.

[0328] Addendum 162. A method according to any one of Addendums 136 to 161, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 150.

[0329] Addendum 163. The method of Addendum 162, wherein the spectral matrix associated with the instrument identifier includes a number of fluorescent dye identifiers in the range of 2 to 100.

[0330] Addendum 164. The method of any one of Addendums 136 to 163, wherein the processor is configured to estimate the maximum number of fluorescent dyes that may be used simultaneously in the fluorescent dye panel.

[0331] Addendum 165. The method of Addendum 164, wherein an optimized numerical stability index is determined for each of multiple values ​​of N when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel.

[0332] The method of claim 165, wherein the plurality of values ​​of N includes i values, where i is the number of fluorescent dye identifiers in the spectral matrix associated with the instrument identifier.

[0333] Addendum 167. The method of Addendum 165 or 166, in which, when evaluating the maximum number of fluorescent dyes that may be used simultaneously in a fluorescent dye panel, a maximum value of N is determined that does not result in an optimized numerical stability index exceeding a threshold.

[0334] Addendum 168. The method of any one of Addendums 136 to 167, wherein the spectral matrix associated with the device identifier includes autofluorescence spectral features.

[0335] Although the foregoing invention has been described in some detail by way of illustration and example for ease of understanding, it will be readily apparent to those skilled in the art, in view of the teachings of the invention, that certain changes and modifications may be made thereto without departing from the spirit and scope of the appended claims.

[0336] Thus, the foregoing merely illustrates the nature of the invention. It is clear that those skilled in the art may devise various configurations that embody the nature of the invention and are within the spirit and scope of the invention, although not expressly described or shown herein. Furthermore, all examples and conditional language described herein are intended essentially to aid the reader in understanding the nature of the invention and the concepts given by the inventors to advance the art, and should not be construed as being limited to the examples and conditions specifically described as such. Furthermore, all descriptions herein that describe the nature, aspects, and embodiments of the invention, as well as specific examples of the invention, are intended to include both structural and functional equivalents of the invention. In addition, such equivalents are intended to include both currently known equivalents and equivalents developed in the future, i.e., all elements developed that perform the same function, regardless of structure. Furthermore, the descriptions disclosed herein are not intended to be public only, regardless of whether such disclosure is expressly set forth in the claims.

[0337] Therefore, it is not intended that the scope of the present invention be limited to the exemplary embodiments shown and described herein. Rather, the scope and spirit of the present invention is embodied in the appended claims. With respect to claims, 35 U.S.C. 112(f) or 35 U.S.C. 112(6) is expressly provided to be invoked with respect to a limitation in a claim only when the precise phrase "means for" or "step for" is recited at the beginning of such limitation in the claim, and if such precise phrase is not used in the limitation of the claim, 35 U.S.C. 112(f) or 35 U.S.C. 112(6) is not invoked.

[0338] CROSS-REFERENCE TO RELATED APPLICATIONS This application claims priority under 35 U.S.C. §119(e) of the filing date of U.S. Provisional Application No. 63 / 305,010, filed January 31, 2022, the entire disclosure of which is incorporated herein by reference.

Claims

1. 1. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2; selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with said instrument identifier, each fluorochrome identifier in said spectral matrix being associated with a spectral feature; A method for identifying the fluorescent dye panel from two or more subsets of N fluorescent dye identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorescent dye identifiers using the spectral features associated with the N fluorescent dye identifiers of the two or more subsets.

2. The method of claim 1 , wherein the numerical stability index is a condition number.

3. 3. The method of claim 2, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a singular value of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

4. The method of claim 3 , further comprising determining a ratio of a maximum singular value to a minimum singular value when calculating the condition number.

5. 3. The method of claim 2, wherein when calculating the condition number, for each of two or more subsets of the N fluorescent dye identifiers, a matrix norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers is determined.

6. When calculating the condition number, for each of two or more subsets of the N fluorochrome identifiers: the Frobenius norm of a spectral matrix associated with a subset of the N fluorochrome identifiers; and the Frobenius norm of the pseudoinverse of each of the spectral matrices associated with two or more subsets of the N fluorochrome identifiers The method of claim 5, further comprising determining

7. 7. The method of claim 6, further comprising determining, for each of two or more subsets of the N fluorescent dye identifiers, the product of the Frobenius norm of a spectral matrix associated with the subset of the N fluorescent dye identifiers and the Frobenius norm of a pseudoinverse of the spectral matrix associated with the subset of the N fluorescent dye identifiers.

8. The method of any one of claims 1 to 7, wherein the fluorescent dye panel comprises an optimized numerical stability index.

9. The method according to any one of claims 1 to 8, wherein said spectral features are spillover values.

10. The method according to any one of claims 1 to 9, wherein said spectral features are obtained from experimental or simulated data.

11. The method according to any one of claims 1 to 10, performed without evaluation of antigen data.

12. The method according to any one of claims 1 to 11, comprising assessing the maximum number of fluorochromes that may be used simultaneously in the fluorochrome panel.

13. 1. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: receiving a request for N fluorochrome identifiers, where N is an integer of at least 2; A method comprising identifying a fluorochrome panel comprising a subset of N fluorochrome identifiers, and performing said method without evaluating antigen data.

35. receiving a device identifier; 35. The method of claim 34, further comprising selecting two or more subsets of N fluorescent dye identifiers from a set of fluorescent dye identifiers in a spectral matrix associated with the instrument identifier, each fluorescent dye identifier in the spectral matrix being associated with a spectral feature.

14. receiving a request for an instrument identifier and N fluorochrome identifiers, where N is an integer of at least 2; selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with said instrument identifier, each fluorochrome identifier in said spectral matrix being associated with a spectral feature; identifying a fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using the spectral features associated with the N fluorochrome identifiers of the two or more subsets; 16. A system comprising: a processor configured to:

15. 1. A method for identifying a fluorochrome panel suitable for use in a flow cytometry protocol, comprising: (a) inputting into a processor an instrument identifier and a request for N fluorochrome identifiers, where N is an integer of at least 2; The processor, selecting two or more subsets of N fluorochrome identifiers from a set of fluorochrome identifiers in a spectral matrix associated with said instrument identifier, each fluorochrome identifier in said spectral matrix being associated with a spectral feature; identifying the fluorochrome panel from the two or more subsets of N fluorochrome identifiers based on a numerical stability index calculated for each of the two or more subsets of N fluorochrome identifiers using spectral features associated with the N fluorochrome identifiers of the two or more subsets; It is structured as follows: (b) receiving the identified fluorochrome panel from the processor.

Images