Data-driven fine-grained region discovery in the mouse brain with transformers.

Technologies such as spatial transcriptomics offer unique opportunities to define the spatial organization of the mouse brain. We developed an unsupervised training scheme and novel transformer-based deep learning architecture to detect spatial domains in mouse whole-brain spatial transcriptomics data. Our model learns local representations of molecular and cellular statistical patterns. These local representations can be clustered to identify spatial domains within the brain from coarse to fine-grained. Discovered domains are spatially regular, even with several hundreds of spatial clusters. They are also consistent with existing anatomical ontologies such as the Allen Mouse Brain Common Coordinate Framework version 3 (CCFv31) and can be visually interpreted at the cell type or transcript level. We demonstrate our method can be used to identify previously uncatalogued subregions, such as in the midbrain, where we uncover gradients of inhibitory neuron complexity and abundance. We apply our method to a separate multi-animal whole-brain spatial transcriptomic dataset and observe that inclusion of both sagittal and coronal tissue slices in region identification improves correspondence of spatial domains to CCF.

Metabolic reprogramming of tumor-associated macrophages using glutamine antagonist JHU083 drives tumor immunity in myeloid-rich prostate and bladder cancer tumors.

Glutamine metabolism in tumor microenvironments critically regulates anti-tumor immunity. Using glutamine-antagonist prodrug JHU083, we report potent tumor growth inhibition in urologic tumors by JHU083-reprogrammed tumor-associated macrophages (TAMs) and tumor-infiltrating monocytes (TIMs). We show JHU083-mediated glutamine antagonism in tumor microenvironments induces TNF, pro-inflammatory, and mTORC1 signaling in intratumoral TAM clusters. JHU083-reprogrammed TAMs also exhibit increased tumor cell phagocytosis and diminished pro-angiogenic capacities. In vivo inhibition of TAM glutamine consumption resulted in increased glycolysis, a broken TCA cycle, and purine metabolism disruption. Although the anti-tumor effect of glutamine antagonism on tumor-infiltrating T cells was moderate, JHU083 promoted a stem cell-like phenotype in CD8+ T cells and decreased Treg abundance. Finally, JHU083 caused a ubiquitous shutdown in glutamine utilizing metabolic pathways in tumor cells, leading to reduced HIF-1alpha, c-MYC phosphorylation, and induction of tumor cell apoptosis, all key anti-tumor features.

Building an ab initio solvated DNA model using Euclidean neural networks.

Accurately modeling large biomolecules such as DNA from first principles is fundamentally challenging due to the steep computational scaling of ab initio quantum chemistry methods. This limitation becomes even more prominent when modeling biomolecules in solution due to the need to include large numbers of solvent molecules. We present a machine-learned electron density model based on a Euclidean neural network framework that includes a built-in understanding of equivariance to model explicitly solvated double-stranded DNA. By training the machine learning model using molecular fragments that sample the key DNA and solvent interactions, we show that the model predicts electron densities of arbitrary systems of solvated DNA accurately, resolves polarization effects that are neglected by classical force fields, and captures the physics of the DNA-solvent interaction at the ab initio level.

Staphylococcus aureus proteases trigger eosinophil-mediated skin inflammation.

Staphylococcus aureus skin colonization and eosinophil infiltration are associated with many inflammatory skin disorders, including atopic dermatitis, bullous pemphigoid, Netherton's syndrome, and prurigo nodularis. However, whether there is a relationship between S. aureus and eosinophils and how this interaction influences skin inflammation is largely undefined. We show in a preclinical mouse model that S. aureus epicutaneous exposure induced eosinophil-recruiting chemokines and eosinophil infiltration into the skin. Remarkably, we found that eosinophils had a comparable contribution to the skin inflammation as T cells, in a manner dependent on eosinophil-derived IL-17A and IL-17F production. Importantly, IL-36R signaling induced CCL7-mediated eosinophil recruitment to the inflamed skin. Last, S. aureus proteases induced IL-36α expression in keratinocytes, which promoted infiltration of IL-17-producing eosinophils. Collectively, we uncovered a mechanism for S. aureus proteases to trigger eosinophil-mediated skin inflammation, which has implications in the pathogenesis of inflammatory skin diseases.

Functional network collapse in neurodegenerative disease.

Cognitive and behavioral deficits in Alzheimer's disease (AD) and frontotemporal dementia (FTD) result from brain atrophy and altered functional connectivity. However, it is unclear how atrophy relates to functional connectivity disruptions across dementia subtypes and stages. We addressed this question using structural and functional MRI from 221 patients with AD (n=82), behavioral variant FTD (n=41), corticobasal syndrome (n=27), nonfluent (n=34) and semantic (n=37) variant primary progressive aphasia, and 100 cognitively normal individuals. Using partial least squares regression, we identified three principal structure-function components. The first component showed overall atrophy correlating with primary cortical hypo-connectivity and subcortical/association cortical hyper-connectivity. Components two and three linked focal syndrome-specific atrophy to peri-lesional hypo-connectivity and distal hyper-connectivity. Structural and functional component scores predicted global and domain-specific cognitive deficits. Anatomically, functional connectivity changes reflected alterations in specific brain activity gradients. Eigenmode analysis identified temporal phase and amplitude collapse as an explanation for atrophy-driven functional connectivity changes.

FTLD targets brain regions expressing recently evolved genes.

In frontotemporal lobar degeneration (FTLD), pathological protein aggregation is associated with a decline in human-specialized social-emotional and language functions. Most disease protein aggregates contain either TDP-43 (FTLD-TDP) or tau (FTLD-tau). Here, we explored whether FTLD targets brain regions that express genes containing human accelerated regions (HARs), conserved sequences that have undergone positive selection during recent human evolution. To this end, we used structural neuroimaging from patients with FTLD and normative human regional transcriptomic data to identify genes expressed in FTLD-targeted brain regions. We then integrated primate comparative genomic data to test our hypothesis that FTLD targets brain regions expressing recently evolved genes. In addition, we asked whether genes expressed in FTLD-targeted brain regions are enriched for genes that undergo cryptic splicing when TDP-43 function is impaired. We found that FTLD-TDP and FTLD-tau subtypes target brain regions that express overlapping and distinct genes, including many linked to neuromodulatory functions. Genes whose normative brain regional expression pattern correlated with FTLD cortical atrophy were strongly associated with HARs. Atrophy-correlated genes in FTLD-TDP showed greater overlap with TDP-43 cryptic splicing genes compared with atrophy-correlated genes in FTLD-tau. Cryptic splicing genes were enriched for HAR genes, and vice versa, but this effect was due to the confounding influence of gene length. Analyses performed at the individual-patient level revealed that the expression of HAR genes and cryptically spliced genes within putative regions of disease onset differed across FTLD-TDP subtypes. Overall, our findings suggest that FTLD targets brain regions that have undergone recent evolutionary specialization and provide intriguing potential leads regarding the transcriptomic basis for selective vulnerability in distinct FTLD molecular-anatomical subtypes.

Network Connectivity Alterations across the MAPT Mutation Clinical Spectrum.

OBJECTIVE: Microtubule-associated protein tau (MAPT) mutations cause frontotemporal lobar degeneration, and novel biomarkers are urgently needed for early disease detection. We used task-free functional magnetic resonance imaging (fMRI) mapping, a promising biomarker, to analyze network connectivity in symptomatic and presymptomatic MAPT mutation carriers. METHODS: We compared cross-sectional fMRI data between 17 symptomatic and 39 presymptomatic carriers and 81 controls with (1) seed-based analyses to examine connectivity within networks associated with the 4 most common MAPT-associated clinical syndromes (ie, salience, corticobasal syndrome, progressive supranuclear palsy syndrome, and default mode networks) and (2) whole-brain connectivity analyses. We applied K-means clustering to explore connectivity heterogeneity in presymptomatic carriers at baseline. Neuropsychological measures, plasma neurofilament light chain, and gray matter volume were compared at baseline and longitudinally between the presymptomatic subgroups defined by their baseline whole-brain connectivity profiles. RESULTS: Symptomatic and presymptomatic carriers had connectivity disruptions within MAPT-syndromic networks. Compared to controls, presymptomatic carriers showed regions of connectivity alterations with age. Two presymptomatic subgroups were identified by clustering analysis, exhibiting predominantly either whole-brain hypoconnectivity or hyperconnectivity at baseline. At baseline, these two presymptomatic subgroups did not differ in neuropsychological measures, although the hypoconnectivity subgroup had greater plasma neurofilament light chain levels than controls. Longitudinally, both subgroups showed visual memory decline (vs controls), yet the subgroup with baseline hypoconnectivity also had worsening verbal memory and neuropsychiatric symptoms, and extensive bilateral mesial temporal gray matter decline. INTERPRETATION: Network connectivity alterations arise as early as the presymptomatic phase. Future studies will determine whether presymptomatic carriers' baseline connectivity profiles predict symptomatic conversion. ANN NEUROL 2023;94:632-646.

Machine Learning for Uncovering Biological Insights in Spatial Transcriptomics Data.

Development and homeostasis in multicellular systems both require exquisite control over spatial molecular pattern formation and maintenance. Advances in spatially-resolved and high-throughput molecular imaging methods such as multiplexed immunofluorescence and spatial transcriptomics (ST) provide exciting new opportunities to augment our fundamental understanding of these processes in health and disease. The large and complex datasets resulting from these techniques, particularly ST, have led to rapid development of innovative machine learning (ML) tools primarily based on deep learning techniques. These ML tools are now increasingly featured in integrated experimental and computational workflows to disentangle signals from noise in complex biological systems. However, it can be difficult to understand and balance the different implicit assumptions and methodologies of a rapidly expanding toolbox of analytical tools in ST. To address this, we summarize major ST analysis goals that ML can help address and current analysis trends. We also describe four major data science concepts and related heuristics that can help guide practitioners in their choices of the right tools for the right biological questions.

Dysfunctional Cortical Gradient Topography in Treatment-Resistant Major Depressive Disorder.

BACKGROUND: Treatment-resistant depression (TRD) refers to patients with major depressive disorder who do not remit after 2 or more antidepressant trials. TRD is common and highly debilitating, but its neurobiological basis remains poorly understood. Recent neuroimaging studies have revealed cortical connectivity gradients that dissociate primary sensorimotor areas from higher-order associative cortices. This fundamental topography determines cortical information flow and is affected by psychiatric disorders. We examined how TRD impacts gradient-based hierarchical cortical organization. METHODS: In this secondary study, we analyzed resting-state functional magnetic resonance imaging data from a mindfulness-based intervention enrolling 56 patients with TRD and 28 healthy control subjects. Using gradient extraction tools, baseline measures of cortical gradient dispersion within and between functional brain networks were derived, compared across groups, and associated with graph theoretical measures of network topology. In patients, correlation analyses were used to associate measures of cortical gradient dispersion with clinical measures of anxiety, depression, and mindfulness at baseline and following the intervention. RESULTS: Cortical gradient dispersion was reduced within major intrinsic brain networks in patients with TRD. Reduced cortical gradient dispersion correlated with increased network degree assessed through graph theory-based measures of network topology. Lower dispersion among default mode, control, and limbic network nodes related to baseline levels of trait anxiety, depression, and mindfulness. Patients' baseline limbic network dispersion predicted trait anxiety scores 24 weeks after the intervention. CONCLUSIONS: Our findings provide preliminary support for widespread alterations in cortical gradient architecture in TRD, implicating a significant role for transmodal and limbic networks in mediating depression, anxiety, and lower mindfulness in patients with TRD.

Accurate Hellmann-Feynman forces from density functional calculations with augmented Gaussian basis sets.

The Hellmann-Feynman (HF) theorem provides a way to compute forces directly from the electron density, enabling efficient force calculations for large systems through machine learning (ML) models for the electron density. The main issue holding back the general acceptance of the HF approach for atom-centered basis sets is the well-known Pulay force which, if naively discarded, typically constitutes an error upward of 10 eV/Å in forces. In this work, we demonstrate that if a suitably augmented Gaussian basis set is used for density functional calculations, the Pulay force can be suppressed, and HF forces can be computed as accurately as analytical forces with state-of-the-art basis sets, allowing geometry optimization and molecular dynamics to be reliably performed with HF forces. Our results pave a clear path forward for the accurate and efficient simulation of large systems using ML densities and the HF theorem.

Increasing empathic concern relates to salience network hyperconnectivity in cognitively healthy older adults with elevated amyloid-ß burden.

Enhanced emotional empathy, the ability to share others' affective experiences, can be a feature of Alzheimer's disease (AD), but whether emotional empathy increases in the preclinical phase of the disease is unknown. We measured emotional empathy over time (range = 0 - 7.3 years, mean = 2.4 years) in 86 older adults during a period in which they were cognitively healthy, functionally normal, and free of dementia symptoms. For each participant, we computed longitudinal trajectories for empathic concern (i.e., an other-oriented form of emotional empathy that promotes prosocial actions) and emotional contagion (i.e., a self-focused form of emotional empathy often accompanied by feelings of distress) from informant ratings of participants' empathy on the Interpersonal Reactivity Index. Amyloid-ß (Aß) positron emission tomography (PET) scans were used to classify participants as either Aß positive (Aß+, n = 23) or negative (Aß-, n = 63) based on Aß-PET cortical binding. Participants also underwent structural and task-free functional magnetic resonance imaging approximately two years on average after their last empathy assessment, at which time most participants remained cognitively healthy. Results indicated that empathic concern, but not emotional contagion, increased more over time in Aß+ participants than in Aß- participants despite no initial group difference at the first measurement. Higher connectivity between certain salience network node-pairs (i.e., pregenual anterior cingulate cortex and periaqueductal gray) predicted longitudinal increases in empathic concern in the Aß+ group but not in the Aß- group. The Aß+ participants also had higher overall salience network connectivity than Aß- participants despite no differences in gray matter volume. These results suggest gains in empathic concern may be a very early feature of AD pathophysiology that relates to hyperconnectivity in the salience network, a system that supports emotion generation and interoception. A better understanding of emotional empathy trajectories in the early stages of AD pathophysiology will broaden the lens on preclinical AD changes and help clinicians to identify older adults who should be screened for AD biomarkers.

RAS oncogene signal strength regulates matrisomal gene expression and tumorigenicity of mouse keratinocytes.

Environmental and molecular carcinogenesis are linked by the discovery that chemical carcinogen induced-mutations in the Hras or Kras genes drives tumor development in mouse skin. Importantly, enhanced expression or allele amplification of the mutant Ras gene contributes to selection of initiated cells, tumor persistence, and progression. To explore the consequences of Ras oncogene signal strength, primary keratinocytes were isolated and cultured from the LSL-HrasG12D and LSL-KrasG12D C57BL/6J mouse models and the mutant allele was activated by adeno-Cre recombinase. Keratinocytes expressing one (H) or two (HH) mutant alleles of HrasG12D, one KrasG12D allele (K), or one of each (HK) were studied. All combinations of activated Ras alleles stimulated proliferation and drove transformation marker expression, but only HH and HK formed tumors. HH, HK, and K sustained long-term keratinocyte growth in vitro, while H and WT could not. RNA-Seq yielded two distinct gene expression profiles; HH, HK, and K formed one cluster while H clustered with WT. Weak MAPK activation was seen in H keratinocytes but treatment with a BRAF inhibitor enhanced MAPK signaling and facilitated tumor formation. K keratinocytes became tumorigenic when they were isolated from mice where the LSL-KrasG12D allele was backcrossed from the C57BL/6 onto the FVB/N background. All tumorigenic keratinocytes but not the non-tumorigenic precursors shared a common remodeling of matrisomal gene expression that is associated with tumor formation. Thus, RAS oncogene signal strength determines cell-autonomous changes in initiated cells that are critical for their tumor-forming potential.

Overview of the Nomenclature and Network of Contributors to the Development of Bioreactors for Human Gut Simulation Using Bibliometric Tools: A Fragmented Landscape.

The evolution of complex in vitro models of the human gastrointestinal system to interrogate the biochemical functionality of the gut microbiome has augmented our understanding of its role in human physiology and pathology. With 5718 authors from 52 countries, gut bioreactor research reflects the growing awareness of our need to understand the contribution of the gut microbiome to human health. Although a large body of knowledge has been generated from in vitro models, it is scattered and defined by application-specific terminologies. To better grasp the capacity of bioreactors and further our knowledge of the human gastrointestinal system, we have conducted a cross-field bibliometric search and mapped the evolution of human gastrointestinal in vitro research. We present reference material with the aim of identifying key authors and bioreactor types to enable researchers to make decisions regarding the choice of method for simulating the human gut in the context of microbiome functionality.

Predicting accurate ab initio DNA electron densities with equivariant neural networks.

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

A dynamic gradient architecture generates brain activity states.

The human brain exhibits a diverse yet constrained range of activity states. While these states can be faithfully represented in a low-dimensional latent space, our understanding of the constitutive functional anatomy is still evolving. Here we applied dimensionality reduction to task-free and task fMRI data to address whether latent dimensions reflect intrinsic systems and if so, how these systems may interact to generate different activity states. We find that each dimension represents a dynamic activity gradient, including a primary unipolar sensory-association gradient underlying the global signal. The gradients appear stable across individuals and cognitive states, while recapitulating key functional connectivity properties including anticorrelation, modularity, and regional hubness. We then use dynamical systems modeling to show that gradients causally interact via state-specific coupling parameters to create distinct brain activity patterns. Together, these findings indicate that a set of dynamic, intrinsic spatial gradients interact to determine the repertoire of possible brain activity states.

Next-generation sequencing of cerebrospinal fluid for clinical molecular diagnostics in pediatric, adolescent and young adult brain tumor patients.

BACKGROUND: Safe sampling of central nervous system tumor tissue for diagnostic purposes may be difficult if not impossible, especially in pediatric patients, and an unmet need exists to develop less invasive diagnostic tests. METHODS: We report our clinical experience with minimally invasive molecular diagnostics using a clinically validated assay for sequencing of cerebrospinal fluid (CSF) cell-free DNA (cfDNA). All CSF samples were collected as part of clinical care, and results reported to both clinicians and patients/families. RESULTS: We analyzed 64 CSF samples from 45 pediatric, adolescent and young adult (AYA) patients (pediatric = 25; AYA = 20) with primary and recurrent brain tumors across 12 histopathological subtypes including high-grade glioma (n = 10), medulloblastoma (n = 10), pineoblastoma (n = 5), low-grade glioma (n = 4), diffuse leptomeningeal glioneuronal tumor (DLGNT) (n = 4), retinoblastoma (n = 4), ependymoma (n = 3), and other (n = 5). Somatic alterations were detected in 30/64 samples (46.9%) and in at least one sample per unique patient in 21/45 patients (46.6%). CSF cfDNA positivity was strongly associated with the presence of disseminated disease at the time of collection (81.5% of samples from patients with disseminated disease were positive). No association was seen between CSF cfDNA positivity and the timing of CSF collection during the patient's disease course. CONCLUSIONS: We identified three general categories where CSF cfDNA testing provided additional relevant diagnostic, prognostic, and/or therapeutic information, impacting clinical assessment and decision making: (1) diagnosis and/or identification of actionable alterations; (2) monitor response to therapy; and (3) tracking tumor evolution. Our findings support broader implementation of clinical CSF cfDNA testing in this population to improve care.

First use of grape waste-derived building blocks to yield antimicrobial materials.

Grape marc is an underutilised waste material that poses significant environmental issues. This study offers the first proof-of-concept investigation into the polymerisation of both crude and purified Sauvignon blanc grape marc extracts using the diacyl chlorides terephthaloyl chloride, succinyl chloride, adipoyl chloride, sebacoyl chloride, and the tartaric acid derivative (4R,5R)-2,2-dimethyl-1,3-dioxolane-4,5-dicarbonyl dichloride to obtain new materials, in what to the best of our knowledge is the first reported example of a direct polymerisation of an agricultural waste extract. A total of 26 novel materials were prepared. It has also shown that quercetin, a phenolic monomer found in grape marc extracts, can be polymerised with (4R,5R)-2,2-dimethyl-1,3-dioxolane-4,5-dicarbonyl dichloride to give a polymer that shows activity towards S. aureus.

RAS induced senescence of skin keratinocytes is mediated through Rho-associated protein kinase (ROCK).

Cellular senescence is a well-documented response to oncogene activation in many tissues. Multiple pathways are invoked to achieve senescence indicating its importance to counteract the transforming activities of oncogenic stimulation. We now report that the Rho-associated protein kinase (ROCK) signaling pathway is a critical regulator of oncogene-induced senescence in skin carcinogenesis. Transformation of mouse keratinocytes with oncogenic RAS upregulates ROCK activity and initiates a senescence response characterized by cell enlargement, growth inhibition, upregulation of senescence associated ß-galactosidase (SAßgal) expression, and release of multiple pro-inflammatory factors comprising the senescence-associated secretory phenotype (SASP). The addition of the ROCK inhibitor Y-27632 and others prevents these senescence responses and maintains proliferating confluent RAS transformed keratinocyte cultures indefinitely. Mechanistically, oncogenic RAS transformation is associated with upregulation of cell cycle inhibitors p15Ink4b , p16Ink4a , and p19Arf and downregulation of p-AKT, all of which are reversed by Y-27632. RNA-seq analysis of Y-27632 treated RAS-transformed keratinocytes indicated that the inhibitor reduced growth-inhibitory gene expression profiles and maintained expression of proliferative pathways. Y-27632 also reduced the expression of NF-κB effector genes and the expression of IκBζ downstream mediators. The senescence inhibition from Y-27632 was reversible, and upon its removal, senescence reoccurred in vitro with rapid upregulation of cell cycle inhibitors, SASP expression, and cell detachment. Y-27632 treated cultured RAS-keratinocytes formed tumors in the absence of the inhibitor when placed in skin orthografts suggesting that factors in the tumor microenvironment can overcome the drive to senescence imparted by overactive ROCK activity.

Epipyrone A, a Broad-Spectrum Antifungal Compound Produced by Epicoccum nigrum ICMP 19927.

We have isolated a filamentous fungus that actively secretes a pigmented exudate when growing on agar plates. The fungus was identified as being a strain of Epicoccum nigrum. The fungal exudate presented strong antifungal activity against both yeasts and filamentous fungi, and inhibited the germination of fungal spores. The chemical characterization of the exudate showed that the pigmented molecule presenting antifungal activity is the disalt of epipyrone A-a water-soluble polyene metabolite with a molecular mass of 612.29 and maximal UV-Vis absorbance at 428 nm. This antifungal compound showed excellent stability to different temperatures and neutral to alkaline pH.

Draft Genome Sequence of Piscirickettsia litoralis, Isolated from Seawater.

One species of Piscirickettsia, a pathogen of salmonid fish, has been described. The genome sequence of a putative second and free-living species may provide insights into the evolution of pathogenicity in the genus.