Neurocomputational mechanisms of affected beliefs The feedback people receive on their behavior shapes the process of belief formation and self-efficacy in mastering a given task. However, the neural and computational mechanisms of how the subjective value of these beliefs, and the corresponding affect, bias the learning process remain unclear. We investigated these mechanisms during the learning of self-efficacy beliefs using fMRI, pupillometry, and computational modeling, and by analyzing individual differences in affective experience. Biases in the formation of self-efficacy beliefs were associated with affect, pupil dilation, and neural activity within the anterior insula, amygdala, ventral tegmental area/ substantia nigra, and mPFC. Specifically, neural and pupil responses map the valence of the prediction errors in correspondence with individuals experienced affect and learning bias during belief formation. Together with the functional connectivity dynamics of the anterior insula within this network, our results hint at neural and computational mechanisms that integrate affect into the process of belief formation.
Spike-timing-dependent plasticity can account for connectivity aftereffects of dual-site transcranial alternating current stimulation Transcranial alternating current stimulation (tACS), applied to two brain sites with different phase lags, has been shown to modulate stimulation-outlasting functional EEG connectivity between the targeted regions. Given the lack of knowledge on mechanisms of tACS aftereffects, it is difficult to further enhance effect sizes and reduce variability in experiments. In this computational study, we tested if spike-timing-dependent plasticity (STDP) can explain stimulation-outlasting connectivity modulation by dual-site tACS and explored the effects of tACS parameter choices. Two populations of spiking neurons were coupled with synapses subject to STDP, and results were validated via a re-analysis of EEG data. Our simulations showed stimulation-outlasting connectivity changes between in- and anti-phase tACS, dependent on both tACS frequency and synaptic conduction delays. Importantly, both a simple network entraining to a wide range of tACS frequencies as well as a more realistic network that spontaneously oscillated at alpha frequency predicted that the largest effects would occur for short conduction delays between the stimulated regions. This finding agreed with experimental EEG connectivity modulation by 10 Hz tACS, showing a clear negative correlation of tACS effects with estimated conduction delays between regions. In conclusion, STDP can explain connectivity aftereffects of dual-site tACS. However, not all combinations of tACS frequency and application sites are expected to effectively modulate connectivity via STDP. We therefore suggest using appropriate computational models and/or EEG analysis for planning and interpretation of dual-site tACS studies relying on aftereffects.

HighlightsO_LINetwork model with STDP explains EEG connectivity change after dual-site tACS
C_LIO_LIEffects are predicted to depend on tACS frequency and conduction delays
C_LIO_LIEEG data confirm dependence on conduction delays between regions
C_LIO_LIModel can be used to estimate and maximize experimental effects
C_LI
Multivariate Genomewide Association Analysis by Iterative Hard Thresholding 1In genome-wide association studies (GWAS), analyzing multiple correlated traits is potentially superior to conducting multiple univariate analyses. Standard methods for multivariate GWAS operate marker-by-marker and are computationally intensive. We present a penalized regression algorithm for multivariate GWAS based on iterative hard thresholding (IHT) and implement it in a convenient Julia package MendelIHT.jl (https://github.com/OpenMendel/MendelIHT.jl). In simulation studies with up to 100 traits, IHT exhibits similar true positive rates, smaller false positive rates, and faster execution times than GEMMAs linear mixed models and mv-PLINKs canonical correlation analysis. On UK Biobank data, our IHT software completed a 3-trait joint analysis in 20 hours and an 18-trait joint analysis in 53 hours, requiring up to 80GB of computer memory. In short, our software enables geneticists to fit a single regression model that simultaneously considers the effect of all SNPs and dozens of traits. Naming of human diseases on the wrong side of history BackgroundIn the medical sphere, understanding naming conventions strengthen the integrity of naming human diseases remains nominal rather than substantial yet. Since the current nosology-based standard for human diseases could not offer a one-size-fits-all corrective mechanism, many idiomatic but flawed names frequently appear in scientific literature and news outlets at the cost of sociocultural impacts.

ObjectiveWe attempt to examine the ethical oversights of current naming practices and propose heuristic rationales and approaches to determine a pithy name instead of an inopportune nosology.

MethodsFirst, we examined the compiled global online news volumes and emotional tones on some inopportune nosology like German measles, Middle Eastern Respiratory Syndrome, Spanish flu, Hong Kong flu, and Huntingtons disease in the wake of COVID-19. Second, we prototypically scrutinize the lexical dynamics and pathological differentials of German measles and common synonyms by leveraging the capacity of the Google Books Ngram Corpus. Third, we demonstrated the empirical approaches to curate an exclusive substitute for an anachronistic nosology German measles based on deep learning models and post-hoc explanations.

ResultsThe infodemiological study shows that the public informed the offensive names with extremely negative tones in textual and visual narratives. The findings of the historiographical study indicate that many synonyms of German measles did not survive, while German measles became an anachronistic usage, and rubella has taken the dominant place since 1994. The PubMedBERT model could identify rubella as a potential substitution for German measles with the highest semantic similarity. The results of the semantic drift experiments further indicate that rubella tends to survive during the ebb and flow of semantic drift.

ConclusionsOur findings indicate that the nosological evolution of anachronistic names could result in sociocultural impacts without a corrective mechanism. To mitigate such impacts, we introduce some ethical principles for formulating an improved naming scheme. Based on deep learning models and post-hoc explanations, our illustrated experiments could provide hallmark references to the remedial mechanism of naming practices and pertinent credit allocations. Mathematical Modeling and Analysis of Mitochondrial Retrograde Signaling Dynamics Mitochondria, semi-autonomous eukaryotic organelles, participate in energy production and metabolism, making mitochondrial quality control crucial. As most mitochondrial proteins are encoded by nuclear genes, quality control depends on proper mitochondria-nucleus communication, designated mitochondrial retrograde signaling. Early studies focused on retrograde signaling participants and specific gene knockouts. However, mitochondrial signal modulation remains elusive. Using yeast, we simulated signal propagation following mitochondrial damage and proposed a mathematical model based ordinary differential equations. Mitochondrial retrograde signaling decisions were described using a Boolean model. Dynamics were analyzed through an ordinary differential equation-based model and extended to evaluate the model response to noisy damage signals. Simulation revealed localized protein concentration dynamics, including waveforms, frequency response, and robustness under noise. Retrograde signaling is bistable with three localized steady states, and increased damage compromises robustness. We elucidated mitochondrial retrograde signaling, thus providing a basis for drug design against yeast and fungi.

Author SummaryThe yeast RTG pathway regulates mitochondrial metabolism and mitochondrial quality through passing mitochondrial signal to the nucleus to modulate gene expressions. Using microscopic data of RTG proteins translocation from cytosol to nucleus, the parameters were found by fitting 16 knockout conditions simulated from the proposed differential equation-based model. Further, through dose response, ultrasensitivity, frequency response and noise interference, we demonstrated a switch-like property of RTG pathway activation, a capacity charging property of RTG protein translocation, a low pass filter property of signals in different frequencies, and noise amplification due to mitochondrial signal. The proposed mathematical model enables us to understand the dynamics and mechanisms in mitochondrial retrograde signaling in yeast and provides potential antifungal treatment strategies. Cell type-specific inference of differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr. Thermodynamic constraints on the assembly and diversity of microbial ecosystems are different near to and far from equilibrium Non-equilibrium thermodynamics has long been an area of substantial interest to ecologists because most fundamental biological processes, such as protein synthesis and respiration, are inherently energy-consuming. However, most of this interest has focused on developing coarse ecosystem-level maximisation principles, providing little insight into underlying mechanisms that lead to such emergent constraints. Microbial communities are a natural system to decipher this mechanistic basis because their interactions in the form of substrate consumption, metabolite production, and cross-feeding can be described explicitly in thermodynamic terms. Previous work has considered how thermodynamic constraints impact competition between pairs of species, but restrained from analysing how this manifests in complex dynamical systems. To address this gap, we develop a thermodynamic microbial community model with fully reversible reaction kinetics, which allows direct consideration of free-energy dissipation. This also allows species to interact via products rather than just substrates, increasing the dynamical complexity, and allowing a more nuanced classification of interaction types to emerge. Using this model, we find that community diversity increases with substrate lability, because greater free-energy availability allows for faster generation of niches. Thus, more niches are generated in the time frame of community establishment, leading to higher final species diversity. We also find that allowing species to make use of near-to-equilibrium reactions increases diversity in a low free-energy regime. In such a regime, two new thermodynamic interaction types that we identify here reach comparable strengths to the conventional (competition and facilitation) types, emphasising the key role that thermodynamics plays in community dynamics. Our results suggest that accounting for realistic thermodynamic constraints is vital for understanding the dynamics of real-world microbial communities.

Author summaryThere is a growing interest in microbial communities due to their important role in biogeochemical cycling as well as plant and animal health. Although our understanding of thermodynamic constraints on individual cells is rapidly improving, the impact of these constraints on complex microbial communities remains largely unexplored theoretically and empirically. Here, we develop a new microbial community model which allows thermodynamic efficiency and entropy production to be calculated directly. We find that availability of substrates with greater free-energy allows for a faster rate of niche generation, leading to higher final species diversity. We also show that when the free-energy availability is low, species with reactions close to thermodynamic equilibrium are favoured, leading to more diverse and efficient communities. In addition to the conventional interaction types (competition and facilitation), our model reveals the existence of two novel interaction types mediated by products rather than substrates. Though the conventional interactions are generally the strongest, the novel interaction types are significant when free-energy availability is low. Our results suggest that non-equilibrium thermodynamics need to be considered when studying microbial community dynamics. In vitro evolution of uracil glycosylase towards DnaKJ and GroEL binding evolves different misfolded states Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK+DnaJ+ATP (DnaKJ) or E.coli cell extract (CE). Using binding to the eUNG inhibitor Ugi as probe for native foldedess, the CE selected population was further divided to Ugi binders (+U) (native) or non-binders (-U). We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection, as well as CE+U and natural evolution of eUNG has its own trajectory. While GroEL and CE-U selected for mutants highly sensitive to protease cleavage, DnaKJ selected for partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. CE+U selected for more neutral mutations than natural evolution. In a more general context, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds mutants misfolding-prone species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL and DnaKJ selected populations, while harboring also unique properties, explained by the existence of additional chaperones in CE, such as Tig, HtpG and ClpB. Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. An amino residue that guides the correct photoassembly the water-oxidation complex but is not required for high affinity Mn2+ binding The assembly of the Mn4O5Ca cluster of the photosystem II (PSII) starts from the initial binding and photooxidation of the first Mn2+ at a high affinity site (HAS). Recent cryo-EM apo-PSII structures reveal an altered geometry of amino ligands in this region and suggest the involvement of D1-Glu189 ligand in the formation of the HAS. We now find that Gln and Lys substitution mutants photoactivate with reduced quantum efficiency compared to the wild-type. However, the affinity of Mn2+ at the HAS in D1-E189K was very similar to the wild-type (~2.2 M). Thus, we conclude that D1-E189 does not form the HAS (~2.9 M) and that the reduced quantum efficiency of photoactivation in D1-E189K cannot be ascribed to the initial photooxidation of Mn2+ at the HAS. Besides reduced quantum efficiency, the D1-E189K mutant exhibits a large fraction of centers that fail to recover activity during photoactivation starting early in the assembly phase, becoming recalcitrant to further assembly. Fluorescence relaxation kinetics indicate on the presence of an alternative route for the charge recombination in Mn-depleted samples in all studied mutants and exclude damage to the photochemical reaction center as the cause for the recalcitrant centers failing to assemble and show that dark incubation of cells reverses some of the inactivation. This reversibility would explain the ability of these mutants to accumulate a significant fraction of active PSII during extended periods of cell growth. The failed recovery in the fraction of inactive centers appears to a reversible mis-assembly involving the accumulation of photooxidized, but non-catalytic high valence Mn at the donor side of photosystem II, and that a reductive mechanism exists for restoration of assembly capacity at sites incurring mis-assembly. Given the established role of Ca2+ in preventing misassembled Mn, we conclude that D1-E189K mutant impairs the ligation of Ca2+ at its effector site in all PSII centers that consequently leads to the mis-assembly resulting in accumulation of non-catalytic Mn at the donor side of PSII. Our data indicate that D1-E189 is not functionally involved in Mn2+ oxidation\binding at the HAS but rather involved in Ca2+ ligation and steps following the initial Mn2+ photooxidation. Behavioral effects of ethanol in the Red Swamp Crayfish (Procambarus clarkii) Alcohol abuse and dependence remains one of the primary preventable sources of human mortality in the United States. Model laboratory species can be used to evaluate behavioral, neuropharmacological and other biological changes associated with alcohol exposure and to identify novel treatment modalities. This report describes methods for evaluating the behavioral effects of ethyl alcohol (ethanol; EtOH) in a crayfish model. Crayfish (Procambarus clarkii) were immersed in ethanol baths with concentrations ranging from 0.1 to 1.0 molar, and for durations of 10-30 minutes. Studies evaluated hemolymph alcohol concentration, locomotor behavior in an open field and anxiety-like behavior using a Light/Dark transfer approach. EtOH immersion produced dose dependent increases in hemolymph EtOH concentration (up to 249 mg/dL) and reductions in open field locomotor behavior that depended on EtOH concentration or exposure duration. Under baseline conditions, crayfish exhibit avoidance of the open parts of the locomotor arena and a preference for the covered portion, when available. Acute EtOH immersion decreased time spent in the covered portion of the Light/Dark arena, consistent with a decrease in anxiety-like behavior. EtOH immersion once per day for 5 days did not alter the acute locomotor behavioral response, however increased activity was observed 3 days after the repeated EtOH regimen. Overall, this study shows that this inexpensive, easily maintained species can be used for behavioral pharmacological experiments designed to assess the acute and repeated effects of EtOH. Age of the bone marrow dictates clonality of smooth muscle-derived cells in the atherosclerotic plaque Aging is the predominant risk factor for atherosclerosis, the leading cause of death. Rare smooth muscle cell (SMC) progenitors clonally expand giving rise to up to [~]70% of atherosclerotic plaque cells; however, the effect of age on SMC clonality is not known. Our results indicate that aged bone marrow (BM)-derived cells non-cell autonomously induce SMC polyclonality and worsen atherosclerosis. Indeed, in myeloid cells from aged mice and humans, TET2 levels are reduced which epigenetically silences integrin {beta}3 resulting in increased tumor necrosis factor [TNF]- signaling. In turn, TNF signals through TNF receptor 1 on SMCs to promote proliferation and induces recruitment and expansion of multiple SMC progenitors into the atherosclerotic plaque. Notably, integrin {beta}3 overexpression in aged BM preserves dominance of the lineage of a single SMC progenitor and attenuates plaque burden. Our results demonstrate a molecular mechanism of aged macrophage-induced SMC polyclonality and atherogenesis and suggest novel therapeutic strategies.

Graphical abstractAge of BM-derived monocytes/macrophages determines clonality of SMC lineage in the atherosclerotic plaque. Atherogenesis is depicted in a young (a) or aged (b) host. Aged monocytes/macrophages have decreased levels of the epigenetic regulator TET2, leading to reduction of the 5-hydroxymethylcytosine (5hmC) mark on the Itgb3 promoter. The resulting low integrin {beta}3 levels in aged monocytes/macrophages induces high TNF- levels, facilitating recruitment and expansion of multiple SMC progenitors (polyclonality) in the atherosclerotic plaque and worse disease burden. In contrast, the young control is characterized by mono/oligoclonal SMC expansion in a smaller plaque.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/476756v2_ufig1.gif" ALT="Figure 1">
View larger version (28K):
org.highwire.dtl.DTLVardef@12f7090org.highwire.dtl.DTLVardef@186ff92org.highwire.dtl.DTLVardef@1f4da69org.highwire.dtl.DTLVardef@2f509c_HPS_FORMAT_FIGEXP  M_FIG C_FIG Tracheal tube fusion in Drosophila involves release of exosomes from multivesicular bodies Extracellular vesicles (EVs) comprise diverse types of cell-released membranous structures that are thought to play important roles in intercellular communication. Despite extensive work on the formation and functions of EVs in cultured cells, studies of EVs in vivo have remained scarce. We report here that EVs are present in the developing lumen of tracheal tubes in Drosophila embryos. We defined two distinct EV subpopulations, one of which contains the Munc13-4 homologue Staccato (Stac) and is spatially and temporally associated with tracheal tube fusion events. The formation of Stac-positive luminal EVs depends on the tip-cell-specific GTPase Arl3, which is also required for the formation of Stac-positive multivesicular bodies, suggesting that Stac-EVs derive from fusion of Stac-MVBs with the luminal membrane in tip cells during anastomosis formation. The GTPases Rab27 and Rab35 cooperate downstream of Arl3 to promote Stac-MVB formation and tube fusion. We propose that Stac-MVBs act as membrane reservoirs that facilitate tracheal lumen fusion in a process regulated by Arl3, Rab27, Rab35, and Stac/Munc13-4. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Analysis of Drosophila cardiac hypertrophy by micro-computerized tomography for genetic dissection of heart growth mechanisms Heart failure is often preceded by pathological cardiac hypertrophy, a thickening of the heart musculature driven by complex gene regulatory and signaling processes. The Drosophila heart has great potential as a genetic model for deciphering the underlying mechanisms of cardiac hypertrophy. However, current methods for evaluating hypertrophy of the Drosophila heart are laborious and difficult to carry out reproducibly. Here we demonstrate that micro-computerized tomography (microCT) is an accessible, highly reproducible method for non-destructive, quantitative analysis of Drosophila heart morphology and size. To validate our microCT approach for analyzing Drosophila cardiac hypertrophy, we show that expression of constitutively active Ras (Ras85DV12), previously shown to cause hypertrophy of the fly heart, results in significant thickening of both adult and larval heart walls when measured from microCT images. We then show using microCT analysis that genetic upregulation of store-operated Ca2+ entry (SOCE) driven by expression of constitutively active Stim (StimCA) or Orai (OraiCA) proteins also results in significant hypertrophy of the Drosophila heart, through a process that specifically depends on Orai Ca2+ influx channels. Intravital imaging of heart contractility revealed significantly reduced end diastolic dimensions in StimCA and OraiCA expressing hearts, consistent with the hypertrophic phenotype. These results demonstrate that increased SOCE activity is an important driver of hypertrophic cardiomyocyte growth, and demonstrate how microCT analysis combined with tractable genetic tools in Drosophila can be used to delineate molecular signaling processes that underlie cardiac hypertrophy and heart failure.

NEW AND NOTEWORTHYGenetic analysis of cardiac hypertrophy in Drosophila holds immense potential for the discovery of new therapeutic targets to prevent and treat heart failure. However, this potential has been hindered by a lack of rapid and effective methods for analysis of heart size in flies. Here we demonstrate that analysis of the Drosophila heart with micro-computerized tomography yields accurate and highly reproducible heart size measurements that can be used to efficiently analyze heart growth and cardiac hypertrophy in Drosophila. Feline tooth resorption findings from a questionnaire survey on 8115 Finnish cats and a case-control study for a more defined population drawn from the survey Tooth resorption (TR) is one of the most common dental diseases in cats. It is a painful disease characterised by progressive dental destruction, which eventually results in loss of teeth. The aetiology of the disease remains unclear, but associations with old age, breed, other oral and dental diseases, and certain environmental factors have been suspected. We collected health and environmental information of 8115 Finnish cats from 41 breeds through an online survey for cat owners to study TR in Finland with an observational study. The frequency of TR that the owner reported as veterinarian-diagnosed was 3.9% in the entire data (316/8115) and 15% in cats that the owner reported as having veterinarian-diagnosed oral or dental disease (316/2070). We determined the risk factors and breed variation of feline TR by logistic regression with the entire data and a case-control study in a more defined population drawn from the entire data. In the entire data, we defined those cats that the owners had reported as being diagnosed with TR by a veterinarian as TR-cats and compared them to all other cats. In the case-control study, we included only cats with owner-reported veterinarian-diagnosed oral and dental diseases. We defined cats with TR-diagnosis made under sedation as TR-cases. Other veterinarian-diagnosed TR-cats were excluded. Cats with other oral and dental diseases were used as controls. The model with the entire population included 308 TR-cats compared with 7508 cats, and the case-control study included 198 TR-cats and 1683 controls. Although the diagnoses were reportedly made by a veterinarian, the owner answered the survey and might have understood, remembered, or interpreted the diagnoses in a way that could affect the results. The results indicated an increasing risk of TR with age. Dental calculus, gingivitis, and periodontitis were potential risk factors for TR. These findings and the interaction between dental calculus and gingivitis suggest that inflammatory changes caused by dental calculus increase the risk of TR. Keeping food available constantly was a potential protective factor. We found that Cornish Rex, European, and Ragdoll breeds were at higher risk for TR. The observed differences between breeds highlight a genetic contribution. DeepVelo: Data-driven Single-cell Transcriptomic Deep Velocity Field Learning with Neural Ordinary Differential Equations Recent advances in single-cell RNA sequencing technology have provided unprecedented opportunities to simultaneously measure the gene expression profile and transcriptional velocity of individual cells, enabling us to sample gene regulatory network dynamics along developmental trajectories. However, traditional methods have faced challenges in modeling gene expression dynamics within individual cells due to sparse, non-linear (e.g., obligate heterodimer transcription factors), and high-dimensional measurements. Here, we present DeepVelo, a neural-network-based ordinary differential equation model that can learn non-linear, high-dimensional single-cell transcriptome dynamics and describe continuous gene expression changes within individual cells across time. We applied DeepVelo to multiple published datasets from different technical platforms and demonstrated its utility to 1) formulate transcriptome dynamics on different timescales, 2) measure the instability of cell states, and 3) identify developmental driver genes upstream of a signaling cascade. Benchmarking against state-of-the-art methods shows that DeepVelo can improve velocity field representation accuracy by at least 50% in out-of-sample cells. Further, perturbation studies revealed that single-cell dynamical systems may exhibit properties similar to those of chaotic systems. In summary, DeepVelo allows for the data-driven discovery of differential equations that delineate single-cell transcriptome dynamics.

TeaserEmbedding neural networks into ordinary differential equations to model gene expression changes within single cells across time. ART duration and immunometabolic state determine efficacy of DC-based treatment restoring functional HIV- specific CD8+ T cells in PLWH. Dysfunction of CD8+ T cells in people living with HIV-1 (PLWH) receiving anti-retroviral therapy (ART) has restricted the efficacy of dendritic cell (DC)-based immunotherapies against HIV-1. Heterogeneous immune exhaustion and metabolic states of CD8+ T cells might differentially associate with dysfunction. However, specific parameters associated to functional restoration of CD8+ T cells after DC treatment have not been investigated in detail. Here, we studied the association of ART duration with memory subsets, exhaustion and metabolic profiles of CD8+ T cells from PLWH and improvement of polyfunctional and effector HIV-1 specific responses after stimulation with Gag-adjuvant-primed DC. HIV-1-specific CD8+ T cell responses from a larger proportion PLWH on ART for more than 10 years (LT-ARTp) improved polyfunctionality and capacity to eliminate autologous p24+ infected CD4+ T cells in vitro. In contrast, CD8+ T cells from PLWH on ART for less than a decade (ST-ARTp) were less responsive to DC treatment and functional improvement was limited in this group. This was associated with lower frequencies of central memory CD8+ T cells, increased co-expression of PD1 and TIGIT and reduced mitochondrial respiration and glycolytic induction upon TCR activation. In contrast, CD8+ T cells from LT-ARTp showed increased frequencies of TIM3+PD1-cells and preserved induction of glycolysis. Treatment of dysfunctional CD8+ T cells from ST-ARTp with combined anti-PD1 and anti-TIGIT antibodies plus a glycolysis promoting drug restored their ability to eliminate infected CD4+ T cells. Together, our study identifies specific immunometabolic parameters for different PLWH subgroups potentially useful for future personalized DC-based HIV-1 vaccines. Extracellular Vesicle (EV) Mitochondria and EV/HSP27 Mixtures Protect the Ischemic Brain Endothelium Ischemic stroke causes brain endothelial cell (BEC) death and damages tight junction integrity of the blood-brain barrier (BBB). We harnessed the innate mitochondrial load of endothelial cell-derived extracellular vesicles (EVs) and utilized mixtures of EV/exogenous heat shock protein 27 (HSP27) as a one-two punch strategy to increase BEC survival (via EV mitochondria) and preserve their tight junction integrity (via HSP27 effects). We demonstrated that the medium-to-large (m/lEV) but not small EVs (sEV) transferred their mitochondrial load, which subsequently colocalized with the mitochondrial network of the recipient primary human BECs. BECs treated with m/lEVs increased relative ATP levels and displayed superior mitochondrial function. m/lEVs isolated from rotenone (mitochondrial complex I inhibitor) exposed BECs (RTN-m/lEVs) did not increase BECs ATP levels compared to naive m/lEVs. In contrast, RTN-sEVs functionality was minimally affected compared to naive sEVs. Compared to EV/HSP27 mixtures, PEG-DET/EV/HSP27 ternary mixtures increased HSP27 protein loading. (PEG-DET/HSP27)/EV ternary and EV/HSP27 binary mixtures decreased the paracellular permeability of small and large molecular mass fluorescent tracers in oxygen glucose-deprived primary human BECs. This one-two punch approach to increase BEC metabolic function and tight junction integrity is a promising strategy for BBB protection and prevention of long-term neurological dysfunction post-ischemic stroke.

Highlights[tpltrtarr] Medium-to-large extracellular vesicles (m/lEVs), not small EVs contain mitochondria
[tpltrtarr]m/lEVs increased ATP and mitochondrial function in brain endothelial cells (BECs)
[tpltrtarr]Donor m/lEV mitochondria colocalized with host mitochondria in recipient BECs
[tpltrtarr]m/lEVs from rotenone-exposed BECs did not increase recipient BEC ATP levels
[tpltrtarr]PEG-DET/EV/HSP27 ternary mixtures increase exogenous HSP27 protein loading


Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/466491v4_ufig1.gif" ALT="Figure 1">
View larger version (43K):
org.highwire.dtl.DTLVardef@155f618org.highwire.dtl.DTLVardef@8fa778org.highwire.dtl.DTLVardef@13bd44aorg.highwire.dtl.DTLVardef@ee771e_HPS_FORMAT_FIGEXP  M_FIG C_FIG CRCM5484: A BET- BDII Selective Compound With Differential Anti-Leukemic Drug Modulation Differentially screening the Fr-PPIChem chemical library on the BET BRD4-BDII versus -BDI bromodomains led to the discovery of a BDII selective tetrahydropyridothienopyrimidinone (THPTP)-based compound. Structure-activity relationship (SAR) and hit-to-lead approaches allowed us to develop CRCM5484, a potent inhibitor of BET proteins with a preferential and 475-fold selectivity for the second bromodomain of the BRD3 protein (BRD3-BDII) over its first bromodomain (BRD3-BDI). Its very low activity was demonstrated in various cell-based assays, corresponding with recent data describing other selective BDII compounds. However, screening on a drug sensitivity and resistance-profiling platform revealed its ability to modulate the anti-leukemic activity in combination with various FDA-approved and/or in-development drugs in a cell- and context-dependent differential manner. Altogether, the results confirm the originality of the THPTP molecular mode of action in the BD cavity and its potential as starting scaffold for the development of potent and selective bromodomain inhibitors. Morphology of immatures of the thelytokous ant, Monomorium triviale Wheeler (Formicidae: Myrmicinae: Solenopsidini) with descriptions of the extraordinary last-instar queen larvae The ant genus Monomorium is one of the most species-rich but taxonomically problematic groups in the hyperdiverse subfamily Myrmicinae. An East Asian species, M. triviale Wheeler, produces both reproductive queens and sterile workers via obligate thelytokous parthenogenesis. Here, we describe the immature forms of M. triviale based on light and scanning electron microscopy observations, with a note on the striking caste dimorphism in the last larval instar. The last-instar queen larvae were easily recognized by their large size, "aphaenogastroid" body shape, and rows of doorknob-like tubercles on the lateral and dorsal body surface. This type of queen-specific structure has not been found in ants in general, let alone congeneric species found in Japan. In stark contrast to the queen larvae, worker larvae showed a "pheidoloid" body shape and a body surface similar to other ants. The worker larvae were estimated to have three instars, consistent with previously described congeners. The pupae of both castes had no cocoon, a characteristic commonly described in other Myrmicinae species. In total, the developmental period from egg to adult worker averaged 59 days under 25{degrees}C. We discuss possible functions of the tubercles of queen larvae based on previous studies. A cell type-specific multiomics uncovers a guard cell-specific RAF15-SnRK2.6/OST1 kinase cascade Multicellular organisms such as plants contain different cell types with specialized functions. Analyzing the characteristics of each cell type reveals specific cell functions and enhances understanding of organization and function at the organismal level. Here we report a highly-sensitive and efficient cell type-specific multiomics pipeline, combining simplified flow cytometry-based cell sorting for fluorescent protoplasts and optimized nanoscale proteomics and metabolomics methods, which allow in-depth analysis of the proteomes and metabolomes of a particular cell type. Using this method, we quantitatively compared the proteomes and metabolomes of guard cells and mesophyll cells and revealed that the enrichment of signal transduction-related proteins enables guard cells to respond rapidly to various environmental stimuli. We uncovered a guard cell-specific kinase cascade whereby RAF15 and OST1 mediate ABA-induced stomatal closure. This pipeline can be applied to various cell types in plant or non-plant systems to learn how cells function in highly organized multicellular organisms. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Long term analysis of social structure: evidence of age-based consistent associations in male Alpine ibex Despite its recognized importance for understanding the evolution of animal sociality as well as for conservation, long term analysis of social networks of animal populations is still relatively uncommon. We investigated social network dynamics in males of a gregarious mountain ungulate (Alpine ibex, Capra ibex) over ten years focusing on groups, sub-groups and individuals, exploring the dynamics of sociality over different scales. Despite the social structure changing between seasons, the Alpine ibex population was highly cohesive: fission-fusion dynamics lead almost every male in the population to associate with each other male at least once. Nevertheless, we found that male Alpine ibex showed preferential associations that were maintained across seasons and years. Age seemed to be the most important factor driving preferential associations while other characteristics, such as social status, appeared less crucial. We also found that centrality measures were influenced by age and were also related to individual physical condition. The multi-scale and long-term frame of our study helped us show that ecological constrains, such as resource availability, may play a role in shaping associations in a gregarious species, but they cannot solely explain sociality and preferential association that are likely also to be driven by life-history linked physiological and social needs. Our results highlight the importance of long-term studies based on individually recognizable subjects to help us build on our understanding of the evolution of animal sociality. On-person adaptive evolution of Staphylococcus aureus during treatment for atopic dermatitis Genetic variation among bacterial strains can contribute to heterogeneity in the severity of chronic inflammatory diseases 1,2, but the degree of variation created by de novo mutation during colonization is not well understood. The inflamed skin of people with atopic dermatitis (AD) is frequently colonized with Staphylococcus aureus, an opportunistic pathogen associated with both asymptomatic colonization of nasal passages and invasive disease3-6. While genetic risk and barrier disruption are critical to AD initiation 7,8, S. aureus colonization worsens disease severity by promoting skin damage9. Here we show, from tracking 25 children treated for AD over 9 months, that S. aureus adapts via de novo mutations during colonization. Patients S. aureus populations are typically dominated by a single lineage, with infrequent invasion by distant genotypes. Variants emerge within each lineage with mutation rates similar to S. aureus in other contexts. Some of these variants replace their ancestors across the body within months, with signatures of adaptive, rather than neutral, forces. Most strikingly, the capsule synthesis gene capD obtained four parallel mutations within one patient and was involved in mutational sweeps in multiple patients. We confirm that selection for capD negativity is common in AD, but not in other contexts, via reanalysis of public S. aureus genomes from 276 people. Our finding of disease-specific selection raises the possibility that adaptation of pathobionts during colonization prolongs the positive feedback cycle of inflammation. The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP3 produced by PI3Ks. Excess PIP3 promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP2 and PIP3. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP3. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP3 allosterically unfolding the N-terminal PIP2 binding domain, positioning it favorably for the polybasic motif interaction with PIP2. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.