Tutorials › Statistical Methods for Single-Cell RNA Sequencing › Droplet Encapsulation and Bead Barcoding

Single-Cell Technologies and Count Data · Tutorial 2 of 700

Droplet Encapsulation and Bead Barcoding

You will learn how droplet encapsulation and bead-borne tags turn RNA from many individual cells into a high-throughput, pooled sequencing library.

Beginner 10 min read

What You'll Learn

  • Describe how microfluidic droplets isolate cells and barcoded beads during RNA capture.
  • Explain how bead-bound oligos copy cell and molecule tags into cDNA.
  • Use a Poisson model to estimate empty, single-cell, and multiplet droplets.
  • Calculate the expected yield of droplets containing one cell and one bead.
  • Recognize why finite tag sequence spaces cannot guarantee unique molecule tags.
  • Identify practical conditions that can disrupt droplet loading or barcoding.

Why Use Droplets and Barcoded Beads?

Single-cell RNA sequencing must preserve which cell each captured RNA molecule came from while processing material from many cells efficiently. Droplet protocols do this by making many small aqueous compartments in oil. A cell and a bead carrying barcode-bearing oligos can meet in a droplet; the droplet temporarily provides a separate reaction environment in which RNA is captured and copied into tagged complementary DNA (cDNA). The resulting cDNA from many droplets can then be pooled for library preparation and sequencing.

The key idea is that the bead carries many copies of oligos with the same cell barcode, while the oligos also carry molecule-tag sequences. During reverse transcription, the bead-bound oligos prime cDNA synthesis from captured RNA. Their tag sequences become part of the resulting cDNA. Thus, the identity information is carried forward even after droplets are broken and material is pooled. The bead is a delivery vehicle for many copies of the tags, not a separate measurement of each gene.

Definition: A droplet-based barcoding protocol partitions a cell suspension into many small aqueous droplets, typically using a microfluidic device and an oil phase. In protocols using barcoded beads, bead-borne oligos capture or prime RNA and contribute a cell barcode and a molecule-specific tag to the cDNA made in that droplet.

The workflow depends on successful isolation and molecular processing, not on an assumption that every droplet contains exactly one cell and one bead. Some droplets are empty, some contain a cell but no bead, and some contain multiple cells or beads. Loading concentrations and the degree to which cells and beads arrive independently determine how often these cases occur.

From Cell Suspension to Tagged cDNA

A typical droplet-generation setup brings together an aqueous suspension of cells, a suspension of beads, and oil. As the fluids pass through a microfluidic junction, the aqueous phase is divided into droplets. Surfactants in the oil help stabilize the droplets and reduce their tendency to merge. The droplets act as separate reaction compartments while cells are lysed and RNA is captured by oligos on the beads.

A bead commonly carries many copies of an oligo with a shared cell barcode and a sequence that can capture polyadenylated RNA, often using a poly(dT) region. Oligos also include molecule-specific tags, often called unique molecular identifiers (UMIs). In a droplet, an oligo binds a transcript and primes reverse transcription. The resulting cDNA contains the oligo’s barcode sequences, so later pooling does not erase the droplet-associated identity recorded in those sequences.

A cell barcode is intended to label cDNA associated with one bead and its droplet. A molecule-specific tag is intended to distinguish captured RNA molecules within a cell, but tags can be shared by different molecules because the tag sequence space is finite. These tags are sequence labels, not guarantees of perfect identification. Their detailed sequence design and later interpretation require care, but the essential point here is that both kinds of tag are copied into the cDNA during the droplet reaction.

Once the reaction has proceeded, the emulsion is broken and cDNA from the droplets is pooled. The pooled material can be amplified and converted into a sequencing library while retaining the tag sequences. The protocol therefore combines isolation during chemistry with sequence-based identity after pooling. The individual droplets need not be physically kept separate throughout the entire downstream process.

Loading and the Poisson Occupancy Model

A useful first approximation for droplet loading is a Poisson model. Let \(N\) be the number of cells entering a randomly selected droplet, and let \(\lambda_c\) be the average number of cells per droplet. If cell arrivals are independent, occur at a constant average rate, and droplets have comparable opportunities to receive cells, then \(N\sim\operatorname{Poisson}(\lambda_c)\). The corresponding probability of exactly \(k\) cells is:

$$ P(N=k)=\frac{e^{-\lambda_c}\lambda_c^k}{k!}, \qquad k=0,1,2,\ldots $$

In particular, \(P(N=0)=e^{-\lambda_c}\), \(P(N=1)=\lambda_c e^{-\lambda_c}\), and \(P(N\geq 2)=1-e^{-\lambda_c}(1+\lambda_c)\). These probabilities describe cell occupancy, not the probability that a droplet also contains a bead. If cell arrivals and bead arrivals are independent, bead occupancy can be modeled separately with its own mean, \(\lambda_b\). The chance of exactly one cell and exactly one bead is then the product of the two single-occupancy probabilities.

The model assumes well-mixed suspensions, independent arrivals, and a roughly stable droplet-generation process. It does not require a droplet to contain a cell. In practice, clumps, uneven mixing, variable droplet sizes, or cell loss can violate the assumptions. The calculations are therefore planning approximations, not a substitute for checking the actual loading and library results.

Worked Example: Estimating Single-Cell and Multiplet Droplets

Suppose a hypothetical run generates 100,000 droplets, with an average cell loading of \(\lambda_c=0.10\) cell per droplet. Assume cell arrivals follow the Poisson model. Estimate the numbers of empty droplets, single-cell droplets, and droplets with two or more cells.

The probability of an empty droplet is:

$$ P(N=0)=e^{-0.10}=0.9048 $$

The expected number of empty droplets is \(100{,}000(0.9048)=90{,}480\). For exactly one cell:

$$ P(N=1)=0.10e^{-0.10}=0.09048 $$

This gives \(100{,}000(0.09048)=9{,}048\) single-cell droplets. For two or more cells:

$$ P(N\geq 2)=1-e^{-0.10}(1+0.10)=0.004679 $$

Thus the expected number of multiplet droplets is \(100{,}000(0.004679)=468\). These estimates total approximately 99,996 droplets; the small discrepancy from 100,000 is due to rounding. Among droplets containing at least one cell, the estimated multiplet fraction is \(0.004679/(1-e^{-0.10})=0.0492\), or about 4.92%. This is a model-based estimate, not a claim that a particular run will have exactly these counts.

Worked Yield With Both Cells and Beads

The cell and bead streams can be considered separately when their arrivals are approximately independent. This makes it possible to estimate the fraction of droplets with the desired occupancy, even though a droplet can fail to receive either component. Such a calculation helps explain why the total number of droplets is not the same as the number of usable cell-and-bead reactions.

Worked Example: Droplets With One Cell and One Bead

Consider a hypothetical experiment generating 50,000 droplets. Let the average cell loading be \(\lambda_c=0.08\), and let the average bead loading be \(\lambda_b=0.90\). Assume independent Poisson arrivals for cells and beads, with stable droplet formation. Estimate the number of droplets containing exactly one cell and exactly one bead.

The probability of exactly one cell is:

$$ P(N_c=1)=0.08e^{-0.08}=0.07385 $$

The probability of exactly one bead is:

$$ P(N_b=1)=0.90e^{-0.90}=0.3659 $$

Under the independence assumption, the probability of both events is their product:

$$ P(N_c=1,N_b=1)=(0.07385)(0.3659)=0.02702 $$

The expected number of such droplets is \(50{,}000(0.02702)=1{,}351\), rounded to the nearest whole droplet. This is not the expected total number of droplets containing any cell or bead; it counts only the exact one-cell, one-bead case. If cell and bead arrivals are dependent, or the loading is uneven, multiplying the two probabilities would not be justified.

Why Tags Are Useful—and Why They Are Not Perfect

After droplets are pooled, a cell barcode helps associate cDNA with the droplet-associated cell, while a molecule-specific tag helps distinguish captured molecules within that barcode. The bead carries many copies of its barcode-bearing oligos so that multiple transcripts from the same droplet can receive the same cell label. The molecule tags vary among oligos, but the tag sequence space is finite. Different molecules can therefore receive the same molecule tag by chance; the tag is intended to distinguish molecules, not guaranteed to be unique for every molecule.

Worked Example: A Finite Molecule-Tag Sequence Space

For illustration, suppose molecule tags were generated uniformly from all possible 8-base sequences, with four possible bases at each position. Suppose 100 captured molecules from one cell independently receive tags from this sequence space. How likely is it that at least two molecules receive the same tag?

There are \(K=4^8=65{,}536\) possible tag sequences. The probability that all 100 tags are different is:

$$ P(\text{no collision})= \prod_{j=0}^{99}\left(1-\frac{j}{65{,}536}\right) $$

Therefore:

$$ P(\text{at least one collision})= 1-\prod_{j=0}^{99}\left(1-\frac{j}{65{,}536}\right) \approx 0.073 $$

The approximate collision probability is 7.3% under these hypothetical assumptions. This is the chance of at least one shared tag among the 100 molecules, not the fraction of molecules that must collide. Real protocols can use different tag lengths and designs, and tag assignment need not match this simple uniform model. The example shows why a molecule tag should not be described as a guaranteed unique label.

Diagnostics, Failure Modes, and Alternatives

Several practical problems can make the idealized model a poor description of a run. Cell clumps can enter droplets together, producing more multiplets than independent Poisson loading predicts. Poor mixing can create uneven occupancy, and variable droplet formation can make individual compartments differ in volume. Bead-free droplets cannot receive the intended bead-borne tags; droplets with multiple cells can place RNA from more than one cell under a shared droplet-associated barcode. Droplet merging or instability can also weaken the separation the method relies on.

Failures can arise after encapsulation as well. Incomplete lysis or inefficient RNA capture can reduce the amount of cDNA produced. If tag-bearing oligos or cDNA are not confined as intended, material from separate reactions may become mixed. A final library can therefore contain fewer useful cell-associated molecules than the nominal number of generated droplets suggests. In practice, investigators assess loading and library behavior using protocol-specific controls and data summaries; the Poisson calculation alone cannot diagnose these technical effects or establish that every barcode represents one valid cell.

Droplet protocols are one way to create many parallel reaction compartments. A broad alternative is to process cells in separate wells or other individually addressed compartments, which can offer different capture and library properties but does not rely on the same high-throughput droplet partitioning. The introductory material on cell capture and gene counts provides the broader context for why the identity of the source cell must be preserved. Choosing between platforms requires considering the experimental goal and protocol characteristics; the detailed comparison is outside the scope of this account.

The main limitation of droplet barcoding is that a droplet is a physical partition, not proof of one biological cell. Loading can create empty droplets and multiplets, while sequence tags have finite capacity and can be affected by technical errors. The method provides a scalable way to keep many cell-associated reactions distinguishable, but its output must be interpreted with those limits in mind.

Common Mistakes

  • Assuming every droplet contains one cell and one bead. Random loading produces empty and multiply occupied droplets as well as the desired occupancy. Use occupancy probabilities rather than treating all generated droplets as successful reactions.
  • Treating a cell barcode as a biological identity guaranteed to represent one cell. It labels material associated with a bead and its droplet; a droplet may contain multiple cells, and barcodes are not globally unique identities.
  • Claiming every molecule gets a unique molecule tag. A finite tag sequence space allows collisions. The tag is intended to distinguish molecules, but two different molecules can share a tag by chance.
  • Using Poisson calculations without checking their assumptions. Clumping, uneven suspensions, or variable droplet generation can produce occupancy patterns unlike the model. Treat the result as an estimate and investigate departures.
Key takeaway: Droplet protocols use physical partitioning for high-throughput reactions and bead-borne sequences to carry cell and molecule labels into cDNA. Poisson occupancy calculations clarify expected loading, while finite tag spaces and imperfect partitioning mean that neither droplets nor tags guarantee perfect one-cell, one-molecule identification.

Check Your Understanding

Work through each question, then compare with the solutions below.

  1. Under a Poisson cell-loading model with \(\lambda_c=0.20\), what is the probability of an empty droplet?
  2. For \(\lambda_c=0.20\), calculate the probability of exactly one cell and the probability of two or more cells.
  3. If cell and bead arrivals are independent, what calculation gives the probability of exactly one cell and exactly one bead?
  4. Why can a molecule-specific tag not guarantee that every captured molecule has its own distinct sequence?

Solutions. (1) \(P(N=0)=e^{-0.20}=0.8187\), so about 81.87% of droplets are expected to be empty. (2) \(P(N=1)=0.20e^{-0.20}=0.1637\); \(P(N\geq2)=1-e^{-0.20}(1+0.20)=0.01752\). (3) Multiply the separate probabilities: \(P(N_c=1,N_b=1)=P(N_c=1)P(N_b=1)\), provided the independent-arrival assumption is reasonable. (4) The possible tag sequences form a finite set, so separate molecules can receive the same sequence by chance.

Clinical Biostats