On clustering fMRI time series

Alterations in regional homogeneity in schizophrenia patients comorbid with metabolic syndrome treated with risperidone or clozapine

The neuroimaging mechanisms that arise in patients with schizophrenia and comorbid metabolic syndrome (MetS) remain poorly understood. This study was devised to examine potential alterations in regional homogeneity (ReHo) that arise in schizophrenia patients with comorbid MetS undergoing risperidone or clozapine treatment. In total, 43 schizophrenia patients undergoing risperidone or clozapine treatment were enrolled in this study, of whom 20 had comorbid MetS (SZ-MetS) while 23 did not (SZ-nMetS). In addition, 28 sex- and age-matched healthy controls (HCs) were analyzed. Analyses of covariance (ANCOVA) were used to compare ReHo in this study, utilizing age, sex, and years of education as covariates, with subsequent post hoc testing. Correlations between brain regions showing differences between groups of patients with schizophrenia and MetS-related indicators were also assessed. Relative to HCs, patients in both the SZ-MetS and SZ-nMetS groups exhibited reductions in ReHo in the right postcentral gyrus, left superior parietal gyrus, and left middle occipital gyrus. A decrease in ReHo was also evident in the left calcarine fissure and surrounding cortex and the right angular gyrus of patients in the SZ-MetS group, while the SZ-nMetS group exhibited reductions in ReHo in the left superior occipital gyrus and the right precuneus, together with an increase in ReHo in the right inferior orbitofrontal gyrus. Compared to the SZ-nMetS group, the SZ-MetS group exhibited a reduction in ReHo in the right inferior orbitofrontal gyrus. The ReHo of the right inferior orbitofrontal gyrus was significantly negatively correlated with BMI, waist circumference, and hip circumference in the SZ-MetS and SZ-nMetS groups. A reduction in ReHo in the right inferior orbitofrontal gyrus may thus be related to MetS in schizophrenia patients undergoing treatment. The findings may provide an imaging basis for brain alterations in patients with schizophrenia combined with MetS and provide a new insight into the neuroimaging mechanisms associated with MetS in patients with schizophrenia.

Neuroplastic Responses to Chiropractic Care: Broad Impacts on Pain, Mood, Sleep, and Quality of Life

OBJECTIVES

This study aimed to elucidate the mechanisms of chiropractic care using resting electroencephalography (EEG), somatosensory evoked potentials (SEPs), clinical health assessments (Fitbit), and Patient-reported Outcomes Measurement Information System (PROMIS-29).

METHODS

Seventy-six people with chronic low back pain (mean age ± SD: 45 ± 11 years, 33 female) were randomised into control (n = 38) and chiropractic (n = 38) groups. EEG and SEPs were collected pre and post the first intervention and post 4 weeks of intervention. PROMIS-29 was measured pre and post 4 weeks. Fitbit data were recorded continuously.

RESULTS

Spectral analysis of resting EEG showed a significant increase in Theta, Alpha and Beta, and a significant decrease in Delta power in the chiropractic group post intervention. Source localisation revealed a significant increase in Alpha activity within the Default Mode Network (DMN) post intervention and post 4 weeks. A significant decrease in N30 SEP peak amplitude post intervention and post 4 weeks was found in the chiropractic group. Source localisation demonstrated significant changes in Alpha and Beta power within the DMN post-intervention and post 4 weeks. Significant improvements in light sleep stage were observed in the chiropractic group along with enhanced overall quality of life post 4 weeks, including significant reductions in anxiety, depression, fatigue, and pain.

CONCLUSIONS

These findings indicate that many health benefits of chiropractic care are due to altered brain activity.

Evolutionarily conserved fMRI network dynamics in the mouse, macaque, and human brain

Evolutionarily relevant networks have been previously described in several mammalian species using time-averaged analyses of fMRI time-series. However, fMRI network activity is highly dynamic and continually evolves over timescales of seconds. Whether the dynamic organization of resting-state fMRI network activity is conserved across mammalian species remains unclear. Using frame-wise clustering of fMRI time-series, we find that intrinsic fMRI network dynamics in awake male macaques and humans is characterized by recurrent transitions between a set of 4 dominant, neuroanatomically homologous fMRI coactivation modes (C-modes), three of which are also plausibly represented in the male rodent brain. Importantly, in all species C-modes exhibit species-invariant dynamic features, including preferred occurrence at specific phases of fMRI global signal fluctuations, and a state transition structure compatible with infraslow coupled oscillator dynamics. Moreover, dominant C-mode occurrence reconstitutes the static organization of the fMRI connectome in all species, and is predictive of ranking of corresponding fMRI connectivity gradients. These results reveal a set of species-invariant principles underlying the dynamic organization of fMRI networks in mammalian species, and offer novel opportunities to relate fMRI network findings across the phylogenetic tree.

A Weighted GraphSAGE-Based Context-Aware Approach for Big Data Access Control

Context information is the key element to realizing dynamic access control of big data. However, existing context-aware access control (CAAC) methods do not support automatic context awareness and cannot automatically model and reason about context relationships. To solve these problems, this article proposes a weighted GraphSAGE-based context-aware approach for big data access control. First, graph modeling is performed on the access record data set and transforms the access control context-awareness problem into a graph neural network (GNN) node learning problem. Then, a GNN model WGraphSAGE is proposed to achieve automatic context awareness and automatic generation of CAAC rules. Finally, weighted neighbor sampling and weighted aggregation algorithms are designed for the model to realize automatic modeling and reasoning of node relationships and relationship strengths simultaneously in the graph node learning process. The experiment results show that the proposed method has obvious advantages in context awareness and context relationship reasoning compared with similar GNN models. Meanwhile, it obtains better results in dynamic access control decisions than the existing CAAC models.

Sex-specific differences in brain activity dynamics of youth with a family history of substance use disorder

An individual's risk of substance use disorder (SUD) is shaped by a complex interplay of potent biosocial factors. Current neurodevelopmental models posit vulnerability to SUD in youth is due to an overreactive reward system and reduced inhibitory control. Having a family history of SUD is a particularly strong risk factor, yet few studies have explored its impact on brain function and structure prior to substance exposure. Herein, we utilized a network control theory approach to quantify sex-specific differences in brain activity dynamics in youth with and without a family history of SUD, drawn from a large cohort of substance-naïve youth from the Adolescent Brain Cognitive Development Study. We summarize brain dynamics by calculating transition energy, which probes the ease with which a whole brain, region or network drives the brain towards a specific spatial pattern of activation (i.e., brain state). Our findings reveal that a family history of SUD is associated with alterations in the brain's dynamics wherein: i) independent of sex, certain regions' transition energies are higher in those with a family history of SUD and ii) there exist sex-specific differences in SUD family history groups at multiple levels of transition energy (global, network, and regional). Family history-by-sex effects reveal that energetic demand is increased in females with a family history of SUD and decreased in males with a family history of SUD, compared to their same-sex counterparts with no SUD family history. Specifically, we localize these effects to higher energetic demands of the default mode network in females with a family history of SUD and lower energetic demands of attention networks in males with a family history of SUD. These results suggest a family history of SUD may increase reward saliency in males and decrease efficiency of top-down inhibitory control in females. This work could be used to inform personalized intervention strategies that may target differing cognitive mechanisms that predispose individuals to the development of SUD.

Investigating Whether a Combination of Electro-Encephalography and Gene Expression Profiling Can Predict the Risk of Chronic Pain: A Protocol for an Observational Prospective Cohort Study

Despite most episodes of low back pain (LBP) being short-lasting, some transition into persistent long-lasting problems. Hence, the need for a deeper understanding of the physiological mechanisms of this is pertinent. Therefore, the aims of the present study are (1) to map pain-induced changes in brain activity and blood gene expression associated with persistent LBP, and (2) to explore whether these brain and gene expression signatures show promise as predictive biomarkers for the development of persistent LBP. The participants will be allocated into three different pain groups (no pain, mild short-lasting, or moderate long-term). One in-person visit, where two blood samples will be collected and sent for RNA sequencing, along with resting 64-channel electro-encephalography measurements before, during, and after a cold pressor test, will be conducted. Thereafter, follow-up questionnaires will be distributed at 2 weeks, 3 months, and 6 months. Recruitment will start during the second quarter of 2024, with expected completion by the last quarter of 2024. The results are expected to provide insight into the relationship between central nervous system activity, gene expression profiles, and LBP. If successful, this study has the potential to provide physiological indicators that are sensitive to the transition from mild, short-term LBP to more problematic, long-term LBP.

Hotspots of biogeochemical activity linked to aridity and plant traits across global drylands

Perennial plants create productive and biodiverse hotspots, known as fertile islands, beneath their canopies. These hotspots largely determine the structure and functioning of drylands worldwide. Despite their ubiquity, the factors controlling fertile islands under conditions of contrasting grazing by livestock, the most prevalent land use in drylands, remain virtually unknown. Here we evaluated the relative importance of grazing pressure and herbivore type, climate and plant functional traits on 24 soil physical and chemical attributes that represent proxies of key ecosystem services related to decomposition, soil fertility, and soil and water conservation. To do this, we conducted a standardized global survey of 288 plots at 88 sites in 25 countries worldwide. We show that aridity and plant traits are the major factors associated with the magnitude of plant effects on fertile islands in grazed drylands worldwide. Grazing pressure had little influence on the capacity of plants to support fertile islands. Taller and wider shrubs and grasses supported stronger island effects. Stable and functional soils tended to be linked to species-rich sites with taller plants. Together, our findings dispel the notion that grazing pressure or herbivore type are linked to the formation or intensification of fertile islands in drylands. Rather, our study suggests that changes in aridity, and processes that alter island identity and therefore plant traits, will have marked effects on how perennial plants support and maintain the functioning of drylands in a more arid and grazed world.

Grouped Spherical Data Modeling Through Hierarchical Nonparametric Bayesian Models and Its Application to fMRI Data Analysis

Recently, spherical data (i.e., L2 normalized vectors) modeling has become a promising research topic in various real-world applications (such as gene expression data analysis, document categorization, and gesture recognition). In this work, we propose a hierarchical nonparametric Bayesian model based on von Mises-Fisher (VMF) distributions for modeling spherical data that involve multiple groups, where each observation within a group is sampled from a VMF mixture model with an infinite number of components allowing them to be shared across groups. Our model is formulated by employing a hierarchical nonparametric Bayesian framework known as the hierarchical Pitman-Yor (HPY) process mixture model, which possesses a power-law nature over the distribution of the components and is particularly useful for data distributions with heavy tails and skewness. To learn the proposed HPY process mixture model with VMF distributions, we systematically develop a closed-form optimization algorithm based on variational Bayes (VB). The merits of the proposed hierarchical Bayesian nonparametric model for modeling grouped spherical data are demonstrated through experiments on both synthetic data and a real-world application about resting-state functional magnetic resonance imaging (fMRI) data analysis.

Anatomically distributed neural representations of instincts in the hypothalamus

Artificial activation of anatomically localized, genetically defined hypothalamic neuron populations is known to trigger distinct innate behaviors, suggesting a hypothalamic nucleus-centered organization of behavior control. To assess whether the encoding of behavior is similarly anatomically confined, we performed simultaneous neuron recordings across twenty hypothalamic regions in freely moving animals. Here we show that distinct but anatomically distributed neuron ensembles encode the social and fear behavior classes, primarily through mixed selectivity. While behavior class-encoding ensembles were spatially distributed, individual ensembles exhibited strong localization bias. Encoding models identified that behavior actions, but not motion-related variables, explained a large fraction of hypothalamic neuron activity variance. These results identify unexpected complexity in the hypothalamic encoding of instincts and provide a foundation for understanding the role of distributed neural representations in the expression of behaviors driven by hardwired circuits.

Automatic Clustering of Excited-State Trajectories: Application to Photoexcited Dynamics

We introduce automatic clustering as a computationally efficient tool for classifying and interpreting trajectories from simulations of photo-excited dynamics. Trajectories are treated as time-series data, with the features for clustering selected by variance mapping of normalized data. The L2-norm and dynamic time warping are proposed as suitable similarity measures for calculating the distance matrices, and these are clustered using the unsupervised density-based DBSCAN algorithm. The silhouette coefficient and the number of trajectories classified as noise are used as quality measures for the clustering. The ability of clustering to provide rapid overview of large and complex trajectory data sets, and its utility for extracting chemical and physical insight, is demonstrated on trajectories corresponding to the photochemical ring-opening reaction of 1,3-cyclohexadiene, noting that the clustering can be used to generate reduced dimensionality representations in an unbiased manner.

Progressive trajectories of schizophrenia across symptoms, genes, and the brain

BACKGROUND

Schizophrenia is characterized by complex psychiatric symptoms and unclear pathological mechanisms. Most previous studies have focused on the morphological changes that occur over the development of the disease; however, the corresponding functional trajectories remain unclear. In the present study, we aimed to explore the progressive trajectories of patterns of dysfunction after diagnosis.

METHODS

Eighty-six patients with schizophrenia and 120 healthy controls were recruited as the discovery dataset. Based on multiple functional indicators of resting-state brain functional magnetic resonance imaging, we conducted a duration-sliding dynamic analysis framework to investigate trajectories in association with disease progression. Neuroimaging findings were associated with clinical symptoms and gene expression data from the Allen Human Brain Atlas database. A replication cohort of patients with schizophrenia from the University of California, Los Angeles, was used as the replication dataset for the validation analysis.

RESULTS

Five stage-specific phenotypes were identified. A symptom trajectory was characterized by positive-dominated, negative ascendant, negative-dominated, positive ascendant, and negative surpassed stages. Dysfunctional trajectories from primary and subcortical regions to higher-order cortices were recognized; these are associated with abnormal external sensory gating and a disrupted internal excitation-inhibition equilibrium. From stage 1 to stage 5, the importance of neuroimaging features associated with behaviors gradually shifted from primary to higher-order cortices and subcortical regions. Genetic enrichment analysis identified that neurodevelopmental and neurodegenerative factors may be relevant as schizophrenia progresses and highlighted multiple synaptic systems.

CONCLUSIONS

Our convergent results indicate that progressive symptoms and functional neuroimaging phenotypes are associated with genetic factors in schizophrenia. Furthermore, the identification of functional trajectories complements previous findings of structural abnormalities and provides potential targets for drug and non-drug interventions in different stages of schizophrenia.

Longitudinal investigation of changes in resting-state co-activation patterns and their predictive ability in the zQ175 DN mouse model of Huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder caused by expanded (≥ 40) glutamine-encoding CAG repeats in the huntingtin gene, which leads to dysfunction and death of predominantly striatal and cortical neurons. While the genetic profile and clinical signs and symptoms of the disease are better known, changes in the functional architecture of the brain, especially before the clinical expression becomes apparent, are not fully and consistently characterized. In this study, we sought to uncover functional changes in the brain in the heterozygous (HET) zQ175 delta-neo (DN) mouse model at 3, 6, and 10 months of age, using resting-state functional magnetic resonance imaging (RS-fMRI). This mouse model shows molecular, cellular and circuitry alterations that worsen through age. Motor function disturbances are manifested in this model at 6 and 10 months of age. Specifically, we investigated, longitudinally, changes in co-activation patterns (CAPs) that are the transient states of brain activity constituting the resting-state networks (RSNs). Most robust changes in the temporal properties of CAPs occurred at the 10-months time point; the durations of two anti-correlated CAPs, characterized by simultaneous co-activation of default-mode like network (DMLN) and co-deactivation of lateral-cortical network (LCN) and vice-versa, were reduced in the zQ175 DN HET animals compared to the wild-type mice. Changes in the spatial properties, measured in terms of activation levels of different brain regions, during CAPs were found at all three ages and became progressively more pronounced at 6-, and 10 months of age. We then assessed the cross-validated predictive power of CAP metrics to distinguish HET animals from controls. Spatial properties of CAPs performed significantly better than the chance level at all three ages with 80% classification accuracy at 6 and 10 months of age.

Effect of LSD and music on the time-varying brain dynamics

RATIONALE

Psychedelics are getting closer to being widely used in clinical treatment. Music is known as a key element of psychedelic-assisted therapy due to its psychological effects, specifically on the emotion, meaning-making, and sensory processing. However, there is still a lack of understanding in how psychedelics influence brain activity in experimental settings involving music listening.

OBJECTIVES

The main goal of our research was to investigate the effect of music, as a part of "setting," on the brain states dynamics after lysergic acid diethylamide (LSD) intake.

METHODS

We used an open dataset, where a group of 15 participants underwent two functional MRI scanning sessions under LSD and placebo influence. Every scanning session contained three runs: two resting-state runs separated by one run with music listening. We applied K-Means clustering to identify the repetitive patterns of brain activity, so-called brain states. For further analysis, we calculated states' dwell time, fractional occupancy and transition probability.

RESULTS

The interaction effect of music and psychedelics led to change in the time-varying brain activity of the task-positive state. LSD, regardless of the music, affected the dynamics of the state of combined activity of DMN, SOM, and VIS networks. Crucially, we observed that the music itself could potentially have a long-term influence on the resting-state, in particular on states involving task-positive networks.

CONCLUSIONS

This study indicates that music, as a crucial element of "setting," can potentially have an influence on the subject's resting-state during psychedelic experience. Further studies should replicate these results on a larger sample size.

Grazing and ecosystem service delivery in global drylands

Grazing represents the most extensive use of land worldwide. Yet its impacts on ecosystem services remain uncertain because pervasive interactions between grazing pressure, climate, soil properties, and biodiversity may occur but have never been addressed simultaneously. Using a standardized survey at 98 sites across six continents, we show that interactions between grazing pressure, climate, soil, and biodiversity are critical to explain the delivery of fundamental ecosystem services across drylands worldwide. Increasing grazing pressure reduced ecosystem service delivery in warmer and species-poor drylands, whereas positive effects of grazing were observed in colder and species-rich areas. Considering interactions between grazing and local abiotic and biotic factors is key for understanding the fate of dryland ecosystems under climate change and increasing human pressure.

Navigating the Statistical Minefield of Model Selection and Clustering in Neuroscience

Model selection is often implicit: when performing an ANOVA, one assumes that the normal distribution is a good model of the data; fitting a tuning curve implies that an additive and a multiplicative scaler describes the behavior of the neuron; even calculating an average implicitly assumes that the data were sampled from a distribution that has a finite first statistical moment: the mean. Model selection may be explicit, when the aim is to test whether one model provides a better description of the data than a competing one. As a special case, clustering algorithms identify groups with similar properties within the data. They are widely used from spike sorting to cell type identification to gene expression analysis. We discuss model selection and clustering techniques from a statistician's point of view, revealing the assumptions behind, and the logic that governs the various approaches. We also showcase important neuroscience applications and provide suggestions how neuroscientists could put model selection algorithms to best use as well as what mistakes should be avoided.

Biclustering fMRI time series: a comparative study

BACKGROUND

The effectiveness of biclustering, simultaneous clustering of rows and columns in a data matrix, was shown in gene expression data analysis. Several researchers recognize its potentialities in other research areas. Nevertheless, the last two decades have witnessed the development of a significant number of biclustering algorithms targeting gene expression data analysis and a lack of consistent studies exploring the capacities of biclustering outside this traditional application domain.

RESULTS

This work evaluates the potential use of biclustering in fMRI time series data, targeting the Region × Time dimensions by comparing seven state-in-the-art biclustering and three traditional clustering algorithms on artificial and real data. It further proposes a methodology for biclustering evaluation beyond gene expression data analysis. The results discuss the use of different search strategies in both artificial and real fMRI time series showed the superiority of exhaustive biclustering approaches, obtaining the most homogeneous biclusters. However, their high computational costs are a challenge, and further work is needed for the efficient use of biclustering in fMRI data analysis.

CONCLUSIONS

This work pinpoints avenues for the use of biclustering in spatio-temporal data analysis, in particular neurosciences applications. The proposed evaluation methodology showed evidence of the effectiveness of biclustering in finding local patterns in fMRI time series data. Further work is needed regarding scalability to promote the application in real scenarios.

Automated Detection of Electric Energy Consumption Load Profile Patterns

Load profiles of energy consumption from smart meters are becoming more and more available, and the amount of data to analyse is huge. In order to automate this analysis, the application of state-of-the-art data mining techniques for time series analysis is reviewed. In particular, the use of dynamic clustering techniques to obtain and visualise temporal patterns characterising the users of electrical energy is deeply studied. The performed review can be used as a guide for those interested in the automatic analysis and groups of behaviour detection within load profile databases. Additionally, a selection of dynamic clustering algorithms have been implemented and the performances compared using an available electric energy consumption load profile database. The results allow experts to easily evaluate how users consume energy, to assess trends and to predict future scenarios.

Analysis of fMRI Signals from Working Memory Tasks and Resting-State of Brain: Neutrosophic-Entropy-Based Clustering Algorithm

This study applies a neutrosophic-entropy-based clustering algorithm (NEBCA) to analyze the fMRI signals. We consider the data obtained from four different working memory tasks and the brain's resting state for the experimental purpose. Three non-overlapping clusters of data related to temporal brain activity are determined and statistically analyzed. Moreover, we used the Uniform Manifold Approximation and Projection (UMAP) method to reduce system dimensionality and present the effectiveness of NEBCA. The results show that using NEBCA, we are able to distinguish between different working memory tasks and resting-state and identify subtle differences in the related activity of brain regions. By analyzing the statistical properties of the entropy inside the clusters, the various regions of interest (ROIs), according to Automated Anatomical Labeling (AAL) atlas crucial for clustering procedure, are determined. The inferior occipital gyrus is established as an important brain region in distinguishing the resting state from the tasks. Moreover, the inferior occipital gyrus and superior parietal lobule are identified as necessary to correct the data discrimination related to the different memory tasks. We verified the statistical significance of the results through the two-sample t-test and analysis of surrogates performed by randomization of the cluster elements. The presented methodology is also appropriate to determine the influence of time of day on brain activity patterns. The differences between working memory tasks and resting-state in the morning are related to a lower index of small-worldness and sleep inertia in the first hours after waking. We also compared the performance of NEBCA to two existing algorithms, KMCA and FKMCA. We showed the advantage of the NEBCA over these algorithms that could not effectively accumulate fMRI signals with higher variability.

Deep Time-Series Clustering: A Review

We present a comprehensive, detailed review of time-series data analysis, with emphasis on deep time-series clustering (DTSC), and a case study in the context of movement behavior clustering utilizing the deep clustering method. Specifically, we modified the DCAE architectures to suit time-series data at the time of our prior deep clustering work. Lately, several works have been carried out on deep clustering of time-series data. We also review these works and identify state-of-the-art, as well as present an outlook on this important field of DTSC from five important perspectives.

The Cancer Surfaceome Atlas integrates genomic, functional and drug response data to identify actionable targets

Cell-surface proteins (SPs) are a rich source of immune and targeted therapies. By systematically integrating single-cell and bulk genomics, functional studies and target actionability, in the present study we comprehensively identify and annotate genes encoding SPs (GESPs) pan-cancer. We characterize GESP expression patterns, recurrent genomic alterations, essentiality, receptor-ligand interactions and therapeutic potential. We also find that mRNA expression of GESPs is cancer-type specific and positively correlates with protein expression, and that certain GESP subgroups function as common or specific essential genes for tumor cell growth. We also predict receptor-ligand interactions substantially deregulated in cancer and, using systems biology approaches, we identify cancer-specific GESPs with therapeutic potential. We have made this resource available through the Cancer Surfaceome Atlas ( http://fcgportal.org/TCSA ) within the Functional Cancer Genome data portal.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells

Background

Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown.

Results

The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D.

Conclusion

Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08057-4.

W-TSS: A Wavelet-Based Algorithm for Discovering Time Series Shapelets

Many approaches to time series classification rely on machine learning methods. However, there is growing interest in going beyond black box prediction models to understand discriminatory features of the time series and their associations with outcomes. One promising method is time-series shapelets (TSS), which identifies maximally discriminative subsequences of time series. For example, in environmental health applications TSS could be used to identify short-term patterns in exposure time series (shapelets) associated with adverse health outcomes. Identification of candidate shapelets in TSS is computationally intensive. The original TSS algorithm used exhaustive search. Subsequent algorithms introduced efficiencies by trimming/aggregating the set of candidates or training candidates from initialized values, but these approaches have limitations. In this paper, we introduce Wavelet-TSS (W-TSS) a novel intelligent method for identifying candidate shapelets in TSS using wavelet transformation discovery. We tested W-TSS on two datasets: (1) a synthetic example used in previous TSS studies and (2) a panel study relating exposures from residential air pollution sensors to symptoms in participants with asthma. Compared to previous TSS algorithms, W-TSS was more computationally efficient, more accurate, and was able to discover more discriminative shapelets. W-TSS does not require pre-specification of shapelet length.

Construction of embedded fMRI resting-state functional connectivity networks using manifold learning

We construct embedded functional connectivity networks (FCN) from benchmark resting-state functional magnetic resonance imaging (rsfMRI) data acquired from patients with schizophrenia and healthy controls based on linear and nonlinear manifold learning algorithms, namely, Multidimensional Scaling, Isometric Feature Mapping, Diffusion Maps, Locally Linear Embedding and kernel PCA. Furthermore, based on key global graph-theoretic properties of the embedded FCN, we compare their classification potential using machine learning. We also assess the performance of two metrics that are widely used for the construction of FCN from fMRI, namely the Euclidean distance and the cross correlation metric. We show that diffusion maps with the cross correlation metric outperform the other combinations.

Systematic Characterization of Recurrent Genomic Alterations in Cyclin-Dependent Kinases Reveals Potential Therapeutic Strategies for Cancer Treatment

Recurrent copy-number alterations, mutations, and transcript fusions of the genes encoding CDKs/cyclins are characterized in >10,000 tumors. Genomic alterations of CDKs/cyclins are dominantly driven by copy number aberrations. In contrast to cell-cycle-related CDKs/cyclins, which are globally amplified, transcriptional CDKs/cyclins recurrently lose copy numbers across cancers. Although mutations and transcript fusions are relatively rare events, CDK12 exhibits recurrent mutations in multiple cancers. Among the transcriptional CDKs, CDK7 and CDK12 show the most significant copy number loss and mutation, respectively. Their genomic alterations are correlated with increased sensitivities to DNA-damaging drugs. Inhibition of CDK7 preferentially represses the expression of genes in the DNA-damage-repair pathways and impairs the activity of homologous recombination. Low-dose CDK7 inhibitor treatment sensitizes cancer cells to PARP inhibitor-induced DNA damage and cell death. Our analysis provides genomic information for identification and prioritization of drug targets for CDKs and reveals rationales for treatment strategies.

Systematic characterization of mutations altering protein degradation in human cancers

The ubiquitin-proteasome system (UPS) is the primary route for selective protein degradation in human cells. The UPS is an attractive target for novel cancer therapies, but the precise UPS genes and substrates important for cancer growth are incompletely understood. Leveraging multi-omics data across more than 9,000 human tumors and 33 cancer types, we found that over 19% of all cancer driver genes affect UPS function. We implicate transcription factors as important substrates and show that c-Myc stability is modulated by CUL3. Moreover, we developed a deep learning model (deepDegron) to identify mutations that result in degron loss and experimentally validated the prediction that gain-of-function truncating mutations in GATA3 and PPM1D result in increased protein stability. Last, we identified UPS driver genes associated with prognosis and the tumor microenvironment. This study demonstrates the important role of UPS dysregulation in human cancer and underscores the potential therapeutic utility of targeting the UPS.

Resting-State Co-activation Patterns as Promising Candidates for Prediction of Alzheimer's Disease in Aged Mice

Alzheimer’s disease (AD), a neurodegenerative disorder marked by accumulation of extracellular amyloid-β (Aβ) plaques leads to progressive loss of memory and cognitive function. Resting-state fMRI (RS-fMRI) studies have provided links between these two observations in terms of disruption of default mode and task-positive resting-state networks (RSNs). Important insights underlying these disruptions were recently obtained by investigating dynamic fluctuations in RS-fMRI signals in old TG2576 mice (a mouse model of amyloidosis) using a set of quasi-periodic patterns (QPP). QPPs represent repeating spatiotemporal patterns of neural activity of predefined temporal length. In this article, we used an alternative methodology of co-activation patterns (CAPs) that represent instantaneous and transient brain configurations that are likely contributors to the emergence of commonly observed RSNs and QPPs. We followed a recently published approach for obtaining CAPs that divided all time frames, instead of those corresponding to supra-threshold activations of a seed region as done traditionally, to extract CAPs from RS-fMRI recordings in 10 TG2576 female mice and eight wild type littermates at 18 months of age. Subsequently, we matched the CAPs from the two groups using the Hungarian method and compared the temporal (duration, occurrence rate) and the spatial (lateralization of significantly co-activated and co-deactivated voxels) properties of matched CAPs. We found robust differences in the spatial components of matched CAPs. Finally, we used supervised learning to train a classifier using either the temporal or the spatial component of CAPs to distinguish the transgenic mice from the WT. We found that while duration and occurrence rates of all CAPs performed the classification with significantly higher accuracy than the chance-level, blood oxygen level-dependent (BOLD) signals of significantly activated voxels from individual CAPs turned out to be a significantly better predictive feature demonstrating a near-perfect classification accuracy. Our results demonstrate resting-state co-activation patterns are a promising candidate in the development of a diagnostic, and potentially, prognostic RS-fMRI biomarker of AD.

Organized Resting-state Functional Dysconnectivity of the Prefrontal Cortex in Patients with Schizophrenia

Schizophrenia has prominent functional dysconnectivity, especially in the prefrontal cortex (PFC). However, it is unclear whether in the same group of patients with schizophrenia, PFC functional dysconnectivity appears in an organized manner or is stochastically located in different subregions. By investigating the resting-state functional connectivity (rsFC) of each PFC subregion from the Brainnetome atlas in 40 schizophrenia patients and 40 healthy subjects, we found 24 altered connections in schizophrenia, and the connections were divided into four categories by a clustering analysis: increased connections within the PFC, increased connections between the inferior PFC and the thalamus/striatum, reduced connections between the PFC and the motor control areas, and reduced connections between the orbital PFC and the emotional perception regions. In addition, the four categories of rsFC showed distinct cognitive engagement patterns. Our findings suggest that PFC subregions have specific functional dysconnectivity patterns in schizophrenia and may reflect heterogeneous symptoms and cognitive deficits in schizophrenia.

Improved long-term time-series predictions of total blood use data from England

BACKGROUND

Red blood cells are essential for modern medicine but managing their collection and supply to cope with fluctuating demands represents a major challenge. As deterministic models based on predicted population changes have been problematic, there remains a need for more precise and reliable prediction of use. Here, we develop three new time-series methods to predict red cell use 4 to 52 weeks ahead.

STUDY DESIGN AND METHODS

From daily aggregates of red blood cell (RBC) units issued from 2005 to 2011 from the NHS Blood and Transplant, we generated a new set of non-overlapping weekly data by summing the daily data over 7 days and derived the average blood use per week over 4-week and 52-week periods. We used three new methods for linear prediction of blood use by computing the coefficients using Minimum Mean Squared Error (MMSE) algorithm.

RESULTS

We optimized the time-window size, order of the prediction, and order of the polynomial fit for our data set. By exploiting the annual periodicity of the data, we achieved significant improvements in long-term predictions, as well as modest improvements in short-term predictions. The new methods predicted mean RBC use with a standard deviation of the percentage error of 2.5% for 4 weeks ahead and 3.4% for 52 weeks ahead.

CONCLUSION

This paradigm allows short- and long-term prediction of RBC use and could provide reliable and precise prediction up to 52 weeks ahead to improve the efficiency of blood services and sufficiency of blood supply with reduced costs.

Fine-grain atlases of functional modes for fMRI analysis

Population imaging markedly increased the size of functional-imaging datasets, shedding new light on the neural basis of inter-individual differences. Analyzing these large data entails new scalability challenges, computational and statistical. For this reason, brain images are typically summarized in a few signals, for instance reducing voxel-level measures with brain atlases or functional modes. A good choice of the corresponding brain networks is important, as most data analyses start from these reduced signals. We contribute finely-resolved atlases of functional modes, comprising from 64 to 1024 networks. These dictionaries of functional modes (DiFuMo) are trained on millions of fMRI functional brain volumes of total size 2.4 ​TB, spanned over 27 studies and many research groups. We demonstrate the benefits of extracting reduced signals on our fine-grain atlases for many classic functional data analysis pipelines: stimuli decoding from 12,334 brain responses, standard GLM analysis of fMRI across sessions and individuals, extraction of resting-state functional-connectomes biomarkers for 2500 individuals, data compression and meta-analysis over more than 15,000 statistical maps. In each of these analysis scenarii, we compare the performance of our functional atlases with that of other popular references, and to a simple voxel-level analysis. Results highlight the importance of using high-dimensional "soft" functional atlases, to represent and analyze brain activity while capturing its functional gradients. Analyses on high-dimensional modes achieve similar statistical performance as at the voxel level, but with much reduced computational cost and higher interpretability. In addition to making them available, we provide meaningful names for these modes, based on their anatomical location. It will facilitate reporting of results.

Data Modeling With Polynomial Representations and Autoregressive Time-Series Representations, and Their Connections

Two of the data modelling techniques - polynomial representation and time-series representation - are explored in this paper to establish their connections and differences. All theoretical studies are based on uniformly sampled data in the absence of noise. This paper proves that all data from an underlying polynomial model of finite degree [Formula: see text] can be represented perfectly by an autoregressive time-series model of order [Formula: see text] and a constant term [Formula: see text] as in equation (2). Furthermore, all polynomials of degree [Formula: see text] are shown to give rise to the same set of time-series coefficients of specific forms with the only possible difference being in the constant term [Formula: see text]. It is also demonstrated that time-series with either non-integer coefficients or integer coefficients not of the aforementioned specific forms represent polynomials of infinite degree. Six numerical explorations, with both generated data and real data, including the UK data and US data on the current Covid-19 incidence, are presented to support the theoretical findings. It is shown that all polynomials of degree [Formula: see text] can be represented by an all-pole filter with [Formula: see text] repeated roots (or poles) at [Formula: see text]. Theoretically, all noise-free data representable by a finite order all-pole filter, whether they come from finite degree or infinite degree polynomials, can be described exactly by a finite order AR time-series; if the values of polynomial coefficients are not of special interest in any data modelling, one may use time-series representations for data modelling.

Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands

A diverse set of white matter connections supports seamless transitions between cognitive states. However, it remains unclear how these connections guide the temporal progression of large-scale brain activity patterns in different cognitive states. Here, we analyze the brain's trajectories across a set of single time point activity patterns from functional magnetic resonance imaging data acquired during the resting state and an n-back working memory task. We find that specific temporal sequences of brain activity are modulated by cognitive load, associated with age, and related to task performance. Using diffusion-weighted imaging acquired from the same subjects, we apply tools from network control theory to show that linear spread of activity along white matter connections constrains the probabilities of these sequences at rest, while stimulus-driven visual inputs explain the sequences observed during the n-back task. Overall, these results elucidate the structural underpinnings of cognitively and developmentally relevant spatiotemporal brain dynamics.

Minimax Bridgeness-Based Clustering for Hyperspectral Data

Hyperspectral (HS) imaging has been used extensively in remote sensing applications like agriculture, forestry, geology and marine science. HS pixel classification is an important task to help identify different classes of materials within a scene, such as different types of crops on a farm. However, this task is significantly hindered by the fact that HS pixels typically form high-dimensional clusters of arbitrary sizes and shapes in the feature space spanned by all spectral channels. This is even more of a challenge when ground truth data is difficult to obtain and when there is no reliable prior information about these clusters (e.g., number, typical shape, intrinsic dimensionality). In this letter, we present a new graph-based clustering approach for hyperspectral data mining that does not require ground truth data nor parameter tuning. It is based on the minimax distance, a measure of similarity between vertices on a graph. Using the silhouette index, we demonstrate that the minimax distance is more suitable to identify clusters in raw hyperspectral data than two other graph-based similarity measures: mutual proximity and shared nearest neighbours. We then introduce the minimax bridgeness-based clustering approach, and we demonstrate that it can discover clusters of interest in hyperspectral data better than comparable approaches.

Resolution-based spectral clustering for brain parcellation using functional MRI

BACKGROUND

Brain parcellation is important for exploiting neuroimaging data. Variability in physiology between individuals has led to the need for data-driven approaches to parcellation, with recent research focusing on simultaneously estimating and partitioning the network structure of the brain.

NEW METHOD

We view data preprocessing, parcellation, and parcel validation from the perspective of predictive modeling. The goal is to identify parcels in a way that best generalizes to unseen data. We utilize an uncertainty quantification approach from image science to define parcels as groups of unresolvable variables in the predictive model. Model parameters are chosen via cross-validation. Parcellation results are compared based on both their repeatability as well as their ability to describe held-out data.

RESULTS

The approach provides insight and strategies for open questions such as the choice of evaluation metrics, selection of model order, and the optimal tuning of preprocessing steps for functional imaging data.

COMPARISON WITH EXISTING METHODS

We compare new and established approaches using functional imaging data, where we find the proposed approach produces parcellations which are both more accurate and more repeatable than the current baseline clustering method. The metrics also demonstrate potential problems with overfitting for the baseline method, and with bias for other methods.

CONCLUSIONS

Clustering of resolution provides a principled and robust approach to brain parcellation which improves on the current baseline.

Prediction paradigm involving time series applied to total blood issues data from England

BACKGROUND

Blood products are essential for modern medicine, but managing their collection and supply in the face of fluctuating demands represents a major challenge. As deterministic models based on predicted changes in population have been problematic, there remains a need for more precise and reliable prediction of demands. Here, we propose a paradigm incorporating four different time‐series methods to predict red blood cell (RBC) issues 4 to 24 weeks ahead.

STUDY DESIGN AND METHODS

We used daily aggregates of RBC units issued from 2005 to 2011 from the National Health Service Blood and Transplant. We generated a new set of nonoverlapping weekly data by summing the daily data over 7 days and derived the average blood issues per week over 4‐week periods. We used four methods for linear prediction of blood demand by computing the coefficients with the minimum mean squared error and weighted least squares error algorithms.

RESULTS

We optimized the time‐window size, order of the prediction, and order of the polynomial fit for our data set. The four time‐series methods, essentially using different weightings to data points, gave very similar results and predicted mean RBC issues with a standard deviation of the percentage error of 3.0% for 4 weeks ahead and 4.0% for 24 weeks ahead.

CONCLUSION

This paradigm allows prediction of demand for RBCs and could be developed to provide reliable and precise prediction up to 24 weeks ahead to improve the efficiency of blood services and sufficiency of blood supply with reduced costs.

Global ecosystem thresholds driven by aridity

Aridity, which is increasing worldwide because of climate change, affects the structure and functioning of dryland ecosystems. Whether aridification leads to gradual (versus abrupt) and systemic (versus specific) ecosystem changes is largely unknown. We investigated how 20 structural and functional ecosystem attributes respond to aridity in global drylands. Aridification led to systemic and abrupt changes in multiple ecosystem attributes. These changes occurred sequentially in three phases characterized by abrupt decays in plant productivity, soil fertility, and plant cover and richness at aridity values of 0.54, 0.7, and 0.8, respectively. More than 20% of the terrestrial surface will cross one or several of these thresholds by 2100, which calls for immediate actions to minimize the negative impacts of aridification on essential ecosystem services for the more than 2 billion people living in drylands.

Visual temporal frequency preference shows a distinct cortical architecture using fMRI

Studies of visual temporal frequency preference typically examine frequencies under 20 Hz and measure local activity to evaluate the sensitivity of different cortical areas to variations in temporal frequencies. Most of these studies have not attempted to map preferred temporal frequency within and across visual areas, nor have they explored in detail, stimuli at gamma frequency, which recent research suggests may have potential clinical utility. In this study, we address this gap by using functional magnetic resonance imaging (fMRI) to measure response to flickering visual stimuli varying in frequency from 1 to 40 Hz. We apply stimulation in both a block design to examine task response and a steady-state design to examine functional connectivity. We observed distinct activation patterns between 1 Hz and 40 Hz stimuli. We also found that the correlation between medial thalamus and visual cortex was modulated by the temporal frequency. The modulation functions and tuned frequencies are different for the visual activity and thalamo-visual correlations. Using both fMRI activity and connectivity measurements, we show evidence for a temporal frequency specific organization across the human visual system.

Infraslow State Fluctuations Govern Spontaneous fMRI Network Dynamics

Spontaneous brain activity as assessed with resting-state fMRI exhibits rich spatiotemporal structure. However, the principles by which brain-wide patterns of spontaneous fMRI activity reconfigure and interact with each other remain unclear. We used a framewise clustering approach to map spatiotemporal dynamics of spontaneous fMRI activity with voxel resolution in the resting mouse brain. We show that brain-wide patterns of fMRI co-activation can be reliably mapped at the group and subject level, defining a restricted set of recurring brain states characterized by rich network structure. Importantly, we document that the identified fMRI states exhibit contrasting patterns of functional activity and coupled infraslow network dynamics, with each network state occurring at specific phases of global fMRI signal fluctuations. Finally, we show that autism-associated genetic alterations entail the engagement of atypical functional states and altered infraslow network dynamics. Our results reveal a novel set of fundamental principles guiding the spatiotemporal organization of resting-state fMRI activity and its disruption in brain disorders.

Comparing 3D Genome Organization in Multiple Species Using Phylo-HMRF

Recent whole-genome mapping approaches for the chromatin interactome have offered new insights into 3D genome organization. However, our knowledge of the evolutionary patterns of 3D genome in mammals remains limited. In particular, there are no existing phylogenetic-model-based methods to analyze chromatin interactions as continuous features. Here, we develop phylogenetic hidden Markov random field (Phylo-HMRF) to identify evolutionary patterns of 3D genome based on multi-species Hi-C data by jointly utilizing spatial constraints among genomic loci and continuous-trait evolutionary models. We used Phylo-HMRF to uncover cross-species 3D genome patterns based on Hi-C data from the same cell type in four primate species (human, chimpanzee, bonobo, and gorilla). The identified evolutionary patterns of 3D genome correlate with features of genome structure and function. This work provides a new framework to analyze multi-species continuous genomic features with spatial constraints and has the potential to help reveal the evolutionary principles of 3D genome organization.

A connectivity-based parcellation improved functional representation of the human cerebellum

The cerebellum is traditionally well known for its role in motor learning and coordination. Recently, it is recognized that the function of the cerebellum is highly diverse and extends to non-motor domains, such as working memory, emotion and language. The diversity of the cerebellum can be appreciated by examining its extensive connectivity to the cerebral regions selective for both motor and cognitive functions. Importantly, the pattern of cerebro-cerebellar connectivity is specific and distinct to different cerebellar subregions. Therefore, to understand the cerebellum and the various functions it involves, it is essential to identify and differentiate its subdivisions. However, most studies are still referring the cerebellum as one brain structure or by its gross anatomical subdivisions, which does not necessarily reflect the functional mapping of the cerebellum. We here employed a data-driven method to generate a functional connectivity-based parcellation of the cerebellum. Our results demonstrated that functional connectivity-based atlas is superior to existing atlases in regards to cluster homogeneity, accuracy of functional connectivity representation and individual identification. Furthermore, our functional atlas improves statistical results of task fMRI analyses, as compared to the standard voxel-based approach and existing atlases. Our detailed functional parcellation provides a valuable tool for elucidating the functional diversity and connectivity of the cerebellum as well as its network relationships with the whole brain.

Impaired cerebellar Purkinje cell potentiation generates unstable spatial map orientation and inaccurate navigation

Cerebellar activity supported by PKC-dependent long-term depression in Purkinje cells (PCs) is involved in the stabilization of self-motion based hippocampal representation, but the existence of cerebellar processes underlying integration of allocentric cues remains unclear. Using mutant-mice lacking PP2B in PCs (L7-PP2B mice) we here assess the role of PP2B-dependent PC potentiation in hippocampal representation and spatial navigation. L7-PP2B mice display higher susceptibility to spatial map instability relative to the allocentric cue and impaired allocentric as well as self-motion goal-directed navigation. These results indicate that PP2B-dependent potentiation in PCs contributes to maintain a stable hippocampal representation of a familiar environment in an allocentric reference frame as well as to support optimal trajectory toward a goal during navigation.

It is known that Purkinje cell PKC-dependent depression is involved in the stabilization of self-motion based hippocampal representation. Here the authors describe decreased stability of hippocampal place cells based on allocentric cues in mice lacking Purkinje cell PP2B-dependent potentiation.

Stochastic Rank Aggregation for the Identification of Functional Neuromarkers

The main challenge in analysing functional magnetic resonance imaging (fMRI) data from extended samples of subject (N > 100) is to extract as much relevant information as possible from big amounts of noisy data. When studying neurodegenerative diseases with resting-state fMRI, one of the objectives is to determine regions with abnormal background activity with respect to a healthy brain and this is often attained with comparative statistical models applied to single voxels or brain parcels within one or several functional networks. In this work, we propose a novel approach based on clustering and stochastic rank aggregation to identify parcels that exhibit a coherent behaviour in groups of subjects affected by the same disorder and apply it to default-mode network independent component maps from resting-state fMRI data sets. Brain voxels are partitioned into parcels through k-means clustering, then solutions are enhanced by means of consensus techniques. For each subject, clusters are ranked according to their median value and a stochastic rank aggregation method, TopKLists, is applied to combine the individual rankings within each class of subjects. For comparison, the same approach was tested on an anatomical parcellation. We found parcels for which the rankings were different among control subjects and subjects affected by Parkinson's disease and amyotrophic lateral sclerosis and found evidence in literature for the relevance of top ranked regions in default-mode brain activity. The proposed framework represents a valid method for the identification of functional neuromarkers from resting-state fMRI data, and it might therefore constitute a step forward in the development of fully automated data-driven techniques to support early diagnoses of neurodegenerative diseases.

On the spectral signature of melanoma: a non-parametric classification framework for cancer detection in hyperspectral imaging of melanocytic lesions

Early detection and diagnosis is a must in secondary prevention of melanoma and other cancerous lesions of the skin. In this work, we present an online, reservoir-based, non-parametric estimation and classification model that allows for this functionality on pigmented lesions, such that detection thresholding can be tuned to maximize accuracy and/or minimize overall false negative rates. This system has been tested in a dataset consisting of 116 patients and a total of 124 hyperspectral images of nevi, raised nevi and melanomas, detecting up to 100% of the suspicious lesions at the expense of some false positives.

Context-Dependent Risk Aversion: A Model-Based Approach

Most research on risk aversion in behavioral science with human subjects has focused on a component of risk aversion that does not adapt itself to context. More recently, studies have explored risk aversion adaptation to changing circumstances in sequential decision-making tasks. It is an open question whether one can identify evidence, at the single subject level, for such risk aversion adaptation. We conducted a behavioral experiment on human subjects, using a sequential decision making task. We developed a model-based approach for estimating the adaptation of risk-taking behavior with single-trial resolution by modeling a subject's goals and internal representation of task contingencies. Using this model-based approach, we estimated the subject-specific adaptation of risk aversion depending on the current task context. We found striking inter-subject variations in the adaptation of risk-taking behavior. We show that these differences can be explained by differences in subjects' internal representations of task contingencies and goals. We discuss that the proposed approach can be adapted to a wide range of experimental paradigms and be used to analyze behavioral measures other than risk aversion.