Therapy-Induced Evolution of Human Lung Cancer Revealed by Single-Cell RNA Sequencing.

Lung cancer, the leading cause of cancer mortality, exhibits heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) of metastatic lung cancer was performed using 49 clinical biopsies obtained from 30 patients before and during targeted therapy. Over 20,000 cancer and tumor microenvironment (TME) single-cell profiles exposed a rich and dynamic tumor ecosystem. scRNA-seq of cancer cells illuminated targetable oncogenes beyond those detected clinically. Cancer cells surviving therapy as residual disease (RD) expressed an alveolar-regenerative cell signature suggesting a therapy-induced primitive cell-state transition, whereas those present at on-therapy progressive disease (PD) upregulated kynurenine, plasminogen, and gap-junction pathways. Active T-lymphocytes and decreased macrophages were present at RD and immunosuppressive cell states characterized PD. Biological features revealed by scRNA-seq were biomarkers of clinical outcomes in independent cohorts. This study highlights how therapy-induced adaptation of the multi-cellular ecosystem of metastatic cancer shapes clinical outcomes.

SpatialData: an open and universal data framework for spatial omics.

Spatially resolved omics technologies are transforming our understanding of biological tissues. However, the handling of uni- and multimodal spatial omics datasets remains a challenge owing to large data volumes, heterogeneity of data types and the lack of flexible, spatially aware data structures. Here we introduce SpatialData, a framework that establishes a unified and extensible multiplatform file-format, lazy representation of larger-than-memory data, transformations and alignment to common coordinate systems. SpatialData facilitates spatial annotations and cross-modal aggregation and analysis, the utility of which is illustrated in the context of multiple vignettes, including integrative analysis on a multimodal Xenium and Visium breast cancer study.

Validation and tuning of in situ transcriptomics image processing workflows with crowdsourced annotations.

Recent advancements in in situ methods, such as multiplexed in situ RNA hybridization and in situ RNA sequencing, have deepened our understanding of the way biological processes are spatially organized in tissues. Automated image processing and spot-calling algorithms for analyzing in situ transcriptomics images have many parameters which need to be tuned for optimal detection. Having ground truth datasets (images where there is very high confidence on the accuracy of the detected spots) is essential for evaluating these algorithms and tuning their parameters. We present a first-in-kind open-source toolkit and framework for in situ transcriptomics image analysis that incorporates crowdsourced annotations, alongside expert annotations, as a source of ground truth for the analysis of in situ transcriptomics images. The kit includes tools for preparing images for crowdsourcing annotation to optimize crowdsourced workers' ability to annotate these images reliably, performing quality control (QC) on worker annotations, extracting candidate parameters for spot-calling algorithms from sample images, tuning parameters for spot-calling algorithms, and evaluating spot-calling algorithms and worker performance. These tools are wrapped in a modular pipeline with a flexible structure that allows users to take advantage of crowdsourced annotations from any source of their choice. We tested the pipeline using real and synthetic in situ transcriptomics images and annotations from the Amazon Mechanical Turk system obtained via Quanti.us. Using real images from in situ experiments and simulated images produced by one of the tools in the kit, we studied worker sensitivity to spot characteristics and established rules for annotation QC. We explored and demonstrated the use of ground truth generated in this way for validating spot-calling algorithms and tuning their parameters, and confirmed that consensus crowdsourced annotations are a viable substitute for expert-generated ground truth for these purposes.

Laterally Aggregated Polyacrylamide Gels for Immunoprobed Isoelectric Focusing.

Immunoprobed isoelectric focusing (IEF) resolves proteins based on differences in isoelectric point (pI) and then identifies protein targets through immunoprobing of IEF-separated proteins that have been immobilized onto a gel scaffold. During the IEF stage, the gel functions as an anti-convective medium and not as a molecular sieving matrix. During the immunoprobing stage, the gel acts as an immobilization scaffold for IEF-focused proteins via photoactive moieties. Here, we characterized the effect of gel pore size on IEF separation and in-gel immunoassay performance. We modulated polyacrylamide (PA) gel pore size via lateral chain aggregation initiated by PEG monomers. During IEF, the 2% PEG highly porous PA gel formulation offered higher resolution (minimum pI difference ∼0.07 ± 0.02) than unmodified 6%T, 3.3%C (benchmark) and 6%T, 8%C (negative control) PA gels. The highly porous gels supported a pH gradient with slope and linearity comparable to benchmark gels. The partition coefficient for antibodies into the highly porous gels (K = 0.35 ± 0.02) was greater than the benchmark (3×) and negative control (1.75×) gels. The highly porous gels also had lower immunoassay background signal than the benchmark (2×) and negative control (3×) gels. Taken together, lateral aggregation creates PA gels that are suitable for both IEF and subsequent in-gel immunoprobing by mitigating immunoprobe exclusion from the gels while facilitating removal of unbound immunoprobe.

Controlling Dispersion during Single-Cell Polyacrylamide-Gel Electrophoresis in Open Microfluidic Devices.

New tools for measuring protein expression in individual cells complement single-cell genomics and transcriptomics. To characterize a population of individual mammalian cells, hundreds to thousands of microwells are arrayed on a polyacrylamide-gel-coated glass microscope slide. In this "open" fluidic device format, we explore the feasibility of mitigating diffusional losses during lysis and polyacrylamide-gel electrophoresis (PAGE) through spatial control of the pore-size of the gel layer. To reduce in-plane diffusion-driven dilution of each single-cell lysate during in-microwell chemical lysis, we photopattern and characterize microwells with small-pore-size sidewalls ringing the microwell except at the injection region. To reduce out-of-plane-diffusion-driven-dilution-caused signal loss during both lysis and single-cell PAGE, we scrutinize a selectively permeable agarose lid layer. To reduce injection dispersion, we photopattern and study a stacking-gel feature at the head of each <1 mm separation axis. Lastly, we explore a semienclosed device design that reduces the cross-sectional area of the chip, thus reducing Joule-heating-induced dispersion during single-cell PAGE. As a result, we observed a 3-fold increase in separation resolution during a 30 s separation and a >2-fold enhancement of the signal-to-noise ratio. We present well-integrated strategies for enhancing overall single-cell-PAGE performance.

Arrayed isoelectric focusing using photopatterned multi-domain hydrogels.

Isoelectric focusing (IEF) is a powerful separation method, useful for resolving subtle changes in the isoelectric point of unlabeled proteins. While microfluidic IEF has reduced the separation times from hours in traditional benchtop IEF to minutes, the enclosed devices hinder post-separation access to the sample for downstream analysis. The two-layer open IEF device presented here comprises a photopatterned hydrogel lid layer containing the chemistries required for IEF and a thin polyacrylamide bottom layer in which the analytes are separated. The open IEF device produces comparable minimum resolvable difference in isoelectric point and gradient stability to enclosed microfluidic devices while providing post-separation sample access by simple removal of the lid layer. Further, using simulations, we determine that the material properties and the length of the separation lanes are the primary factors that affect the electric field magnitude in the separation region. Finally, we demonstrate self-indexed photomasks for alignment-free fabrication of multi-domain hydrogels. We leverage this approach to generate arrayed pH gradients with a total of 80 concurrent separation lanes, which to our knowledge is the first demonstration of multiple IEF separations in series addressed by a single pair of electrodes.

Fabrication of an Open Microfluidic Device for Immunoblotting.

Given the wide adoption of polydimethylsiloxane (PDMS) for the rapid fabrication of microfluidic networks and the utility of polyacrylamide gel electrophoresis (PAGE), we develop a technique for fabrication of PAGE molecular sieving gels in PDMS microchannel networks. In developing the fabrication protocol, we trade-off constraints on materials properties of these two polymer materials: PDMS is permeable to O2 and the presence of O2 inhibits the polymerization of polyacrylamide. We present a fabrication method compatible with performing PAGE protein separations in a composite PDMS-glass microdevice, that toggles from an "enclosed" microchannel for PAGE and blotting to an "open" PA gel lane for immunoprobing and readout. To overcome the inhibitory effects of O2, we coat the PDMS channel with a 10% benzophenone solution, which quenches the inhibiting effect of O2 when exposed to UV, resulting in a PAGE-in-PDMS device. We then characterize the PAGE separation performance. Using a ladder of small-to-mid mass proteins (Trypsin Inhibitor (TI); Ovalbumin (OVA); Bovine Serum Albumin (BSA)), we observe resolution of the markers in <60 s, with separation resolution exceeding 1.0 and CVs of 8.4% for BSA-OVA and 2.4% for OVA-TI, with comparable reproducibility to glass microdevice PAGE. We show that benzophenone groups incorporated into the gel through methacrylamide can be UV-activated multiple times to photocapture protein. PDMS microchannel network is reversibly bonded to a glass slide allowing direct access to separated proteins and subsequent in situ diffusion-driven immunoprobing and total protein Sypro red staining. We see this PAGE-in-PDMS fabrication technique as expanding the application and use of microfluidic PAGE without the need for a glass microfabrication infrastructure.

Detection of Isoforms Differing by a Single Charge Unit in Individual Cells.

To measure protein isoforms in individual mammalian cells, we report single-cell resolution isoelectric focusing (scIEF) and high-selectivity immunoprobing. Microfluidic design and photoactivatable materials establish the tunable pH gradients required by IEF and precisely control the transport and handling of each 17-pL cell lysate during analysis. The scIEF assay resolves protein isoforms with resolution down to a single-charge unit, including both endogenous cytoplasmic and nuclear proteins from individual mammalian cells.

Summit: Automated Analysis of Arrayed Single-Cell Gel Electrophoresis.

New pipelines are required to automate the quantitation of emerging high-throughput electrophoretic (EP) assessment of DNA damage, or proteoform expression in single cells. EP cytometry consists of thousands of Western blots performed on a microscope slide-sized gel microwell array for single cells. Thus, EP cytometry images pose an analysis challenge that blends requirements for accurate and reproducible analysis encountered for both standard Western blots and protein microarrays. Here, we introduce the Summit algorithm to automate array segmentation, peak background subtraction, and Gaussian fitting for EP cytometry. The data structure storage of parameters allows users to perform quality control on identically processed data, yielding a ~6.5% difference in coefficient of quartile variation (CQV) of protein peak area under the curve (AUC) distributions measured by four users. Further, inspired by investigations of background subtraction methods to reduce technical variation in protein microarray measurements, we aimed to understand the trade-offs between EP cytometry analysis throughput and variation. We found an 11%-50% increase in protein peaks that passed quality control with a subtraction method similar to microarray "average on-boundary" versus an axial subtraction method. The background subtraction method only mildly influences AUC CQV, which varies between 1% and 4.5%. Finally, we determined that the narrow confidence interval for peak location and peak width parameters from Gaussian fitting yield minimal uncertainty in protein sizing. The AUC CQV differed by only ~1%-2% when summed over the peak width bounds versus the 95% peak width confidence interval. We expect Summit to be broadly applicable to other arrayed EP separations, or traditional Western blot analysis.

ReScan, a Multiplex Diagnostic Pipeline, Pans Human Sera for SARS-CoV-2 Antigens.

Comprehensive understanding of the serological response to SARS-CoV-2 infection is important for both pathophysiologic insight and diagnostic development. Here, we generate a pan-human coronavirus programmable phage display assay to perform proteome-wide profiling of coronavirus antigens enriched by 98 COVID-19 patient sera. Next, we use ReScan, a method to efficiently sequester phage expressing the most immunogenic peptides and print them onto paper-based microarrays using acoustic liquid handling, which isolates and identifies nine candidate antigens, eight of which are derived from the two proteins used for SARS-CoV-2 serologic assays: spike and nucleocapsid proteins. After deployment in a high-throughput assay amenable to clinical lab settings, these antigens show improved specificity over a whole protein panel. This proof-of-concept study demonstrates that ReScan will have broad applicability for other emerging infectious diseases or autoimmune diseases that lack a valid biomarker, enabling a seamless pipeline from antigen discovery to diagnostic using one recombinant protein source.

Subcellular western blotting of single cells.

Although immunoassays are the de facto standard for determining subcellular protein localization in individual cells, antibody probe cross-reactivity and fixation artifacts remain confounding factors. To enhance selectivity while providing single-cell resolution, we introduce a subcellular western blotting technique capable of separately assaying proteins in the 14 pL cytoplasm and 2 pL nucleus of individual cells. To confer precision fluidic control, we describe a passive multilayer microdevice that leverages the rapid transport times afforded by miniaturization. After isolating single cells in microwells, we apply single-cell differential detergent fractionation to lyse and western blot the cytoplasmic lysate, whereas the nucleus remains intact in the microwell. Subsequently, we lyse the intact nucleus and western blot the nuclear lysate. To index each protein analysis to the originating subcellular compartment, we utilize bi-directional electrophoresis, a multidimensional separation that assays the lysate from each compartment in a distinct region of the separation axis. Single-cell bi-directional electrophoresis eliminates the need for semi-subjective image segmentation algorithms required in immunocytochemistry. The subcellular, single-cell western blot is demonstrated for six targets per cell, and successfully localizes spliceosome-associated proteins solubilized from large protein and RNA complexes, even for closely sized proteins (a 7 kDa difference). Measurement of NF-κB translocation dynamics in unfixed cells at 15-min intervals demonstrates reduced technical variance compared with immunofluorescence. This chemical cytometry assay directly measures the nucleocytoplasmic protein distribution in individual unfixed cells, thus providing insight into protein signaling in heterogeneous cell populations.

Profiling protein expression in circulating tumour cells using microfluidic western blotting.

Circulating tumour cells (CTCs) are rare tumour cells found in the circulatory system of certain cancer patients. The clinical and functional significance of CTCs is still under investigation. Protein profiling of CTCs would complement the recent advances in enumeration, transcriptomic and genomic characterization of these rare cells and help define their characteristics. Here we describe a microfluidic western blot for an eight-plex protein panel for individual CTCs derived from estrogen receptor-positive (ER+) breast cancer patients. The precision handling and analysis reveals a capacity to assay sparingly available patient-derived CTCs, a biophysical CTC phenotype more lysis-resistant than breast cancer cell lines, a capacity to report protein expression on a per CTC basis and two statistically distinct GAPDH subpopulations within the patient-derived CTCs. Targeted single-CTC proteomics with the capacity for archivable, multiplexed protein analysis offers a unique, complementary taxonomy for understanding CTC biology and ascertaining clinical impact.

Single cell-resolution western blotting.

This protocol describes how to perform western blotting on individual cells to measure cell-to-cell variation in protein expression levels and protein state. Like conventional western blotting, single-cell western blotting (scWB) is particularly useful for protein targets that lack selective antibodies (e.g., isoforms) and in cases in which background signal from intact cells is confounding. scWB is performed on a microdevice that comprises an array of microwells molded in a thin layer of a polyacrylamide gel (PAG). The gel layer functions as both a molecular sieving matrix during PAGE and a blotting scaffold during immunoprobing. scWB involves five main stages: (i) gravity settling of cells into microwells; (ii) chemical lysis of cells in each microwell; (iii) PAGE of each single-cell lysate; (iv) exposure of the gel to UV light to blot (immobilize) proteins to the gel matrix; and (v) in-gel immunoprobing of immobilized proteins. Multiplexing can be achieved by probing with antibody cocktails and using antibody stripping/reprobing techniques, enabling detection of 10+ proteins in each cell. We also describe microdevice fabrication for both uniform and pore-gradient microgels. To extend in-gel immunoprobing to gels of small pore size, we describe an optional gel de-cross-linking protocol for more effective introduction of antibodies into the gel layer. Once the microdevice has been fabricated, the assay can be completed in 4-6 h by microfluidic novices and it generates high-selectivity, multiplexed data from single cells. The technique is relevant when direct measurement of proteins in single cells is needed, with applications spanning the fundamental biosciences to applied biomedicine.

Integrating qPLM and biomechanical test data with an anisotropic fiber distribution model and predictions of TGF-ß1 and IGF-1 regulation of articular cartilage fiber modulus.

A continuum mixture model with distinct collagen (COL) and glycosaminoglycan elastic constituents was developed for the solid matrix of immature bovine articular cartilage. A continuous COL fiber volume fraction distribution function and a true COL fiber elastic modulus ([Formula: see text] were used. Quantitative polarized light microscopy (qPLM) methods were developed to account for the relatively high cell density of immature articular cartilage and used with a novel algorithm that constructs a 3D distribution function from 2D qPLM data. For specimens untreated and cultured in vitro, most model parameters were specified from qPLM analysis and biochemical assay results; consequently, [Formula: see text] was predicted using an optimization to measured mechanical properties in uniaxial tension and unconfined compression. Analysis of qPLM data revealed a highly anisotropic fiber distribution, with principal fiber orientation parallel to the surface layer. For untreated samples, predicted [Formula: see text] values were 175 and 422 MPa for superficial (S) and middle (M) zone layers, respectively. TGF-[Formula: see text]1 treatment was predicted to increase and decrease [Formula: see text] values for the S and M layers to 281 and 309 MPa, respectively. IGF-1 treatment was predicted to decrease [Formula: see text] values for the S and M layers to 22 and 26 MPa, respectively. A novel finding was that distinct native depth-dependent fiber modulus properties were modulated to nearly homogeneous values by TGF-[Formula: see text]1 and IGF-1 treatments, with modulated values strongly dependent on treatment.