Population activity measurement by calcium imaging can be combined with cellular resolution optogenetic activity perturbations to enable the mapping of neural connectivity in vivo. This requires accurate inference of perturbed and unperturbed neural activity from calcium imaging measurements, which are noisy and indirect, and can also be contaminated by photostimulation artifacts. We have developed a new fully Bayesian approach to jointly inferring spiking activity and neural connectivity from in vivo all-optical perturbation experiments. In contrast to standard approaches that perform spike inference and analysis in two separate maximum-likelihood phases, our joint model is able to propagate uncertainty in spike inference to the inference of connectivity and vice versa. We use the framework of variational autoencoders to model spiking activity using discrete latent variables, low-dimensional latent common input, and sparse spike-and-slab generalized linear coupling between neurons. Additionally, we model two properties of the optogenetic perturbation: off-target photostimulation and photostimulation transients. Our joint model includes at least two sets of discrete random variables; to avoid the dramatic slowdown typically caused by being unable to differentiate such variables, we introduce two strategies that have not, to our knowledge, been used with variational autoencoders. Using this model, we were able to fit models on 30 minutes of data in just 10 minutes. We performed an all-optical circuit mapping experiment in primary visual cortex of the awake mouse, and use our approach to predict neural connectivity between excitatory neurons in layer 2/3. Predicted connectivity is sparse and consistent with known correlations with stimulus tuning, spontaneous correlation and distance.

Understanding learning through synaptic plasticity rules in the brain is a grand challenge for neuroscience. Here we introduce a novel computational framework for inferring plasticity rules from experimental data on neural activity trajectories and behavioral learning dynamics. Our methodology parameterizes the plasticity function to provide theoretical interpretability and facilitate gradient-based optimization. For instance, we use Taylor series expansions or multilayer perceptrons to approximate plasticity rules, and we adjust their parameters via gradient descent over entire trajectories to closely match observed neural activity and behavioral data. Notably, our approach can learn intricate rules that induce long nonlinear time-dependencies, such as those incorporating postsynaptic activity and current synaptic weights. We validate our method through simulations, accurately recovering established rules, like Oja’s, as well as more complex hypothetical rules incorporating reward-modulated terms. We assess the resilience of our technique to noise and, as a tangible application, apply it to behavioral data from Drosophila during a probabilistic reward-learning experiment. Remarkably, we identify an active forgetting component of reward learning in flies that enhances the predictive accuracy of previous models. Overall, our modeling framework provides an exciting new avenue to elucidate the computational principles governing synaptic plasticity and learning in the brain.

The spatial association between fluorescently tagged biomolecules in situ provides valuable insight into their biological relationship. Within the limits of diffraction, such association can be measured using either Pearson's Correlation Coefficient (PCC) or Spearman's Rank Coefficient (SRC), which are designed to measure linear and monotonic correlations, respectively. However, the relationship between real biological signals is often more complex than these measures assume, rendering their results difficult to interpret. Here, we have adapted methods from the field of information theory to measure the association between two probes' concentrations based on their statistical dependence. Our approach is mathematically more general than PCC or SRC, making no assumptions about the type of relationship between the probes. We show that when applied to biological images, our measures provide more intuitive results that are also more robust to outliers and the presence of multiple relationships than PCC or SRC. We also devise a display technique to highlight regions in the input images where the probes' association is higher versus lower. We expect that our methods will allow biologists to more accurately and robustly quantify and visualize the association between two probes in a pair of fluorescence images. © 2018 International Society for Advancement of Cytometry.

A remaining challenge in protein modeling is to predict structures for sequences with no sequence similarity to any experimentally solved structure. Based on earlier observations, the library of protein backbone supersecondary structure motifs (Smotifs) saturated about a decade ago. Therefore, it should be possible to build any structure from a combination of existing Smotifs with the help of limited experimental data that are sufficient to relate the backbone conformations of Smotifs between target proteins and known structures. Here, we present a hybrid modeling algorithm that relies on an exhaustive Smotif library and on nuclear magnetic resonance chemical shift patterns without any input of primary sequence information. In a test of 102 proteins, the algorithm delivered 90 homology-model-quality models, among them 24 high-quality ones, and a topologically correct solution for almost all cases. The current approach opens a venue to address the modeling of larger protein structures for which chemical shifts are available.

The weak pixel counts surrounding the Bragg spots in a diffraction image are important for establishing a model of the background underneath the peak and estimating the reliability of the integrated intensities. Under certain circumstances, particularly with equipment not optimized for low-intensity measurements, these pixel values may be corrupted by corrections applied to the raw image. This can lead to truncation of low pixel counts, resulting in anomalies in the integrated Bragg intensities, such as systematically higher signal-to-noise ratios. A correction for this effect can be approximated by a three-parameter lognormal distribution fitted to the weakly positive-valued pixels at similar scattering angles. The procedure is validated by the improved refinement of an atomic model against structure factor amplitudes derived from corrected micro-electron diffraction (MicroED) images.

Primary aldosteronism (PA) is the most frequent form of secondary hypertension. Over the past two decades, major advances have been made in our understanding of the disease with the identification of germline or somatic mutations in ion channels and pumps. These mutations enhance calcium signaling, the main trigger of aldosterone biosynthesis.

In developing brains, axons exhibit remarkable precision in selecting synaptic partners among many non-partner cells. Evolutionarily conserved teneurins are transmembrane proteins that instruct synaptic partner matching. However, how intracellular signaling pathways execute teneurins' functions is unclear. Here, we use in situ proximity labeling to obtain the intracellular interactome of a teneurin (Ten-m) in the Drosophila brain. Genetic interaction studies using quantitative partner matching assays in both olfactory receptor neurons (ORNs) and projection neurons (PNs) reveal a common pathway: Ten-m binds to and negatively regulates a RhoGAP, thus activating the Rac1 small GTPases to promote synaptic partner matching. Developmental analyses with single-axon resolution identify the cellular mechanism of synaptic partner matching: Ten-m signaling promotes local F-actin levels and stabilizes ORN axon branches that contact partner PN dendrites. Combining spatial proteomics and high-resolution phenotypic analyses, this study advanced our understanding of both cellular and molecular mechanisms of synaptic partner matching.

The ATP-dependent chromatin-remodeling complex SWR1 exchanges a variant histone H2A.Z/H2B dimer for a canonical H2A/H2B dimer at nucleosomes flanking histone-depleted regions, such as promoters. This localization of H2A.Z is conserved throughout eukaryotes. SWR1 is a 1 megadalton complex containing 14 different polypeptides, including the AAA+ ATPases Rvb1 and Rvb2. Using electron microscopy, we obtained the three-dimensional structure of SWR1 and mapped its major functional components. Our data show that SWR1 contains a single heterohexameric Rvb1/Rvb2 ring that, together with the catalytic subunit Swr1, brackets two independently assembled multisubunit modules. We also show that SWR1 undergoes a large conformational change upon engaging a limited region of the nucleosome core particle. Our work suggests an important structural role for the Rvbs and a distinct substrate-handling mode by SWR1, thereby providing a structural framework for understanding the complex dimer-exchange reaction.