Improving Multicolor Colocalization in Single-Vesicle Flow Cytometry with Vesicle Transit Time.

Immunophenotyping of vesicles, such as extracellular vesicles (EVs), is essential to understanding their origin and biological role. We previously described a custom-built flow analyzer that utilizes a gravity-driven flow, high numerical aperture objective, and micrometer-sized flow channels to reach the sensitivity needed for fast multidimensional analysis of the surface proteins of EVs, even down to the smallest EVs (e.g., ∼30-40 nm). It is difficult to flow focus small EVs, and thus, the transiting EVs exhibit a distribution in particle velocities due to the laminar flow. This distribution of vesicle velocities leads to potentially incorrect results when immunophenotyping nanometer-sized vesicles using cross-correlation analysis (Xcorr), as the order of appearance of the vesicles might not be the same at different spatially offset laser excitation regions. Here, we describe an alternative cross-correlation analysis strategy (Scorr), which uses information on particle transit time across the laser excitation beam width to improve multicolor colocalization in single-vesicle immunoprofiling. We tested the performance of the algorithm for colocalization analysis of multicolor nanobeads and EVs experimentally and via simulations and found that Scorr improved both the efficiency and accuracy of colocalization versus Xcorr. As shown from Monte Carlo simulations, Scorr provided an ∼1.2-4.7-fold increase in the number of colocalized peaks and ensured negligible colocalization of peaks. In silico results were in good agreement with experimental data, which showed an increase in colocalized peaks of ∼1.3-2.5-fold and ∼1.2-2-fold for multicolor beads and EVs, respectively.

Sizing Extracellular Vesicles Using Membrane Dyes and a Single Molecule-Sensitive Flow Analyzer.

Extracellular vesicles (EVs) are membranous particles released by most cells in our body, which are involved in many cell-to-cell signaling processes. Given the nanometer sizes and heterogeneity of EVs, highly sensitive methods with single-molecule resolution are fundamental to investigating their biophysical properties. Here, we demonstrate the sizing of EVs using a fluorescence-based flow analyzer with single-molecule sensitivity. Using a dye that selectively partitions into the vesicle's membrane, we show that the fluorescence intensity of a vesicle is proportional to its diameter. We discuss the constraints in sample preparation which are inherent to sizing nanoscale vesicles with a fluorescent membrane dye and propose several guidelines to improve data consistency. After optimizing staining conditions, we were able to measure the size of vesicles in the range ∼35-300 nm, covering the spectrum of EV sizes. Lastly, we developed a method to correct the signal intensity from each vesicle based on its traveling speed inside the microfluidic channel, by operating at a high sampling rate (10 kHz) and measuring the time required for the particle to cross the laser beam. Using this correction, we obtained a threefold greater accuracy in EV sizing, with a precision of ±15-25%.

Detection of 14 High-Risk Human Papillomaviruses Using Digital LAMP Assays on a Self-Digitization Chip.

Cervical cancer is the fourth-leading cause of cancer deaths among women worldwide and most cases occur in developing countries. Detection of high-risk (HR) HPV, the etiologic agent of cervical cancer, is a primary screening method for cervical cancer. However, the current gold standard for HPV detection, real-time PCR, is expensive, time-consuming, and instrumentation-intensive. A rapid, low-cost HPV detection method is needed for cervical cancer screening in low-resource settings. We previously developed a digital loop-mediated isothermal amplification (dLAMP) assay for rapid, quantitative detection of nucleic acids without the need for thermocycling. This assay employs a microfluidic self-digitization chip to automatically digitize a sample into an array of nanoliter wells in a simple assay format. Here we evaluate the dLAMP assay and self-digitization chip for detection of the commonly tested 14 high-risk HPVs in clinical samples. The dLAMP platform provided reliable genotyping and quantitative detection of the 14 high-risk HPVs with high sensitivity, demonstrating its potential for simple, rapid, and low-cost diagnosis of HPV infection.

High-Throughput Counting and Superresolution Mapping of Tetraspanins on Exosomes Using a Single-Molecule Sensitive Flow Technique and Transistor-like Semiconducting Polymer Dots.

A method for high-throughput counting and superresolution mapping of surface proteins on exosomes is described. The method combines a single-molecule sensitive flow technique and an adaptive superresolution imaging method. Exosomes stained with membrane dye and dye-conjugated antibodies were analyzed using a microfluidic platform at a flow rate of 100 exosome s-1 to determine size and protein copy number. Superresolution mapping was performed with exosomes labeled with novel transistor-like, semiconducting polymer dots (Pdots), which exhibit spontaneous blinking with <5 nm localization error and a broad range of optical-adjustable duty cycles. Based on the copy numbers extracted from the flow analysis, the switch-on frequency of the Pdots were finely adjusted so that structures of hundreds of exosomes were obtained within five minutes. The high throughput and high sensitivity of this method offer clear advantages for characterization of exosomes and similar biological vesicles.

Statistical Analysis of Nonuniform Volume Distributions for Droplet-Based Digital PCR Assays.

We present a method to determine the concentration of nucleic acids in a sample by partitioning it into droplets with a nonuniform volume distribution. This digital PCR method requires no special equipment for partitioning, unlike other methods that require nearly identical volumes. Droplets are generated by vortexing a sample in an immiscible oil to create an emulsion. PCR is performed, and droplets in the emulsion are imaged. Droplets with one or more copies of a nucleic acid are identified, and the nucleic acid concentration of the sample is determined. Numerical simulations of droplet distributions were used to estimate measurement error and dynamic range and to examine the effects of the total volume of droplets imaged and the shape of the droplet size distribution on measurement accuracy. The ability of the method to resolve 1.5- and 3-fold differences in concentration was assessed by using simulations of statistical power. The method was validated experimentally; droplet shrinkage and fusion during amplification were also assessed experimentally and showed negligible effects on measured concentration.

Isolating Rare Cells and Circulating Tumor Cells with High Purity by Sequential eDAR.

Isolation and analysis of circulating tumor cells (CTCs) from the blood of patients at risk of metastatic cancers is a promising approach to improving cancer treatment. However, CTC isolation is difficult due to low CTC abundance and heterogeneity. Previously, we reported an ensemble-decision aliquot ranking (eDAR) platform for the rare cell and CTC isolation with high throughput, greater than 90% recovery, and high sensitivity, allowing detection of low surface antigen-expressing cells linked to metastasis. Here we demonstrate a sequential eDAR platform capable of isolating rare cells from whole blood with high purity. This improvement in purity is achieved by using a sequential sorting and flow stretching design in which whole blood is sorted and fluid elements are stretched using herringbone features and the parabolic flow profile being sorted a second time. This platform can be used to collect single CTCs in a multiwell plate for downstream analysis.

Single-Molecule Flow Platform for the Quantification of Biomolecules Attached to Single Nanoparticles.

We describe here a flow platform for quantifying the number of biomolecules on individual fluorescent nanoparticles. The platform combines line-confocal fluorescence detection with near nanoscale channels (1-2 µm in width and height) to achieve high single-molecule detection sensitivity and throughput. The number of biomolecules present on each nanoparticle was determined by deconvolving the fluorescence intensity distribution of single-nanoparticle-biomolecule complexes with the intensity distribution of single biomolecules. We demonstrate this approach by quantifying the number of streptavidins on individual semiconducting polymer dots (Pdots); streptavidin was rendered fluorescent using biotin-Alexa647. This flow platform has high-throughput (hundreds to thousands of nanoparticles detected per second) and requires minute amounts of sample (∼5 µL at a dilute concentration of 10 pM). This measurement method is an additional tool for characterizing synthetic or biological nanoparticles.

Optically Encoded Semiconducting Polymer Dots with Single-Wavelength Excitation for Barcoding and Tracking of Single Cells.

Multiplexed optical encoding is emerging as a powerful technique for high-throughput cellular analysis and molecular assays. Most of the developed optical barcodes, however, either suffer from large particle size or are incompatible with most commercial optical instruments. Here, a new type of nanoscale fluorescent barcode (Pdot barcodes) was prepared from semiconducting polymers. The Pdot barcodes possess the merits of small size (∼20 nm in diameter), narrow emission bands (full-width-at-half-maximum (fwhm) of 30-40 nm), three-color emissions (blue, green, and red) under single-wavelength excitation, a high brightness, good pH and thermal stability, and efficient cellular uptake. The Pdot barcodes were prepared using a three-color and six-intensity encoding strategy; for ratiometric readout of the barcodes, one of the colors might be used as an internal reference. We used the Pdot barcodes to label 20 sets of cancer cells and then distinguished and identified each set based on the Pdot barcodes using flow cytometry. We also monitored and tracked single cells labeled with different Pdot barcodes, even through rounds of cell division. These results suggest Pdot barcodes are strong candidates for discriminating different labeled cell and for long-term cell tracking.

Squaraine-based polymer dots with narrow, bright near-infrared fluorescence for biological applications.

This article describes the design and development of squaraine-based semiconducting polymer dots (Pdots) that show large Stokes shifts and narrow-band emissions in the near-infrared (NIR) region. Fluorescent copolymers containing fluorene and squaraine units were synthesized and used as precursors for preparing the Pdots, where exciton diffusion and likely through-bond energy transfer led to highly bright and narrow-band NIR emissions. The resulting Pdots exhibit the emission full width at half-maximum of ∼36 nm, which is ∼2 times narrower than those of inorganic quantum dots in the same wavelength region (∼66 nm for Qdot705). The squaraine-based Pdots show a high fluorescence quantum yield (QY) of 0.30 and a large Stokes shift of ∼340 nm. Single-particle analysis indicates that the average per-particle brightness of the Pdots is ∼6 times higher than that of Qdot705. We demonstrate bioconjugation of the squaraine Pdots and employ the Pdot bioconjugates in flow cytometry and cellular imaging applications. Our results suggest that the narrow bandwidth, high QY, and large Stokes shift are promising for multiplexed biological detections.

Comparison of EV characterization by commercial high-sensitivity flow cytometers and a custom single-molecule flow cytometer.

High-sensitivity flow cytometers have been developed for multi-parameter characterization of single extracellular vesicles (EVs), but performance varies among instruments and calibration methods. Here we compare the characterization of identical (split) EV samples derived from human colorectal cancer (DiFi) cells by three high-sensitivity flow cytometers, two commercial instruments, CytoFLEX/CellStream, and a custom single-molecule flow cytometer (SMFC). DiFi EVs were stained with the membrane dye di-8-ANEPPS and with PE-conjugated anti-EGFR or anti-tetraspanin (CD9/CD63/CD81) antibodies for estimation of EV size and surface protein copy numbers. The limits of detection (LODs) for immunofluorescence and vesicle size based on calibration using cross-calibrated, hard-dyed beads were ∼10 PE/∼80 nm EV diameter for CytoFLEX and ∼10 PEs/∼67 nm for CellStream. For the SMFC, the LOD for immunofluorescence was 1 PE and ≤ 35 nm for size. The population of EVs detected by each system (di-8-ANEPPS+/PE+ particles) differed widely depending on the LOD of the system; for example, CellStream/CytoFLEX detected only 5.7% and 1.5% of the tetraspanin-labelled EVs detected by SMFC, respectively, and median EV diameter and antibody copy numbers were much larger for CellStream/CytoFLEX than for SMFC as measured and validated using super-resolution/single-molecule TIRF microscopy. To obtain a dataset representing a common EV population analysed by all three platforms, we filtered out SMFC and CellStream measurements for EVs below the CytoFLEX LODs as determined by bead calibration (10 PE/80 nm). The inter-platform agreement using this filtered dataset was significantly better than for the unfiltered dataset, but even better concordance between results was obtained by applying higher cutoffs (21 PE/120 nm) determined by threshold analysis using the SMFC data. The results demonstrate the impact of specifying LODs to define the EV population analysed on inter-instrument reproducibility in EV flow cytometry studies, and the utility of threshold analysis of SMFC data for providing semi-quantitative LOD values for other flow cytometers.

Supercoiled pseudocircular domains in single-twisted DNAs under tension: Elastic constants and unwinding dynamics in complexes with Topo I.

Extension versus twist data of Koster et al. (Nature 2005, 434, 671-674) are analyzed to obtain C for the main-chain segments and the twist energy parameter (ET ) for the supercoiled pseudocircular (sp) domain(s) from which C is estimated via simulations. The torsional rigidity in the tension-free sp domain(s) (C = 163 fJ fm) is typical of the unstrained DNA and is less than half the value in the main-chain segments under tension (C = 350-410 fJ fm). Tension is suggested to induce a structural transition to a torsionally stiffer state. Data of Koster et al. for the rate of extension owing to unwinding of a covalent complex of DNA with human Topoisomerase Ib (H Topo I) are analyzed to determine the torque and rate of rotation from which an effective friction coefficient is obtained. A Langevin equation for the unwinding motion in a supercoiled DNA:H Topo I complex is solved to obtain the temporal trajectory of the average winding angle and the time-dependent distribution of winding angles. The mean rate constant for the religation reaction is estimated from the measured probability of reaction per turn. We predict that unwinding proceeds rather far during a single-cleavage and religation cycle, and is effectively completely equilibrated during the 3.2 cleavage and religation cycles that occur during each noncovalent binding and dissociation event. H Topo I is predicted to be completely processive as in accord with observations on calf-thymus Topo I (Brewood et al., Biochemistry 2010, 49, 3367-3380).

Semiconducting polymer dots doped with europium complexes showing ultranarrow emission and long luminescence lifetime for time-gated cellular imaging.

Bright dots: Semiconducting polymer dots (Pdots) doped with europium complexes possess line-like fluorescence emission, high quantum yield, and long fluorescence lifetime. The Pdots successfully labeled receptors on cells. The long fluorescence lifetime of the Pdots was used to distinguish them from other red fluorescence emitting nanoparticles, and improve the signal-to-noise ratio for time-gated cellular imaging. PVK=poly(9-vinylcarbazole).

Protein quantification at the single vesicle level reveals that a subset of synaptic vesicle proteins are trafficked with high precision.

Protein sorting represents a potential point of regulation in neurotransmission because it dictates the protein composition of synaptic vesicles, the organelle that mediates transmitter release. Although the average number of most vesicle proteins has been estimated using bulk biochemical approaches (Takamori et al., 2006), no information exists on the intervesicle variability of protein number, and thus on the precision with which proteins are sorted to vesicles. To address this, we adapted a single molecule quantification approach (Mutch et al., 2007) and used it to quantify both the average number and variance of seven integral membrane proteins in brain synaptic vesicles. We report that four vesicle proteins, SV2, the proton ATPase, Vglut1, and synaptotagmin 1, showed little intervesicle variation in number, indicating they are sorted to vesicles with high precision. In contrast, the apparent number of VAMP2/synaptobrevin 2, synaptophysin, and synaptogyrin demonstrated significant intervesicle variability. These findings place constraints on models of protein function at the synapse and raise the possibility that changes in vesicle protein expression affect vesicle composition and functioning.

Single-molecule fluorescence quantification with a photobleached internal standard.

In cellular and molecular biology, fluorophores are employed to aid in tracking and quantifying molecules involved in cellular function. We previously developed a sensitive single-molecule quantification technique to count the number of proteins and the variation of the protein number over the population of individual subcellular organelles. However, environmental effects on the fluorescent intensity of fluorophores can make it difficult to accurately quantify proteins using these sensitive techniques. In this letter, we demonstrate the use of photobleaching to extract an accurate single-molecule calibration intensity distribution from the sample directly to avoid any differences in environment that may alter the count. Using this technique, we were able to show that goat antimouse IgG antibody labeled with Alexa Fluor 488, an environmentally insensitive fluorophore, exhibited an average fluorescence equivalent to 4.6 single fluorophores. SynaptopHluorin vesicles, which contain the environmentally sensitive green fluorescent protein, exhibited an average of 4.4 single green fluorescent proteins per vesicle.

Measurements of the acidification kinetics of single SynaptopHluorin vesicles.

Uptake of neurotransmitters into synaptic vesicles is driven by the proton gradient established across the vesicle membrane. The acidification of synaptic vesicles, therefore, is a crucial component of vesicle function. Here we present measurements of acidification rate constants from isolated, single synaptic vesicles. Vesicles were purified from mice expressing a fusion protein termed SynaptopHluorin created by the fusion of VAMP/synaptobrevin to the pH-sensitive super-ecliptic green fluorescent protein. We calibrated SynaptopHluorin fluorescence to determine the relationship between fluorescence intensity and internal vesicle pH, and used these values to measure the rate constant of vesicle acidification. We also measured the effects of ATP, glutamate, and chloride on acidification. We report acidification time constants of 500 ms to 1 s. The rate of acidification increased with increasing extravesicular concentrations of ATP and glutamate. These data provide an upper and a lower bound for vesicle acidification and indicate that vesicle readiness can be regulated by changes in energy and transmitter availability.

Large structural change in isolated synaptic vesicles upon loading with neurotransmitter.

The size of a synaptic vesicle (SV) is generally thought to be determined by the amount of lipid and membrane protein it contains. Once formed, it is thought to remain constant in size. Using fluorescence correlation spectroscopy and cryogenic electron microscopy, we show that glutamatergic vesicles reversibly increase their size upon filling with glutamate. The increase ( approximately 25% in diameter) corresponds to an increase in surface area of approximately 50% and in volume of approximately 100%. This large size increase implies a large structural change in the SV upon loading with neurotransmitters. Vesicles lacking SV protein 2A (SV2A) did not manifest a change in size after loading with glutamate, indicating that SV2A is required for this phenomenon.

Sizing subcellular organelles and nanoparticles confined within aqueous droplets.

This article describes two complementary techniques, single-particle tracking and correlation spectroscopy, for accurately sizing nanoparticles confined within picoliter volume aqueous droplets. Single-particle tracking works well with bright particles that can be continuously illuminated and imaged, and we demonstrated this approach for sizing single fluorescent beads. Fluorescence correlation spectroscopy detects small intensity bursts from particles or molecules diffusing through the confocal probe volume, which works well with dim and rapidly diffusing particles or molecules; we demonstrated FCS for sizing synaptic vesicles confined in aqueous droplets. In combination with recent advances in droplet manipulations and analysis, we anticipate this capability to size single nanoparticles and molecules in free solution will complement existing tools for probing cellular systems, subcellular organelles, and nanoparticles.

Estimation of the persistence length of DNA from the torsion elastic constant and supercoiling free energy: effect of ethylene glycol.

Two different methods are proposed to estimate the persistence length ( P) of DNA from the measured torsion elastic constant (alpha) and the twist energy parameter ( E T ) that governs the supercoiling free energy. The first method involves Monte Carlo simulations and reversible-work calculations of E T for model DNAs that possess the measured alpha and selected trial values of P. Comparison of the computed E T values with the experimental value allows estimation of P (or equivalently the bending elastic constant (kappa beta)) by interpolation. A far simpler, though less accurate, alternative is to solve a previously conjectured analytical relation connecting E T , alpha, kappa beta (or P), and an unknown "constant" ( B). The present simulations are used to ascertain the optimum value of B and to assess the validity and accuracy of that relation. Within the simulation errors, P values obtained from the measured alpha and E T via this analytical expression agree with those determined from the simulations and E T values reckoned from the input alpha and kappa beta by this analytical expression agree with the corresponding simulated values. Although B is found to be insensitive to variation in alpha, it appears to decline slightly with increasing kappa beta. The original analytical expression is modified to take this apparent variation of B with kappa beta into account. By using this modified analytical relation to estimate P (from the measured alpha and E T ) or E T (from the input alpha and kappa beta), much closer agreement is obtained respectively with the values of P or E T obtained from the simulations. As specific examples, these methods are applied to determine P in 0 and 20 w/v % ethylene glycol, which has been shown to induce a structural transition in duplex DNA.

Self-digitization chip for single-cell genotyping of cancer-related mutations.

Cancer is a heterogeneous disease, and patient-level genetic assessments can guide therapy choice and impact prognosis. However, little is known about the impact of genetic variability within a tumor, intratumoral heterogeneity (ITH), on disease progression or outcome. Current approaches using bulk tumor specimens can suggest the presence of ITH, but only single-cell genetic methods have the resolution to describe the underlying clonal structures themselves. Current techniques tend to be labor and resource intensive and challenging to characterize with respect to sources of biological and technical variability. We have developed a platform using a microfluidic self-digitization chip to partition cells in stationary volumes for cell imaging and allele-specific PCR. Genotyping data from only confirmed single-cell volumes is obtained and subject to a variety of relevant quality control assessments such as allele dropout, false positive, and false negative rates. We demonstrate single-cell genotyping of the NPM1 type A mutation, an important prognostic indicator in acute myeloid leukemia, on single cells of the cell line OCI-AML3, describing a more complex zygosity distribution than would be predicted via bulk analysis.

Optical gradient flow focusing.

This paper describes a new method for carrying out flow cytometry, which employs optical gradient forces to guide and focus particles in the fluid flow. An elliptically shaped Gaussian beam was focused at the center of a microchannel to exert radiation pressure on suspended nanoparticles that are passing through the channel, such that these particles are guided to the center of the channel for efficient detection and sorting. To verify the efficiency of this optical-gradient-flow-focusing method, we present numerical simulations of the trajectories of the nanoparticles in both electroosmotic flow (EOF) and pressure-driven flow (PDF).