Design and Evaluation of a Robust CRISPR Kinetic Assay for Hot-Spot Genotyping.

Next-generation sequencing offers highly multiplexed and accurate detection of nucleic acid sequences but at the expense of complex workflows and high input requirements. The ease of use of CRISPR-Cas12 assays is attractive and may enable highly accurate detection of sequences implicated in, for example, cancer pathogenic variants. CRISPR assays often employ end-point measurements of Cas12 trans-cleavage activity after Cas12 activation by the target; however, end point-based methods can be limited in accuracy and robustness by arbitrary experimental choices. To overcome such limitations, we develop and demonstrate here an accurate assay targeting a mutation of the epidermal growth factor gene implicated in lung cancer (exon 19 deletion). The assay is based on characterizing the kinetics of Cas12 trans-cleavage to discriminate the mutant from wild-type targets. We performed extensive experiments (780 reactions) to calibrate key assay design parameters, including the guide RNA sequence, reporter sequence, reporter concentration, enzyme concentration, and DNA target type. Interestingly, we observed a competitive reaction between the target and reporter molecules that has important consequences for the design of CRISPR assays, which use preamplification to improve sensitivity. Finally, we demonstrate the assay on 18 tumor-extracted amplicons and 100 training iterations with 99% accuracy and discuss discrimination parameters and models to improve wild type versus mutant classification.

Analytical solutions for viscoelectric effects in electrokinetic nanochannels.

Understanding electrokinetic transport in nanochannels and nanopores is essential for emerging biological and electrochemical applications. The viscoelectric effect is an important mechanism implicated in the increase of local viscosity due to the polarization of a solvent under a strong electric field. However, most analyses of the viscoelectric effect have been limited to numerical analyses. In this work, we present a set of analytical solutions applicable to the physical description of viscoelectric effects in nanochannel electrokinetic systems. To achieve such closed-form solutions, we employ the Debye-Hückel approximation of small diffuse charge layer potentials compared to the thermal potential. We analyze critical parameters, including electroosmotic flow profiles, electroosmotic mobility, flow rate, and channel conductance. We compare and benchmark our analytical solutions with published predictions from numerical models. Importantly, we leverage these analytical solutions to identify essential thermophysical and nondimensional parameters that govern the behavior of these systems. We identify scaling parameters and relations among surface charge density, ionic strength, and nanochannel height.

Isotachophoresis: Theory and Microfluidic Applications.

Isotachophoresis (ITP) is a versatile electrophoretic technique that can be used for sample preconcentration, separation, purification, and mixing, and to control and accelerate chemical reactions. Although the basic technique is nearly a century old and widely used, there is a persistent need for an easily approachable, succinct, and rigorous review of ITP theory and analysis. This is important because the interest and adoption of the technique has grown over the last two decades, especially with its implementation in microfluidics and integration with on-chip chemical and biochemical assays. We here provide a review of ITP theory starting from physicochemical first-principles, including conservation of species, conservation of current, approximation of charge neutrality, pH equilibrium of weak electrolytes, and so-called regulating functions that govern transport dynamics, with a strong emphasis on steady and unsteady transport. We combine these generally applicable (to all types of ITP) theoretical discussions with applications of ITP in the field of microfluidic systems, particularly on-chip biochemical analyses. Our discussion includes principles that govern the ITP focusing of weak and strong electrolytes; ITP dynamics in peak and plateau modes; a review of simulation tools, experimental tools, and detection methods; applications of ITP for on-chip separations and trace analyte manipulation; and design considerations and challenges for microfluidic ITP systems. We conclude with remarks on possible future research directions. The intent of this review is to help make ITP analysis and design principles more accessible to the scientific and engineering communities and to provide a rigorous basis for the increased adoption of ITP in microfluidics.

Electrochemical Methods for Water Purification, Ion Separations, and Energy Conversion.

Agricultural development, extensive industrialization, and rapid growth of the global population have inadvertently been accompanied by environmental pollution. Water pollution is exacerbated by the decreasing ability of traditional treatment methods to comply with tightening environmental standards. This review provides a comprehensive description of the principles and applications of electrochemical methods for water purification, ion separations, and energy conversion. Electrochemical methods have attractive features such as compact size, chemical selectivity, broad applicability, and reduced generation of secondary waste. Perhaps the greatest advantage of electrochemical methods, however, is that they remove contaminants directly from the water, while other technologies extract the water from the contaminants, which enables efficient removal of trace pollutants. The review begins with an overview of conventional electrochemical methods, which drive chemical or physical transformations via Faradaic reactions at electrodes, and proceeds to a detailed examination of the two primary mechanisms by which contaminants are separated in nondestructive electrochemical processes, namely electrokinetics and electrosorption. In these sections, special attention is given to emerging methods, such as shock electrodialysis and Faradaic electrosorption. Given the importance of generating clean, renewable energy, which may sometimes be combined with water purification, the review also discusses inverse methods of electrochemical energy conversion based on reverse electrosorption, electrowetting, and electrokinetic phenomena. The review concludes with a discussion of technology comparisons, remaining challenges, and potential innovations for the field such as process intensification and technoeconomic optimization.

A neural network model for rapid prediction of analyte focusing in isotachophoresis.

We present the development and demonstration of a neural network (NN) model for fast and accurate prediction of whether or not a chosen analyte is focused by an isotachophoresis (ITP) buffer system. The NN model is useful in the rapid evaluation of possible ITP chemistries applicable to analytes of interest. We trained and tested the NN model for univalent species based on extensive data sets of over 10,000 anionic and 10,000 cationic ITP simulations. The NN model uses as inputs the mobilities and the acid dissociation constants of leading electrolyte ion, trailing electrolyte ion, counterion, and a single analyte as well as the leading-to-counterion concentration ratio of the leading zone. The output then indicates whether the chosen electrolyte system yields stable ITP focusing of the analyte. The prediction accuracy of the NN model is over 97.7%. We demonstrate the applicability of the NN by validating its predictions with reported experimental data for anionic and cationic ITP. We have packaged the NN model in a free, web-based application named IONN (isotachophoresis on neural network), which can be used to rapidly screen ITP electrolyte systems.

Electric field-driven microfluidics for rapid CRISPR-based diagnostics and its application to detection of SARS-CoV-2.

The rapid spread of COVID-19 across the world has revealed major gaps in our ability to respond to new virulent pathogens. Rapid, accurate, and easily configurable molecular diagnostic tests are imperative to prevent global spread of new diseases. CRISPR-based diagnostic approaches are proving to be useful as field-deployable solutions. In one basic form of this assay, the CRISPR-Cas12 enzyme complexes with a synthetic guide RNA (gRNA). This complex becomes activated only when it specifically binds to target DNA and cleaves it. The activated complex thereafter nonspecifically cleaves single-stranded DNA reporter probes labeled with a fluorophore-quencher pair. We discovered that electric field gradients can be used to control and accelerate this CRISPR assay by cofocusing Cas12-gRNA, reporters, and target within a microfluidic chip. We achieve an appropriate electric field gradient using a selective ionic focusing technique known as isotachophoresis (ITP) implemented on a microfluidic chip. Unlike previous CRISPR diagnostic assays, we also use ITP for automated purification of target RNA from raw nasopharyngeal swab samples. We here combine this ITP purification with loop-mediated isothermal amplification and the ITP-enhanced CRISPR assay to achieve detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA (from raw sample to result) in about 35 min for both contrived and clinical nasopharyngeal swab samples. This electric field control enables an alternate modality for a suite of microfluidic CRISPR-based diagnostic assays.

Detection and Discrimination of Single Nucleotide Polymorphisms by Quantification of CRISPR-Cas Catalytic Efficiency.

The specificity of CRISPR-Cas12 assays is attractive for the detection of single nucleotide polymorphisms (SNPs) implicated in, e.g., cancer and SARS-CoV-2 variants. Such assays often employ endpoint measurements of SNP or wild type (WT) activated Cas12 trans-cleavage activity; however, the fundamental kinetic effects of SNP versus WT activation remain unknown. We here show that endpoint-based assays are limited by arbitrary experimental choices (like used reporter concentration and assay duration) and work best for known target concentrations. More importantly, we show that SNP (versus WT) activation results in measurable kinetic shifts in the Cas12 trans-cleavage substrate affinity (KM) and apparent catalytic efficiency (kcat*/KM). To address endpoint-based assay limitations, we then develop an assay based on the quantification of Michaelis-Menten parameters and apply this assay to a 20 base pair WT target of the SARS-CoV-2 E gene. We find that the kcat*/KM measured for WT is 130-fold greater than the lowest kcat*/KM among all 60 measured SNPs (compared to a 4.8-fold for endpoint fluorescence of the same SNP). KM also offers a strong ability to distinguish SNPs, varies 27-fold over all the cases, and, importantly, is insensitive to the target concentration. Last, we point out trends among kinetic rates and SNP base and location within the CRISPR-Cas12 targeted region.

Enzyme Kinetics and Detector Sensitivity Determine Limits of Detection of Amplification-Free CRISPR-Cas12 and CRISPR-Cas13 Diagnostics.

Interest in CRISPR-Cas12 and CRISPR-Cas13 detection continues to increase as these detection schemes enable the specific recognition of nucleic acids. The fundamental sensitivity limits of these schemes (and their applicability in amplification-free assays) are governed by kinetic rates. However, these kinetic rates remain poorly understood, and their reporting has been inconsistent. We quantify kinetic parameters for several enzymes (LbCas12a, AsCas12a, AapCas12b, LwaCas13a, and LbuCas13a) and their corresponding limits of detection (LoD). Collectively, we present quantification of enzyme kinetics for 14 guide RNAs (gRNAs) and nucleic acid targets for a total of 50 sets of kinetic rate parameters and 25 LoDs. We validate the self-consistency of our measurements by comparing trends and limiting behaviors with a Michaelis-Menten trans-cleavage reaction kinetics model. For our assay conditions, activated Cas12 and Cas13 enzymes exhibit trans-cleavage catalytic efficiencies between order 105 and 106 M-1 s-1. For assays that use fluorescent reporter molecules (ssDNA and ssRNA) for target detection, the kinetic rates at the current assay conditions result in an amplification-free LoD in the picomolar range. The results suggest that successful detection of target requires cleavage (by an activated CRISPR enzyme) of the order of at least 0.1% of the fluorescent reporter molecules. This fraction of reporters cleaved is required to differentiate the signal from the background, and we hypothesize that this required fraction is largely independent of the detection method (e.g., endpoint vs reaction velocity) and detector sensitivity. Our results demonstrate the fundamental nature by which kinetic rates and background signal limit LoDs and thus highlight areas of improvement for the emerging field of CRISPR diagnostics.

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer. Corrigendum.

A figure in the article by Huyke et al. [(2021), J. Synchrotron Rad. 28, 1100-1113] is corrected.

Liquid Heterostructures: Generation of Liquid-Liquid Interfaces in Free-Flowing Liquid Sheets.

Chemical reactions and biological processes are frequently governed by the structure and dynamics of the interface between two liquid phases, but these interfaces are often difficult to study due to the relative abundance of the bulk liquids. Here, we demonstrate a method for generating multilayer thin film stacks of liquids, which we call liquid heterostructures. These free-flowing layered liquid sheets are produced with a microfluidic nozzle that impinges two converging jets of one liquid onto opposite sides of a third jet of another liquid. The resulting sheet consists of two layers of the first liquid enveloping an inner layer of the second liquid. Infrared microscopy, white light reflectivity, and imaging ellipsometry measurements demonstrate that the buried liquid layer has a tunable thickness and displays well-defined liquid-liquid interfaces and that this inner layer can be only tens of nanometers thick. The demonstrated multilayer liquid sheets minimize the amount of bulk liquid relative to their buried interfaces, which makes them ideal targets for spectroscopy and scattering experiments.

Lossless Medical Image Compression by Using Difference Transform.

This paper introduces a new method of compressing digital images by using the Difference Transform applied in medical imaging. The Difference Transform algorithm performs the decorrelation process of image data, and in this way improves the encoding process, achieving a file with a smaller size than the original. The proposed method proves to be competitive and in many cases better than the standards used for medical images such as TIFF or PNG. In addition, the Difference Transform can replace other transforms like Cosine or Wavelet.

Uncertainty Quantification of Michaelis-Menten Kinetic Rates and Its Application to the Analysis of CRISPR-Based Diagnostics.

Michaelis-Menten kinetics is an essential model to rationalize enzyme reactions. The quantification of Michaelis-Menten parameters can be very challenging as it is sensitive to even small experimental errors. We here present a quantification of the uncertainty inherent to the experimental determination of kinetic rate parameters for enzymatic reactions. We study the influence of several sources of uncertainty and bias, including the inner filter effect, pipetting errors, number of points in the Michaelis-Menten curve, and flat-field correction. Using Monte Carlo simulations and analyses of experimental data, we compute typical uncertainties of k c a t ${{k}_{cat}}$ , K M ${{K}_{M}}$ , and catalytic efficiency k c a t / K M ${{k}_{cat}/{K}_{M}}$ . As a salient example, we analyze the extraction of such parameters for CRISPR-Cas systems. CRISPR diagnostics have recently attracted much interest and yet reports of these enzymatic kinetic rates have been highly unreliable and inconsistent.

CRISPR Enzyme Kinetics for Molecular Diagnostics.

CRISPR-diagnostic assays have gained significant interest in the last few years. This interest has grown rapidly during the current COVID-19 pandemic, where CRISPR-diagnostics have been frontline contenders for rapid testing solutions. This surge in CRISPR-diagnostic research prompts the following question: what exactly are the achievable limits of detection and associated assay times enabled by the kinetics of enzymes such as Cas12 and Cas13? To explore this question, we here present a model based on Michaelis-Menten enzyme kinetics theory applied to CRISPR enzymes. We use the model to develop analytical solutions for reaction kinetics and develop back-of-the-envelope criteria to validate and check for consistency in reported enzyme kinetic parameters. We applied our analyses to all studies known to us, which report Michaelis-Menten-type kinetic data for CRISPR-associated enzymes. These studies include all subtypes of Cas12 and Cas13 and orthologs. We found all but one study clearly violate at least two of our three rules and therefore present data that violate basic physical limits. We performed an experimental study of reaction kinetics of LbCas12a with both ssDNA and dsDNA activators and use these data to validate our model and its predicted scaling. The validated model is used to explore CRISPR reaction time scales and the degree of reaction completion for practically relevant target concentrations applicable to CRISPR-diagnostic assays. The results have broad implications for achievable limits of detection and assay times of emerging, amplification-free CRISPR-detection methods.

Species Abundance and Reaction Off-Rate Regulate Product Formation in Reactions Accelerated Using Isotachophoresis.

We present a model for second-order and pseudo-first-order reversible chemical reactions accelerated using peak-mode isotachophoresis (ITP). In such systems, ITP preconcentrates and co-locates the reactants between the leading and trailing electrolyte zones, and this significantly accelerates chemical reactions. Our model quantifies the effects of reaction rate constants and species abundance on product formation rate. We identify two key non-dimensional parameters, which are specific groupings of reaction rate constants, species concentrations, and influx rates. We then use a regular perturbation to study the effects of reverse reaction rate and relative species abundance (and relative rates of species accumulation) on production rate. We also use this perturbation method to derive an analytical expression for the quasi-steady-state production rate achievable by ITP. Our analytical models and numerical solutions are generally applicable to a wide range of systems, which use ITP to enhance reactions. The model is also an interesting case study of the complex coupling of electric field-driven species transport and reaction kinetics.

Web-Based Open-Source Tool for Isotachophoresis.

We present the development of a client-side web-based simulator for complex electrophoresis phenomena, including isotachophoresis. The simulation tool is called Client-based Application for Fast Electrophoresis Simulation (CAFES). CAFES uses the broad cross-browser compatibility of JavaScript to provide a rapid and easy-to-use tool for coupled unsteady electromigration, diffusion, and equilibrium electrolyte reactions among multiple weak electrolytes. The code uses a stationary grid (for simplicity) and an adaptive time step to provide reliable estimates of ion concentration dynamics (including pH profile evolution), requiring no prior installation nor compilation. CAFES also offers a large database of commonly used species and their relevant physicochemical properties. We present a validation of predictions from CAFES by comparing them to experimental data of peak- and plateau-mode isotachophoresis experiments. The code yields accurate estimates of interface velocity, plateau length and relative intensity, and pH variations while significantly reducing the computation time compared to existing codes. The tool is open-source and available for free at https://microfluidics.stanford.edu/cafes.

Millisecond timescale reactions observed via X-ray spectroscopy in a 3D microfabricated fused silica mixer.

Determination of electronic structures during chemical reactions remains challenging in studies which involve reactions in the millisecond timescale, toxic chemicals, and/or anaerobic conditions. In this study, a three-dimensionally (3D) microfabricated microfluidic mixer platform that is compatible with time-resolved X-ray absorption and emission spectroscopy (XAS and XES, respectively) is presented. This platform, to initiate reactions and study their progression, mixes a high flow rate (0.50-1.5 ml min-1) sheath stream with a low-flow-rate (5-90 µl min-1) sample stream within a monolithic fused silica chip. The chip geometry enables hydrodynamic focusing of the sample stream in 3D and sample widths as small as 5 µm. The chip is also connected to a polyimide capillary downstream to enable sample stream deceleration, expansion, and X-ray detection. In this capillary, sample widths of 50 µm are demonstrated. Further, convection-diffusion-reaction models of the mixer are presented. The models are experimentally validated using confocal epifluorescence microscopy and XAS/XES measurements of a ferricyanide and ascorbic acid reaction. The models additionally enable prediction of the residence time and residence time uncertainty of reactive species as well as mixing times. Residence times (from initiation of mixing to the point of X-ray detection) during sample stream expansion as small as 2.1 ± 0.3 ms are also demonstrated. Importantly, an exploration of the mixer operational space reveals a theoretical minimum mixing time of 0.91 ms. The proposed platform is applicable to the determination of the electronic structure of conventionally inaccessible reaction intermediates.

Nondestructive nanostraw intracellular sampling for longitudinal cell monitoring.

Here, we report a method for time-resolved, longitudinal extraction and quantitative measurement of intracellular proteins and mRNA from a variety of cell types. Cytosolic contents were repeatedly sampled from the same cell or population of cells for more than 5 d through a cell-culture substrate, incorporating hollow 150-nm-diameter nanostraws (NS) within a defined sampling region. Once extracted, the cellular contents were analyzed with conventional methods, including fluorescence, enzymatic assays (ELISA), and quantitative real-time PCR. This process was nondestructive with >95% cell viability after sampling, enabling long-term analysis. It is important to note that the measured quantities from the cell extract were found to constitute a statistically significant representation of the actual contents within the cells. Of 48 mRNA sequences analyzed from a population of cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs), 41 were accurately quantified. The NS platform samples from a select subpopulation of cells within a larger culture, allowing native cell-to-cell contact and communication even during vigorous activity such as cardiomyocyte beating. This platform was applied both to cell lines and to primary cells, including CHO cells, hiPSC-CMs, and human astrocytes derived in 3D cortical spheroids. By tracking the same cell or group of cells over time, this method offers an avenue to understand dynamic cell behavior, including processes such as induced pluripotency and differentiation.

Using Ultramicroporous Carbon for the Selective Removal of Nitrate with Capacitive Deionization.

The contamination of water resources with nitrate is a growing and significant problem. Here we report the use of ultramicroporous carbon as a capacitive deionization (CDI) electrode for selectively removing nitrate from an anion mixture. Through moderate activation, we achieve a micropore-size distribution consisting almost exclusively of narrow (<1 nm) pores that are well suited for adsorbing the planar, weakly hydrated nitrate molecule. Cyclic voltammetry measurements reveal an enhanced capacitance for nitrate when compared to chloride as well as significant ion sieving effects when sulfate is the only anion present. We measure high selectivities (S) of both nitrate over sulfate (SNO3/SO4 = 17.8 ± 3.6 at 0.6 V) and nitrate over chloride (SNO3/Cl = 6.1 ± 0.4 at 0.6 V) when performing a constant voltage CDI separation on 3.33 mM/3.33 mM/1.67 mM Cl/NO3/SO4 feedwater. These results are particularly encouraging considering that a divalent interferant was present in the feed. Using molecular dynamics simulations, we examine the solvation characteristics of these ions to better understand why nitrate is preferentially electrosorbed over sulfate and chloride.

Efficient Production of On-Target Reads for Small RNA Sequencing of Single Cells Using Modified Adapters.

Although single-cell mRNA sequencing has been a powerful tool to explore cellular heterogeneity, the sequencing of small RNA at the single-cell level (sc-sRNA-seq) remains a challenge, as these have no consensus sequence, are relatively low abundant, and are difficult to amplify in a bias-free fashion. We present two methods of single-cell-lysis that enable sc-sRNA-seq. The first method is a chemical-based technique with overnight freezing while the second method leverages on-chip electrical lysis of plasma membrane and physical extraction and separation of cytoplasmic RNA via isotachophoresis. We coupled these two methods with off-chip small RNA library preparation using CleanTag modified adapters to prevent the formation of adapter dimers. We then demonstrated sc-sRNA-seq with single K562 human leukemic cells. Our approaches offer a relatively short hands-on time of 6 h and efficient generation of on-target reads. The sc-sRNA-seq with our approaches showed detection of miRNA with various abundances ranging from 16 000 copies/cell to about 10 copies/cell. We anticipate this approach will create a new opportunity to explore cellular heterogeneity through small RNA expression.

Thermodynamics of Ion Separation by Electrosorption.

We present a simple, top-down approach for the calculation of minimum energy consumption of electrosorptive ion separation using variational form of the (Gibbs) free energy. We focus and expand on the case of electrostatic capacitive deionization (CDI). The theoretical framework is independent of details of the double-layer charge distribution and is applicable to any thermodynamically consistent model, such as the Gouy-Chapman-Stern and modified Donnan models. We demonstrate that, under certain assumptions, the minimum required electric work energy is indeed equivalent to the free energy of separation. Using the theory, we define the thermodynamic efficiency of CDI. We show that the thermodynamic efficiency of current experimental CDI systems is currently very low, around 1% for most existing systems. We applied this knowledge and constructed and operated a CDI cell to show that judicious selection of the materials, geometry, and process parameters can lead to a 9% thermodynamic efficiency and 4.6 kT per removed ion energy cost. This relatively high thermodynamic efficiency is, to our knowledge, by far the highest thermodynamic efficiency ever demonstrated for traditional CDI. We hypothesize that efficiency can be further improved by further reduction of CDI cell series resistances and optimization of operational parameters.