In the last decade, single-cell RNA sequencing routinely performed on large numbers of single cells has greatly advanced our understanding of the underlying heterogeneity of complex biological systems. Technological advances have also enabled protein measurements, further contributing to the elucidation of cell types and states present in complex tissues. Recently, there have been independent advances in mass spectrometric techniques bringing us one step closer to characterizing single-cell proteomes. Here we discuss the challenges of detecting proteins in single cells by both mass spectrometry and sequencing-based methods. We review the state of the art for these techniques and propose that there is a space for technological advancements and complementary approaches that maximize the advantages of both classes of technologies.

Capture array-based spatial transcriptomics methods have been widely used to resolve gene expression in tissues; however, their spatial resolution is limited by the density of the array. Here we present expansion spatial transcriptomics to overcome this limitation by clearing and expanding tissue prior to capturing the entire polyadenylated transcriptome with an enhanced protocol. This approach enables us to achieve higher spatial resolution while retaining high library quality, which we demonstrate using mouse brain samples.

We report the evolution of mScarlet3, a cysteine-free monomeric red fluorescent protein with fast and complete maturation, as well as record brightness, quantum yield (75%) and fluorescence lifetime (4.0 ns). The mScarlet3 crystal structure reveals a barrel rigidified at one of its heads by a large hydrophobic patch of internal residues. mScarlet3 behaves well as a fusion tag, displays no apparent cytotoxicity and it surpasses existing red fluorescent proteins as a Förster resonance energy transfer acceptor and as a reporter in transient expression systems.

As a miniature RNA-guided endonuclease, IscB is presumed to be the ancestor of Cas9 and to share similar functions. IscB is less than half the size of Cas9 and thus more suitable for in vivo delivery. However, the poor editing efficiency of IscB in eukaryotic cells limits its in vivo applications. Here we describe the engineering of OgeuIscB and its corresponding ωRNA to develop an IscB system that is highly efficient in mammalian systems, named enIscB. By fusing enIscB with T5 exonuclease (T5E), we found enIscB-T5E exhibited comparable targeting efficiency to SpG Cas9 while showing reduced chromosome translocation effects in human cells. Furthermore, by fusing cytosine or adenosine deaminase with enIscB nickase, we generated miniature IscB-derived base editors (miBEs), exhibiting robust editing efficiency (up to 92%) to induce DNA base conversions. Overall, our work establishes enIscB-T5E and miBEs as versatile tools for genome editing.

High-throughput profiling methods (such as genomics or imaging) have accelerated basic research and made deep molecular characterization of patient samples routine. These approaches provide a rich portrait of genes, molecular pathways and cell types involved in disease phenotypes. Machine learning (ML) can be a useful tool for extracting disease-relevant patterns from high-dimensional datasets. However, depending upon the complexity of the biological question, machine learning often requires many samples to identify recurrent and biologically meaningful patterns. Rare diseases are inherently limited in clinical cases, leading to few samples to study. In this Perspective, we outline the challenges and emerging solutions for using ML for small sample sets, specifically in rare diseases. Advances in ML methods for rare diseases are likely to be informative for applications beyond rare diseases for which few samples exist with high-dimensional data. We propose that the method community prioritize the development of ML techniques for rare disease research.

The discovery of biomolecular condensates transformed our understanding of intracellular compartmentalization of molecules. To integrate interdisciplinary scientific knowledge about the function and composition of biomolecular condensates, we developed the crowdsourcing condensate database and encyclopedia ( cd-code.org ). CD-CODE is a community-editable platform, which includes a database of biomolecular condensates based on the literature, an encyclopedia of relevant scientific terms and a crowdsourcing web application. Our platform will accelerate the discovery and validation of biomolecular condensates, and facilitate efforts to understand their role in disease and as therapeutic targets.

Understanding the routing of neuronal information requires the functional characterization of connections. Neuronal projections recruit large postsynaptic ensembles with distinct postsynaptic response types (PRTs). PRT is typically probed by low-throughput whole-cell electrophysiology and is not a selection criterion for single-cell RNA-sequencing (scRNA-seq). To overcome these limitations and target neurons based on specific PRTs for soma harvesting and subsequent scRNA-seq, we created Voltage-Seq. We established all-optical voltage imaging and recorded the PRT of 8,347 neurons in the mouse periaqueductal gray (PAG) evoked by the optogenetic activation of ventromedial hypothalamic (VMH) terminals. PRTs were classified and spatially resolved in the entire VMH-PAG connectome. We built an onsite analysis tool named VoltView to navigate soma harvesting towards target PRTs guided by a classifier that used the VMH-PAG connectome database as a reference. We demonstrated Voltage-seq by locating VMH-driven γ-aminobutyric acid-ergic neurons in the PAG, guided solely by the onsite classification in VoltView.

Genome-wide measurements of RNA structure can be obtained using reagents that react with unpaired bases, leading to adducts that can be identified by mutational profiling on next-generation sequencing machines. One drawback of these experiments is that short sequencing reads can rarely be mapped to specific transcript isoforms. Consequently, information is acquired as a population average in regions that are shared between transcripts, thus blurring the underlying structural landscape. Here, we present nanopore dimethylsulfate mutational profiling (Nano-DMS-MaP)-a method that exploits long-read sequencing to provide isoform-resolved structural information of highly similar RNA molecules. We demonstrate the value of Nano-DMS-MaP by resolving the complex structural landscape of human immunodeficiency virus-1 transcripts in infected cells. We show that unspliced and spliced transcripts have distinct structures at the packaging site within the common 5' untranslated region, likely explaining why spliced viral RNAs are excluded from viral particles. Thus, Nano-DMS-MaP is a straightforward method to resolve biologically important transcript-specific RNA structures that were previously hidden in short-read ensemble analyses.

Analyzing proteins from single cells by tandem mass spectrometry (MS) has recently become technically feasible. While such analysis has the potential to accurately quantify thousands of proteins across thousands of single cells, the accuracy and reproducibility of the results may be undermined by numerous factors affecting experimental design, sample preparation, data acquisition and data analysis. We expect that broadly accepted community guidelines and standardized metrics will enhance rigor, data quality and alignment between laboratories. Here we propose best practices, quality controls and data-reporting recommendations to assist in the broad adoption of reliable quantitative workflows for single-cell proteomics. Resources and discussion forums are available at https://single-cell.net/guidelines .

The count table, a numeric matrix of genes × cells, is the basic input data structure in the analysis of single-cell RNA-sequencing data. A common preprocessing step is to adjust the counts for variable sampling efficiency and to transform them so that the variance is similar across the dynamic range. These steps are intended to make subsequent application of generic statistical methods more palatable. Here, we describe four transformation approaches based on the delta method, model residuals, inferred latent expression state and factor analysis. We compare their strengths and weaknesses and find that the latter three have appealing theoretical properties; however, in benchmarks using simulated and real-world data, it turns out that a rather simple approach, namely, the logarithm with a pseudo-count followed by principal-component analysis, performs as well or better than the more sophisticated alternatives. This result highlights limitations of current theoretical analysis as assessed by bottom-line performance benchmarks.

An unambiguous description of an experiment, and the subsequent biological observation, is vital for accurate data interpretation. Minimum information guidelines define the fundamental complement of data that can support an unambiguous conclusion based on experimental observations. We present the Minimum Information About Disorder Experiments (MIADE) guidelines to define the parameters required for the wider scientific community to understand the findings of an experiment studying the structural properties of intrinsically disordered regions (IDRs). MIADE guidelines provide recommendations for data producers to describe the results of their experiments at source, for curators to annotate experimental data to community resources and for database developers maintaining community resources to disseminate the data. The MIADE guidelines will improve the interpretability of experimental results for data consumers, facilitate direct data submission, simplify data curation, improve data exchange among repositories and standardize the dissemination of the key metadata on an IDR experiment by IDR data sources.

We propose two new measures of resolution anisotropy for cryogenic electron microscopy maps: Fourier shell occupancy (FSO), and the Bingham test (BT). FSO varies from 1 to 0, with 1 representing perfect isotropy, and lower values indicating increasing anisotropy. The threshold FSO = 0.5 occurs at Fourier shell correlation resolution. BT is a hypothesis test that complements the FSO to ensure the existence of anisotropy. FSO and BT allow visualization of resolution anisotropy. We illustrate their use with different experimental cryogenic electron microscopy maps.

As long-read sequencing technologies are becoming increasingly popular, a number of methods have been developed for the discovery and analysis of structural variants (SVs) from long reads. Long reads enable detection of SVs that could not be previously detected from short-read sequencing, but computational methods must adapt to the unique challenges and opportunities presented by long-read sequencing. Here, we summarize over 50 long-read-based methods for SV detection, genotyping and visualization, and discuss how new telomere-to-telomere genome assemblies and pangenome efforts can improve the accuracy and drive the development of SV callers in the future.

Neuronal computation is achieved through connections of individual neurons into a larger network. To expand the repertoire of endogenous cellular communication, we developed a synthetic, photon-assisted synaptic transmission (PhAST) system. PhAST is based on luciferases and channelrhodopsins that enable the transmission of a neuronal state across space, using photons as neurotransmitters. PhAST overcomes synaptic barriers and rescues the behavioral deficit of a glutamate mutant with conditional, calcium-triggered photon emission between two neurons of the Caenorhabditis elegans nociceptive avoidance circuit. To demonstrate versatility and flexibility, we generated de novo synaptic transmission between two unconnected cells in a sexually dimorphic neuronal circuit, suppressed endogenous nocifensive response through activation of an anion channelrhodopsin and switched attractive to aversive behavior in an olfactory circuit. Finally, we applied PhAST to dissect the calcium dynamics of the temporal pattern generator in a motor circuit for ovipositioning. In summary, we established photon-based synaptic transmission that facilitates the modification of animal behavior.

Cryogenic electron microscopy and data processing enable the determination of structures of isolated macromolecules to near-atomic resolution. However, these data do not provide structural information in the cellular environment where macromolecules perform their native functions, and vital molecular interactions can be lost during the isolation process. Cryogenic focused ion beam (FIB) fabrication generates thin lamellae of cellular samples and tissues, enabling structural studies on the near-native cellular interior and its surroundings by cryogenic electron tomography (cryo-ET). Cellular cryo-ET benefits from the technological developments in electron microscopes, detectors and data processing, and more in situ structures are being obtained and at increasingly higher resolution. In this Review, we discuss recent studies employing cryo-ET on FIB-generated lamellae and the technological developments in ultrarapid sample freezing, FIB fabrication of lamellae, tomography, data processing and correlative light and electron microscopy that have enabled these studies. Finally, we explore the future of cryo-ET in terms of both methods development and biological application.

Even though the recent advances in 'complete genomics' revealed the previously inaccessible genomic regions, analysis of variations in centromeres and other extra-long tandem repeats (ETRs) faces an algorithmic challenge since there are currently no tools for accurate sequence comparison of ETRs. Counterintuitively, the classical alignment approaches, such as the Smith-Waterman algorithm, fail to construct biologically adequate alignments of ETRs. We present UniAligner-the parameter-free sequence alignment algorithm with sequence-dependent alignment scoring that automatically changes for any pair of compared sequences. UniAligner prioritizes matches of rare substrings that are more likely to be relevant to the evolutionary relationship between two sequences. We apply UniAligner to estimate the mutation rates in human centromeres, and quantify the extremely high rate of large duplications and deletions in centromeres. This high rate suggests that centromeres may represent some of the most rapidly evolving regions of the human genome with respect to their structural organization.

Homology-directed repair (HDR), a method for repair of DNA double-stranded breaks can be leveraged for the precise introduction of mutations supplied by synthetic DNA donors, but remains limited by low efficiency and off-target effects. In this study, we report HDRobust, a high-precision method that, via the combined transient inhibition of nonhomologous end joining and microhomology-mediated end joining, resulted in the induction of point mutations by HDR in up to 93% (median 60%, s.e.m. 3) of chromosomes in populations of cells. We found that, using this method, insertions, deletions and rearrangements at the target site, as well as unintended changes at other genomic sites, were largely abolished. We validated this approach for 58 different target sites and showed that it allows efficient correction of pathogenic mutations in cells derived from patients suffering from anemia, sickle cell disease and thrombophilia.

Genetically encoded indicators engineered from G-protein-coupled receptors are important tools that enable high-resolution in vivo neuromodulator imaging. Here, we introduce a family of sensitive multicolor norepinephrine (NE) indicators, which includes nLightG (green) and nLightR (red). These tools report endogenous NE release in vitro, ex vivo and in vivo with improved sensitivity, ligand selectivity and kinetics, as well as a distinct pharmacological profile compared with previous state-of-the-art GRABNE indicators. Using in vivo multisite fiber photometry recordings of nLightG, we could simultaneously monitor optogenetically evoked NE release in the mouse locus coeruleus and hippocampus. Two-photon imaging of nLightG revealed locomotion and reward-related NE transients in the dorsal CA1 area of the hippocampus. Thus, the sensitive NE indicators introduced here represent an important addition to the current repertoire of indicators and provide the means for a thorough investigation of the NE system.

Unsupervised clustering of single-cell RNA-sequencing data enables the identification of distinct cell populations. However, the most widely used clustering algorithms are heuristic and do not formally account for statistical uncertainty. We find that not addressing known sources of variability in a statistically rigorous manner can lead to overconfidence in the discovery of novel cell types. Here we extend a previous method, significance of hierarchical clustering, to propose a model-based hypothesis testing approach that incorporates significance analysis into the clustering algorithm and permits statistical evaluation of clusters as distinct cell populations. We also adapt this approach to permit statistical assessment on the clusters reported by any algorithm. Finally, we extend these approaches to account for batch structure. We benchmarked our approach against popular clustering workflows, demonstrating improved performance. To show practical utility, we applied our approach to the Human Lung Cell Atlas and an atlas of the mouse cerebellar cortex, identifying several cases of over-clustering and recapitulating experimentally validated cell type definitions.

Secreted proteins play critical roles in cellular communication. Methods enabling concurrent measurement of cellular protein secretion, phenotypes and transcriptomes are still unavailable. Here we describe time-resolved assessment of protein secretion from single cells by sequencing (TRAPS-seq). Released proteins are trapped onto the cell surface and probed by oligonucleotide-barcoded antibodies before being simultaneously sequenced with transcriptomes in single cells. We demonstrate that TRAPS-seq helps unravel the phenotypic and transcriptional determinants of the secretion of pleiotropic TH1 cytokines (IFNγ, IL-2 and TNF) in activated T cells. In addition, we show that TRAPS-seq can be used to track the secretion of multiple cytokines over time, uncovering unique molecular signatures that govern the dynamics of single-cell cytokine secretions. Our results revealed that early central memory T cells with CD45RA expression (TCMRA) are important in both the production and maintenance of polyfunctional cytokines. TRAPS-seq presents a unique tool for seamless integration of secretomics measurements with multi-omics profiling in single cells.

Efficient methods for the generation of specific mutations enable the study of functional variations in natural populations and lead to advances in genetic engineering applications. Here, we present a new approach, mutagenesis by template-guided amplicon assembly (MEGAA), for the rapid construction of kilobase-sized DNA variants. With this method, many mutations can be generated at a time to a DNA template at more than 90% efficiency per target in a predictable manner. We devised a robust and iterative protocol for an open-source laboratory automation robot that enables desktop production and long-read sequencing validation of variants. Using this system, we demonstrated the construction of 31 natural SARS-CoV2 spike gene variants and 10 recoded Escherichia coli genome fragments, with each 4 kb region containing up to 150 mutations. Furthermore, 125 defined combinatorial adeno-associated virus-2 cap gene variants were easily built using the system, which exhibited viral packaging enhancements of up to 10-fold compared with wild type. Thus, the MEGAA platform enables generation of multi-site sequence variants quickly, cheaply, and in a scalable manner for diverse applications in biotechnology.

A substantial fraction of metabolic features remains undetermined in mass spectrometry (MS)-based metabolomics, and molecular formula annotation is the starting point for unraveling their chemical identities. Here we present bottom-up tandem MS (MS/MS) interrogation, a method for de novo formula annotation. Our approach prioritizes MS/MS-explainable formula candidates, implements machine-learned ranking and offers false discovery rate estimation. Compared with the mathematically exhaustive formula enumeration, our approach shrinks the formula candidate space by 42.8% on average. Method benchmarking on annotation accuracy was systematically carried out on reference MS/MS libraries and real metabolomics datasets. Applied on 155,321 recurrent unidentified spectra, our approach confidently annotated >5,000 novel molecular formulae absent from chemical databases. Beyond the level of individual metabolic features, we combined bottom-up MS/MS interrogation with global optimization to refine formula annotations while revealing peak interrelationships. This approach allowed the systematic annotation of 37 fatty acid amide molecules in human fecal data. All bioinformatics pipelines are available in a standalone software, BUDDY ( https://github.com/HuanLab/BUDDY ).

Spatiotemporal regulation of the cellular transcriptome is crucial for proper protein expression and cellular function. However, the intricate subcellular dynamics of RNA remain obscured due to the limitations of existing transcriptomics methods. Here, we report TEMPOmap-a method that uncovers subcellular RNA profiles across time and space at the single-cell level. TEMPOmap integrates pulse-chase metabolic labeling with highly multiplexed three-dimensional in situ sequencing to simultaneously profile the age and location of individual RNA molecules. Using TEMPOmap, we constructed the subcellular RNA kinetic landscape in various human cells from transcription and translocation to degradation. Clustering analysis of RNA kinetic parameters across single cells revealed 'kinetic gene clusters' whose expression patterns were shaped by multistep kinetic sculpting. Importantly, these kinetic gene clusters are functionally segregated, suggesting that subcellular RNA kinetics are differentially regulated in a cell-state- and cell-type-dependent manner. Spatiotemporally resolved transcriptomics provides a gateway to uncovering new spatiotemporal gene regulation principles.

Single-molecule localization microscopy in a typical wide-field setup has been widely used for investigating subcellular structures with super resolution; however, field-dependent aberrations restrict the field of view (FOV) to only tens of micrometers. Here, we present a deep-learning method for precise localization of spatially variant point emitters (FD-DeepLoc) over a large FOV covering the full chip of a modern sCMOS camera. Using a graphic processing unit-based vectorial point spread function (PSF) fitter, we can fast and accurately model the spatially variant PSF of a high numerical aperture objective in the entire FOV. Combined with deformable mirror-based optimal PSF engineering, we demonstrate high-accuracy three-dimensional single-molecule localization microscopy over a volume of ~180 × 180 × 5 μm3, allowing us to image mitochondria and nuclear pore complexes in entire cells in a single imaging cycle without hardware scanning; a 100-fold increase in throughput compared to the state of the art.

Spatially resolved genomic technologies have allowed us to study the physical organization of cells and tissues, and promise an understanding of local interactions between cells. However, it remains difficult to precisely align spatial observations across slices, samples, scales, individuals and technologies. Here, we propose a probabilistic model that aligns spatially-resolved samples onto a known or unknown common coordinate system (CCS) with respect to phenotypic readouts (for example, gene expression). Our method, Gaussian Process Spatial Alignment (GPSA), consists of a two-layer Gaussian process: the first layer maps observed samples' spatial locations onto a CCS, and the second layer maps from the CCS to the observed readouts. Our approach enables complex downstream spatially aware analyses that are impossible or inaccurate with unaligned data, including an analysis of variance, creation of a dense three-dimensional (3D) atlas from sparse two-dimensional (2D) slices or association tests across data modalities.

Jointly profiling the transcriptome, chromatin accessibility and other molecular properties of single cells offers a powerful way to study cellular diversity. Here we present MultiVI, a probabilistic model to analyze such multiomic data and leverage it to enhance single-modality datasets. MultiVI creates a joint representation that allows an analysis of all modalities included in the multiomic input data, even for cells for which one or more modalities are missing. It is available at scvi-tools.org .

Here we present a fluorescence microscope light path that enables imaging, during free behavior, of thousands of neurons in mice and hundreds of neurons in juvenile songbirds. The light path eliminates traditional illumination optics, allowing for head-mounted microscopes that have both a lower weight and a larger field of view (FOV) than previously possible. Using this light path, we designed two microscopes: one optimized for FOV (~4 mm FOV; 1.4 g), and the other optimized for weight (1.0 mm FOV; 1.0 g).

Phylogenetic models of molecular evolution are central to numerous biological applications spanning diverse timescales, from hundreds of millions of years involving orthologous proteins to just tens of days relating to single cells within an organism. A fundamental problem in these applications is estimating model parameters, for which maximum likelihood estimation is typically employed. Unfortunately, maximum likelihood estimation is a computationally expensive task, in some cases prohibitively so. To address this challenge, we here introduce CherryML, a broadly applicable method that achieves several orders of magnitude speedup by using a quantized composite likelihood over cherries in the trees. The massive speedup offered by our method should enable researchers to consider more complex and biologically realistic models than previously possible. Here we demonstrate CherryML's utility by applying it to estimate a general 400 × 400 rate matrix for residue-residue coevolution at contact sites in three-dimensional protein structures; we estimate that using current state-of-the-art methods such as the expectation-maximization algorithm for the same task would take >100,000 times longer.

In biological systems, spatial organization and function are interconnected. Here we present photoselective sequencing, a new method for genomic and epigenomic profiling within morphologically distinct regions. Starting with an intact biological specimen, photoselective sequencing uses targeted illumination to selectively unblock a photocaged fragment library, restricting the sequencing-based readout to microscopically identified spatial regions. We validate photoselective sequencing by measuring the chromatin accessibility profiles of fluorescently labeled cell types within the mouse brain and comparing with published data. Furthermore, by combining photoselective sequencing with a computational strategy for decomposing bulk accessibility profiles, we find that the oligodendrocyte-lineage-cell population is relatively enriched for oligodendrocyte-progenitor cells in the cortex versus the corpus callosum. Finally, we leverage photoselective sequencing at the subcellular scale to identify features of chromatin that are correlated with positioning at the nuclear periphery. These results collectively demonstrate that photoselective sequencing is a flexible and generalizable platform for exploring the interplay of spatial structures with genomic and epigenomic properties.

High-speed three-dimensional (3D) intravital imaging in animals is useful for studying transient subcellular interactions and functions in health and disease. Light-field microscopy (LFM) provides a computational solution for snapshot 3D imaging with low phototoxicity but is restricted by low resolution and reconstruction artifacts induced by optical aberrations, motion and noise. Here, we propose virtual-scanning LFM (VsLFM), a physics-based deep learning framework to increase the resolution of LFM up to the diffraction limit within a snapshot. By constructing a 40 GB high-resolution scanning LFM dataset across different species, we exploit physical priors between phase-correlated angular views to address the frequency aliasing problem. This enables us to bypass hardware scanning and associated motion artifacts. Here, we show that VsLFM achieves ultrafast 3D imaging of diverse processes such as the beating heart in embryonic zebrafish, voltage activity in Drosophila brains and neutrophil migration in the mouse liver at up to 500 volumes per second.

Most approaches to transcript quantification rely on fixed reference annotations; however, the transcriptome is dynamic and depending on the context, such static annotations contain inactive isoforms for some genes, whereas they are incomplete for others. Here we present Bambu, a method that performs machine-learning-based transcript discovery to enable quantification specific to the context of interest using long-read RNA-sequencing. To identify novel transcripts, Bambu estimates the novel discovery rate, which replaces arbitrary per-sample thresholds with a single, interpretable, precision-calibrated parameter. Bambu retains the full-length and unique read counts, enabling accurate quantification in presence of inactive isoforms. Compared to existing methods for transcript discovery, Bambu achieves greater precision without sacrificing sensitivity. We show that context-aware annotations improve quantification for both novel and known transcripts. We apply Bambu to quantify isoforms from repetitive HERVH-LTR7 retrotransposons in human embryonic stem cells, demonstrating the ability for context-specific transcript expression analysis.

Optical microscopy methods such as calcium and voltage imaging enable fast activity readout of large neuronal populations using light. However, the lack of corresponding advances in online algorithms has slowed progress in retrieving information about neural activity during or shortly after an experiment. This gap not only prevents the execution of real-time closed-loop experiments, but also hampers fast experiment-analysis-theory turnover for high-throughput imaging modalities. Reliable extraction of neural activity from fluorescence imaging frames at speeds compatible with indicator dynamics and imaging modalities poses a challenge. We therefore developed FIOLA, a framework for fluorescence imaging online analysis that extracts neuronal activity from calcium and voltage imaging movies at speeds one order of magnitude faster than state-of-the-art methods. FIOLA exploits algorithms optimized for parallel processing on GPUs and CPUs. We demonstrate reliable and scalable performance of FIOLA on both simulated and real calcium and voltage imaging datasets. Finally, we present an online experimental scenario to provide guidance in setting FIOLA parameters and to highlight the trade-offs of our approach.

The development of transgenic mouse models that express genes of interest in specific cell types has transformed our understanding of basic biology and disease. However, generating these models is time- and resource-intensive. Here we describe a model system, SELective Expression and Controlled Transduction In Vivo (SELECTIV), that enables efficient and specific expression of transgenes by coupling adeno-associated virus (AAV) vectors with Cre-inducible overexpression of the multi-serotype AAV receptor, AAVR. We demonstrate that transgenic AAVR overexpression greatly increases the efficiency of transduction of many diverse cell types, including muscle stem cells, which are normally refractory to AAV transduction. Superior specificity is achieved by combining Cre-mediated AAVR overexpression with whole-body knockout of endogenous Aavr, which is demonstrated in heart cardiomyocytes, liver hepatocytes and cholinergic neurons. The enhanced efficacy and exquisite specificity of SELECTIV has broad utility in development of new mouse model systems and expands the use of AAV for gene delivery in vivo.

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve 'active' regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain 'inactive' regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Metagenomic binning has revolutionized the study of uncultured microorganisms. Here we compare single- and multi-coverage binning on the same set of samples, and demonstrate that multi-coverage binning produces better results than single-coverage binning and identifies contaminant contigs and chimeric bins that other approaches miss. While resource expensive, multi-coverage binning is a superior approach and should always be performed over single-coverage binning.

Multiplexed antibody-based imaging enables the detailed characterization of molecular and cellular organization in tissues. Advances in the field now allow high-parameter data collection (>60 targets); however, considerable expertise and capital are needed to construct the antibody panels employed by these methods. Organ mapping antibody panels are community-validated resources that save time and money, increase reproducibility, accelerate discovery and support the construction of a Human Reference Atlas.

The Cell Tracking Challenge is an ongoing benchmarking initiative that has become a reference in cell segmentation and tracking algorithm development. Here, we present a significant number of improvements introduced in the challenge since our 2017 report. These include the creation of a new segmentation-only benchmark, the enrichment of the dataset repository with new datasets that increase its diversity and complexity, and the creation of a silver standard reference corpus based on the most competitive results, which will be of particular interest for data-hungry deep learning-based strategies. Furthermore, we present the up-to-date cell segmentation and tracking leaderboards, an in-depth analysis of the relationship between the performance of the state-of-the-art methods and the properties of the datasets and annotations, and two novel, insightful studies about the generalizability and the reusability of top-performing methods. These studies provide critical practical conclusions for both developers and users of traditional and machine learning-based cell segmentation and tracking algorithms.

Single-cell datasets are routinely collected to investigate changes in cellular state between control cells and the corresponding cells in a treatment condition, such as exposure to a drug or infection by a pathogen. To better understand heterogeneity in treatment response, it is desirable to deconvolve variations enriched in treated cells from those shared with controls. However, standard computational models of single-cell data are not designed to explicitly separate these variations. Here, we introduce contrastive variational inference (contrastiveVI; https://github.com/suinleelab/contrastiveVI ), a framework for deconvolving variations in treatment-control single-cell RNA sequencing (scRNA-seq) datasets into shared and treatment-specific latent variables. Using three treatment-control scRNA-seq datasets, we apply contrastiveVI to perform a variety of analysis tasks, including visualization, clustering and differential expression testing. We find that contrastiveVI consistently achieves results that agree with known ground truths and often highlights subtle phenomena that may be difficult to ascertain with standard workflows. We conclude by generalizing contrastiveVI to accommodate joint transcriptome and surface protein measurements.

Joint profiling of chromatin accessibility and gene expression in individual cells provides an opportunity to decipher enhancer-driven gene regulatory networks (GRNs). Here we present a method for the inference of enhancer-driven GRNs, called SCENIC+. SCENIC+ predicts genomic enhancers along with candidate upstream transcription factors (TFs) and links these enhancers to candidate target genes. To improve both recall and precision of TF identification, we curated and clustered a motif collection with more than 30,000 motifs. We benchmarked SCENIC+ on diverse datasets from different species, including human peripheral blood mononuclear cells, ENCODE cell lines, melanoma cell states and Drosophila retinal development. Next, we exploit SCENIC+ predictions to study conserved TFs, enhancers and GRNs between human and mouse cell types in the cerebral cortex. Finally, we use SCENIC+ to study the dynamics of gene regulation along differentiation trajectories and the effect of TF perturbations on cell state. SCENIC+ is available at scenicplus.readthedocs.io .

The fluorescent glutamate indicator iGluSnFR enables imaging of neurotransmission with genetic and molecular specificity. However, existing iGluSnFR variants exhibit low in vivo signal-to-noise ratios, saturating activation kinetics and exclusion from postsynaptic densities. Using a multiassay screen in bacteria, soluble protein and cultured neurons, we generated variants with improved signal-to-noise ratios and kinetics. We developed surface display constructs that improve iGluSnFR's nanoscopic localization to postsynapses. The resulting indicator iGluSnFR3 exhibits rapid nonsaturating activation kinetics and reports synaptic glutamate release with decreased saturation and increased specificity versus extrasynaptic signals in cultured neurons. Simultaneous imaging and electrophysiology at individual boutons in mouse visual cortex showed that iGluSnFR3 transients report single action potentials with high specificity. In vibrissal sensory cortex layer 4, we used iGluSnFR3 to characterize distinct patterns of touch-evoked feedforward input from thalamocortical boutons and both feedforward and recurrent input onto L4 cortical neuron dendritic spines.

Highly multiplexed imaging holds enormous promise for understanding how spatial context shapes the activity of the genome and its products at multiple length scales. Here, we introduce a deep learning framework called CAMPA (Conditional Autoencoder for Multiplexed Pixel Analysis), which uses a conditional variational autoencoder to learn representations of molecular pixel profiles that are consistent across heterogeneous cell populations and experimental perturbations. Clustering these pixel-level representations identifies consistent subcellular landmarks, which can be quantitatively compared in terms of their size, shape, molecular composition and relative spatial organization. Using high-resolution multiplexed immunofluorescence, this reveals how subcellular organization changes upon perturbation of RNA synthesis, RNA processing or cell size, and uncovers links between the molecular composition of membraneless organelles and cell-to-cell variability in bulk RNA synthesis rates. By capturing interpretable cellular phenotypes, we anticipate that CAMPA will greatly accelerate the systematic mapping of multiscale atlases of biological organization to identify the rules by which context shapes physiology and disease.

Discovery of off-target CRISPR-Cas activity in patient-derived cells and animal models is crucial for genome editing applications, but currently exhibits low sensitivity. We demonstrate that inhibition of DNA-dependent protein kinase catalytic subunit accumulates the repair protein MRE11 at CRISPR-Cas-targeted sites, enabling high-sensitivity mapping of off-target sites to positions of MRE11 binding using chromatin immunoprecipitation followed by sequencing. This technique, termed DISCOVER-Seq+, discovered up to fivefold more CRISPR off-target sites in immortalized cell lines, primary human cells and mice compared with previous methods. We demonstrate applicability to ex vivo knock-in of a cancer-directed transgenic T cell receptor in primary human T cells and in vivo adenovirus knock-out of cardiovascular risk gene PCSK9 in mice. Thus, DISCOVER-Seq+ is, to our knowledge, the most sensitive method to-date for discovering off-target genome editing in vivo.

Single-particle cryogenic electron microscopy (cryo-EM) allows reconstruction of high-resolution structures of proteins in different conformations. Protein function often involves transient functional conformations, which can be resolved using time-resolved cryo-EM (trEM). In trEM, reactions are arrested after a defined delay time by rapid vitrification of protein solution on the EM grid. Despite the increasing interest in trEM among the cryo-EM community, making trEM samples with a time resolution below 100 ms remains challenging. Here we report the design and the realization of a time-resolved cryo-plunger that combines a droplet-based microfluidic mixer with a laser-induced generator of microjets that allows rapid reaction initiation and plunge-freezing of cryo-EM grids. Using this approach, a time resolution of 5 ms was achieved and the protein density map was reconstructed to a resolution of 2.1 Å. trEM experiments on GroEL:GroES chaperonin complex resolved the kinetics of the complex formation and visualized putative short-lived conformations of GroEL-ATP complex.

Genetically encoded voltage indicators (GEVIs) enable optical recording of electrical signals in the brain, providing subthreshold sensitivity and temporal resolution not possible with calcium indicators. However, one- and two-photon voltage imaging over prolonged periods with the same GEVI has not yet been demonstrated. Here, we report engineering of ASAP family GEVIs to enhance photostability by inversion of the fluorescence-voltage relationship. Two of the resulting GEVIs, ASAP4b and ASAP4e, respond to 100-mV depolarizations with ≥180% fluorescence increases, compared with the 50% fluorescence decrease of the parental ASAP3. With standard microscopy equipment, ASAP4e enables single-trial detection of spikes in mice over the course of minutes. Unlike GEVIs previously used for one-photon voltage recordings, ASAP4b and ASAP4e also perform well under two-photon illumination. By imaging voltage and calcium simultaneously, we show that ASAP4b and ASAP4e can identify place cells and detect voltage spikes with better temporal resolution than commonly used calcium indicators. Thus, ASAP4b and ASAP4e extend the capabilities of voltage imaging to standard one- and two-photon microscopes while improving the duration of voltage recordings.

The brain is a complex tissue whose function relies on coordinated anatomical and molecular features. However, the molecular annotation of the spatial organization of the brain is currently insufficient. Here, we describe microfluidic indexing-based spatial assay for transposase-accessible chromatin and RNA-sequencing (MISAR-seq), a method for spatially resolved joint profiling of chromatin accessibility and gene expression. By applying MISAR-seq to the developing mouse brain, we study tissue organization and spatiotemporal regulatory logics during mouse brain development.

Rapid, highly multiplexed, nondestructive imaging that spans the molecular to the supra-cellular scale would be a powerful tool for tissue analysis. However, the physical constraints of established imaging methods limit the simultaneous improvement of these parameters. Whole-organism to atomic-level imaging is possible with tissue-penetrant, picometer-wavelength X-rays. To enable highly multiplexed X-ray imaging, we developed multielement Z-tag X-ray fluorescence (MEZ-XRF) that can operate at kHz speeds when combined with signal amplification by exchange reaction (SABER)-amplified Z-tag reagents. We demonstrated parallel imaging of 20 Z-tag or SABER Z-tag reagents at subcellular resolution in cell lines and multiple human tissues. We benchmarked MEZ-XRF against imaging mass cytometry and demonstrated the nondestructive multiscale repeat imaging capabilities of MEZ-XRF with rapid tissue overview scans, followed by slower, more sensitive imaging of low-abundance markers such as immune checkpoint proteins. The unique multiscale, nondestructive nature of MEZ-XRF, combined with SABER Z-tags for high sensitivity or enhanced speed, enables highly multiplexed bioimaging across biological scales.

Evolutionarily conserved signaling pathways are essential for early embryogenesis, and reducing or abolishing their activity leads to characteristic developmental defects. Classification of phenotypic defects can identify the underlying signaling mechanisms, but this requires expert knowledge and the classification schemes have not been standardized. Here we use a machine learning approach for automated phenotyping to train a deep convolutional neural network, EmbryoNet, to accurately identify zebrafish signaling mutants in an unbiased manner. Combined with a model of time-dependent developmental trajectories, this approach identifies and classifies with high precision phenotypic defects caused by loss of function of the seven major signaling pathways relevant for vertebrate development. Our classification algorithms have wide applications in developmental biology and robustly identify signaling defects in evolutionarily distant species. Furthermore, using automated phenotyping in high-throughput drug screens, we show that EmbryoNet can resolve the mechanism of action of pharmaceutical substances. As part of this work, we freely provide more than 2 million images that were used to train and test EmbryoNet.

Gene regulatory networks (GRNs) are key determinants of cell function and identity and are dynamically rewired during development and disease. Despite decades of advancement, challenges remain in GRN inference, including dynamic rewiring, causal inference, feedback loop modeling and context specificity. To address these challenges, we develop Dictys, a dynamic GRN inference and analysis method that leverages multiomic single-cell assays of chromatin accessibility and gene expression, context-specific transcription factor footprinting, stochastic process network and efficient probabilistic modeling of single-cell RNA-sequencing read counts. Dictys improves GRN reconstruction accuracy and reproducibility and enables the inference and comparative analysis of context-specific and dynamic GRNs across developmental contexts. Dictys' network analyses recover unique insights in human blood and mouse skin development with cell-type-specific and dynamic GRNs. Its dynamic network visualizations enable time-resolved discovery and investigation of developmental driver transcription factors and their regulated targets. Dictys is available as a free, open-source and user-friendly Python package.

Genetically encoded calcium indicators (GECIs) are indispensable tools for real-time monitoring of intracellular calcium signals and cellular activities in living organisms. Current GECIs face the challenge of suboptimal peak signal-to-baseline ratio (SBR) with limited resolution for reporting subtle calcium transients. We report herein the development of a suite of calcium sensors, designated NEMO, with fast kinetics and wide dynamic ranges (>100-fold). NEMO indicators report Ca2+ transients with peak SBRs around 20-fold larger than the top-of-the-range GCaMP6 series. NEMO sensors further enable the quantification of absolution calcium concentration with ratiometric or photochromic imaging. Compared with GCaMP6s, NEMOs could detect single action potentials in neurons with a peak SBR two times higher and a median peak SBR four times larger in vivo, thereby outperforming most existing state-of-the-art GECIs. Given their high sensitivity and resolution to report intracellular Ca2+ signals, NEMO sensors may find broad applications in monitoring neuronal activities and other Ca2+-modulated physiological processes in both mammals and plants.

Learning is thought to involve changes in glutamate receptors at synapses, submicron structures that mediate communication between neurons in the central nervous system. Due to their small size and high density, synapses are difficult to resolve in vivo, limiting our ability to directly relate receptor dynamics to animal behavior. Here we developed a combination of computational and biological methods to overcome these challenges. First, we trained a deep-learning image-restoration algorithm that combines the advantages of ex vivo super-resolution and in vivo imaging modalities to overcome limitations specific to each optical system. When applied to in vivo images from transgenic mice expressing fluorescently labeled glutamate receptors, this restoration algorithm super-resolved synapses, enabling the tracking of behavior-associated synaptic plasticity with high spatial resolution. This method demonstrates the capabilities of image enhancement to learn from ex vivo data and imaging techniques to improve in vivo imaging resolution.

Cryogenic-electron tomography enables the visualization of cellular environments in extreme detail, however, tools to analyze the full amount of information contained within these densely packed volumes are still needed. Detailed analysis of macromolecules through subtomogram averaging requires particles to first be localized within the tomogram volume, a task complicated by several factors including a low signal to noise ratio and crowding of the cellular space. Available methods for this task suffer either from being error prone or requiring manual annotation of training data. To assist in this crucial particle picking step, we present TomoTwin: an open source general picking model for cryogenic-electron tomograms based on deep metric learning. By embedding tomograms in an information-rich, high-dimensional space that separates macromolecules according to their three-dimensional structure, TomoTwin allows users to identify proteins in tomograms de novo without manually creating training data or retraining the network to locate new proteins.

Droplet-based single-cell assays, including single-cell RNA sequencing (scRNA-seq), single-nucleus RNA sequencing (snRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq), generate considerable background noise counts, the hallmark of which is nonzero counts in cell-free droplets and off-target gene expression in unexpected cell types. Such systematic background noise can lead to batch effects and spurious differential gene expression results. Here we develop a deep generative model based on the phenomenology of noise generation in droplet-based assays. The proposed model accurately distinguishes cell-containing droplets from cell-free droplets, learns the background noise profile and provides noise-free quantification in an end-to-end fashion. We implement this approach in the scalable and robust open-source software package CellBender. Analysis of simulated data demonstrates that CellBender operates near the theoretically optimal denoising limit. Extensive evaluations using real datasets and experimental benchmarks highlight enhanced concordance between droplet-based single-cell data and established gene expression patterns, while the learned background noise profile provides evidence of degraded or uncaptured cell types.

Widefield microscopy can provide optical access to multi-millimeter fields of view and thousands of neurons in mammalian brains at video rate. However, tissue scattering and background contamination results in signal deterioration, making the extraction of neuronal activity challenging, laborious and time consuming. Here we present our deep-learning-based widefield neuron finder (DeepWonder), which is trained by simulated functional recordings and effectively works on experimental data to achieve high-fidelity neuronal extraction. Equipped with systematic background contribution priors, DeepWonder conducts neuronal inference with an order-of-magnitude-faster speed and improved accuracy compared with alternative approaches. DeepWonder removes background contaminations and is computationally efficient. Specifically, DeepWonder accomplishes 50-fold signal-to-background ratio enhancement when processing terabytes-scale cortex-wide functional recordings, with over 14,000 neurons extracted in 17 h.

Major aims of single-cell proteomics include increasing the consistency, sensitivity and depth of protein quantification, especially for proteins and modifications of biological interest. Here, to simultaneously advance all these aims, we developed prioritized Single-Cell ProtEomics (pSCoPE). pSCoPE consistently analyzes thousands of prioritized peptides across all single cells (thus increasing data completeness) while maximizing instrument time spent analyzing identifiable peptides, thus increasing proteome depth. These strategies increased the sensitivity, data completeness and proteome coverage over twofold. The gains enabled quantifying protein variation in untreated and lipopolysaccharide-treated primary macrophages. Within each condition, proteins covaried within functional sets, including phagosome maturation and proton transport, similarly across both treatment conditions. This covariation is coupled to phenotypic variability in endocytic activity. pSCoPE also enabled quantifying proteolytic products, suggesting a gradient of cathepsin activities within a treatment condition. pSCoPE is freely available and widely applicable, especially for analyzing proteins of interest without sacrificing proteome coverage. Support for pSCoPE is available at http://scp.slavovlab.net/pSCoPE .

Modeling flexible macromolecules is one of the foremost challenges in single-particle cryogenic-electron microscopy (cryo-EM), with the potential to illuminate fundamental questions in structural biology. We introduce Three-Dimensional Flexible Refinement (3DFlex), a motion-based neural network model for continuous molecular heterogeneity for cryo-EM data. 3DFlex exploits knowledge that conformational variability of a protein is often the result of physical processes that transport density over space and tend to preserve local geometry. From two-dimensional image data, 3DFlex enables the determination of high-resolution 3D density, and provides an explicit model of a flexible protein's motion over its conformational landscape. Experimentally, for large molecular machines (tri-snRNP spliceosome complex, translocating ribosome) and small flexible proteins (TRPV1 ion channel, αVβ8 integrin, SARS-CoV-2 spike), 3DFlex learns nonrigid molecular motions while resolving details of moving secondary structure elements. 3DFlex can improve 3D density resolution beyond the limits of existing methods because particle images contribute coherent signal over the conformational landscape.

Imaging mass cytometry (IMC) is a highly multiplexed, antibody-based imaging method that captures heterogeneous spatial protein expression patterns at subcellular resolution. Here we report the extension of IMC to low-abundance markers through incorporation of the DNA-based signal amplification by exchange reaction, immuno-SABER. We applied SABER-IMC to image the tumor immune microenvironment in human melanoma by simultaneous imaging of 18 markers with immuno-SABER and 20 markers without amplification. SABER-IMC enabled the identification of immune cell phenotypic markers, such as T cell co-receptors and their ligands, that are not detectable with IMC.

Connectomes of human cortical gray matter require high-contrast homogeneously stained samples sized at least 2 mm on a side, and a mouse whole-brain connectome requires samples sized at least 5-10 mm on a side. Here we report en bloc staining and embedding protocols for these and other applications, removing a key obstacle for connectomic analyses at the mammalian whole-brain level.

Recently, various small Cas9 orthologs and variants have been reported for use in in vivo delivery applications. Although small Cas9s are particularly suited for this purpose, selecting the most optimal small Cas9 for use at a specific target sequence continues to be challenging. Here, to this end, we have systematically compared the activities of 17 small Cas9s for thousands of target sequences. For each small Cas9, we have characterized the protospacer adjacent motif and determined optimal single guide RNA expression formats and scaffold sequence. High-throughput comparative analyses revealed distinct high- and low-activity groups of small Cas9s. We also developed DeepSmallCas9, a set of computational models predicting the activities of the small Cas9s at matched and mismatched target sequences. Together, this analysis and these computational models provide a useful guide for researchers to select the most suitable small Cas9 for specific applications.

Our understanding of nerve regeneration can be enhanced by delineating its underlying molecular activities at single-neuron resolution in model organisms such as Caenorhabditis elegans. Existing cell isolation techniques cannot isolate neurons with specific regeneration phenotypes from C. elegans. We present femtosecond laser microdissection (fs-LM), a single-cell isolation method that dissects specific cells directly from living tissue by leveraging the micrometer-scale precision of fs-laser ablation. We show that fs-LM facilitates sensitive and specific gene expression profiling by single-cell RNA sequencing (scRNA-seq), while mitigating the stress-related transcriptional artifacts induced by tissue dissociation. scRNA-seq of fs-LM isolated regenerating neurons revealed transcriptional programs that are correlated with either successful or failed regeneration in wild-type and dlk-1 (0) animals, respectively. This method also allowed studying heterogeneity displayed by the same type of neuron and found gene modules with expression patterns correlated with axon regrowth rate. Our results establish fs-LM as a spatially resolved single-cell isolation method for phenotype-to-genotype mapping.

Photoactivatable drugs and peptides can drive quantitative studies into receptor signaling with high spatiotemporal precision, yet few are compatible with behavioral studies in mammals. We developed CNV-Y-DAMGO-a caged derivative of the mu opioid receptor-selective peptide agonist DAMGO. Photoactivation in the mouse ventral tegmental area produced an opioid-dependent increase in locomotion within seconds of illumination. These results demonstrate the power of in vivo photopharmacology for dynamic studies into animal behavior.

The ability to quantify structural changes of the endoplasmic reticulum (ER) is crucial for understanding the structure and function of this organelle. However, the rapid movement and complex topology of ER networks make this challenging. Here, we construct a state-of-the-art semantic segmentation method that we call ERnet for the automatic classification of sheet and tubular ER domains inside individual cells. Data are skeletonized and represented by connectivity graphs, enabling precise and efficient quantification of network connectivity. ERnet generates metrics on topology and integrity of ER structures and quantifies structural change in response to genetic or metabolic manipulation. We validate ERnet using data obtained by various ER-imaging methods from different cell types as well as ground truth images of synthetic ER structures. ERnet can be deployed in an automatic high-throughput and unbiased fashion and identifies subtle changes in ER phenotypes that may inform on disease progression and response to therapy.

Advancements in sequencing technologies and assembly methods enable the regular production of high-quality genome assemblies characterizing complex regions. However, challenges remain in efficiently interpreting variation at various scales, from smaller tandem repeats to megabase rearrangements, across many human genomes. We present a PanGenome Research Tool Kit (PGR-TK) enabling analyses of complex pangenome structural and haplotype variation at multiple scales. We apply the graph decomposition methods in PGR-TK to the class II major histocompatibility complex demonstrating the importance of the human pangenome for analyzing complicated regions. Moreover, we investigate the Y-chromosome genes, DAZ1/DAZ2/DAZ3/DAZ4, of which structural variants have been linked to male infertility, and X-chromosome genes OPN1LW and OPN1MW linked to eye disorders. We further showcase PGR-TK across 395 complex repetitive medically important genes. This highlights the power of PGR-TK to resolve complex variation in regions of the genome that were previously too complex to analyze.

Single-molecule Förster-resonance energy transfer (smFRET) experiments allow the study of biomolecular structure and dynamics in vitro and in vivo. We performed an international blind study involving 19 laboratories to assess the uncertainty of FRET experiments for proteins with respect to the measured FRET efficiency histograms, determination of distances, and the detection and quantification of structural dynamics. Using two protein systems with distinct conformational changes and dynamics, we obtained an uncertainty of the FRET efficiency ≤0.06, corresponding to an interdye distance precision of ≤2 Å and accuracy of ≤5 Å. We further discuss the limits for detecting fluctuations in this distance range and how to identify dye perturbations. Our work demonstrates the ability of smFRET experiments to simultaneously measure distances and avoid the averaging of conformational dynamics for realistic protein systems, highlighting its importance in the expanding toolbox of integrative structural biology.

Characterizing multifaceted individual differences in brain function using neuroimaging is central to biomarker discovery in neuroscience. We provide an integrative toolbox, Reliability eXplorer (ReX), to facilitate the examination of individual variation and reliability as well as the effective direction for optimization of measuring individual differences in biomarker discovery. We also illustrate gradient flows, a two-dimensional field map-based approach to identifying and representing the most effective direction for optimization when measuring individual differences, which is implemented in ReX.

Prime editors have a broad range of potential research and clinical applications. However, methods to delineate their genome-wide editing activities have generally relied on indirect genome-wide editing assessments or the computational prediction of near-cognate sequences. Here we describe a genome-wide approach for the identification of potential prime editor off-target sites, which we call PE-tag. This method relies on the attachment or insertion of an amplification tag at sites of prime editor activity to allow their identification. PE-tag enables genome-wide profiling of off-target sites in vitro using extracted genomic DNA, in mammalian cell lines and in the adult mouse liver. PE-tag components can be delivered in a variety of formats for off-target site detection. Our studies are consistent with the high specificity previously described for prime editor systems, but we find that off-target editing rates are influenced by prime editing guide RNA design. PE-tag represents an accessible, rapid and sensitive approach for the genome-wide identification of prime editor activity and the evaluation of prime editor safety.

Interferometric scattering (iSCAT) microscopy is a label-free optical method capable of detecting single proteins, localizing their binding positions with nanometer precision, and measuring their mass. In the ideal case, iSCAT is limited by shot noise such that collection of more photons should extend its detection sensitivity to biomolecules of arbitrarily low mass. However, a number of technical noise sources combined with speckle-like background fluctuations have restricted the detection limit in iSCAT. Here, we show that an unsupervised machine learning isolation forest algorithm for anomaly detection pushes the mass sensitivity limit by a factor of 4 to below 10 kDa. We implement this scheme both with a user-defined feature matrix and a self-supervised FastDVDNet and validate our results with correlative fluorescence images recorded in total internal reflection mode. Our work opens the door to optical investigations of small traces of biomolecules and disease markers such as α-synuclein, chemokines and cytokines.

Hierarchical organization of integral membrane proteins (IMP) and lipids at the membrane is essential for regulating myriad downstream signaling. A quantitative understanding of these processes requires both detections of oligomeric organization of IMPs and lipids directly from intact membranes and determination of key membrane components and properties that regulate them. Addressing this, we have developed a platform that enables native mass spectrometry (nMS) analysis of IMP-lipid complexes directly from intact and customizable lipid membranes. Both the lipid composition and membrane properties (such as curvature, tension, and fluidity) of these bilayers can be precisely customized to a target membrane. Subsequent direct nMS analysis of these intact proteolipid vesicles can yield the oligomeric states of the embedded IMPs, identify bound lipids, and determine the membrane properties that can regulate the observed IMP-lipid organization. Applying this method, we show how lipid binding regulates neurotransmitter release and how membrane composition regulates the functional oligomeric state of a transporter.

BigNeuron is an open community bench-testing platform with the goal of setting open standards for accurate and fast automatic neuron tracing. We gathered a diverse set of image volumes across several species that is representative of the data obtained in many neuroscience laboratories interested in neuron tracing. Here, we report generated gold standard manual annotations for a subset of the available imaging datasets and quantified tracing quality for 35 automatic tracing algorithms. The goal of generating such a hand-curated diverse dataset is to advance the development of tracing algorithms and enable generalizable benchmarking. Together with image quality features, we pooled the data in an interactive web application that enables users and developers to perform principal component analysis, t-distributed stochastic neighbor embedding, correlation and clustering, visualization of imaging and tracing data, and benchmarking of automatic tracing algorithms in user-defined data subsets. The image quality metrics explain most of the variance in the data, followed by neuromorphological features related to neuron size. We observed that diverse algorithms can provide complementary information to obtain accurate results and developed a method to iteratively combine methods and generate consensus reconstructions. The consensus trees obtained provide estimates of the neuron structure ground truth that typically outperform single algorithms in noisy datasets. However, specific algorithms may outperform the consensus tree strategy in specific imaging conditions. Finally, to aid users in predicting the most accurate automatic tracing results without manual annotations for comparison, we used support vector machine regression to predict reconstruction quality given an image volume and a set of automatic tracings.