Ultrafast population coding and axo-somatic compartmentalization

A model of task-level human stepping regulation yields semistable walking

A simple lateral dynamic walker, with swing leg dynamics and three adjustable input parameters, is used to study how motor regulation affects frontal plane stepping. Motivated by experimental observations and phenomenological models, we imposed task-level multiobjective regulation targeting the walker's optimal lateral foot placement at each step. The regulator prioritizes achieving step width and lateral body position goals to varying degrees by choosing a mixture parameter. Our model thus integrates a lateral mechanical template, which captures fundamental mechanics of frontal-plane walking, with a lateral motor regulation template, an empirically verified model of how humans manipulate lateral foot placements in a goal-directed manner. The model captures experimentally observed stepping fluctuation statistics and demonstrates how linear empirical models of stepping dynamics can emerge from first-principles nonlinear mechanics. We find that task-level regulation gives rise to a goal equivalent manifold in the system's extended state space of mechanical states and inputs, a subset of which contains a continuum of period-1 gaits forming a semistable set: perturbations off of any of its gaits result in transients that return to the set, though typically to different gaits.

How neuronal morphology impacts the synchronisation state of neuronal networks

The biophysical properties of neurons not only affect how information is processed within cells, they can also impact the dynamical states of the network. Specifically, the cellular dynamics of action-potential generation have shown relevance for setting the (de)synchronisation state of the network. The dynamics of tonically spiking neurons typically fall into one of three qualitatively distinct types that arise from distinct mathematical bifurcations of voltage dynamics at the onset of spiking. Accordingly, changes in ion channel composition or even external factors, like temperature, have been demonstrated to switch network behaviour via changes in the spike onset bifurcation and hence its associated dynamical type. A thus far less addressed modulator of neuronal dynamics is cellular morphology. Based on simplified and anatomically realistic mathematical neuron models, we show here that the extent of dendritic arborisation has an influence on the neuronal dynamical spiking type and therefore on the (de)synchronisation state of the network. Specifically, larger dendritic trees prime neuronal dynamics for in-phase-synchronised or splayed-out activity in weakly coupled networks, in contrast to cells with otherwise identical properties yet smaller dendrites. Our biophysical insights hold for generic multicompartmental classes of spiking neuron models (from ball-and-stick-type to anatomically reconstructed models) and establish a connection between neuronal morphology and the susceptibility of neural tissue to synchronisation in health and disease.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing non-genetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupts acquired resistance in GBM.

Methionine availability influences essential H3K36me3 dynamics during cell differentiation

Histone modifications are integral to epigenetics through their influence on gene expression and cellular status. While it's established that metabolism, including methionine metabolism, can impact histone methylation, the direct influence of methionine availability on crucial histone marks that determine the epigenomic process remains poorly understood. In this study, we demonstrate that methionine, through its metabolic product, S-adenosylmethionine (SAM), dynamically regulates H3K36me3, a cancer-associated histone modification known to influence cellular status, and myogenic differentiation of mouse myoblast cells. We further demonstrate that the methionine-dependent effects on differentiation are mediated in part through the histone methyltransferase SETD2. Methionine restriction leads to preferential decreases in H3K36me3 abundance and genome accessibility of genes involved in myogenic differentiation. Importantly, the effects of methionine restriction on differentiation and chromatin accessibility can be phenocopied by the deletion of Setd2. Collectively, this study demonstrates that methionine metabolism through its ability to be sensed by chromatin modifying enzymes can have a direct role in influencing cell fate determination.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

A Manifold-Based Framework for Studying the Dynamics of the Vaginal Microbiome

The vaginal bacterial community plays a crucial role in preventing infections. The composition of this community can be classified into five main groups, termed community state types (CSTs). Four of these CSTs, which are primarily consisted of Lactobacillus species, are considered healthy, while the fifth, which is composed of non-Lactobacillus populations, is considered less protective. This latter CST is often considered to represent a state termed Bacterial vaginosis (BV) - a common disease condition associated with unpleasant symptoms and increased susceptibility to sexually transmitted diseases. However, the exact mechanisms underlying BV development are not yet fully understood, including specifically, the dynamics of the vaginal microbiome in BV, and the possible routes it may take from a healthy to a BV state. This study aims to identify the progression from healthy Lactobacillus-dominant populations to symptomatic BV by analyzing 8,026 vaginal samples and using a manifold-detection framework. This approach is inspired by single-cell analysis and aims to identify low-dimensional trajectories in the high-dimensional composition space. This framework further order samples along these trajectories and assign a score (pseudo-time) to each sample based on its proximity to the BV state. Our results reveal distinct routes of progression between healthy and BV state for each CST, with pseudo-time scores correlating with community diversity and quantifying the health state of each sample. BV indicators, including Nugent score, positive Amsel's test, and several Amsel's criteria, can also be successfully predicted based on pseudo-time scores. Additionally, Gardnerella vaginalis can be identified as a key taxon in BV development using this approach, with increased abundance in samples with high pseudo-time, indicating an unhealthier state across all BV-development routes on the manifold. Taken together, these findings demonstrate how manifold detection can be used to successfully characterizes the progression from healthy Lactobacillus-dominant populations to BV and to accurately quantify the health condition of new samples along the route of BV development.

Microglia mediate contact-independent neuronal pruning via secreted Neuraminidase-3 associated with extracellular vesicles

Neurons communicate with each other through electrochemical transmission at synapses. Microglia, the resident immune cells of the central nervous system, can prune these synapses through a variety of contact-dependent and -independent means. Microglial secretion of active sialidase enzymes upon exposure to inflammatory stimuli is one unexplored mechanism of pruning. Recent work from our lab showed that treatment of neurons with bacterial sialidases disrupts neuronal network connectivity. Here, we find that activated microglia secrete Neuraminidase-3 (Neu3) associated with fusogenic extracellular vesicles. Furthermore, we show Neu3 mediates contact-independent pruning of neurons and subsequent disruption of neuronal networks through neuronal glycocalyx remodeling. We observe that NEU3 is transcriptionally upregulated upon exposure to inflammatory stimuli, and that a genetic knock-out of NEU3 abrogates the sialidase activity of inflammatory microglial secretions. Moreover, we demonstrate that Neu3 is associated with a subpopulation of extracellular vesicles, possibly exosomes, that are secreted by microglia upon inflammatory insult. Finally, we demonstrate that Neu3 is both necessary and sufficient to both desialylate neurons and decrease neuronal network connectivity. These results implicate Neu3 in remodeling of the glycocalyx leading to aberrant network-level activity of neurons, with implications in neuroinflammatory diseases such as Parkinson's disease and Alzheimer's disease.

Low dimensional criticality embedded in high dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. Here we offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality - a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescale agrees with established theory of critical dynamics. Our results suggest that cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Cortical Patterns Shift from Sequence Feature Separation during Planning to Integration during Motor Execution

Performing sequences of movements from memory and adapting them to changing task demands is a hallmark of skilled human behavior, from handwriting to playing a musical instrument. Prior studies showed a fine-grained tuning of cortical primary motor, premotor, and parietal regions to motor sequences: from the low-level specification of individual movements to high-level sequence features, such as sequence order and timing. However, it is not known how tuning in these regions unfolds dynamically across planning and execution. To address this, we trained 24 healthy right-handed human participants (14 females, 10 males) to produce four five-element finger press sequences with a particular finger order and timing structure in a delayed sequence production paradigm entirely from memory. Local cortical fMRI patterns during preparation and production phases were extracted from separate No-Go and Go trials, respectively, to tease out activity related to these perimovement phases. During sequence planning, premotor and parietal areas increased tuning to movement order or timing, regardless of their combinations. In contrast, patterns reflecting the unique integration of sequence features emerged in these regions during execution only, alongside timing-specific tuning in the ventral premotor, supplementary motor, and superior parietal areas. This was in line with the participants' behavioral transfer of trained timing, but not of order to new sequence feature combinations. Our findings suggest a general informational state shift from high-level feature separation to low-level feature integration within cortical regions for movement execution. Recompiling sequence features trial-by-trial during planning may enable flexible last-minute adjustment before movement initiation.SIGNIFICANCE STATEMENT Musicians and athletes can modify the timing and order of movements in a sequence trial-by-trial, allowing for a vast repertoire of flexible behaviors. How does the brain put together these high-level sequence features into an integrated whole? We found that, trial-by-trial, the control of sequence features undergoes a state shift from separation during planning to integration during execution across a network of motor-related cortical areas. These findings have implications for understanding the hierarchical control of skilled movement sequences, as well as how information in brain areas unfolds across planning and execution.