A gain and dynamic range independent index to quantify spillover spread to aid panel design in flow cytometry.

In conventional flowcytometry one detector (primary) is dedicated for one fluorochrome. However, photons usually end up in other detectors too (fluorescence spillover). 'Compensation' is a process that corrects the spillover signal from all detectors except the primary detector. Post 'compensation', the photon counting error of spillover signals become evident as spreading of the data. The spreading induced by spillover impairs the ability to resolve stained cell population from the unstained one, potentially reducing or completely losing cell populations. For successful multi-color panel design, it is important to know the expected spillover to maximize the data resolution. The Spillover Spreading Matrix (SSM) can be used to estimate the spread, but the outcome is dependent on detector sensitivity. Simply, the same single stained sample produces different spillover spread values when detector(s) sensitivity is altered. Many researchers mistakenly use this artifact to "reduce" the spread by decreasing detector sensitivity. This can result in diminished capacity to resolve dimly expressing cell populations. Here, we introduce SQI (Spread Quantification Index), that can quantify the spillover spread independent of detector sensitivity and independent of dynamic range. This allows users to compare spillover spread between instruments having different types of detectors, which is not possible using SSM.

Introduction

Flow Cytometry has improved greatly since it was first introduced, evolving into a powerful technique to monitor and study huge number of cells within a few minutes. Modern machines have excellent sensitivity both on Fluorescence and scatter[1,2]. In conventional flow cytometers, each detector should ideally receive photons from one specific fluorochrome as a labeled cell travels through the laser interrogation point. However, fluorescent emission usually ends up in different detectors too[3]. The extent of unwanted signals depends on the characteristics of the band pass (BP) associated with the secondary detectors and the shape/nature of the emission spectra. The details and rules of rectifying this unwanted signal in secondary detectors, compensation, can be found extensively explained by others[4-8]. This unwanted signal is spillover. Spillover reduces the sensitivity, the lowest signal that can be detected over accumulation of signals coming from all other sources. In Flow Cytometers, spillover causes a non-linear spreading of properly compensated data (appears as straight line in log–log plot)[4]. Higher spillover causes higher spread in the secondary detector, this spread can mask dim signals. For Immunophenotypic panels, spreading can decrease the detection sensitivity for low/dimly expressing antigens.

Detectors like Photo Multiplier Tubes (PMTs) have been a stable choice for many manufacturers over many years. These detectors are highly efficient at converting the light/incoming photons into photo electrons for detection. New types of detectors, such as Avalanche Photo Diodes (APDs) and Silicon Photo Multipliers (SiPMs), are now gaining popularity among manufactures, as they are as good as PMTs, but more efficient in converting longer wavelength photons than PMTs[9,10]. The photo electron generation process can be controlled by ‘PMT voltage’ for PMTs or by gain for APDs and SiPMs. Increasing PMT voltage or gain may increases the detectors sensitivity, which may mean better detection of dim signals, but can also increase detector noise if increased too much, which negatively impacts data resolution[11]. The best way to avoid this problem is to find a combination of bright fluorochromes with dimly expressing antigens.

With passing time, the complexity of using multiple fluorochromes has increased a lot. Recently, scientists have shown 31 fluorochromes can be used together for conventional systems[12] and over 40 fluorochromes can be used in full spectral flow cytometers[3,13]. With higher number of fluorochromes combined, predicting, and bypassing the possible spillover spread becomes more important than ever. Full spectral cytometers allow researchers to assemble bigger panels much cheaper than before, but they are not immune to spillover spreading issue. For designing a good panel, it is key to know how much spread can be expected in advance.

Useful tools/methods have been developed to aid investigators in this process. The most well-known and widely used is the Spillover Spreading Matrix (SSM), which provides a metric for signal variance introduced by spillover fluorescence for every fluorochrome-detector combination in a given panel[14]. The SSM can provide valuable insight into which combinations of fluorochromes are expected to work well together in a panel. Each value of the SSM is called Intrinsic Spillover Spread (ISS). However, care should be taken when using this method. We have found that ISS values are dependent on the detector sensitivity (detector voltage/gain). For example, if users decrease the ‘PMT Voltage’ or gain of all the detectors, then ISS values drop for all detectors. This phenomenon can be used by the researcher to ‘reduce spread’ by lowering the voltage/gain (instead of redesigning the panel) while lower sensitivity (due to lowering the detector sensitivity) can lead to reduced resolution between positive and negative populations. The effect of detector voltage/gain is intrinsically present in SSM and cannot be removed. It is also possible to use SSM to compare spread between machines (only between same make and models), but the process is tedious[14].

The output of the detectors is displayed using a scale, which can vary depending on the vendor and model. Aria Fusion from BD has 262,144 scale units (the scale is also known as dynamic range), CytoFLEX from Beckman Coulter has 16,777,215, Aurora from Cytek has 4,194,304 and Quanteon from ACEA has 107.2 scale units respectively. In this article, we have introduced a new method to quantify spread, named Spread Quantification Index (SQI). For a given set of samples, SQI is independent of detector voltage/gain as well as independent of the dynamic range of the detector. It is easy to implement; researchers only need to carefully perform compensation by following all the rules to get a reproducible SQI value. SQI works correctly when the data is within the detector’s linearity range. Spreads from different machines independent of make/model/dynamic range or vendor can be compared directly using SQI.

Results

Concept of SQI

The SQI value is defined as the difference in median fluorescent intensity (MFI) between the 99th percentile and the 50th percentile in detector B, of the positive population of fluorochrome X (Fig. 1). Detector B is the recipient detector where spread from fluorochrome X is being registered. This value is a measure of total spread including spillover fluorescence in a secondary detector. Some examples of SQI can be found in Fig. 2. The following abbreviations are used to describe the calculation procedure as explained in Fig. 1.

Figure 1

Concept of the SQI: The figure shows the concept of the SQI value. Spillover fluorescence from primary fluorochrome X results in spread in detector B, and vice versa. The spread is quantified as the difference of the 99th percentile and 50th percentile in dimension B. Subsequently, this value is multiplied by a normalization factor to make the SQI value independent of detector gain settings. The 99th percentile is marked by a Green line.

Figure 2

Effect of PMT voltage on ISS: In (A), the primary detector (V605) voltage was changed from 600 to 400 V in steps of 100 V. This did not change the rSD and SQI values measuring the spread in the secondary detector (V677). However, the ISS values changed noticeably. In (A–C), voltages of both detectors were changed simultaneously to 600 V, 500 V or 400 V respectively. The SQI values remained the same for both detectors in all cases, but the ISS values decreased considerably. In (D) only voltage of V677 was lowered which resulted in a huge increase of ISS but SQI remain same.

XB50 = 50th percentile of compensated positive population of fluorochrome X in detector B.

XB99 = 99th percentile of compensated positive population of fluorochrome X in detector B.

UnB50 = 50th percentile of compensated unstained population of in detector B.

YStB50 = 50th percentile of compensated stained population of fluorochrome Y in detector B.

NFYB = Normalization factor of fluorochrome Y in detector B (Normalization Factor, in Fig. 1).

(YStB50 – UnB50) = The difference between the positive and negative populations (scale units) (Difference 2, D in Fig. 1).

If P (any constant number, in this article we have used P = 9000) is the number of photons received by the detector B and (YStB50 – UnB50) are the number of scale units of detector B that represent P.

Then (YStB50–UnB50) ≡ P, but not equal.

1 scale unit in detector B = P/(YStB50 – UnB50) = NFYB

If SQI value XB = Spread of fluorochrome X in detector B,

Then SQI XB = (XB99ber of photons received by the detecto – XB50) × NFYB.

By applying the stated normalization factor, the SQI value is made independent from voltage/gain settings of the detectors, and dynamic range. The normalization factor is founded on two assumptions. First, the number of incoming photons remains the same when settings of the secondary detector are changed. This means that changing the detector voltage/gain will not change the excitation of the fluorochrome by the laser and the number of emitting photons will remain same when they pass through the same filters and mirrors. Second, differences in signal intensities are linearly proportional to changes in detector settings. This means that the conversion of photons into photo electrons is constant. In one gain setting, 1 photon converts into 3 photoelectron and 5 photons into 15 photoelectrons, while in another gain setting, if 1 photon converts to 5 photoelectrons, then 5 photons must convert into 25 photoelectrons.

The normalization factor allows for the SQI value to be determined without prior optimization of detector settings, provided that all obtained data is within the linear range and negative values are outside the detector noise. P is the difference in number of photons that are coming from positive and negative population in detector B. This number will stay practically constant when saturatedly stained samples will be maximally excited. The lasers of almost all of the presently available cytometers are designed to get maximal excitation. For this reason, we can get the same P values in all commercial cytometers. Maximal excitation makes P practically a constant. The constant nature of P in these terms allows for direct comparison of SQI values between instruments with different dynamic ranges. SQI showed a high degree of reproducibility in an average of three runs on the Fusion (Supplementary Fig. S1). In the discussion, we elaborate on the conditions that must be met. Importantly, when multiple fluorochromes are associated with the same detector, the same amount of spread results in a higher SQI in secondary detectors when associated with dim fluorochromes (Supplementary Table S1).

Classification of SQI values

Since there are no clear consensus or characteristics that divide spread in categories[14,15], we created a different, semi-quantitative method to categorize SQI outcome as described below:

For example, 262,144 is the total number of scale units for BD machines (Fusion and Symphony). SQI between 0 to 120, 121 to 220, 221 to 300 and 300 + matches nicely with 0 to 0.5%, 0.5% to 1%, 1% to 1.5% and 1.5% + of scale unit occupancy respectively. Comparing percent of scale units occupied with SQI will provide the user with a better understanding of spread between different types of machines (Table 1) and between multiple machines of the same type (Table 2). We have color coded SQI values in 4 categories (Tables 1 and 2): 1–120 (Green); 121–220 (Yellow); 221–300 (Orange); and 300 + (Red). While Green and yellow both provide acceptable discriminative power in most cases, Green should be reserved for markers with low or unknown expression levels. Values with an Orange or Red SQI value should be used for markers on distinctly separated cell populations or for markers with exclusive expression on specific cell subsets. For example, if BV421 is combined with Alexa Fluor 647 (AF647), the SQI value is 20 (Green). Note that the SQI value is not zero due to measurement error of the autofluorescent signal and biological variance. In comparison, if PE-Cy5 is combined with AF647 the SQI value is 611 (Red). At these different levels of spread, 15% and 5% double positive events were identified respectively (Supplementary Fig. S2). Double positive events are those which are outside the spreading error, received by the secondary/recipient detector (R670). The lower number of double positive events registered for the PE Cy5 compare to BV421 because of the need to account the higher spreading error coming from the PE Cy5.

Table 1

Comparison of spread in four different flow cytometers: ten different single stained bead sets were run on four different instruments.

SQI was calculated for every combination. For BV711 there was no data available for the Cytoflex S. The SQI values are categorized as follows: 1–120 (Green); 121–199 (Yellow); 200–299 (Orange); and 300 + (Red).

Table 2

Comparison of all four Fusion setup: the SQI values for 10 fluorochromes run on three different Fusions are shown.

For Fusion 2, “2_High” denotes higher laser power as for the other two instruments. The SQI values are categorized as follows: 1–120 (Green); 121–199 (Yellow); 200–299 (Orange); and 300 + (Red).

Intrinsic Spillover Spread values are dependent on detector voltage/gain settings

The ISS value of BV605 changed in the recipient detector when voltage of the primary detector changed, because the intensity difference (ΔF) between the positive and negative populations incorrectly describes a change in photonic input in the secondary detector (Fig. 2A). Importantly, the variance as represented by the rSD and SQI value remained the same. When voltages of both the detector were reduced, the rSD and ISS values were lower too (Fig. 2A–C). Conversely, the SQI value was unchanged. In Fig. 2D we only reduced the PMT voltage of the secondary detector which decreased the ΔF of the secondary detector resulting in increasing the related ISS.

SQI values are independent of detector gain settings

Our main goal is to quantify spread independent of detector voltage/gain. The normalization factor was introduced for this reason. We devised this simple experiment to test the accuracy of the normalization factor. Single stained beads were measured with two different PMT voltage sets. The voltage was manipulated to position the positive population to ~ 110,000 MFI or 50,000 MFI in all detectors. The SQI values were the same in both voltage sets, but the ISS values were quite different (Table 3). This observation is possible because SQI does not include any component from the primary detector. When the voltage of detector B was decreased, the spread from fluorochrome X in detector B was reduced, but the intensity of fluorochrome Y, which was used to calculate normalization factor, reduced proportionally, rendering SQI values independent of detector voltage. In contrast, when the same data were used to obtain the conventional Spillover Spread Matrices, the lower voltage set showed lower ISS values (Table 3). In another scenario, a series of voltages were applied in which the voltage of the secondary detector was kept constant, but the voltage of the primary detector was changed gradually (Fig. 3A,B). Alternatively, the voltage of the primary detector was kept constant, and the voltage of the secondary detector was changed (Fig. 3C,D). In both cases the SQI values remained the same, however, the ISS values changed considerably.

Table 3

PMT voltage independent measure of spread: two different sets of PMT voltages were applied to set the MFI of the positive population to 110,000 (110 K) or 50,000 (50 K).

ISS values are shown in Italics.

Six single stained bead samples were run using these two voltage sets, and the ISS and SQI values were calculated. The SQI values were similar at the different voltages, whereas the ISS values were voltage different. ISS values are shown in Italics.

Figure 3

Effect of detector settings on SQI and ISS: In (A, B) we have only changed the PMT voltage for the primary detector. SQI remains same but ISS dropped with increase of voltage. In (C, D), we only changed PMT voltage for secondary detector. SQI again remain same but ISS increased with voltage.

Application of SQI: compare instruments with different dynamic range

To investigate which instrument is the best for running a given panel in terms of spread, we compared the SQI values between a Fusion, Symphony, Quanteon and a CytoFLEX (Table 1). The CytoFLEX showed the lowest number of SQI values above 120. The other instruments showed a higher number of values above 120. The Quanteon and CytoFLEX had the lowest number of SQI values over 300. The Quanteon had the highest number of SQI values between 120 to 199.

Application of SQI: compare instrument performance

While all three Fusions had similar optical configuration and laser power, Fusion 3 performed best with only 6 instances of 300 + range. The performance of The Fusion 1 was considered worst with 12 occurrences of Range values over 300. The findings are summarized in Table 2. The results suggest that Fusion 3 is more suitable for the most sensitive experiments.

Application of SQI: investigate the effect of instrument characteristics on spread

SQI was used to investigate the effect of high and low laser powers in Fusion 2 which showed some interesting findings. At an output power of 140 mW, the Red laser resulted in lower spillover spread from APC and APC-R700 compared to lower output power of 100 mW (Table 2). In other cases, higher laser power showed minimal effect on spread. For instance, at 100 mW, the 405 nm laser showed virtually the same spread in all three violet detectors as compared to measurement at 85 mW (Supplementary Table 2). Our data indicates that in our systems, higher laser powers other than for the red laser, doesn’t lower the data spread (Table 2).

Validation assay

To test the usefulness of SQI, we ran two small panels separately. The experiment was performed on Cytek Aurora, a full spectral cytometer. Antibodies in both panels have identical clones and the exact same amount of antibody was used. In Panel 2, we replaced CD3-BV605, CD45-Alexa Fluor 488 with CD3-APC-R700 and CD45-BV650. We did this change to introduce more spread in the APC and PerCP-eFluor 710 detectors. We also choose BUV395, a dim fluorochrome with CD25 (low expression) to show clearly the adverse effect of lowering the gain. Both panels were run four times separately using four different gain settings (details are in material and methods). When analyzing for specific populations (Supplementary Fig. S3), ISS and SQI calculated using beads to show with change of gain, SQI remains the same but ISS drops (Supplementary Fig. S4). Further, reducing gain resulted in the loss of data resolution of the specific populations, marked by the black arrows in panel 1, and spread-induced masking of the specific populations, marked by the red arrows in panel 2 (Supplementary Fig. S4). These data, combined with the single color SQI with these same antibodies (Supplementary Table S4), demonstrate that BV650 and APC-R700 would produce much larger spread in APC and PerCP-eFluor 710 detectors resulting in masking of dimly expressing cells. This data confirms that SQI works nicely with full spectral cytometers as well as conventional cytometers.

Discussion

To show how ISS values are dependent on detector voltage/gain, we explored two scenarios. First, the voltage of the primary detector was changed (Fig. 2A). The spread remained the same as was evident by the unchanged rSD, but the ISS values changed because the value of √ΔF changed. More specifically, the lower voltage on the primary detector increased the ISS value, because the same ∆σC [Eq. (1) of Ref.[14]] was divided by a smaller number (smaller ΔF). In the Fig. 3, we have shown SQI remains unchanged if we change primary or secondary voltages. However, the voltage changes do affect the ISS values considerably. In the second scenario (Fig. 2A–C), voltages of the secondary detectors were also changed. In this example, lower voltages resulted in lower ISS values. Unfortunately, sometimes users do “optimized” the PMT voltages to get low ISS values. The motivation behind this is to reduce the spread so that they will be able to identify the double positive population better. This, however, may result in loss of resolution and obscuration of cells that have low expression of markers. Rather, users should adjust their panel to a combinations that have low spread, increasing the chances of getting a good separation. Optimization of detector voltage/gain should be achieved by monitoring the separation of signal over the detector noise[11,16]. SQI is not designed nor recommended for voltage/gain optimization. The SQI value is independent of detector voltage/gain (bias induced by voltage is prevented). The inherit relation of ISS and detector voltage/gain will always make the ISS values subjective. Even when users are not comparing different voltage sets, the influence of detector voltage/gain will always be there. SQI is, to the best of our knowledge, currently the only way to remove this subjectivity.

Furthermore, because variance and MFI are proportional to amplification settings in equal fashion, the MFI can be used to correct for the effect of changing secondary detector settings on spread. Additionally, the MFI of the unstained population (UnB50) is subtracted from the MFI of the Stained population (YStB50) to normalize for differences in autofluorescence between channels. As shown in Table 3, lowering both voltages lowered the ISS values, with same samples. As SQI is independent of voltage, the spread quantitation will remain the same. Ashhurst et al. also showed that spread do not change with voltage alteration[8]. In Fig. 2D we showed this by altering the PMT voltage of secondary detector which resulted in change of ISS but SQI remain same.

We presented a special scenario where the SQI value of one fluorochrome was calculated for two different fluorochromes, used separately in the secondary/recipient detector (Supplementary Table S1). For example, the SQI value of APC in detector B695 is different when different secondary fluorochromes (BB700 or PerCP-Cy5.5) are considered. Because the SQI value is relative to the intensity of the secondary fluorochrome, the same spread from APC results in a relatively higher occupancy in scale units when combined with PerCP-Cy5.5 as compared to BB700. It is helpful in the panel design phase to know that APC with PerCP-Cy5.5 may be a more problematic combination for the identification of a double positive population. The calculation for ISS does not consider the MFI of the positive population from BB700 or PerCP-Cy5.5. This makes the ISS values independent of the fluorochromes associated with the secondary detector. This will result in the spread in terms of ISS remains the same for both BB700 and PerCP-Cy5.5. While this specific fluorochrome was not used in the validation panels, the concept can be seen in that panel. Specifically, when looking at CD3 vs TCRδγ, where the double positive population is much easier to identify with less spread in panel 1, when compared to panel 2. In addition, when gain is increased compared to the manufacturer suggested values, there is a substantial increase in a single positive population that is not reported in the literature (Supplementary Fig. S4). These interactions are also important to other fluorochromes within that staining panel, even with ideally designed panels. For example, even with panels designed to have little to no spread, gain values are still important, as decreasing the gain did not result in large changes in the above double positive population, but the data revealed a loss in a rare CD25+ CD127− population (Supplementary Fig. S4), indicating that dim fluorochromes on antibodies for low expressing surface proteins are more sensitive to changes in detector sensitivity. The SQI and ISS values PE in BUV395 detector is not available for lower gain setups because the XB99 and XB50 values are outside the detector linear range which makes the calculation invalid. We deliberately introduced BV650 and APC-R700 in panel 2 to introduce spread in APC and PerCP-eFluor-710 detectors (Supplementary Table S4). In panel 1, CD3− CD4+ dim population lost the separation in lower gain value but because of the introduction of spread from BV650 and APC-R700 in APC channel becomes impossible to separate in panel 2. CD3+ TCRγδ++ population is visible in both panels but due to the added spread CD3 + TCRγδ + is not clearly visible in panel 2 (Red arrow) but in panel 1 (Black arrow). We have found an interesting population in panel 1 (CD3− TCRγδ +) which gradually disappears with decrease in gain. To best of our knowledge, this is an unreported population. Because of spread introduced in panel 2, this population become unidentifiable. Importantly, these data (all SQI data for various machines in Tables 1, 2, 3 and Supplementary Table S4 and Fig. S4) demonstrate that SQI is a platform independent method that will aid in better panel design.

It is not easy to compare ISS values between systems with different dynamic range, because instruments need to be calibrated first[14]. The ΔF accounts for the dependency of spread on the number of fluorescent photons, the ∆F needs to be numerically similar on both instruments to have a meaningful comparison of ISS values. This calibration is a tedious process. However, because the intensity values in the primary detector don’t influence the SQI values, this is not necessary in our method. Furthermore, due to the normalization factor, the scaling on every secondary detector is comparable. This makes the SQI value well suited to compare panel performance on different instruments.

As proof of concept, we compared the SQI values for a given set of single stains on different instruments with different dynamic range. Due to higher quantum efficiency for longer wavelength fluorescent light, we anticipated to find less spread in APD and SiPM based instruments in the red emission region as compared to PMT based instruments[9,10]. Indeed, both APD and SiPM detector-based instruments showed lower SQI values in detectors for APC and APC-R700 (Table 1).

We matched the other instrument characteristics as best as possible, but these were not the same. For example, laser output powers of the Quanteon and the Symphony were matched, but there was no option to alter the laser power for the CytoFLEX. Furthermore, the filter configuration (characteristics and positioning) between the instruments was also not the same (Supplementary Table S2). However, this example shows that it is possible to use SQI values for the purpose of comparing instrument performance with regards to spread for a given set of fluorochromes.

SQI values are empirical and very robust because it’s not a variance-based calculation of spread. For a Gaussian distribution, the 99th percentile can be calculated as 2.325 times the standard deviation (SD) of the mean. However, flow cytometric data follows the Poisson distribution and is known for positively skewed distributions[17]. For this reason, the factor of rSD required to calculate the 99th percentile for each detector is highly variable. The median of 99th percentile was obtained using FlowJo statistics tools. However, since the determination of the 99th percentile is effectively based on just 1% of the data, it could potentially be distorted by outliers. To address this issue, careful singlet selection was performed, the median was used instead of the mean as a measure of central tendency, and a minimum of 35,000 single beads were recorded. This practice has been shown to provide highly robust results with less than 5% variation between tests (Supplementary Fig. S1).

While the SQI value is easy to use, there are some important things to keep in mind. First, the precision of the SQI value depends on correct ‘compensation’. For that reason, it is very important that all data is within the linear range of detection. Also, a thorough washing of the beads is an important criterion, as any non-specific binding of the antibody to the negative beads will change the photon difference (P) and introduce error in the normalization factor. Lastly, the SQI is invalid if the values are outside the linearity range.

Another important consideration to ensure reproducibility is the amount of bead bound antibody. This needs to be the same between repeats to ensure the difference in number of photons (P) between the positive and negative populations remain equal between runs. Using antibodies at saturating concentrations is therefore required. At the same time, this will also ensure maximum spread. On most of the commercially available machines, manufactures use a set of excitation optics to excite the fluorochromes maximally for practical use. The maximal excitation will ensure that almost the same number of photons will be emitted for both positive and negative populations. This is especially important with kits containing many different dyes, including tandem dyes. CD4 kits consist of many different tandem dyes. It has long been established that within different lots, the brightness and emission spectra of Tandem dyes differs[4,5]. To get reproducible data, it is required to use CD4 kits from same lot. This way, we can maintain the consistency of the emission spectra from the tandem dyes.

Nguyen et al. already proved that monitoring spread via ISS values is an excellent way to identify poor laser alignment, failing of the laser, other optical problems and to reveal detectors that may need optimization. As such, SQI can also be used to track systems over time. However, fluorochrome stability, of tandem dyes specifically, is of major concern since tandem dye replacement or time dependent degradation will affect the outcome. A set of hard dyed plastic beads, which mimic single stains for all available channels would prove useful for long-term tracking.

Conclusion

The main take home messages can be summarized as follows:

SQI values are unaffected by detector voltage/gain settings.

SQI values are independent of detector type or dynamic range. Therefore, SQI values can be directly compared between any instruments.

The SQI method is a useful tool to aid in panel design.

We strongly believe that this article helps novice users to understand data spread and how to deal with the effects better. Finally, we are certain that the voltage and dynamic range independent nature of SQI values will be instrumental in preventing manipulation of detector gain settings to minimize spread, rather than to use settings that provide the best separation and sensitivity.

Materials and methods

Reagents

5 ml Falcon polystyrene round-bottom tubes (Corning Life Sciences, cat. no. 352054).

5 ml Falcon polystyrene round-bottom tube, with cell strainer snap cap (Corning Life Sciences, cat. no. 352235).

PBS as staining buffer.

Peripheral blood mononuclear cells (PBMCs) from a young healthy male adult was isolated from a leukocyte reduction chamber (LRC) purchased from the Oklahoma Blood Institute, with Lymphoprep (StemCell Technologies). Oklahoma Blood Institute performed relevant informed consents for blood product use for research purposes. Because the chambers are a byproduct of platelet donation and do not cause additional risk to the donor, the ECU Institutional Review Board does not require board review of protocols. Protocols were approved by East Carolina University Institutional Biosafety Committee protocol #01–19. All the methods were carried out in accordance with relevant guidelines and regulations.

Instruments

We have used 3 BD FACSAria Fusion cell sorters (Fusion) and a FACSSymphony A5 (Symphony) from BD Biosciences, a NovoCyte Quanteon from ACEA (Quanteon), and a CytoFLEX S (CytoFLEX) from Beckman Coulter Life Sciences. Fusion 1, 2 and 3 were used with identical optical configuration and laser output power unless stated otherwise. For one of the experiments, data was also generated at higher laser output powers on the Fusion 2 (Fusion 2_High). The output of all lasers of the Symphony and the Quanteon were set to 100mW. A detailed overview of the optical configurations (laser and filters) is given in the Supplementary Table S2. Description of Aurora from Cytek (full spectral cytometer) can be found in Ref.[13].

Antibody-fluorochrome conjugates

All experiments were done using human CD4 evaluation kit (BD Biosciences, cat no 566352). The following fluorochromes were used: BB700, PerCP-Cy5.5, BV421, BV605, BV711, PE, PE-CF594, PE-Cy5, PE-Cy7, Alexa Fluor 647, APC, APC-R700 and Alexa Fluor 700. Clone matched (SK3) Alexa Fluor 647 conjugated CD4 antibody (cat no 344635), and anti-human Alexa Fluor 647 conjugated CD45RA antibody, clone HI100 (cat no 304153) were obtained from BioLegend. All reagents were stored as recommended by the manufacturer.

Software

BD FACSDiva Software (version 8.0.1), NovoExpress and CytExpert (version 2) were used only for data acquisition. Generation of compensation matrices and other analyses were done using FlowJo version 10.6.1. Only in Aurora Spectroflo (v2.2) was used for data acquisition and unmixing.

The script

Compensated single stain fcs files are needed for the script to work. SQI values are computed using R (R 4.0.3). To do so, the population of stained beads, and unstained beads for fluorochrome X as well as the population of stained beads for fluorochrome Y need to be defined. This is done via statistical analysis using an analysis of mixture distribution. A first statistical model, using EM algorithm, is performed on the absolute value of the intensities measured by the primary detector corresponding to fluorochrome X. This analysis, performed with the R function normalmixEM from R package mixtools, decomposes the intensities measured on the primary detector as two distinct normal distributions, i.e., the unstained control beads and the stained beads populations. By taking the first and the third quartile of these distributions, low and high boundaries of both unstained and stained populations can be derived, which allows for gating of the data. The same approach is performed to gate the population of the fluorochrome Y stained beads. Once the three different population of beads are defined, the relevant parameters for SQI value computations are obtained. The script can be found in the webpage: github.com/RenTissier/RangeComputation.

Data collection

At least 35,000 single beads or 100,000 single cells were recorded. We made sure that the population intensities stay within the linear dynamic range of the detector. Manufacturer recommended quality controls were performed before every experiment.

Sample preparation

Staining of antibody capture beads

Antibody capture beads (both blank and positive) were vortex for ~ 3 s.

One drop of each capture bead was mixed in a tube.

The saturating antibody concentration was determined in titration series for every antibody-fluorochrome conjugate used.

0.25 µg of antibody was added to all samples (at saturation). 0.25 µg is well above the saturation amount for any antibody we have used.

The samples were stained at 4 °C for 30 min in the dark.

Beads were washed twice in 3 ml of staining buffer, 5 min centrifugation at 500 g (acceleration 9, deceleration 5).

Re-suspended in 200 µl of PBS.

Staining of PBMC

The cells were thawed in a 37 °C water bath.

Filtered gently.

Centrifuged at 500×g for 10 min.

The supernatant was discarded.

1 × 106 cells were re-suspended in 100 µl of staining buffer.

The saturating antibody concentration was determined in titration series for every antibody-fluorochrome conjugate used.

0.06 µg of CD4 antibody was added (at saturation).

The samples were stained at 4 °C for 30 min in the dark.

Beads were washed twice in 3 ml of staining buffer, 5 min centrifugation at 500×g (acceleration 9, deceleration 5).

Re-suspended in 200 µl of staining buffer.

Intrinsic spillover spread values are dependent on detector gain settings

Antibody capture beads were stained as mentioned above. Single stained beads were run only changing the primary PMT voltage from 600 to 400 V in steps of 100 V. A mixture of these single stains were run using different PMT voltages. ISS values were calculated and compared to the Robust Standard Deviation (rSD).

Saturated single stains of six different fluorochromes (BB700, BV605, BV711, PE-CF594, PE-Cy5 and APC-R700) were prepared as described above.

Two sets of PMT voltages were used.

Voltages in group one was set to place the positive population at ~ 110,000 MFI.

Voltages in the second group were set to place the positive population at ~ 50,000 MFI.

Range was calculated for all spillover signals.

In a separate experiment single stained beads were run for a series of PMT voltages. In one set only the voltage for the primary detector was changed, and in the second set only the voltage of secondary detector was changed. The ISS and SQI values were calculated for both sets.

Application of SQI to compare instrument performance

We have three Fusions. The comparison was done using same laser power for all. Only in case of the Fusion 2 another higher laser power setup was used to test the effect of laser power on spread.

Single stains were run in all three fusions.

Each experiment in each machine was performed in triplicates.

SQI value from all the repeat experiments were calculated and compared to check intra assay variation.

SQI, as dynamic range independent measure of spread

Fusion, Symphony, Quanteon and CytoFLEX were used for this comparison. Fusion and Symphony has a 5 decade scale. CytoFLEX and Quanteon has 7 and 7.2 decade scale respectively.

For Quanteon and CytoFLEX, positive signals were set to approximately 250,000 and 1,000,000 MFI. For the Symphony and Fusion these were set to 110,000 MFI. Spillover from BV711 was not analyzed in the CytoFLEX, because it lacked the appropriate detector.

Saturated single stain antibody capture beads were run in each system in triplicates.

Range was calculated for all spillover signals.

Single stained beads were run at different laser output powers. The SQI values were compared (same Fusion) to investigate the effect of excitation intensity on spread.

Usefulness of SQI

Human PBMCs were used as samples

7 X 106 Cells were thawed in a 37 °C water bath.

The cells were then filtered gently and centrifuged at 500×g for 10 min.

The supernatant was discarded, and the cells were re-suspended in 700 µl of staining buffer.

100 µl of the cell suspension was kept aside as complete unstained.

600 µl cells were stained 1st with BV421 or PE-Cy5 or BB700 conjugated CD4 antibody (0.36 µg) for 15 min in dark at 4 °C.

This was divided in 6 equal portions.

Five different amounts (0.01, 0.025, 0.05, 0.5 and 1.25 µg) of CD45RA antibody were added.

The samples were stained for another 30 min in dark at 4 °C,

Cells were washed twice in 3 ml of staining buffer, 5 min centrifugation at 500×g (acceleration 9, deceleration 5).

Re-suspended in 200 µl of staining buffer.

Single stains of BV421, PE-Cy5, BB700 and Alexa Fluor 647 were prepared using BD compensation beads as described previously for automated compensation.

Staining for the validation assay

Cells were stained with different antibody as described above. Before adding the antibodies 5 µl of Brilliant stain buffer (BD, Catalog No 563794) and 5 µl True Stain Monocyte blocker (BioLegend, Catalog no 426101) was added in every tube. As recommended in Ref.[13], cells were 1st stained with TCRγδ antibody for 15 min then stained with the master mix of the rest of the antibodies for another 15 min. Cells were then washed twice with 3 ml of staining buffer and finally resuspended in 250 µl of PBS, kept in ice until analyzed using a Cytek Aurora (full spectral cytometer). Final staining volume was 200 µl. We have used two panels, 1 and 2. All the details are in Supplementary Table S3.

One drop of compensation bead (UltraComp eBeads Compensation Beads: Thermo Fisher, catalog 01-2222-41) was stained with half of the antibody used to stain cell for 20 min and washed twice with PBS for single stained control.

Flow cytometer setup for the validation assay

Both panels were run using exactly same gain settings.

Four sets of gain were used.

Cytek Assay Settings or CAS was used as standard gain setup.

One setup was created by increasing the gain of all fluorescence channel by 80%. Other two setups were created by decreasing the gain of all fluorescence channel by 50 and 70% respectively.

For every gain set antibody capture beads were run to perform the unmixing, followed by the fully stained cells.

Supplementary Information.