Showing papers in "Nature Methods in 2018"

Accurate detection of complex structural variations using single-molecule sequencing.

Abstract: Structural variations are the greatest source of genetic variation, but they remain poorly understood because of technological limitations. Single-molecule long-read sequencing has the potential to dramatically advance the field, although high error rates are a challenge with existing methods. Addressing this need, we introduce open-source methods for long-read alignment (NGMLR; https://github.com/philres/ngmlr ) and structural variant identification (Sniffles; https://github.com/fritzsedlazeck/Sniffles ) that provide unprecedented sensitivity and precision for variant detection, even in repeat-rich regions and for complex nested events that can have substantial effects on human health. In several long-read datasets, including healthy and cancerous human genomes, we discovered thousands of novel variants and categorized systematic errors in short-read approaches. NGMLR and Sniffles can automatically filter false events and operate on low-coverage data, thereby reducing the high costs that have hindered the application of long reads in clinical and research settings.

Deep generative modeling for single-cell transcriptomics.

Abstract: Single-cell transcriptome measurements can reveal unexplored biological diversity, but they suffer from technical noise and bias that must be modeled to account for the resulting uncertainty in downstream analyses. Here we introduce single-cell variational inference (scVI), a ready-to-use scalable framework for the probabilistic representation and analysis of gene expression in single cells ( https://github.com/YosefLab/scVI ). scVI uses stochastic optimization and deep neural networks to aggregate information across similar cells and genes and to approximate the distributions that underlie observed expression values, while accounting for batch effects and limited sensitivity. We used scVI for a range of fundamental analysis tasks including batch correction, visualization, clustering, and differential expression, and achieved high accuracy for each task.

Species-level functional profiling of metagenomes and metatranscriptomes.

Abstract: Functional profiles of microbial communities are typically generated using comprehensive metagenomic or metatranscriptomic sequence read searches, which are time-consuming, prone to spurious mapping, and often limited to community-level quantification. We developed HUMAnN2, a tiered search strategy that enables fast, accurate, and species-resolved functional profiling of host-associated and environmental communities. HUMAnN2 identifies a community's known species, aligns reads to their pangenomes, performs translated search on unclassified reads, and finally quantifies gene families and pathways. Relative to pure translated search, HUMAnN2 is faster and produces more accurate gene family profiles. We applied HUMAnN2 to study clinal variation in marine metabolism, ecological contribution patterns among human microbiome pathways, variation in species' genomic versus transcriptional contributions, and strain profiling. Further, we introduce 'contributional diversity' to explain patterns of ecological assembly across different microbial community types.

Strelka2: fast and accurate calling of germline and somatic variants.

Abstract: We describe Strelka2 ( https://github.com/Illumina/strelka ), an open-source small-variant-calling method for research and clinical germline and somatic sequencing applications. Strelka2 introduces a novel mixture-model-based estimation of insertion/deletion error parameters from each sample, an efficient tiered haplotype-modeling strategy, and a normal sample contamination model to improve liquid tumor analysis. For both germline and somatic calling, Strelka2 substantially outperformed the current leading tools in terms of both variant-calling accuracy and computing cost.

Highly parallel direct RNA sequencing on an array of nanopores.

Abstract: Direct sequencing of RNA molecules in real time using nanopores allows for the detection of splice variants and hold promises for profiling RNA modifications.

Statistics versus machine learning

Abstract: Statistics draws population inferences from a sample, and machine learning finds generalizable predictive patterns.

Bioconda: sustainable and comprehensive software distribution for the life sciences.

Abstract: We present Bioconda (https://bioconda.github.io), a distribution of bioinformatics software for the lightweight, multi-platform and language-agnostic package manager Conda. Currently, Bioconda offers a collection of over 3000 software packages, which is continuously maintained, updated, and extended by a growing global community of more than 200 contributors. Bioconda improves analysis reproducibility by allowing users to define isolated environments with defined software versions, all of which are easily installed and managed without administrative privileges.

Content-aware image restoration: pushing the limits of fluorescence microscopy.

Abstract: Fluorescence microscopy is a key driver of discoveries in the life sciences, with observable phenomena being limited by the optics of the microscope, the chemistry of the fluorophores, and the maximum photon exposure tolerated by the sample. These limits necessitate trade-offs between imaging speed, spatial resolution, light exposure, and imaging depth. In this work we show how content-aware image restoration based on deep learning extends the range of biological phenomena observable by microscopy. We demonstrate on eight concrete examples how microscopy images can be restored even if 60-fold fewer photons are used during acquisition, how near isotropic resolution can be achieved with up to tenfold under-sampling along the axial direction, and how tubular and granular structures smaller than the diffraction limit can be resolved at 20-times-higher frame rates compared to state-of-the-art methods. All developed image restoration methods are freely available as open source software in Python, FIJI, and KNIME. Content-aware image restoration (CARE) uses deep learning to improve microscopy images. CARE bypasses the trade-offs between imaging speed, resolution, and maximal light exposure that limit fluorescence imaging to enable discovery.

SAVER: gene expression recovery for single-cell RNA sequencing.

Abstract: In single-cell RNA sequencing (scRNA-seq) studies, only a small fraction of the transcripts present in each cell are sequenced. This leads to unreliable quantification of genes with low or moderate expression, which hinders downstream analysis. To address this challenge, we developed SAVER (single-cell analysis via expression recovery), an expression recovery method for unique molecule index (UMI)-based scRNA-seq data that borrows information across genes and cells to provide accurate expression estimates for all genes.

Bias, robustness and scalability in single-cell differential expression analysis

Abstract: Many methods have been used to determine differential gene expression from single-cell RNA (scRNA)-seq data. We evaluated 36 approaches using experimental and synthetic data and found considerable differences in the number and characteristics of the genes that are called differentially expressed. Prefiltering of lowly expressed genes has important effects, particularly for some of the methods developed for bulk RNA-seq data analysis. However, we found that bulk RNA-seq analysis methods do not generally perform worse than those developed specifically for scRNA-seq. We also present conquer, a repository of consistently processed, analysis-ready public scRNA-seq data sets that is aimed at simplifying method evaluation and reanalysis of published results. Each data set provides abundance estimates for both genes and transcripts, as well as quality control and exploratory analysis reports.

scmap: projection of single-cell RNA-seq data across data sets

Abstract: Single-cell RNA-seq (scRNA-seq) allows researchers to define cell types on the basis of unsupervised clustering of the transcriptome. However, differences in experimental methods and computational analyses make it challenging to compare data across experiments. Here we present scmap (http://bioconductor.org/packages/scmap; web version at http://www.sanger.ac.uk/science/tools/scmap), a method for projecting cells from an scRNA-seq data set onto cell types or individual cells from other experiments.

Inferring single-trial neural population dynamics using sequential auto-encoders.

Abstract: Neuroscience is experiencing a revolution in which simultaneous recording of thousands of neurons is revealing population dynamics that are not apparent from single-neuron responses. This structure is typically extracted from data averaged across many trials, but deeper understanding requires studying phenomena detected in single trials, which is challenging due to incomplete sampling of the neural population, trial-to-trial variability, and fluctuations in action potential timing. We introduce latent factor analysis via dynamical systems, a deep learning method to infer latent dynamics from single-trial neural spiking data. When applied to a variety of macaque and human motor cortical datasets, latent factor analysis via dynamical systems accurately predicts observed behavioral variables, extracts precise firing rate estimates of neural dynamics on single trials, infers perturbations to those dynamics that correlate with behavioral choices, and combines data from non-overlapping recording sessions spanning months to improve inference of underlying dynamics.

Acoustic tweezers for the life sciences

Abstract: Acoustic tweezers are a versatile set of tools that use sound waves to manipulate bioparticles ranging from nanometer-sized extracellular vesicles to millimeter-sized multicellular organisms. Over the past several decades, the capabilities of acoustic tweezers have expanded from simplistic particle trapping to precise rotation and translation of cells and organisms in three dimensions. Recent advances have led to reconfigured acoustic tweezers that are capable of separating, enriching, and patterning bioparticles in complex solutions. Here, we review the history and fundamentals of acoustic-tweezer technology and summarize recent breakthroughs. This review discusses and contrasts different acoustic-tweezer technologies and their applications in biology.

A comparison of methods to assess cell mechanical properties.

Abstract: The mechanical properties of cells influence their cellular and subcellular functions, including cell adhesion, migration, polarization, and differentiation, as well as organelle organization and trafficking inside the cytoplasm. Yet reported values of cell stiffness and viscosity vary substantially, which suggests differences in how the results of different methods are obtained or analyzed by different groups. To address this issue and illustrate the complementarity of certain approaches, here we present, analyze, and critically compare measurements obtained by means of some of the most widely used methods for cell mechanics: atomic force microscopy, magnetic twisting cytometry, particle-tracking microrheology, parallel-plate rheometry, cell monolayer rheology, and optical stretching. These measurements highlight how elastic and viscous moduli of MCF-7 breast cancer cells can vary 1,000-fold and 100-fold, respectively. We discuss the sources of these variations, including the level of applied mechanical stress, the rate of deformation, the geometry of the probe, the location probed in the cell, and the extracellular microenvironment.

STED super-resolved microscopy.

Abstract: Stimulated emission depletion (STED) microscopy provides subdiffraction resolution while preserving useful aspects of fluorescence microscopy, such as optical sectioning, and molecular specificity and sensitivity However, sophisticated microscopy architectures and high illumination intensities have limited STED microscopy's widespread use in the past Here we summarize the progress that is mitigating these problems and giving substantial momentum to STED microscopy applications We discuss the future of this method in regard to spatiotemporal limits, live-cell imaging and combination with spectroscopy Advances in these areas may elevate STED microscopy to a standard method for imaging in the life sciences

Qiita: rapid, web-enabled microbiome meta-analysis.

Abstract: Multi-omic insights into microbiome function and composition typically advance one study at a time. However, in order for relationships across studies to be fully understood, data must be aggregated into meta-analyses. This makes it possible to generate new hypotheses by finding features that are reproducible across biospecimens and data layers. Qiita dramatically accelerates such integration tasks in a web-based microbiome-comparison platform, which we demonstrate with Human Microbiome Project and Integrative Human Microbiome Project (iHMP) data.

Deep generative models of genetic variation capture the effects of mutations.

Abstract: The functions of proteins and RNAs are defined by the collective interactions of many residues, and yet most statistical models of biological sequences consider sites nearly independently Recent approaches have demonstrated benefits of including interactions to capture pairwise covariation, but leave higher-order dependencies out of reach Here we show how it is possible to capture higher-order, context-dependent constraints in biological sequences via latent variable models with nonlinear dependencies We found that DeepSequence ( https://githubcom/debbiemarkslab/DeepSequence ), a probabilistic model for sequence families, predicted the effects of mutations across a variety of deep mutational scanning experiments substantially better than existing methods based on the same evolutionary data The model, learned in an unsupervised manner solely on the basis of sequence information, is grounded with biologically motivated priors, reveals the latent organization of sequence families, and can be used to explore new parts of sequence space

Spatial organization of the somatosensory cortex revealed by osmFISH

Abstract: Global efforts to create a molecular census of the brain using single-cell transcriptomics are producing a large catalog of molecularly defined cell types. However, spatial information is lacking and new methods are needed to map a large number of cell type-specific markers simultaneously on large tissue areas. Here, we describe a cyclic single-molecule fluorescence in situ hybridization methodology and define the cellular organization of the somatosensory cortex.

Label-free prediction of three-dimensional fluorescence images from transmitted-light microscopy.

Abstract: Understanding cells as integrated systems is central to modern biology. Although fluorescence microscopy can resolve subcellular structure in living cells, it is expensive, is slow, and can damage cells. We present a label-free method for predicting three-dimensional fluorescence directly from transmitted-light images and demonstrate that it can be used to generate multi-structure, integrated images. The method can also predict immunofluorescence (IF) from electron micrograph (EM) inputs, extending the potential applications. Convolutional neural networks enable prediction of fluorescently labeled structures from three-dimensional time-lapse transmitted-light images. Applications include multiplexed long time-lapse imaging and prediction of fluorescence in electron micrographs.

SpatialDE: identification of spatially variable genes

Abstract: Technological advances have made it possible to measure spatially resolved gene expression at high throughput. However, methods to analyze these data are not established. Here we describe SpatialDE, a statistical test to identify genes with spatial patterns of expression variation from multiplexed imaging or spatial RNA-sequencing data. SpatialDE also implements 'automatic expression histology', a spatial gene-clustering approach that enables expression-based tissue histology.

Identifying metabolites by integrating metabolome databases with mass spectrometry cheminformatics

Abstract: Novel metabolites distinct from canonical pathways can be identified through the integration of three cheminformatics tools: BinVestigate, which queries the BinBase gas chromatography-mass spectrometry (GC-MS) metabolome database to match unknowns with biological metadata across over 110,000 samples; MS-DIAL 2.0, a software tool for chromatographic deconvolution of high-resolution GC-MS or liquid chromatography-mass spectrometry (LC-MS); and MS-FINDER 2.0, a structure-elucidation program that uses a combination of 14 metabolome databases in addition to an enzyme promiscuity library. We showcase our workflow by annotating N-methyl-uridine monophosphate (UMP), lysomonogalactosyl-monopalmitin, N-methylalanine, and two propofol derivatives.

Precision and accuracy of single-molecule FRET measurements—a multi-laboratory benchmark study

Abstract: Single-molecule Forster resonance energy transfer (smFRET) is increasingly being used to determine distances, structures, and dynamics of biomolecules in vitro and in vivo. However, generalized protocols and FRET standards to ensure the reproducibility and accuracy of measurements of FRET efficiencies are currently lacking. Here we report the results of a comparative blind study in which 20 labs determined the FRET efficiencies (E) of several dye-labeled DNA duplexes. Using a unified, straightforward method, we obtained FRET efficiencies with s.d. between ±0.02 and ±0.05. We suggest experimental and computational procedures for converting FRET efficiencies into accurate distances, and discuss potential uncertainties in the experiment and the modeling. Our quantitative assessment of the reproducibility of intensity-based smFRET measurements and a unified correction procedure represents an important step toward the validation of distance networks, with the ultimate aim of achieving reliable structural models of biomolecular systems by smFRET-based hybrid methods.

An enhanced CRISPR repressor for targeted mammalian gene regulation

Abstract: The RNA-guided endonuclease Cas9 can be converted into a programmable transcriptional repressor, but inefficiencies in target-gene silencing have limited its utility. Here we describe an improved Cas9 repressor based on the C-terminal fusion of a rationally designed bipartite repressor domain, KRAB-MeCP2, to nuclease-dead Cas9. We demonstrate the system's superiority in silencing coding and noncoding genes, simultaneously repressing a series of target genes, improving the results of single and dual guide RNA library screens, and enabling new architectures of synthetic genetic circuits.

Systematic characterization of maturation time of fluorescent proteins in living cells

Abstract: The slow maturation time of fluorescent proteins (FPs) limits the temporal accuracy of measurements of rapid processes such as gene expression dynamics and effectively reduces fluorescence signal in growing cells. We used high-precision time-lapse microscopy to characterize the maturation kinetics of 50 FPs that span the visible spectrum at two different temperatures in Escherichia coli cells. We identified fast-maturing FPs from this set that yielded the highest signal-to-noise ratio and temporal resolution in individual growing cells.

Stability, affinity, and chromatic variants of the glutamate sensor iGluSnFR.

Abstract: Single-wavelength fluorescent reporters allow visualization of specific neurotransmitters with high spatial and temporal resolution. We report variants of intensity-based glutamate-sensing fluorescent reporter (iGluSnFR) that are functionally brighter; detect submicromolar to millimolar amounts of glutamate; and have blue, cyan, green, or yellow emission profiles. These variants could be imaged in vivo in cases where original iGluSnFR was too dim, resolved glutamate transients in dendritic spines and axonal boutons, and allowed imaging at kilohertz rates.

Using deep learning to model the hierarchical structure and function of a cell.

Abstract: Although artificial neural networks are powerful classifiers, their internal structures are hard to interpret. In the life sciences, extensive knowledge of cell biology provides an opportunity to design visible neural networks (VNNs) that couple the model's inner workings to those of real systems. Here we develop DCell, a VNN embedded in the hierarchical structure of 2,526 subsystems comprising a eukaryotic cell (http://d-cell.ucsd.edu/). Trained on several million genotypes, DCell simulates cellular growth nearly as accurately as laboratory observations. During simulation, genotypes induce patterns of subsystem activities, enabling in silico investigations of the molecular mechanisms underlying genotype-phenotype associations. These mechanisms can be validated, and many are unexpected; some are governed by Boolean logic. Cumulatively, 80% of the importance for growth prediction is captured by 484 subsystems (21%), reflecting the emergence of a complex phenotype. DCell provides a foundation for decoding the genetics of disease, drug resistance and synthetic life.

BoxCar acquisition method enables single-shot proteomics at a depth of 10,000 proteins in 100 minutes

Abstract: Great advances have been made in sensitivity and acquisition speed on the Orbitrap mass analyzer, enabling increasingly deep proteome coverage. However, these advances have been mainly limited to the MS2 level, whereas ion beam sampling for the MS1 scans remains extremely inefficient. Here we report a data-acquisition method, termed BoxCar, in which filling multiple narrow mass-to-charge segments increases the mean ion injection time more than tenfold as compared to that of a standard full scan. In 1-h analyses, the method provided MS1-level evidence for more than 90% of the proteome of a human cancer cell line that had previously been identified in 24 fractions, and it quantified more than 6,200 proteins in ten of ten replicates. In mouse brain tissue, we detected more than 10,000 proteins in only 100 min, and sensitivity extended into the low-attomolar range.

High-precision automated reconstruction of neurons with flood-filling networks.

Abstract: Reconstruction of neural circuits from volume electron microscopy data requires the tracing of cells in their entirety, including all their neurites. Automated approaches have been developed for tracing, but their error rates are too high to generate reliable circuit diagrams without extensive human proofreading. We present flood-filling networks, a method for automated segmentation that, similar to most previous efforts, uses convolutional neural networks, but contains in addition a recurrent pathway that allows the iterative optimization and extension of individual neuronal processes. We used flood-filling networks to trace neurons in a dataset obtained by serial block-face electron microscopy of a zebra finch brain. Using our method, we achieved a mean error-free neurite path length of 1.1 mm, and we observed only four mergers in a test set with a path length of 97 mm. The performance of flood-filling networks was an order of magnitude better than that of previous approaches applied to this dataset, although with substantially increased computational costs.

Burden-driven feedback control of gene expression

Abstract: In this CRISPR-based feedback control system, sgRNA expression is triggered by the burden of protein overexpression, and the sgRNA directs repression of the exogenous gene promoter to reduce burdensome expression and restore growth of the cell. Cells use feedback regulation to ensure robust growth despite fluctuating demands for resources and differing environmental conditions. However, the expression of foreign proteins from engineered constructs is an unnatural burden that cells are not adapted for. Here we combined RNA-seq with an in vivo assay to identify the major transcriptional changes that occur in Escherichia coli when inducible synthetic constructs are expressed. We observed that native promoters related to the heat-shock response activated expression rapidly in response to synthetic expression, regardless of the construct. Using these promoters, we built a dCas9-based feedback-regulation system that automatically adjusts the expression of a synthetic construct in response to burden. Cells equipped with this general-use controller maintained their capacity for native gene expression to ensure robust growth and thus outperformed unregulated cells in terms of protein yield in batch production. This engineered feedback is to our knowledge the first example of a universal, burden-based biomolecular control system and is modular, tunable and portable.

Genetically engineered cerebral organoids model brain tumor formation.

Abstract: Brain tumors are among the most lethal and devastating cancers Their study is limited by genetic heterogeneity and the incompleteness of available laboratory models Three-dimensional organoid culture models offer innovative possibilities for the modeling of human disease Here we establish a 3D in vitro model called a neoplastic cerebral organoid (neoCOR), in which we recapitulate brain tumorigenesis by introducing oncogenic mutations in cerebral organoids via transposon- and CRISPR-Cas9-mediated mutagenesis By screening clinically relevant mutations identified in cancer genome projects, we defined mutation combinations that result in glioblastoma-like and central nervous system primitive neuroectodermal tumor (CNS-PNET)-like neoplasms We demonstrate that neoCORs are suitable for use in investigations of aspects of tumor biology such as invasiveness, and for evaluation of drug effects in the context of specific DNA aberrations NeoCORs will provide a valuable complement to the current basic and preclinical models used to study brain tumor biology