Solar flares, email exchanges, and many natural or social systems exhibit bursty dynamics, with periods of intense activity separated by long inactivity. These patterns often follow power- law distributions in inter-event intervals or event rates. Existing models typically capture only one of these features and rely on non-local memory, which complicates analysis and mechanistic interpretation. We introduce a novel self-reinforcing point process whose event rates are governed by local, Markovian nonlinear dynamics and post-event resets. The model generates power-law tails for both inter-event intervals and event rates over a broad range of exponents observed empirically across natural and human phenomena. Compared to non-local models such as Hawkes processes, our approach is mechanistically simpler, highly analytically tractable, and also easier to simulate. We provide methods for model fitting and validation, establishing this framework as a versatile foundation for the study of bursty phenomena.

Signaling pathways induce stereotyped transcriptional changes as stem cells progress into mature cell types during embryogenesis. Signaling perturbations are necessary to discover which genes are responsive or insensitive to pathway activity. However, gene regulation is additionally dependent on cell state-specific factors like chromatin modifications or transcription factor binding. Thus, transcriptional profiles need to be assayed in single cells to identify potentially multiple, distinct perturbation responses among heterogeneous cell states in an embryo. In perturbation studies, comparing heterogeneous transcriptional states among experimental conditions often requires samples to be collected over multiple independent experiments, which can introduce confounding batch effects. We present Design-Aware Integration of Single Cell ExpEriments (DAISEE), a new algorithm that models perturbation responses in single-cell datasets collected according to complex experimental designs. We demonstrate that DAISEE improves upon a previously available integrative nonnegative matrix factorization framework, more efficiently separating perturbation responses from confounding variation. We use DAISEE to integrate newly collected single-cell RNA sequencing datasets from 5-h-old zebrafish embryos expressing optimized photoswitchable MEK (psMEK), which globally activates the extracellular signal-regulated kinase (ERK), a signaling molecule involved in many cell specification events. psMEK drives some cells that are normally not exposed to ERK signals toward other wild type states and induces novel states expressing early-acting endothelial genes. Overactive signaling is therefore capable of producing unexpected gene expression states in developing embryos.

bioRxiv preprint: https://www.doi.org/10.1101/2024.09.05.610903

Optical nanoscopy of intact biological specimens has been transformed by recent advancements in hydrogel-based tissue clearing and expansion, enabling the imaging of cellular and subcellular structures with molecular contrast. However, existing high-resolution fluorescence microscopes are physically limited by objective-to-specimen distance, which prevents the study of whole-mount specimens without physical sectioning. To address this challenge, we developed a photochemical strategy for spatially precise sectioning of specimens. By combining serial photochemical sectioning with lattice light-sheet imaging and petabyte-scale computation, we imaged and reconstructed axons and myelin sheaths across entire mouse olfactory bulbs at nanoscale resolution. An olfactory bulb–wide analysis of myelinated and unmyelinated axons revealed distinctive patterns of axon degeneration and de-/dysmyelination in the neurodegenerative brain, highlighting the potential for peta- to exabyte-scale super-resolution studies using this approach. High-resolution microscopes have a short working distance, making it difficult to see deep within large biological samples such as an intact brain. Slicing the tissue with a blade can reach deeper, but this often distorts or destroys the fine structures that scientists want to study. By embedding a sample in a light-sensitive hydrogel, Wang et al. demonstrated a gentler approach using a precise ray or sheet of light to dissolve or cut away tissue layer by layer. After each layer is removed, the newly exposed surface is imaged, allowing for a complete, high-resolution, three-dimensional reconstruction without damaging physical contact. 

bioRxiv preprint: https://www.biorxiv.org/content/10.1101/2024.08.01.605857v1

Neural mechanisms underlying sexually dimorphic social behaviors remain enigmatic in most species. In Drosophila, sexually dimorphic P1/pC1x neurons have been described as a site of sensory integration that regulates mating and aggressive behaviors. We show that the male P1/pC1x population forms a highly intertwined network with male-specific mAL and aSP-a neurons that is poised to regulate male behavior. The 48 P1/pC1x cell types exhibit heterogeneous synaptic connections with a subset receiving strong input from identified sensory pathways. We also describe circuit motifs by which P1 and sexually dimorphic aIPg neurons co-regulate social behaviors. Genetic driver lines for these cell types were generated and used to discover distinct roles for P1/pC1x cell types in promoting social acoustic signaling and male-male interactions. Our results reveal unexpected diversity in the connectivity and behavioral roles of the P1/pC1x cell types and provide essential genetic tools for interrogating their neurophysiological and behavioral functions.

Power-law scaling in coarse-grained data suggests critical dynamics, but the true source of this scaling often remains unclear. Here, we analyze neural activity recorded during spatial navigation, reproducing power-law scaling under a phenomenological renormalization group (PRG) procedure that clusters units by activity similarity. Such scaling was previously linked to criticality. Here, we show that the iterative nature of the procedure itself leads to the emergence of power laws when applied to heterogeneous, non-interacting units obeying spatially structured activity without requiring critical interactions. Furthermore, the scaling exponents produced by heteregeneous non-interacting units match the observed exponents in recorded neural data. A simplified version of the PRG further reveals how heterogeneity smooths transitions across scales, mimicking critical behavior. The resulting exponents depend systematically on system and population size, predictions confirmed by subsampling the data.

Researchers have long noted the differences in synapse count between different EM reconstructions of similar circuitry. In this paper we attempt to determine the portion of these differences that may be due to different sample preparation and imaging techniques, in particular serial-section transmission imaging (SS-TEM) compared to focused ion beam with scanning electron microscopy (FIB-SEM). To do this, we compare synapse detection in the major Drosophila EM reconstructions - FANC, MANC, FAFB (with original and new synapses), male CNS, BANC, and HemiBrain, plus several smaller reconstructions. We look at raw synapse counts to avoid any dependence on proofreading, and compensate insofar as possible for the confounds of sample sizes differences and different software detection efficiency. The result are estimates, per compartment and for the sample as a whole, of the number of synapses that would be visible to a skilled human observer. These are then compared across all samples, using regions which are reconstructed in common for each sample pair. We find that in almost all known cases where a volume has been reconstructed by both techniques, isotropic FIB-SEM reconstructions show more human-visible synapses than microtome sliced reconstructions, typically by more than 40%. This strongly suggests, but does not conclusively prove, that synapses are easier to see in isotropic FIB-SEM data.

fMRI signals were traditionally seen as slow and sampled in the order of seconds, but recent technological advances have enabled much faster sampling rates. We hypothesized that high-frequency fMRI signals can capture spontaneous neural activity that index brain states. Using fast fMRI (TR=378ms) and simultaneous EEG in 27 humans drifting between sleep and wakefulness, we found that fMRI spectral power increased during NREM sleep (compared to wakefulness) across several frequency ranges as fast as 1Hz. This fast fMRI power was correlated with canonical arousal-linked EEG rhythms (alpha and delta), with spatiotemporal correlation patterns for each rhythm reflecting a combination of shared arousal dynamics and rhythm-specific neural signatures. Using machine learning, we found that alpha and delta EEG rhythms can be decoded from fast fMRI signals, in subjects held-out from the training set, showing that fMRI as fast as 0.9Hz (alpha) and 0.7Hz (delta) contains reliable neurally-coupled information that generalizes across individuals. Finally, we demonstrate that this fast fMRI acquisition allows for EEG rhythms to be decoded from 3.8s windows of fMRI data. These results reveal that high-frequency fMRI signals are coupled to dynamically varying brain states, and that fast fMRI sampling allows for more temporally precise quantification of spontaneous neural activity than previously thought possible.

Matriglycan is a linear glycan (xylose-β1,3-glucuronate), which binds proteins in the extracellular matrix that contain laminin-globular domains and Lassa Fever Virus. It is indispensable for neuromuscular function. Matriglycan of insufficient length can cause muscular dystrophy with abnormal brain and eye development. LARGE1 (Like-acetylglucosaminyltransferase-1) uniquely synthesizes matriglycan on dystroglycan. The mechanism of matriglycan synthesis is not obvious from cryo-EM reconstructions of LARGE1. However, by reconstituting activity in vitro on recombinant prodystroglycan we show that the presence of the dystroglycan N-terminal domain (DGN), phosphorylated core M3, and a xylose-glucuronate primer are necessary for matriglycan polymerization by LARGE1. By introducing active site mutations, we demonstrate that LARGE1 processively polymerizes matriglycan on prodystroglycan, with its length regulated by the dystroglycan prodomain, DGN. Our enzymatic analysis of LARGE1 uncovers the mechanism of matriglycan synthesis on dystroglycan, which can form the basis for therapeutic strategies to treat matriglycan-deficient neuromuscular disorders and arenaviral infections.

Two simple models—vaulting over stiff legs and rebounding over compliant legs—are employed to describe the mechanics of legged locomotion. It is agreed that compliant legs are necessary for describing running, and that leg compliance is also present during walking. Stiff legs continue to be employed to model walking under the assumption that the compliance of the leg during walking is low enough to be considered stiff. Here we study gait choice and walk-to-run transition in a biped with compliance and show that the principles underlying gait choice and transition are completely different from stiff legs. Two findings underpin our conclusions: First, at the same speed, step length and stance duration, multiple gaits that differ in the number of times the leg expands and contracts during a single stance are possible. Among them, humans and other animals choose the (normal) gait with M-shaped vertical ground reaction forces (vGRF) not just because of energy considerations but also constraints from forces. Second, the transition from walking to running occurs because of three factors: vGRF minimum at mid-stance characteristic of normal walking, synchronization of horizontal and vertical motions during single support, and velocity redirection during the double support. The insight above required an analytical approximation of the double spring-loaded pendulum (DSLIP) model describing the intricate oscillatory dynamics that relate single and double support phases. Additionally, we also examined DSLIP as a quantitative model for locomotion and conclude that DSLIP speed range is limited. However, insights gleaned from the analytical treatment of DSLIP are general and will inform the construction of more accurate models of walking.

bioRxiv preprint: https://doi.org/10.1101/2024.09.23.612940

We have identified a Drosophila species in which males exhibit spontaneous, elaborate, and robust intermale sexual behavior. Males of D. santomea, a West African island endemic, distinguish conspecific sexes but court males and females promiscuously and seldom attack. Elevated intermale courtship derives from at least three changes in two separate pheromone systems. In males, the sexually monomorphic cuticular pheromone 7-tricosene promotes rather than inhibits courtship and the courtship-inhibiting olfactory pheromone cVA is reduced 84-92% compared to close relatives, including the sibling species D. yakuba. The third change is surprisingly in D. santomea females, where cVA suppresses rather than promotes sexual receptivity. The female cVA switch and male cVA reduction may have co-evolved to maintain efficient intraspecific mating in D. santomea but prevent sympatric hybridization with D. yakuba, or to reduce intraspecific aggression. We find that high intermale courtship and low cVA also co-occur and appear selectively derived in a distant monomorphic species D. persimilis, implying pheromonal and behavioral convergence in the two recently speciated taxa. The data suggest that sequential changes in the behavioral valence and levels of pheromones explain the recent evolutionary emergence of intermale sexual behavior in Drosophila.