ABSTRACT
The advent of increasingly sophisticated imaging platforms has allowed for the visualization of the murine nervous system at single-cell resolution. However, current experimental approaches have not yet produced whole-brain maps of a comprehensive set of neuronal and nonneuronal types that approaches the cellular diversity of the mammalian cortex. Here, we aim to fill in this gap in knowledge with an open-source computational pipeline, Matrix Inversion and Subset Selection (MISS), that can infer quantitatively validated distributions of diverse collections of neural cell types at 200-µm resolution using a combination of single-cell RNA sequencing (RNAseq) and in situ hybridization datasets. We rigorously demonstrate the accuracy of MISS against literature expectations. Importantly, we show that gene subset selection, a procedure by which we filter out low-information genes prior to performing deconvolution, is a critical preprocessing step that distinguishes MISS from its predecessors and facilitates the production of cell-type maps with significantly higher accuracy. We also show that MISS is generalizable by generating high-quality cell-type maps from a second independently curated single-cell RNAseq dataset. Together, our results illustrate the viability of computational approaches for determining the spatial distributions of a wide variety of cell types from genetic data alone.
Subject(s)
Brain Mapping , Brain , Neurons , Animals , Brain/cytology , Brain Mapping/methods , Mice , Neurons/classification , Neurons/metabolism , RNA-Seq , Single-Cell AnalysisABSTRACT
The neurodegenerative disorder amyotrophic lateral sclerosis (ALS) is characterized by the progressive loss of upper and lower motor neurons, with pathological involvement of cerebral motor and extra-motor areas in a clinicopathological spectrum with frontotemporal dementia (FTD). A key unresolved issue is how the non-random distribution of pathology in ALS reflects differential network vulnerability, including molecular factors such as regional gene expression, or preferential spread of pathology via anatomical connections. A system of histopathological staging of ALS based on the regional burden of TDP-43 pathology observed in postmortem brains has been supported to some extent by analysis of distribution of in vivo structural MRI changes. In this paper, computational modeling using a Network Diffusion Model (NDM) was used to investigate whether a process of focal pathological 'seeding' followed by structural network-based spread recapitulated postmortem histopathological staging and, secondly, whether this had any correlation to the pattern of expression of a panel of genes implicated in ALS across the healthy brain. Regionally parcellated T1-weighted MRI data from ALS patients (baseline n=79) was studied in relation to a healthy control structural connectome and a database of associated regional cerebral gene expression. The NDM provided strong support for a structural network-based basis for regional pathological spread in ALS, but no simple relationship to the spatial distribution of ALS-related genes in the healthy brain. Interestingly, OPTN gene was identified as a significant but a weaker non-NDM contributor within the network-gene interaction model (LASSO). Intriguingly, the critical seed regions for spread within the model were not within the primary motor cortex but basal ganglia, thalamus and insula, where NDM recapitulated aspects of the postmortem histopathological staging system. Within the ALS-FTD clinicopathological spectrum, non-primary motor structures may be among the earliest sites of cerebral pathology.
Subject(s)
Amyotrophic Lateral Sclerosis , Connectome , Frontotemporal Dementia , Amyotrophic Lateral Sclerosis/diagnostic imaging , Amyotrophic Lateral Sclerosis/genetics , Amyotrophic Lateral Sclerosis/metabolism , Brain/metabolism , Frontotemporal Dementia/pathology , Humans , Motor NeuronsABSTRACT
Defects in axonal transport may partly underpin the differences between the observed pathophysiology of Alzheimer's disease (AD) and that of other non-amyloidogenic tauopathies. Particularly, pathological tau variants may have molecular properties that dysregulate motor proteins responsible for the anterograde-directed transport of tau in a disease-specific fashion. Here we develop the first computational model of tau-modified axonal transport that produces directional biases in the spread of tau pathology. We simulated the spatiotemporal profiles of soluble and insoluble tau species in a multicompartment, two-neuron system using biologically plausible parameters and time scales. Changes in the balance of tau transport feedback parameters can elicit anterograde and retrograde biases in the distributions of soluble and insoluble tau between compartments in the system. Aggregation and fragmentation parameters can also perturb this balance, suggesting a complex interplay between these distinct molecular processes. Critically, we show that the model faithfully recreates the characteristic network spread biases in both AD-like and non-AD-like mouse tauopathy models. Tau transport feedback may therefore help link microscopic differences in tau conformational states and the resulting variety in clinical presentations.
Subject(s)
Axonal Transport/physiology , tau Proteins/metabolism , Alzheimer Disease/metabolism , Animals , Computational Biology , Computer Simulation , Dendrites/metabolism , Disease Models, Animal , Feedback, Physiological , Humans , Mice , Models, Neurological , Neurodegenerative Diseases/metabolism , Protein Conformation , Protein Folding , Solubility , Spatio-Temporal Analysis , Tauopathies/metabolism , tau Proteins/chemistryABSTRACT
We introduce a computational model for the cellular level effects of firing rate filtering due to the major forms of neuronal injury, including demyelination and axonal swellings. Based upon experimental and computational observations, we posit simple phenomenological input/output rules describing spike train distortions and demonstrate that slow-gamma frequencies in the 38-41 Hz range emerge as the most robust to injury. Our signal-processing model allows us to derive firing rate filters at the cellular level for impaired neural activity with minimal assumptions. Specifically, we model eight experimentally observed spike train transformations by discrete-time filters, including those associated with increasing refractoriness and intermittent blockage. Continuous counterparts for the filters are also obtained by approximating neuronal firing rates from spike trains convolved with causal and Gaussian kernels. The proposed signal processing framework, which is robust to model parameter calibration, is an abstraction of the major cellular-level pathologies associated with neurodegenerative diseases and traumatic brain injuries that affect spike train propagation and impair neuronal network functionality. Our filters are well aligned with the spectrum of dynamic memory fields including working memory, visual consciousness, and other higher cognitive functions that operate in a frequency band that is - at a single cell level - optimally guarded against common types of pathological effects. In contrast, higher-frequency neural encoding, such as is observed with short-term memory, are susceptible to neurodegeneration and injury.
Subject(s)
Brain Injuries, Traumatic/physiopathology , Computer Simulation , Gamma Rhythm/physiology , Models, Neurological , Neurodegenerative Diseases/physiopathology , Action Potentials , Animals , Awareness/physiology , Axons/physiology , Cognition Disorders/physiopathology , Consciousness/physiology , Forecasting , Hippocampus/physiopathology , Humans , Memory, Short-Term/physiology , Rats , Synaptic Transmission/physiology , Vision, Ocular/physiologyABSTRACT
Alzheimer's disease, the most common form of dementia, is characterized by the emergence and spread of senile plaques and neurofibrillary tangles, causing widespread neurodegeneration. Though the progression of Alzheimer's disease is considered to be stereotyped, the significant variability within clinical populations obscures this interpretation on the individual level. Of particular clinical importance is understanding where exactly pathology, e.g. tau, emerges in each patient and how the incipient atrophy pattern relates to future spread of disease. Here we demonstrate a newly developed graph theoretical method of inferring prior disease states in patients with Alzheimer's disease and mild cognitive impairment using an established network diffusion model and an L1-penalized optimization algorithm. Although the 'seeds' of origin using our inference method successfully reproduce known trends in Alzheimer's disease staging on a population level, we observed that the high degree of heterogeneity between patients at baseline is also reflected in their seeds. Additionally, the individualized seeds are significantly more predictive of future atrophy than a single seed placed at the hippocampus. Our findings illustrate that understanding where disease originates in individuals is critical to determining how it progresses and that our method allows us to infer early stages of disease from atrophy patterns observed at diagnosis.
Subject(s)
Alzheimer Disease/diagnostic imaging , Cognitive Dysfunction/diagnosis , White Matter/pathology , Aged , Aged, 80 and over , Alzheimer Disease/cerebrospinal fluid , Amyloid beta-Peptides/cerebrospinal fluid , Atrophy/etiology , Cognitive Dysfunction/cerebrospinal fluid , Cohort Studies , Connectome , Correlation of Data , Disease Progression , Female , Functional Laterality , Humans , Magnetic Resonance Imaging , Male , Neurofibrillary Tangles/pathology , Psychiatric Status Rating Scales , Reference Values , White Matter/diagnostic imagingABSTRACT
Using a model for the dynamics of the full somatic nervous system of the nematode C. elegans, we address how biological network architectures and their functionality are degraded in the presence of focal axonal swellings (FAS) arising from neurodegenerative disease and/or traumatic brain injury. Using biophysically measured FAS distributions and swelling sizes, we are able to simulate the effects of injuries on the neural dynamics of C. elegans, showing how damaging the network degrades its low-dimensional dynamical responses. We visualize these injured neural dynamics by mapping them onto the worm's low-dimensional postures, i.e. eigenworm modes. We show that a diversity of functional deficits arise from the same level of injury on a connectomic network. Functional deficits are quantified using a statistical shape analysis, a procrustes analysis, for deformations of the limit cycles that characterize key behaviors such as forward crawling. This procrustes metric carries information on the functional outcome of injuries in the model. Furthermore, we apply classification trees to relate injury structure to the behavioral outcome. This makes testable predictions for the structure of an injury given a defined functional deficit. More critically, this study demonstrates the potential role of computational simulation studies in understanding how neuronal networks process biological signals, and how this processing is impacted by network injury.
Subject(s)
Caenorhabditis elegans/physiology , Connectome , Models, Neurological , Nerve Net/injuries , Nerve Net/physiopathology , Animals , Computational BiologyABSTRACT
The accurate diagnosis and assessment of neurodegenerative disease and traumatic brain injuries (TBI) remain open challenges. Both cause cognitive and functional deficits due to focal axonal swellings (FAS), but it is difficult to deliver a prognosis due to our limited ability to assess damaged neurons at a cellular level in vivo. We simulate the effects of neurodegenerative disease and TBI using convolutional neural networks (CNNs) as our model of cognition. We utilize biophysically relevant statistical data on FAS to damage the connections in CNNs in a functionally relevant way. We incorporate energy constraints on the brain by pruning the CNNs to be less over-engineered. Qualitatively, we demonstrate that damage leads to human-like mistakes. Our experiments also provide quantitative assessments of how accuracy is affected by various types and levels of damage. The deficit resulting from a fixed amount of damage greatly depends on which connections are randomly injured, providing intuition for why it is difficult to predict impairments. There is a large degree of subjectivity when it comes to interpreting cognitive deficits from complex systems such as the human brain. However, we provide important insight and a quantitative framework for disorders in which FAS are implicated.
Subject(s)
Brain Injuries, Traumatic/physiopathology , Brain/physiopathology , Cognition Disorders/physiopathology , Models, Neurological , Neural Networks, Computer , Neurodegenerative Diseases/physiopathology , Brain Injuries, Traumatic/complications , Cognition/physiology , Cognition Disorders/etiology , Humans , Neurodegenerative Diseases/complications , Neurons/physiologyABSTRACT
The presence of diffuse Focal Axonal Swellings (FAS) is a hallmark cellular feature in many neurological diseases and traumatic brain injury. Among other things, the FAS have a significant impact on spike-train encodings that propagate through the affected neurons, leading to compromised signal processing on a neuronal network level. This work merges, for the first time, three fields of study: (i) signal processing in excitatory-inhibitory (EI) networks of neurons via population codes, (ii) decision-making theory driven by the production of evidence from stimulus, and (iii) compromised spike-train propagation through FAS. As such, we demonstrate a mathematical architecture capable of characterizing compromised decision-making driven by cellular mechanisms. The computational model also leads to several novel predictions and diagnostics for understanding injury level and cognitive deficits, including a key finding that decision-making reaction times, rather than accuracy, are indicative of network level damage. The results have a number of translational implications, including that the level of network damage can be characterized by the reaction times in simple cognitive and motor tests.
Subject(s)
Brain Injuries, Traumatic/diagnosis , Decision Making , Models, Neurological , Reaction Time , Brain Injuries , Humans , NeuronsABSTRACT
Morphological reconstructions of axon segments reveal the abundance of geometrical ultrastructures that can dramatically affect the propagation of Action Potentials (AP). Moreover, deformations and swellings in axons resulting from brain traumas are associated to many neural dysfunctions and disorders. Our aim is to develop a computational framework to distinguish between geometrical enlargements that lead to minor changes in propagation from those that result in critical phenomenon such as reflection or blockage of the original traveling spike. We use a few geometrical parameters to model a prototypical shaft enlargement and explore the parameter space characterizing all possible propagation regimes and dynamics in an unmylienated AP model. Contrary to earlier notions that large diameter increases mostly lead to blocking, we demonstrate transmission is stable provided the geometrical changes occur in a slow manner. Our method also identifies a narrow range of parameters leading to a reflection regime. The distinction between these three regimes can be evaluated by a simple function of the geometrical parameters inferred through numerical simulations. We suggest that evaluating this function along axon segments can detect regions most susceptible to (i) transmission failure due to perturbations, (ii) structural plasticity, (iii) critical swellings caused by brain traumas and/or (iv) neurological disorders associated with the break down of spike train propagation.
Subject(s)
Action Potentials/physiology , Axons/physiology , Models, Neurological , Neurons/cytology , Neurons/physiology , Animals , Computer Simulation , Fourier Analysis , Humans , Numerical Analysis, Computer-AssistedABSTRACT
Axonal swellings are almost universal in neurodegenerative diseases of the central nervous system, including Alzheimer's and Parkinson's disease. Concussions and traumatic brain injuries can also produce cognitive and behavioral deficits by compromising neuronal morphology. Using a spike metric analysis, we characterize computationally the effects of such axonal varicosities on spike train propagation by comparing Poisson spike train classes before and after propagation through a prototypical axonal enlargement, or focused axonal swelling. Misclassification of spike train classes and low-pass filtering of firing rate activity increases with more pronounced axonal injury. We show that confusion matrices and a calculation of the loss of transmitted information provide a very practical way to characterize how injured neurons compromise the signal processing and faithful conductance of spike trains. The method demonstrates that (i) neural codes encoded with low firing rates are more robust to injury than those encoded with high firing rates, (ii) classification depends upon the length of the spike train used to encode information, and (iii) axonal injuries reduce the variance of spike trains within a given stimulus class. The work introduces a novel theoretical and computational framework to quantify the interplay between electrophysiological dynamics with focused axonal swellings generated by injury or other neurodegenerative processes. It further suggests how pharmacology and plasticity may play a role in recovery of neural computation. Ultimately, the work bridges vast experimental observations of in vitro morphological pathologies with post-traumatic cognitive and behavioral dysfunction.
Subject(s)
Axons/physiology , Brain Injuries/physiopathology , Models, Neurological , Neurodegenerative Diseases/physiopathology , Neurons/physiology , Brain Injuries/pathology , Computer Simulation , Humans , Neurodegenerative Diseases/pathology , Time FactorsABSTRACT
Brain regions in Alzheimer's (AD) exhibit distinct vulnerability to the disease's hallmark pathology, with the entorhinal cortex and hippocampus succumbing early to tau tangles while others like primary sensory cortices remain resilient. The quest to understand how local/regional genetic factors, pathogenesis, and network-mediated spread of pathology together govern this selective vulnerability (SV) or resilience (SR) is ongoing. Although many risk genes in AD are known from gene association and transgenic studies, it is still not known whether and how their baseline expression signatures confer SV or SR to brain structures. Prior analyses have yielded conflicting results, pointing to a disconnect between the location of genetic risk factors and downstream tau pathology. We hypothesize that a full accounting of genes' role in mediating SV/SR would require the modeling of network-based vulnerability, whereby tau misfolds, aggregates, and propagates along fiber projections. We therefore employed an extended network diffusion model (eNDM) and tested it on tau pathology PET data from 196 AD patients from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Thus the fitted eNDM model becomes a reference process from which to assess the role of innate genetic factors. Using the residual (observed - model-predicted) tau as a novel target outcome, we obtained its association with 100 top AD risk-genes, whose baseline spatial transcriptional profiles were obtained from the Allen Human Brain Atlas (AHBA). We found that while many risk genes at baseline showed a strong association with regional tau, many more showed a stronger association with residual tau. This suggests that both direct vulnerability, related to the network, as well as network-independent vulnerability, are conferred by risk genes. We then classified risk genes into four classes: network-related SV (SV-NR), network-independent SV (SV-NI), network-related SR (SR-NR), and network-independent SR (SR-NI). Each class has a distinct spatial signature and associated vulnerability to tau. Remarkably, we found from gene-ontology analyses, that genes in these classes were enriched in distinct functional processes and encompassed different functional networks. These findings offer new insights into the factors governing innate vulnerability or resilience in AD pathophysiology and may prove helpful in identifying potential intervention targets.
ABSTRACT
We accurately reconstruct the Local Field Potential time series obtained from anesthetized and awake rats, both before and during CO 2 euthanasia. We apply the Eigensystem Realization Algorithm to identify an underlying linear dynamical system capable of generating the observed data. Time series exhibiting more intricate dynamics typically lead to systems of higher dimensions, offering a means to assess the complexity of the brain throughout various phases of the experiment. Our results indicate that anesthetized brains possess complexity levels similar to awake brains before CO 2 administration. This resemblance undergoes significant changes following euthanization, as signals from the awake brain display a more resilient complexity profile, implying a state of heightened neuronal activity or a last fight response during the euthanasia process. In contrast, anesthetized brains seem to enter a more subdued state early on. Our data-driven techniques can likely be applied to a broader range of electrophysiological recording modalities.
Subject(s)
Algorithms , Brain , Animals , Brain/physiology , Rats , Wakefulness/physiology , Euthanasia , Male , Euthanasia, Animal/methods , Carbon DioxideABSTRACT
Neurodegenerative diseases such as Alzheimer's disease (AD) exhibit pathological changes in the brain that proceed in a stereotyped and regionally specific fashion, but the cellular and molecular underpinnings of regional vulnerability are currently poorly understood. Recent work has identified certain subpopulations of neurons in a few focal regions of interest, such as the entorhinal cortex, that are selectively vulnerable to tau pathology in AD. However, the cellular underpinnings of regional susceptibility to tau pathology are currently unknown, primarily because whole-brain maps of a comprehensive collection of cell types have been inaccessible. Here, we deployed a recent cell-type mapping pipeline, Matrix Inversion and Subset Selection (MISS), to determine the brain-wide distributions of pan-hippocampal and neocortical neuronal and non-neuronal cells in the mouse using recently available single-cell RNA sequencing (scRNAseq) data. We then performed a robust set of analyses to identify general principles of cell-type-based selective vulnerability using these cell-type distributions, utilizing 5 transgenic mouse studies that quantified regional tau in 12 distinct PS19 mouse models. Using our approach, which constitutes the broadest exploration of whole-brain selective vulnerability to date, we were able to discover cell types and cell-type classes that conferred vulnerability and resilience to tau pathology. Hippocampal glutamatergic neurons as a whole were strongly positively associated with regional tau deposition, suggesting vulnerability, while cortical glutamatergic and GABAergic neurons were negatively associated. Among glia, we identified oligodendrocytes as the single-most strongly negatively associated cell type, whereas microglia were consistently positively correlated. Strikingly, we found that there was no association between the gene expression relationships between cell types and their vulnerability or resilience to tau pathology. When we looked at the explanatory power of cell types versus GWAS-identified AD risk genes, cell type distributions were consistently more predictive of end-timepoint tau pathology than regional gene expression. To understand the functional enrichment patterns of the genes that were markers of the identified vulnerable or resilient cell types, we performed gene ontology analysis. We found that the genes that are directly correlated to tau pathology are functionally distinct from those that constitutively embody the vulnerable cells. In short, we have demonstrated that regional cell-type composition is a compelling explanation for the selective vulnerability observed in tauopathic diseases at a whole-brain level and is distinct from that conferred by risk genes. These findings may have implications in identifying cell-type-based therapeutic targets.
ABSTRACT
The prion-like transsynaptic propagation of misfolded tau along the brain's connectome has previously been modeled using connectome-based network diffusion models. In addition to the connectome, interactions between the general neurological "milieu" in the neurodegenerative brain and proteinopathic species can also contribute to pathology propagation. Such a molecular nexopathy framework posits that the distinct characteristics of neurodegenerative disorders stem from interactions between the network and surrounding molecular players. However, the effects of these modulators remain unquantified. Here, we present Nexopathy in silico ("Nexis"), a quantitative model of tau progression augmenting earlier models by including parameters of pathology propagation defined by the molecular modulators of connectome-based spread. Our Nexis:microglia model provides the first quantitative characterization of this effect on the whole brain by expanding previous models of neuropathology progression by incorporating microglial influence. We show that Trem2, but not microglial homeostasis genes, significantly improved the model's predictive power. Trem2 appears to reduce tau accumulation rate while increasing its interregional spread from the hippocampal seed area, causing higher tau burden in the striatum, pallidum, and contralateral hippocampus. Nexis provides an improved understanding and quantification of microglial contribution to tau propagation and can be flexibly modified to include other modulators of progressive neurodegeneration.
Subject(s)
NeuropathologyABSTRACT
Computational modeling of the neural activity in the human spinal cord may help elucidate the underlying mechanisms involved in the complex processing of painful stimuli. In this study, we use a biologically-plausible model of the dorsal horn circuitry as a platform to simulate pain processing under healthy and pathological conditions. Specifically, we distort signals in the receptor fibers akin to what is observed in axonal damage and monitor the corresponding changes in five quantitative markers associated with the pain response. Axonal damage may lead to spike-train delays, evoked potentials, an increase in the refractoriness of the system, and intermittent blockage of spikes. We demonstrate how such effects applied to mechanoreceptor and nociceptor fibers in the pain processing circuit can give rise to dramatically distinct responses at the network/population level. The computational modeling of damaged neuronal assemblies may help unravel the myriad of responses observed in painful neuropathies and improve diagnostics and treatment protocols.
ABSTRACT
Most organisms suffer neuronal damage throughout their lives, which can impair performance of core behaviors. Their neural circuits need to maintain function despite injury, which in particular requires preserving key system outputs. In this work, we explore whether and how certain structural and functional neuronal network motifs act as injury mitigation mechanisms. Specifically, we examine how (i) Hebbian learning, (ii) high levels of noise, and (iii) parallel inhibitory and excitatory connections contribute to the robustness of the olfactory system in the Manduca sexta moth. We simulate injuries on a detailed computational model of the moth olfactory network calibrated to data. The injuries are modeled on focal axonal swellings, a ubiquitous form of axonal pathology observed in traumatic brain injuries and other brain disorders. Axonal swellings effectively compromise spike train propagation along the axon, reducing the effective neural firing rate delivered to downstream neurons. All three of the network motifs examined significantly mitigate the effects of injury on readout neurons, either by reducing injury's impact on readout neuron responses or by restoring these responses to pre-injury levels. These motifs may thus be partially explained by their value as adaptive mechanisms to minimize the functional effects of neural injury. More generally, robustness to injury is a vital design principle to consider when analyzing neural systems.
ABSTRACT
BACKGROUND: The release of a broad, longitudinal anatomical dataset by the Parkinson's Progression Markers Initiative promoted a surge of machine-learning studies aimed at predicting disease onset and progression. However, the excessive number of features used in these models often conceals their relationship to the Parkinsonian symptomatology. OBJECTIVES: The aim of this study is two-fold: (i) to predict future motor and cognitive impairments up to four years from brain features acquired at baseline; and (ii) to interpret the role of pivotal brain regions responsible for different symptoms from a neurological viewpoint. METHODS: We test several deep-learning neural network configurations, and report our best results obtained with an autoencoder deep-learning model, run on a 5-fold cross-validation set. Comparison with Existing Methods: Our approach improves upon results from standard regression and others. It also includes neuroimaging biomarkers as features. RESULTS: The relative contributions of pivotal brain regions to each impairment change over time, suggesting a dynamical reordering of culprits as the disease progresses. Specifically, the Putamen is initially the most critical region accounting for the overall cognitive state, only being surpassed by the Substantia Nigra in later years. The Pallidum is the first region to influence motor scores, followed by the parahippocampal and ambient gyri, and the anterior orbital gyrus. CONCLUSIONS: While the causal link between regional brain atrophy and Parkinson symptomatology is poorly understood, our methods demonstrate that the contributions of pivotal regions to cognitive and motor impairments are more dynamical than generally appreciated.
ABSTRACT
There is enormous clinical value in inferring the brain regions initially atrophied in Parkinson disease for individual patients and understanding its relationship with clinical and genetic risk factors. The aim of this study is to leverage a new seed-inference algorithm demonstrated for Alzheimer's disease to the Parkinsonian context and to cluster patients in meaningful subgroups based on these incipient atrophy patterns. Instead of testing brain regions separately as the likely initiation site for each patient, we solve an L1-penalized optimization problem that can return a more predictive heterogeneous, multi-locus seed patterns. A cluster analysis of the individual seed patterns reveals two distinct subgroups (S1 versus S2). The S1 subgroup is characterized by the involvement of the brainstem and ventral nuclei, and S2 by cortex and striatum. Post hoc analysis in features not included in the clustering shows significant differences between subgroups regarding age of onset and local transcriptional patterns of Parkinson-related genes. Top genes associated with regional microglial abundance are strongly associated with subgroup S1 but not with S2. Our results suggest two distinct aetiological mechanisms operative in Parkinson disease. The interplay between immune-related genes, lysosomal genes, microglial abundance and atrophy initiation sites may explain why the age of onset for patients in S1 is on average 4.5 years later than for those in S2. We highlight and compare the most prominently affected brain regions for both subgroups. Altogether, our findings may improve current screening strategies for early Parkinson onsetters.
ABSTRACT
Local climate conditions play a major role in the biology of the Aedes aegypti mosquito, the main vector responsible for transmitting dengue, zika, chikungunya and yellow fever in urban centers. For this reason, a detailed assessment of periods in which changes in climate conditions affect the number of human cases may improve the timing of vector-control efforts. In this work, we develop new machine-learning algorithms to analyze climate time series and their connection to the occurrence of dengue epidemic years for seven Brazilian state capitals. Our method explores the impact of two key variables-frequency of precipitation and average temperature-during a wide range of time windows in the annual cycle. Our results indicate that each Brazilian state capital considered has its own climate signatures that correlate with the overall number of human dengue-cases. However, for most of the studied cities, the winter preceding an epidemic year shows a strong predictive power. Understanding such climate contributions to the vector's biology could lead to more accurate prediction models and early warning systems.
Subject(s)
Dengue/epidemiology , Forecasting/methods , Aedes/metabolism , Aedes/pathogenicity , Algorithms , Animals , Brazil/epidemiology , Chikungunya Fever/epidemiology , Chikungunya Fever/transmission , Cities/epidemiology , Climate , Dengue/transmission , Dengue Virus , Environment , Humans , Insect Vectors , Machine Learning , Mosquito Vectors , Rain , Seasons , Temperature , Yellow Fever/epidemiology , Yellow Fever/transmission , Zika Virus , Zika Virus Infection/epidemiology , Zika Virus Infection/transmissionABSTRACT
Neurodegenerative diseases and traumatic brain injuries (TBI) are among the main causes of cognitive dysfunction in humans. At a neuronal network level, they both extensively exhibit focal axonal swellings (FAS), which in turn, compromise the information encoded in spike trains and lead to potentially severe functional deficits. There are currently no satisfactory quantitative predictors of decline in memory-encoding neuronal networks based on the impact and statistics of FAS. Some of the challenges of this translational approach include our inability to access small scale injuries with non-invasive methods, the overall complexity of neuronal pathologies, and our limited knowledge of how networks process biological signals. The purpose of this computational study is three-fold: (i) to extend Hopfield's model for associative memory to account for the effects of FAS, (ii) to calibrate FAS parameters from biophysical observations of their statistical distribution and size, and (iii) to systematically evaluate deterioration rates for different memory-recall tasks as a function of FAS injury. We calculate deterioration rates for a face-recognition task to account for highly correlated memories and also for a discrimination task of random, uncorrelated memories with a size at the capacity limit of the Hopfield network. While it is expected that the performance of any injured network should decrease with injury, our results link, for the first time, the memory recall ability to observed FAS statistics. This allows for plausible estimates of cognitive decline for different stages of brain disorders within neuronal networks, bridging experimental observations following neurodegeneration and TBI with compromised memory recall. The work lends new insights to help close the gap between theory and experiment on how biological signals are processed in damaged, high-dimensional functional networks, and towards positing new diagnostic tools to measure cognitive deficits.