Overloading And unpacKing (OAK) - droplet-based combinatorial indexing for ultra-high throughput single-cell multiomic profiling

Bing Wu,Hayley M. Bennett,Xin Ye,Akshayalakshmi Sridhar,Celine Eidenschenk,Christine Everett,Evgeniya V. Nazarova,Hsu-Hsin Chen,Ivana K. Kim,Margaret Deangelis,Leah A. Owen,Cynthia Chen,Julia Lau,Minyi Shi,Jessica M. Lund,Ana Xavier-Magalhaes,Neha Patel,Yuxin Liang,Zora Modrusan,Spyros Darmanis

DOI: https://doi.org/10.1101/2024.01.23.576918

2024-01-24

Abstract:Multiomic profiling of single cells by sequencing is a powerful technique for investigating cellular diversity in complex biological systems. Although the existing droplet-based microfluidic methods have advanced single-cell sequencing, they produce a plethora of cell-free droplets and underutilize barcoding capacities due to their low cell concentration prerequisites. Meanwhile, combinatorial indexing on microplates can index cells in a more effective way; however, it requires time-consuming and laborious protocols involving multiple splitting and pooling steps. Addressing these constraints, we have developed “Overloading And unpacKing” (OAK). With reduced labor intensity, OAK can provide cost-effective multiomic profiling for hundreds of thousands of cells, offering detection sensitivity on par with commercial droplet-based methods. To demonstrate OAK’s versatility, we conducted single-cell RNA sequencing (scRNA-Seq) as well as joint single-nucleus RNA sequencing (snRNA-Seq) and single-nucleus Assay for Transposase Accessible Chromatin with sequencing (snATAC-Seq) using cell lines. We further showcased OAK’s performance on more complex samples, including differentiated bronchial epithelial cells and primary retinal tissues. Finally, we examined transcriptomic responses of 408,000 melanoma cells across around 1,000 starting lineages over a 90-day treatment with a RAF inhibitor, belvarafenib. We discovered a rare cell population (0.12%) that underwent a sequence of transcriptomic changes, resulting in belvarafenib resistance. Ultra-high throughput, broad compatibility with diverse molecular modalities, high detection sensitivity, and simplified experimental procedures distinguish OAK from previous methods, and render OAK a powerful tool for large-scale analysis of molecular signatures, even for rare cells.

Molecular Biology

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is the limitations of existing single - cell multi - omics analysis techniques in terms of high - throughput and cost - effectiveness. Specifically:

1. **Problems with existing techniques**:
- **Droplet - based microfluidic methods**: Although these methods can achieve large - scale parallel analysis, they require a low cell concentration to reduce the number of multi - cell droplets, resulting in the generation of a large number of cell - free droplets, thus wasting barcode capacity and reagents.
- **Microplate - based combinatorial indexing methods**: This method can index cells more effectively, but it requires complex operational steps, including multiple splitting and merging of cells, which is time - consuming and labor - intensive.

2. **Proposal of the OAK (Overloading And unpacKing) method**:
- **Increasing throughput**: By overloading the microfluidic chip, the number of cell - free droplets is reduced, while the number of single - cell and multi - cell droplets is increased.
- **Simplifying the experimental process**: By breaking the droplets to release all the encapsulated cells and then randomly distributing them into multiple aliquots for secondary indexing, efficient cell separation and indexing are achieved.
- **High detection sensitivity**: OAK can maintain detection sensitivity comparable to commercial droplet methods in large - scale cell samples.

3. **Application scenarios**:
- **Single - cell RNA sequencing (scRNA - Seq)**: It shows the efficiency and accuracy of OAK in processing large - scale cell samples.
- **Combined single - nucleus RNA sequencing (snRNA - Seq) and single - nucleus ATAC sequencing (snATAC - Seq)**: It verifies the applicability of OAK in the combined analysis of multiple molecular modalities.
- **Cell lineage tracing**: By conducting 90 - day RAF inhibitor (belvarafenib) treatment in melanoma cells, it shows the application of OAK in tracing rare cell lineage changes and drug resistance mechanisms.

In summary, this paper aims to overcome the limitations of existing single - cell multi - omics analysis techniques by developing the OAK method, providing a new tool with high - throughput, high - sensitivity and simple operation, especially suitable for the analysis of large - scale cell samples and the detection of rare cells.