Profile hidden Markov models (profile-HMMs) are sensitive tools for remote protein homology detection, but the main scoring algorithms, Viterbi or Forward, require considerable time to search large sequence databases.

Dopaminergic neurons with distinct projection patterns and physiological properties compose memory subsystems in a brain. However, it is poorly understood whether or how they interact during complex learning. Here, we identify a feedforward circuit formed between dopamine subsystems and show that it is essential for second-order conditioning, an ethologically important form of higher-order associative learning. The Drosophila mushroom body comprises a series of dopaminergic compartments, each of which exhibits distinct memory dynamics. We find that a slow and stable memory compartment can serve as an effective “teacher” by instructing other faster and transient memory compartments via a single key interneuron, which we identify by connectome analysis and neurotransmitter prediction. This excitatory interneuron acquires enhanced response to reward-predicting odor after first-order conditioning and, upon activation, evokes dopamine release in the “student” compartments. These hierarchical connections between dopamine subsystems explain distinct properties of first- and second-order memory long known by behavioral psychologists.

The anterodorsal projection neuron lineage of Drosophila melanogaster produces 40 neuronal types in a stereotypic order. Here we take advantage of this complete lineage sequence to examine the role of known temporal fating factors, including Chinmo and the Hb/Kr/Pdm/Cas transcriptional cascade, within this diverse central brain lineage. Kr mutation affects the temporal fate of the neuroblast (NB) itself, causing a single fate to be skipped, whereas Chinmo null only elicits fate transformation of NB progeny without altering cell counts. Notably, Chinmo operates in two separate windows to prevent fate transformation (into the subsequent Chinmo-indenpendent fate) within each window. By contrast, Hb/Pdm/Cas play no detectable role, indicating that Kr either acts outside of the cascade identified in the ventral nerve cord or that redundancy exists at the level of fating factors. Therefore, hierarchical fating mechanisms operate within the lineage to generate neuronal diversity in an unprecedented fashion.

Understanding how neural signals control muscle activity during behavior is a key challenge in motor neuroscience. To this end, recent advances in intramuscular multielectrode arrays have enabled high-quality multichannel recordings of many motor unit action potentials (MUAPs) in freely moving subjects. However, identifying individual MUAP events within multichannel recordings is a significant challenge for existing spike sorting methods, which are typically optimized for identifying action potentials from neurons in the brain. To overcome this challenge, we developed the Enhanced Motor Unit sorter (EMUsort), an extension of Kilosort4 (KS4) that achieves high-performance MUAP spike sorting. We applied EMUsort to high-resolution intramuscular recordings from rat forelimb during locomotion and monkey forelimb during a reaching task. EMUsort improves upon prior methods by addressing key challenges encountered with MUAP datasets, including: 1) long time delays across electrodes due to propagation along muscle fibers, 2) more complex waveform shapes compared to neuronal action potentials, and 3) a high degree of MUAP overlap due to cumulative motor unit recruitment. We compared EMUsort to existing spike sorting methods quantitatively using simulated datasets that closely emulated the rat and monkey datasets we recorded. EMUsort provided median error rate reductions of 67.5% and 49.9% during periods of high motor unit activation for the rat and monkey datasets, respectively. In sum, EMUsort provides a substantial improvement to MUAP spike sorter accuracy, especially during regions of high MUAP overlap, in an easy-to-use software package.

The challenge of recovering the topology of massive neuronal circuits can potentially be met by high throughput Electron Microscopy (EM) imagery. Segmenting a 3-dimensional stack of EM images into the individual neurons is difficult, due to the low depth-resolution in existing high-throughput EM technology, such as serial section Transmission EM (ssTEM). In this paper we propose methods for detecting the high resolution locations of membranes from low depth-resolution images. We approach this problem using both a method that learns a discriminative, over-complete dictionary and a kernel SVM. We test this approach on tomographic sections produced in simulations from high resolution Focused Ion Beam (FIB) images and on low depth-resolution images acquired with ssTEM and evaluate our results by comparing it to manual labeling of this data.

Random scattering and aberrations severely limit the imaging depth in optical microscopy. We introduce a rapid, parallel wavefront compensation technique that efficiently compensates even highly complex phase distortions. Using coherence gated backscattered light as a feedback signal, we focus light deep inside highly scattering brain tissue. We demonstrate that the same wavefront optimization technique can also be used to compensate spectral phase distortions in ultrashort laser pulses using nonlinear iterative feedback. We can restore transform limited pulse durations at any selected target location and compensate for dispersion that has occurred in the optical train and within the sample.

We describe a procedure for designing proteins with backbones produced by varying the parameters in the Crick coiled coil-generating equations. Combinatorial design calculations identify low-energy sequences for alternative helix supercoil arrangements, and the helices in the lowest-energy arrangements are connected by loop building. We design an antiparallel monomeric untwisted three-helix bundle with 80-residue helices, an antiparallel monomeric right-handed four-helix bundle, and a pentameric parallel left-handed five-helix bundle. The designed proteins are extremely stable (extrapolated ΔGfold > 60 kilocalories per mole), and their crystal structures are close to those of the design models with nearly identical core packing between the helices. The approach enables the custom design of hyperstable proteins with fine-tuned geometries for a wide range of applications.

Two-photon microscopy together with fluorescent proteins and fluorescent protein-based biosensors are commonly used tools in neuroscience. To enhance their experimental scope, it is important to optimize fluorescent proteins for two-photon excitation. Directed evolution of fluorescent proteins under one-photon excitation is common, but many one-photon properties do not correlate with two-photon properties. A simple system for expressing fluorescent protein mutants is colonies on an agar plate. The small focal volume of two-photon excitation makes creating a high throughput screen in this system a challenge for a conventional point-scanning approach. We present an instrument and accompanying software that solves this challenge by selectively scanning each colony based on a colony map captured under one-photon excitation. This instrument, called the GIZMO, can measure the two-photon excited fluorescence of 10,000 colonies in 7 hours. We show that the GIZMO can be used to evolve a fluorescent protein under two-photon excitation.

Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,478 SLE cases and 12,656 controls from six East Asian cohorts to identify new SLE susceptibility loci and better localize known loci. We identified ten new loci and confirmed 20 known loci with genome-wide significance. Among the new loci, the most significant locus was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta = 3.75 × 10(-117), odds ratio (OR) = 2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We identified the most likely functional variants at each locus by analyzing epigenetic marks and gene expression data. Ten candidate variants are known to alter gene expression in cis or in trans. Enrichment analysis highlights the importance of these loci in B cell and T cell biology. The new loci, together with previously known loci, increase the explained heritability of SLE to 24%. The new loci share functional and ontological characteristics with previously reported loci and are possible drug targets for SLE therapeutics.

We combined photoactivated localization microscopy (PALM) with live-cell single-particle tracking to create a new method termed sptPALM. We created spatially resolved maps of single-molecule motions by imaging the membrane proteins Gag and VSVG, and obtained several orders of magnitude more trajectories per cell than traditional single-particle tracking enables. By probing distinct subsets of molecules, sptPALM can provide insight into the origins of spatial and temporal heterogeneities in membranes.

Commentary: As a stepping stone to true live cell PALM (see above), our collaborator Jennifer Lippincott-Schwartz suggested using the sparse photoactivation principle of PALM to track the nanoscale motion of thousands of individual molecules within a single living cell. Termed single particle tracking PALM (sptPALM), Jennifer’s postdocs Suliana Manley and Jen Gillette used the method in our PALM rig to create spatially resolved maps of diffusion rates in the plasma membrane of live cells. sptPALM is a powerful tool to study the active cytoskeletal or passive diffusional transport of individual molecules with far more measurements per cell than is possible without sparse photoactivation.