MultiMatch: geometry-informed colocalization in multi-color super-resolution microscopy.

With recent advances in multi-color super-resolution light microscopy, it is possible to simultaneously visualize multiple subunits within biological structures at nanometer resolution. To optimally evaluate and interpret spatial proximity of stainings on such an image, colocalization analysis tools have to be able to integrate prior knowledge on the local geometry of the recorded biological complex. We present MultiMatch to analyze the abundance and location of chain-like particle arrangements in multi-color microscopy based on multi-marginal optimal unbalanced transport methodology. Our object-based colocalization model statistically addresses the effect of incomplete labeling efficiencies enabling inference on existent, but not fully observable particle chains. We showcase that MultiMatch is able to consistently recover existing chain structures in three-color STED images of DNA origami nanorulers and outperforms geometry-uninformed triplet colocalization methods in this task. MultiMatch generalizes to an arbitrary number of color channels and is provided as a user-friendly Python package comprising colocalization visualizations.

Toward quantitative super-resolution microscopy: molecular maps with statistical guarantees.

Quantifying the number of molecules from fluorescence microscopy measurements is an important topic in cell biology and medical research. In this work, we present a consecutive algorithm for super-resolution (stimulated emission depletion (STED)) scanning microscopy that provides molecule counts in automatically generated image segments and offers statistical guarantees in form of asymptotic confidence intervals. To this end, we first apply a multiscale scanning procedure on STED microscopy measurements of the sample to obtain a system of significant regions, each of which contains at least one molecule with prescribed uniform probability. This system of regions will typically be highly redundant and consists of rectangular building blocks. To choose an informative but non-redundant subset of more naturally shaped regions, we hybridize our system with the result of a generic segmentation algorithm. The diameter of the segments can be of the order of the resolution of the microscope. Using multiple photon coincidence measurements of the same sample in confocal mode, we are then able to estimate the brightness and number of molecules and give uniform confidence intervals on the molecule counts for each previously constructed segment. In other words, we establish a so-called molecular map with uniform error control. The performance of the algorithm is investigated on simulated and real data.

Towards Unbiased Fluorophore Counting in Superresolution Fluorescence Microscopy.

With the advent of fluorescence superresolution microscopy, nano-sized structures can be imaged with a previously unprecedented accuracy. Therefore, it is rapidly gaining importance as an analytical tool in the life sciences and beyond. However, the images obtained so far lack an absolute scale in terms of fluorophore numbers. Here, we use, for the first time, a detailed statistical model of the temporal imaging process which relies on a hidden Markov model operating on two timescales. This allows us to extract this information from the raw data without additional calibration measurements. We show this on the basis of added data from experiments on single Alexa 647 molecules as well as GSDIM/dSTORM measurements on DNA origami structures with a known number of labeling positions.

Analysis of patchclamp recordings: model-free multiscale methods and software.

Analysis of patchclamp recordings is often a challenging issue. We give practical guidance how such recordings can be analyzed using the model-free multiscale idealization methodology JSMURF, JULES, and HILDE. We provide an operational manual how to use the accompanying software available as an R-package and as a graphical user interface. This includes selection of the right approach and tuning of parameters. We also discuss advantages and disadvantages of model-free approaches in comparison to hidden Markov model approaches and explain how they complement each other.

An antibiotic-resistance conferring mutation in a neisserial porin: Structure, ion flux, and ampicillin binding.

Gram-negative bacteria cause the majority of highly drug-resistant bacterial infections. To cross the outer membrane of the complex Gram-negative cell envelope, antibiotics permeate through porins, trimeric channel proteins that enable the exchange of small polar molecules. Mutations in porins contribute to the development of drug-resistant phenotypes. In this work, we show that a single point mutation in the porin PorB from Neisseria meningitidis, the causative agent of bacterial meningitis, can strongly affect the binding and permeation of beta-lactam antibiotics. Using X-ray crystallography, high-resolution electrophysiology, atomistic biomolecular simulation, and liposome swelling experiments, we demonstrate differences in drug binding affinity, ion selectivity and drug permeability of PorB. Our work further reveals distinct interactions between the transversal electric field in the porin eyelet and the zwitterionic drugs, which manifest themselves under applied electric fields in electrophysiology and are altered by the mutation. These observations may apply more broadly to drug-porin interactions in other channels. Our results improve the molecular understanding of porin-based drug-resistance in Gram-negative bacteria.

Colocalization for super-resolution microscopy via optimal transport.

Super-resolution fluorescence microscopy is a widely used technique in cell biology. Stimulated emission depletion (STED) microscopy enables the recording of multiple-color images with subdiffraction resolution. The enhanced resolution leads to new challenges regarding colocalization analysis of macromolecule distributions. We demonstrate that well-established methods for the analysis of colocalization in diffraction-limited datasets and for coordinate-stochastic nanoscopy are not equally well suited for the analysis of high-resolution STED images. We propose optimal transport colocalization, which measures the minimal transporting cost below a given spatial scale to match two protein intensity distributions. Its validity on simulated data as well as on dual-color STED recordings of yeast and mammalian cells is demonstrated. We also extend the optimal transport colocalization methodology to coordinate-stochastic nanoscopy.

Heterogeneous Idealization of Ion Channel Recordings - Open Channel Noise.

We propose a new model-free segmentation method for idealizing ion channel recordings. This method is designed to deal with heterogeneity of measurement errors. This in particular applies to open channel noise which, in general, is particularly difficult to cope with for model-free approaches. Our methodology is able to deal with lowpass filtered data which provides a further computational challenge. To this end we propose a multiresolution testing approach, combined with local deconvolution to resolve the lowpass filter. Simulations and statistical theory confirm that the proposed idealization recovers the underlying signal very accurately at presence of heterogeneous noise, even when events are shorter than the filter length. The method is compared to existing approaches in computer experiments and on real data. We find that it is the only one which allows to identify openings of the PorB porine at two different temporal scales. An implementation is available as an R package.

Multiple haplotype reconstruction from allele frequency data.

Because haplotype information is of widespread interest in biomedical applications, effort has been put into their reconstruction. Here, we propose an efficient method, called haploSep, that is able to accurately infer major haplotypes and their frequencies just from multiple samples of allele frequency data. Even the accuracy of experimentally obtained allele frequencies can be improved by re-estimating them from our reconstructed haplotypes. From a methodological point of view, we model our problem as a multivariate regression problem where both the design matrix and the coefficient matrix are unknown. Compared to other methods, haploSep is very fast, with linear computational complexity in the haplotype length. We illustrate our method on simulated and real data focusing on experimental evolution and microbial data.

Testing for dependence on tree structures.

Tree structures, showing hierarchical relationships and the latent structures between samples, are ubiquitous in genomic and biomedical sciences. A common question in many studies is whether there is an association between a response variable measured on each sample and the latent group structure represented by some given tree. Currently, this is addressed on an ad hoc basis, usually requiring the user to decide on an appropriate number of clusters to prune out of the tree to be tested against the response variable. Here, we present a statistical method with statistical guarantees that tests for association between the response variable and a fixed tree structure across all levels of the tree hierarchy with high power while accounting for the overall false positive error rate. This enhances the robustness and reproducibility of such findings.

Molecular contribution function in RESOLFT nanoscopy.

The ultimate objective of a microscope of the highest resolution is to map the molecules of interest in the sample. Traditionally, linear imaging systems are characterized by their spatial frequency transfer function, which is given, in real space, by the point spread function (PSF). By extending the concept of the PSF towards the molecular contribution function (MCF), that quantifies the average contribution of a single fluorophore to the image, a straightforward concept for counting fluorophores is obtained. Using reversible saturable optical fluorescence transitions (RESOLFT), fluorophores are effectively activated only in a small, subdiffraction-sized volume before they are read out. During readout the signal exhibits an increased variance due to the stochastic nature of prior activation, which scales quadratically with the brightness of the active fluorophores while the mean of the signal scales only linearly with it. Using a two-state Markov model for the activation, showing comparable behavior to the switching kinetics of the switchable fluorescent protein rsEGFP2, we can approximate quantitatively the MCF of RESOLFT nanoscopy allowing to count the number of fluorophores within a subdiffraction-sized region of the sample. The method is validated on measurements of tubulin structures in Drosophila melagonaster larvae. Modeling and estimation of the MCF is a promising approach to quantitative microscopy.

High-resolution experimental and computational electrophysiology reveals weak ß-lactam binding events in the porin PorB.

The permeation of most antibiotics through the outer membrane of Gram-negative bacteria occurs through porin channels. To design drugs with increased activity against Gram-negative bacteria in the face of the antibiotic resistance crisis, the strict constraints on the physicochemical properties of the permeants imposed by these channels must be better understood. Here we show that a combination of high-resolution electrophysiology, new noise-filtering analysis protocols and atomistic biomolecular simulations reveals weak binding events between the ß-lactam antibiotic ampicillin and the porin PorB from the pathogenic bacterium Neisseria meningitidis. In particular, an asymmetry often seen in the electrophysiological characteristics of ligand-bound channels is utilised to characterise the binding site and molecular interactions in detail, based on the principles of electro-osmotic flow through the channel. Our results provide a rationale for the determinants that govern the binding and permeation of zwitterionic antibiotics in porin channels.

Fully Automatic Multiresolution Idealization for Filtered Ion Channel Recordings: Flickering Event Detection.

We propose a new model-free segmentation method, JULES, which combines recent statistical multiresolution techniques with local deconvolution for idealization of ion channel recordings. The multiresolution criterion takes into account scales down to the sampling rate enabling the detection of flickering events, i.e., events on small temporal scales, even below the filter frequency. For such small scales the deconvolution step allows for a precise determination of dwell times and, in particular, of amplitude levels, a task which is not possible with common thresholding methods. This is confirmed theoretically and in a comprehensive simulation study. In addition, JULES can be applied as a preprocessing method for a refined hidden Markov analysis. Our new methodology allows us to show that gramicidin A flickering events have the same amplitude as the slow gating events. JULES is available as an R function jules in the package clampSeg.

A studentized permutation test for three-arm trials in the 'gold standard' design.

The 'gold standard' design for three-arm trials refers to trials with an active control and a placebo control in addition to the experimental treatment group. This trial design is recommended when being ethically justifiable and it allows the simultaneous comparison of experimental treatment, active control, and placebo. Parametric testing methods have been studied plentifully over the past years. However, these methods often tend to be liberal or conservative when distributional assumptions are not met particularly with small sample sizes. In this article, we introduce a studentized permutation test for testing non-inferiority and superiority of the experimental treatment compared with the active control in three-arm trials in the 'gold standard' design. The performance of the studentized permutation test for finite sample sizes is assessed in a Monte Carlo simulation study under various parameter constellations. Emphasis is put on whether the studentized permutation test meets the target significance level. For comparison purposes, commonly used Wald-type tests, which do not make any distributional assumptions, are included in the simulation study. The simulation study shows that the presented studentized permutation test for assessing non-inferiority in three-arm trials in the 'gold standard' design outperforms its competitors, for instance the test based on a quasi-Poisson model, for count data. The methods discussed in this paper are implemented in the R package ThreeArmedTrials which is available on the comprehensive R archive network (CRAN). Copyright © 2016 John Wiley & Sons, Ltd.

Partial least squares for dependent data.

We consider the partial least squares algorithm for dependent data and study the consequences of ignoring the dependence both theoretically and numerically. Ignoring nonstationary dependence structures can lead to inconsistent estimation, but a simple modification yields consistent estimation. A protein dynamics example illustrates the superior predictive power of the proposed method.

Design and analysis of three-arm trials with negative binomially distributed endpoints.

A three-arm clinical trial design with an experimental treatment, an active control, and a placebo control, commonly referred to as the gold standard design, enables testing of non-inferiority or superiority of the experimental treatment compared with the active control. In this paper, we propose methods for designing and analyzing three-arm trials with negative binomially distributed endpoints. In particular, we develop a Wald-type test with a restricted maximum-likelihood variance estimator for testing non-inferiority or superiority. For this test, sample size and power formulas as well as optimal sample size allocations will be derived. The performance of the proposed test will be assessed in an extensive simulation study with regard to type I error rate, power, sample size, and sample size allocation. For the purpose of comparison, Wald-type statistics with a sample variance estimator and an unrestricted maximum-likelihood estimator are included in the simulation study. We found that the proposed Wald-type test with a restricted variance estimator performed well across the considered scenarios and is therefore recommended for application in clinical trials. The methods proposed are motivated and illustrated by a recent clinical trial in multiple sclerosis. The R package ThreeArmedTrials, which implements the methods discussed in this paper, is available on CRAN.

Mapping molecules in scanning far-field fluorescence nanoscopy.

In fluorescence microscopy, the distribution of the emitting molecule number in space is usually obtained by dividing the measured fluorescence by that of a single emitter. However, the brightness of individual emitters may vary strongly in the sample or be inaccessible. Moreover, with increasing (super-) resolution, fewer molecules are found per pixel, making this approach unreliable. Here we map the distribution of molecules by exploiting the fact that a single molecule emits only a single photon at a time. Thus, by analysing the simultaneous arrival of multiple photons during confocal imaging, we can establish the number and local brightness of typically up to 20 molecules per confocal (diffraction sized) recording volume. Subsequent recording by stimulated emission depletion microscopy provides the distribution of the number of molecules with subdiffraction resolution. The method is applied to mapping the three-dimensional nanoscale organization of internalized transferrin receptors on human HEK293 cells.

Multiscale DNA partitioning: statistical evidence for segments.

MOTIVATION: DNA segmentation, i.e. the partitioning of DNA in compositionally homogeneous segments, is a basic task in bioinformatics. Different algorithms have been proposed for various partitioning criteria such as Guanine/Cytosine (GC) content, local ancestry in population genetics or copy number variation. A critical component of any such method is the choice of an appropriate number of segments. Some methods use model selection criteria and do not provide a suitable error control. Other methods that are based on simulating a statistic under a null model provide suitable error control only if the correct null model is chosen. RESULTS: Here, we focus on partitioning with respect to GC content and propose a new approach that provides statistical error control: as in statistical hypothesis testing, it guarantees with a user-specified probability [Formula: see text] that the number of identified segments does not exceed the number of actually present segments. The method is based on a statistical multiscale criterion, rendering this as a segmentation method that searches segments of any length (on all scales) simultaneously. It is also accurate in localizing segments: under benchmark scenarios, our approach leads to a segmentation that is more accurate than the approaches discussed in the comparative review of Elhaik et al. In our real data examples, we find segments that often correspond well to features taken from standard University of California at Santa Cruz (UCSC) genome annotation tracks. AVAILABILITY AND IMPLEMENTATION: Our method is implemented in function smuceR of the R-package stepR available at http://www.stochastik.math.uni-goettingen.de/smuce.

Fluorescence nanoscopy by polarization modulation and polarization angle narrowing.

When excited with rotating linear polarized light, differently oriented fluorescent dyes emit periodic signals peaking at different times. We show that measurement of the average orientation of fluorescent dyes attached to rigid sample structures mapped to regularly defined (50 nm)(2) image nanoareas can provide subdiffraction resolution (super resolution by polarization demodulation, SPoD). Because the polarization angle range for effective excitation of an oriented molecule is rather broad and unspecific, we narrowed this range by simultaneous irradiation with a second, de-excitation, beam possessing a polarization perpendicular to the excitation beam (excitation polarization angle narrowing, ExPAN). This shortened the periodic emission flashes, allowing better discrimination between molecules or nanoareas. Our method requires neither the generation of nanometric interference structures nor the use of switchable or blinking fluorescent probes. We applied the method to standard wide-field microscopy with camera detection and to two-photon scanning microscopy, imaging the fine structural details of neuronal spines.

Aggregated motion estimation for real-time MRI reconstruction.

PURPOSE: In real-time MRI serial images are generally reconstructed from highly undersampled datasets as the iterative solutions of an inverse problem. While practical realizations based on regularized nonlinear inversion (NLINV) have hitherto been surprisingly successful, strong assumptions about the continuity of image features may affect the temporal fidelity of the estimated reconstructions. THEORY AND METHODS: The proposed method for real-time image reconstruction integrates the deformations between nearby frames into the data consistency term of the inverse problem. The aggregated motion estimation (AME) is not required to be affine or rigid and does not need additional measurements. Moreover, it handles multi-channel MRI data by simultaneously determining the image and its coil sensitivity profiles in a nonlinear formulation which also adapts to non-Cartesian (e.g., radial) sampling schemes. The new method was evaluated for real-time MRI studies using highly undersampled radial gradient-echo sequences. RESULTS: AME reconstructions for a motion phantom with controlled speed as well as for measurements of human heart and tongue movements demonstrate improved temporal fidelity and reduced residual undersampling artifacts when compared with NLINV reconstructions without motion estimation. CONCLUSION: Nonlinear inverse reconstructions with aggregated motion estimation offer improved image quality and temporal acuity for visualizing rapid dynamic processes by real-time MRI.