The Daily Dongle

Three Rules for Compensation Controls

First and foremost, there must be a single stained control for every parameter in the experiment!

In addition, there are three rules for “good” compensation controls:

1) Controls need to be at least as bright or brighter than any sample the          compensation will be applied to

2) Background fluorescence should be the same for the positive and negative control

3) Compensation controls MUST match the exact experimental fluorochrome

 

1) Controls need to be at least as bright or brighter than any sample the compensation will be applied to

An important consideration is to select the sample with the brightest fluorescence of the experiment. “Dimness” is relatively irrelevant. Only brightest matters, and that is so that low spillovers can be accurately estimated. For example, if a spillover is so low that a MFI of 10,000 doesn't cause enough spillover to be above autofluorescence, then the system assumes no compensation is necessary. At a MFI of 100,000, the spillover becomes apparent and then compensation value can be accurately assessed. Compensation is only about estimating the slope. The bottom line is that because the compensation coefficients are computed based on the RATIO of the DIFFERENCE in MFI's (of the spillover channel and the primary channel), so small absolute errors in the position of the negative control become irrelevant as the positive controls become brighter. The error in the compensation coefficient is the sum of the absolute errors in the MFI's of both the negative and the positive control; the latter has an inherently much larger absolute error than the former.

2) Background fluorescence should be the same for the positive and negative control

Any carrier for binding fluorochromes can be used for single stain compensation controls, such as cells or particles. However, the positive and negative carrier of a parameter must have the same autofluorescence. This is because compensation is a subtraction algorithm. It is imperative to NOT include autofluorescence in the compensation calculation, so if the positive and negative have the same autofluorescence, then the autofluorescence contribution to the compensation spillover calculation will be zero. If this is met, one can apply the compensation matrix to any population. For example, one can compensate on particles and apply that to cells.

3) Compensation control fluorochromes MUST match the exact experimental fluorochrome

Each flurochrome has a unique emission profile.  Therefore, the amount of spillover will be different, even for fluorochromes that emit light at about the same wavelength (e.g. FITC and Alexa Fluor 488) 

This rule is even more restrictive when applied to tandem dyes.  Each lot of tandem dye (PE-TR, PE-Cy5, PerCP-Cy5.5, APC-Cy7, etc.) should be considered unique and require its own single stain control. If a user is using two different lots of PE-Cy7 in an experiment, then they need to have two PE-Cy7 compensation single stain controls, one from each lot. Different lots will have different conjugation ratios, i.e. more Cy7 conjugates to PE or less.

One final note

Finally, compensation controls must be treated in the same manner as experimental samples. This is because exposure to light and treatments like fixation/permeabilization may alter the fluorochrome, particularly the tandem conjugation ratio, i.e. lose some Cy7 on each PE molecule.

Compensation particles versus cells for single stain controls

There is no difference in the accuracy of the two approaches for compensation. However, compensation particles do have numerous benefits over using cells. First and foremost, precious sample does not need to be wasted on single stain controls. All the cells can then be used for the experimental samples. In addition, compensation particles typically provide the brightest signal possible for any given parameter. Compensation particles are also more precise. The reason that particles are more precise for compensation is fairly straightforward. Cells have a large variance in background fluorescence (a high CV), higher than particles. This means that the spillover computation has significant error for compensation coefficients where the measurement of spillover fluorescence on cells is dominated by the error in the autofluorescence measurement. On the other hand, the particles have a much smaller error in the distribution of background fluorescence, meaning that the spillover computation is far more precise.

Of course, particles have some limitations.   

1) Compensation particles cannot be used for dyes like PI, DAPI, or EMA -- but they can be used them with amine-reactive viability dyes.  (Also, some manufacturers are now providing specific dyes preloaded into particles to use as single stain compensation controls).

2)  Particles do not bind all antibody reagents and in some cases they simply are not bright enough.

3) For some experimental conditions using tandems (e.g. permeabilization/fixation), one must ensure that the fluorescence spectrum of the experiment does not alter the emission spectrum of the tandems attached to particles in a different manner than it would the tandems attached to cells.

So, in fact, in many experiments, a user may have one or two cell-based compensation controls for some parameters used together with bead based compensation controls for the other parameters.

- Question of comments can be sent to nick@treestar.com

 

New Class: "FlowJo for Integrators" @ 2011Northwest Regional Cytometry Meeting

No single analysis solution is performant enough to handle all needs.  For many institutions, improving flow cytometry analysis comes from integrating with databases and ever more hardware.
As a contribution to the Northwest Regional meeting, Tree Star is offering a "FlowJo for Integrators" practicum.  Spend a day with FlowJo architects Jay Almarode and Adam Treister to get the straight story on what's inside the software.
This one day workshop covers new or unpublicized features of FlowJo intended for advanced users who wish to:

- work with remote file repositories
- deploy FlowJo in a client/server, grid or cloud architecture
- execute FlowJo templates using command line scripts
- understand the various settings for optimizing the FlowJo calculation engine caches
- integrate results of other applications into FlowJo, or FlowJo results into other applications
- utilize FlowDx capabilities to manage calculate comparative gating statistics

Paid preregistration of $200 is required per institution with a limit of 8.  Up to 3 people per institution.
All proceeds will go to the Friday Night pub crawl.  Send an email to learning@treestar.com to sign up.

For more on NRCM, go here.

 

Cheers,

Adam Treister
President, FlowJo LLC

Why do My Platform Stats Differ from my Workspace Stats?

Users of the Proliferation and Cell Cycle Platforms have undoubtedly noticed that the statistics the model produces do not transfer exactly to the workspace when Gates are created from the model. 

For example, the figure below shows a Watson-Pragmatic cell cycle model for one of our demo files.  The associated statistics are noted on the model.  The stats pertinent to this discussion are the G₀/G₁, S, G₂/M frequencies which in the model are roughly 57% / 24% / 18%.

If you were to press the Create Gates button at the top of the model window, and accept the “Stage” prefix you would get the workspace shown below the figure. Notice that the frequencies are roughly 64% / 14% / 21%.  Why have 7% of the cells “shifted” from the S phase into the G₀/G₁ and another 3% shifted from S into the G₂/M?

The statistics are actually calculated in two different ways.  When using the models the frequencies are calculated using probabilities of membership, based on the model fit.  For cells with intensity measured within a given range, the model can be used to determine the prob

ability of the cell belonging to each of the identified populations plus sub-G₀/G₁ and super G₂/M.  The probabilities of all cells are then added up to computethe phase probabilities for the entire sample.  The summation of these probabilities will be another probability.  To create fractions that add up to ~100% a normalization step is required.  To preserve the relative relationship between phases we do not simply normalize to 100% but rather map the probabilities to those that sum to 100%.  This process can sometimes lead to small negative numbers, particularly in the sub-G0/G1 or super G2/M populations.  In these cases it is reasonable to think of a population with a negative frequency as simply a very small number of cells. 

When you press the Create Gates button FlowJo places hard line gates on the model at the points where it becomes less probable that a cell belongs to one population and more probable that it belongs to another.  Cells that are in the gate are treated as having a probability of membership of 100% for the “winning” phase and 0% for all other. 

You’ll notice in the above example that the difference in the statistics is that the gated frequencies do not take into account that there are probably a decent number of S phase cells that spill into the other two distributions, so if you have the opportunity to use the model statistics, do so.  The stochastic approach models the data better.  However, if you use the model created subpopulations, continue to add layers of gates, and in the end need to come up with fractions that add up to 100%, you will need to use the gate frequencies. 

- Bad Mo'Flow

Send comments and questions to john@treestar.com

 

 

Dongle Transmitted Diseases (D.T.D's) ?

One of the amusing perks of working here at the FJHQ is being witness to the occasional handling of USB dongles which need exchanging. We've had bad batches where the plastic casing comes off, but we've also had good-batch-dongles getting stepped on, or just worn out over the years of use. 

When these arrive, we handle them with gloves and store them up in a big plastic bag in a locked cabinet (obviously not air tight.) 

I'm curious what would be the worst case scenario for handling a dongle that may have been exposed either to unprocessed blood, or live infectious samples? What organisms can jump blood or cell suspension and end up infecting an office worker? 

inquiring minds want to know, comments welcome - also pointers to any published guidelines on handling Infectious Dongles :)

how to check your fluorescence spillover matrices ?

So here's a comp matrix as seen in FlowJo. It has some numbers. How do you check those numbers are right? 

 

Here's one idea:

make a plot for each comp control, showing the control parameter vs. it's biggest spill value- so for the PE-A control (top row), I plot it vs. PE-Cy5, etc. 

Then I overlay the plot with itself - except comped, dressed up with some math. Here's what that looks like, just for the PE control tube: 

 

To get the above figure, make the two separate figs first (PE vs PeCy5 and comp-PE vs comp-PeCy5, in this example). Set both figs background "N" (transparent). Then select both and tap T (top align) and L (left align). 

To get the formulas to evaluate in your textboxes, simply wrap them in {{}} brackets around the expression. You can mix FJML with regular math operations. For example, the formula for the median % error was: 

{{abs(1-<fjml for uncomped median> / <fjml for comped median>) * 100 }}

and for the vertical spread error, 

{{abs (1- (<fjml for SD1> + <fjml for CV1>)/(<fjml for SD2> + <fjml for CV2> ))*100 }}

These stat calculations are the interesting bit which needs some more formulation. The easy check is to line up the medians of the data before and after compensation (shown in the top textbox.)

The harder bit is to check that the vertical spread hasn't changed. In the plot above, the red and black populations should have the same footprint on the Y scale. I tried to quantitate this using SD/CV math but I get much bigger errors.. has anyone done this sort of thing systematically, I'd love to get some feedback!

of interest.. cfr11 part 21

we're getting more and more questions about compliance to this. I intend to grow my own understanding of these regulations and will use this post (static link) to keep track of that.

Most interesting fact about this, which captures some of the definition of what it means to be "compliant":

Q: Can a vendor guarantee compliant software for Part 11?

A: It is not possible, nor correct, for any vendor to offer a turnkey 'part 11 compliant system', as this would involve both the vendor AND the user to comply with procedural controls.

more of such facts here - http://www.21cfrpart11.com/pages/faq/index.htm

also fun to read FDA's warning letters - http://www.fda.gov/foi/warning_letters/s6847c.htm

and wiki entry on this regulation which suggests it's getting a lot of challenges and resistance from the industry:

As of 2007, broad sections of the regulation have been challenged as excessive, and the FDA has stated in guidance that it will exercise enforcement discretion on many parts of the rule. This has led to confusion on exactly what is required, and the rule is being revised. In practice, the requirements on access controls are the only part routinely enforced. The "predicate rules" which required the records to be kept in the first place are still in effect. If electronic records are illegible, inaccessible, or corrupted the manufacturers are still subject to those requirements.

interesting world map projections

http://www.worldmapper.org/

Many of us have already experienced the double-take when we saw the usual (Mercator's) world map vs. Peter's projection - all of the sudden, North Hemisphere countries appear compressed on  log scale, and Africa becomes huge.

The link above has several other takes. Like science research funding:

or poverty:

so check out the link and "see the world as you've never seen it".

awesome - motion magnification

A research group at MIT presented a cool technique called motion magnification back at siggraph05. I just got the memo ;) if this is new to you as well, see webpage here - http://people.csail.mit.edu/celiu/motionmag/motionmag.html

I wonder if there is a way to apply this in bioimaging - this technique would be not unlike biexponential scaling - in that it helps to visualize the data better.

Anyway, if you really wanna be blown away by what these guys are doing, watch the demo movie [link] especially the parts towards the end.

sparklines - glyphs for describing complex phenotypes with simple shapes

I watched an informative presentation by Janet Siebert of cytoanalytics regarding the utilization of Tufte's Sparklines in flow cytometry reports.

Sparklines are glyphs, or "shapes" associated with data distribution - which is defined by the type of data you're working with (fully customizable to any kind of data..). So, you could assign a "shape" to a specific pattern of a phenotypic distribution and use that instead of many words or complex figures. You can see more examples of this in action on this page - here.

Mean, Median, Mode.

Many novice users call and ask us about 'mfi' of their data. Since the acronym is ambiguous, we ask the scientist to specify if they are asking about median, mean, or geometric mean fluorescence intensity.

Mean and Geometric mean are easiest to understand - here are their formulae (from wikipedia):

Mean (aka arithmetic average) is defined as the sum of N things divided by N:

If your data is scaled in Log, you might want to use the log-friendly version of Mean, the Geometric Mean. This expression takes the Nth root of the product of N things:

Median is the number which divides the first half of your data from the second half, if sorted by some variable. Suppose you had 1 event with fluorescence of "262144" and 19 events with fluorescence of "1". The "mean" fluorescence would be (262144+19)/20 = 13108.

The Median fluorescence is 1 - if we sort the events by fluorescence and divide them in two parts, half of the events will have intensity no less than 1, the other half, no more than  1.

Depending on how you look at your data, Median usually gives the representation of a "typical" event. In the example above, no single event had fluorescence of "13108" which is what Mean would have us believe!

Also here is an interesting figure demonstrating the relationship of mode/median/mean in some types of distributions.

Courtesy of Jerrold H. Zar from Biostatistical Analysis, fourth edition, Simon & Schuester 1999