Coordinated multiplexing of information about separate objects in visual cortex Sensory receptive fields are large enough that they can contain more than one perceptible stimulus. How, then, can the brain encode information about each of the stimuli that may be present at a given moment? We recently showed that when more than one stimulus is present, single neurons can fluctuate between coding one vs. the other(s) across some time period, suggesting a form of neural multiplexing of different stimuli (Caruso et al., 2018). Here we investigate (a) whether such coding fluctuations occur in early visual cortical areas; (b) how coding fluctuations are coordinated across the neural population; and (c) how coordinated coding fluctuations depend on the parsing of stimuli into separate vs. fused objects. We found coding fluctuations do occur in V1 but only when the two stimuli form separate objects. Such separate objects evoked a novel pattern of V1 spike count ("noise") correlations involving distinct distributions of positive and negative values. This bimodal correlation pattern was most pronounced among pairs of neurons showing the strongest evidence for coding fluctuations or multiplexing. Whether a given pair of neurons exhibited positive or negative correlations depended on whether the two neurons both responded better to the same object or had different object preferences. Distinct distributions of spike count correlations based on stimulus preferences were also seen in V4 for separate objects but not when two stimuli fused to form one object. These findings suggest multiple objects evoke different response dynamics than those evoked by single stimuli, lending support to the multiplexing hypothesis and suggesting a means by which information about multiple objects can be preserved despite the apparent coarseness of sensory coding.

Significance StatementHow the brain separates information about multiple objects despite overlap in the neurons responsive to each item is not well understood. Here we show that some neurons in V1 exhibit coding fluctuations in response to two objects, and that these coding fluctuations are coordinated at the population level in ways that are not observed for single objects. Broadly similar results were obtained in V4. These response dynamics lend support to the hypothesis that information about individual objects may be multiplexed across the neural population, preserving information about each item despite the coarseness of sensory coding.
The Temporal Pattern of Synaptic Activation Determines the Longevity of Structural Plasticity at Single Dendritic Spines Individual synapses are the points at which information is passed between neurons, yet it is unknown how the diverse patterns of activity that are observed in vivo effect plasticity at the level of single inputs. Here, we aimed to determine what are the structural plasticity consequences of naturalistic patterns of activity at single spines, as these reflect changes in synaptic efficacy. Utilizing two- photon fluorescence imaging and glutamate uncaging, we studied structural plasticity of individual CA1 hippocampal dendritic spines using activation patterns sampled from a Poisson distribution, which resemble endogenous firing patterns from their CA3 inputs. We found that while the majority of inputs initially undergo structural changes, the robustness of this plasticity is determined by the timing structure of the Poisson sampled naturalistic stimulation patterns. Further, we found that structural plasticity elicited by these naturalistic patterns is both NMDAR and protein synthesis dependent, consistent with requirements for other forms of plasticity. Lastly, we found that during the delivery of naturalistic activity patterns, spines underwent rapid and dynamic structural growth that predicted the longevity of plasticity, which was not the case during non-naturalistic stimulation protocols. These data suggest that dendritic spines are able to integrate incoming temporal information and accordingly modulate the longevity of plasticity that is induced.

HighlightsO_LINaturalistic stimulation of single dendritic spines of CA1 hippocampal neurons induces long lasting structural plasticity that depends on the temporal distribution of the synaptic events.
C_LIO_LIStructural plasticity induced by naturalistic stimulation patterns requires NMDA receptor activation and new protein-synthesis.
C_LIO_LIRapid spine structural dynamics during naturalistic activity, but not regular patterns, predict the longevity of subsequent structural plasticity.
C_LI
A distinct brain expression profile for genes within large introgression deserts and under positive selection in Homo sapiens Analyses of ancient DNA from extinct hominins have provided unique insights into the complex evolutionary history of Homo sapiens, intricately related to that of the Neanderthals and the Denisovans as revealed by several instances of admixture events. These analyses have also allowed the identification of introgression deserts: genomic regions in our species that are depleted of  archaic haplotypes. The presence of genes like FOXP2 in these deserts has been taken to be suggestive of brain-related functional differences between Homo species. Here, we seek a deeper characterization of these regions, taking into account signals of positive selection in our lineage. Analyzing publicly available transcriptomic data from the human brain at different developmental stages, we found that structures outside the cerebral neocortex, in particular the cerebellum, the striatum and the mediodorsal nucleus of the thalamus show the most divergent transcriptomic profiles when considering genes within large introgression deserts and under positive selection. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. Robust retrieval of data stored in DNA by de Bruijn graph-based de novo strand assembly DNA data storage is a rapidly developing technology with great potential due to its high density, long-term durability, and low maintenance cost. The major technical challenges include various errors, such as the strand breaks, rearrangements, and indels that frequently arise during DNA synthesis, amplification, sequencing, and preservation. Through a de novo assembly strategy, we developed an algorithm based on the de Bruijn graph and greedy path search (DBGPS) to address these issues. DBGPS shows distinct advantages in handling DNA breaks, rearrangements, and indels. The robustness of DBGPS is demonstrated by accelerated aging, multiple independent data retrievals, deep error-prone PCR, and large data scale simulations. Remarkably, 6.8 MB of data can be retrieved accurately from a seriously corrupted sample that has been treated at 70 {degrees}C for 70 days. With DBGPS, we were able to achieve a logical density of 1.30 bits/cycle and a physical density of 295 PB/g.

One-Sentence SummaryA de Bruijn graph-based de novo assembly algorithm for DNA data storage enables fast and robust data readouts even with DNA samples that have been severely corrupted. GLOBAL LAND-USE AND LAND-COVER DATA: HISTORICAL, CURRENT AND FUTURE SCENARIOS Land-use land-cover (LULC) data are important predictors of species occurrence and biodiversity threat. Although there are LULC datasets available under current conditions, there is a lack of such data under historical and future climatic conditions. This hinders, for example, projecting niche and distribution models under global change scenarios at different time scenarios. The Land Use Harmonization Project (LUH2) is a global terrestrial dataset at 0.25{degrees} spatial resolution that provides LULC data from 850 to 2300 for 12 LULC state classes. The dataset, however, is compressed in a file format (NetCDF) that is incompatible for many analyses and requires layer extractions and transformations that are intractable for most researchers. Here we selected and transformed the LUH2 in a standard GIS format data to make it more user-friendly. We provide LULC for every year from 850 to 2100, and from 2015 on, the LULC dataset is provided under two Shared Socioeconomic Pathways (SSP2-4.5 and SSP5-8.5). We provide two types of files for each year: separate files with continuous values for each of the 12 LULC state classes, and a single categorical file with all state classes combined. To create the categorical layer, we assigned the state with the highest value in a given pixel among the 12 continuous data. LUH2 predicts a pronounced decrease in primary forest, particularly noticeable in the Amazon, the Brazilian Atlantic Forest, the Congo Basin and the boreal forests, an equally pronounced increase in secondary forest and non-forest lands, and in croplands in the Brazilian Atlantic Forest and sub-Saharan Africa. The final dataset provides LULC data for 1251 years that will be of interest for macroecology, ecological niche modeling, global change analysis, and other applications in ecology and conservation. An integrative analysis to predict the Active Ingredients and explore Polypharmacological Mechanisms of Orthosiphon stamineus Benth. for human health BackgroundOrthosiphon stamineus Benth. was a dietary supplement or a traditional folk herb with widespread clinical application, but it was lack of comprehensive understanding of its active ingredients and polypharmacological mechanisms. The aim of this work was to systematically study its natural compounds and the molecular mechanisms of the O. stamineus via network pharmacology.

MethodsCompounds from O. stamineus were collected by literature retrieval and evaluated by SwissADME with the physicochemical property (ADMET model) and the likelihood for a natural medicine. The connection of active ingredients and target genes was built and confirmed by Cytoscape and AutoDock vina. Then, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were separately performed to obtain a more in-depth understanding of O. stamineus. Finally, the relationship among active ingredients, targets, and diseases was built to clarify the polypharmacological mechanisms and found relevant active substances for further drug discovery.

ResultsA total of 159 compounds of O. stamineus were collected, and 22 potentially active ingredients were screened out. The Ingredient-Target Interaction Network was built with 12 flavonoids, 3 diterpenes, 3 phenols, and 4 volatile oils, and 65 targets. The Docking analysis indicated that the ingredient-target interaction network was reliable; most ligand-receptor had a strong binding affinity (lowest binding energy: -6.9 kcal/mol). After pathway analysis, 185 significant biological processes and 36 signal pathways were found, and the ingredient-target-disease network of O. stamineus was constructed for polypharmacological mechanisms.

ConclusionOur study clarified the polypharmacological mechanisms via the relationship among active ingredients, targets, and diseases and provided better guidance for subsequent experiments and potential active ingredients for drug discovery or health promotion. Molecular diversity maintained by long-term balancing selection in mating loci defines multiple mating types in fungi Balancing selection, an evolutionary force that retains genetic diversity, has been detected in multiple genes and organisms, such as the sexual mating loci in fungi. However, to quantify the strength of balancing selection and define the mating-related genes require a large number of specimens. In tetrapolar basidiomycete fungi, sexual type is determined by two unlinked loci, MATA and MATB. Genes in both loci defines mating type identity, control successful mating and completion of the life cycle. These loci are usually highly diverse. Previous studies have speculated, based on culture crosses, that species of the non-model genus Trichaptum (Hymenochaetales, Basidiomycota) possess a tetrapolar mating system, with multiple alleles. Here, we sequenced a hundred and eighty specimens of three Trichaptum species. We characterized the chromosomal location of MATA and MATB, the molecular structure of MAT regions and their allelic richness. Our sequencing effort was sufficient to molecularly characterize multiple MAT alleles segregating before the speciation event of Trichaptum species. Our analyses suggested that long-term balancing selection has generated trans-species polymorphisms. Mating sequences were classified in different allelic classes based on an amino acid identity (AAI) threshold supported by phylogenetics. The inferred allelic information mirrored the outcome of in vitro crosses, thus allowing us to support the degree of allelic divergence needed for successful mating. Even with the high amount of divergence, key amino acids in functional domains are conserved. The observed allelic classes could potentially generate 14,560 different mating types. We conclude that the genetic diversity of mating in Trichaptum loci is due to long-term balancing selection, with limited recombination and duplication activity. Our large number of sequenced specimens highlighted the importance of sequencing multiple individuals from different species to detect the mating-related genes, the mechanisms generating diversity and the evolutionary forces maintaining them.

Author summaryFungi have complex mating systems, and basidiomycete fungi can encode thousands of mating types. Individuals with the same mating type cannot mate. This sexual system has evolved to facilitate sexual mating, increasing the chances to recombine into advantageous allelic combination and prune deleterious alleles. We explored the genomes of hundred and eighty specimens, combined with experimental mating studies of selected specimens, from a non-model organism (Trichaptum). We characterized the genomic regions controlling sex. The mating ability of the specimens confirmed the role of the mating alleles observed in the genomic data. The detailed analyses of many specimens allowed us to observe gene duplication and rearrangements within the mating loci, increasing the diversity within these loci. We supported previous suggestions of balancing selection in this region, an evolutionary force that maintains genomic diversity. These results supports that our fungal specimens are prone to outcross, which might facilitate the adaptation to new conditions. Identifying Cross-Cancer Similar Patients via a Semi-Supervised Deep Clustering Approach The treatment decisions for a cancer patient are typically based on the patients diagnosed cancer type. With the characterization of cancer tumors at the molecular level, there have been reports of patients being similar despite being diagnosed with different cancer types. Motivated from these observations, we aim at discovering cross-cancer patients, which we define as patients whose tumors are more similar to patient tumors diagnosed with another cancer type. We develop DeepCrossCancer to identify cross-cancer patients that always co-cluster with the other patient from another cancer type. The input to DeepCrossCancer is the transcriptomic profiles of the patient tumors, the age, and sex of the patient. To solve the clustering problem, we use a semi-supervised deep learning-based clustering method in which the clustering task is supervised by cancer type labels and the survival times of the patients. Applying the method to patient data from nine different cancers, we discover 20 cross-cancer patients that consistently co-cluster. By analyzing the predictive genes of the cross-cancer patients and other genomic information available for the patient such as somatic mutations and copy number variations, we identify striking genomic similarities across these patients providing support. The detection of cross-cancer patients opens up possibilities for transferring clinical decisions across patients at a single patient level. The source code is available at github.com/tastanlab/DeepCrossCancer Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. Glucose dynamics during ozone exposure measured using radiotelemetry: Stress drivers BackgroundStress-related neurobehavioral and metabolic disorders are associated with altered circulating adrenal-derived hormones and hyperglycemia. Temporal assessment of glucose and these hormones is critical for insights on an individuals health.

ObjectivesHere we used implantable-telemetry in rats to assess real-time changes in circulating glucose during and after exposure to the air pollutant ozone, and link responses to circulating neuroendocrine stress and metabolic hormones. We also compared rodent glucose and corticosterone (cortisol in humans) responses to humans exposed to ozone.

MethodsFirst, using a cross-over design, we monitored glucose levels during single or repeated ozone exposures (0.0, 0.2, 0.4 and 0.8-ppm) and non-exposure periods in male Wistar-Kyoto-rats implanted with glucose-telemeters. A second cohort of un-implanted rats was exposed to ozone (0.0, 0.4 or 0.8-ppm) for 30-min, 1-hour, 2-hour, or 4-hour with hormones measured immediately after exposure. Then we assessed glucose metabolism in sham and adrenalectomized rats with or without pharmacological interventions of adrenergic and glucocorticoid receptors. Finally, we assessed glucose and cortisol in serum samples from a clinical study involving exposure of human volunteers to air or 0.3 ppm ozone.

ResultsOzone (0.8-ppm) caused hyperglycemia and hypothermia beginning 90-min into exposure, with reversal of effects 4-6 hours post-exposure. Glucose monitoring during four daily 4-hour ozone exposures revealed duration of hyperglycemia, adaptation, and diurnal variations. Ozone-induced hyperglycemia was preceded by increased adrenocorticotropic hormone, corticosterone, and epinephrine, but depletion of thyroid-stimulating, prolactin, and luteinizing hormones. Hyperglycemia was inhibited in rats that were adrenalectomized and/or treated with glucocorticoid inhibitor. The depletion of cortisol was dampened in humans exposed to ozone during intermittent exercise.

DiscussionWe demonstrate for the first time the temporality of neuroendocrine-stress-mediated biological sequalae responsible for ozone-induced metabolic dysfunction as exposure occurs. Real-time glucose monitoring with stress hormones assessment may be useful in identifying interactions among pollutants and stress-related illnesses. Assigning metabolic rate measurements to torpor and euthermy in heterothermic endotherms: 'torpor', a new package for R. Torpor is a state of controlled reduction of metabolic rate (M) in endotherms. Assigning measurements of M to torpor or euthermy can be challenging, especially when the difference between euthermic M and torpid M is small, in species defending a high minimal body temperature in torpor, in thermolabile species, and slightly below the thermoneutral zone (TNZ). Here, we propose a novel method for distinguishing torpor from euthermy. We use the variation in M measured during euthermic rest and torpor at varying ambient temperatures (Ta) to objectively estimate the lower critical temperature (Tlc) of the TNZ and to assign measurements to torpor, euthermic rest or rest within TNZ. In addition, this method allows the prediction of M during euthermic rest and torpor at varying Ta, including resting M within the TNZ. The present method has shown highly satisfactory results using 28 published sets of metabolic data obtained by respirometry on 26 species of mammals. Ultimately, this novel method aims to facilitate analysis of respirometry data in heterothermic endotherms. Finally, the development of the associated R-package (torpor) will enable widespread use of the method amongst biologists.

Summary statementThe presented method and its associated R-package (torpor) enable the assignment of metabolic rate measurements to torpor or euthermy, ultimately improving the standardization of respirometry analyses in heterotherms. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Multivariate genome-wide association study on tissue-sensitive diffusion metrics identifies key molecular pathways for axonal growth, synaptogenesis, and astrocyte-mediated neuroinflammation The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discovered 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% were validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development were found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. INPP4B drives lysosome biogenesis to restrict leukemic stem cell differentiation and promote leukemogenesis Signaling pathways that control vital features of leukemic stem cells including multipotency, self-renewal, clonal expansion and quiescence remain unclear. Emerging studies illustrate critical roles for lysosomes in hematopoietic and leukemic stem cell fate. By investigating consequences of INPP4B alterations in AML, we have discovered its role in driving leukemic  stemness. We observed that INPP4B is highly expressed leukemic stem cell populations and Inpp4b-deficeint leukemias demonstrate increased disease latency, reduced leukemia initiating potential which is associated with a differentiated leukemic phenotype. Molecular analyses show that Inpp4b-deficient leukemias have compromised lysosomal gene expression, lysosomal content, and lysosomal activity. Our discovery of a novel pathway linking INPP4B, lysosomal biogenesis and leukemic stemness, provides a mechanism to explain the association of high INPP4B expression with poor AML prognosis, and highlights novel patient stratification strategies and LSC-specific leukemic therapies.

Key PointsOur findings highlight a novel pathway linking INPP4B, lysosomal function and leukemic stemness that explains the prognostic role of INPP4B in AML.
Our data reveal the utility of INPP4B as a biomarker of aggressive AML and provide a rationale to explore INPP4B and its associated function in lysosome biology as novel strategies to target LSC and AML Predicting the re-distribution of antibiotic molecules caused byinter-species interactions in microbial communities Microbes associate in nature forming complex communities, but they are often studied in purified form. Here I show that differences in antibiotic sensitivity between co-existing bacterial species can modify the diffusion of antibiotic molecules, and alter the drug efficacy measured in conventional in vitro assays. I developed a model, validated experimentally using an artificial two-species community, showing antibiotic molecules re-distribute evenly across all species if they have similar drug sensitivity. The result is lower drug efficacy with respect to conventional in vitro assays. However, when species have different sensitivities, the re-distribution of antibiotic is uneven and most sensitive species accumulate more drug. The result is increased drug efficacy against these species. Despite this discrepancy with conventional laboratory assays, my data suggests the changes in drug efficacy are predictable and provides insight into how communities can affect the diffusion of other molecules in the environment. Changes in Behavior and Salivary Serum Amyloid A in cows with Subclinical Mastitis The aim of this study was to identify detailed changes in behavior, and in salivary serum amyloid A (SAA), associated with subclinical mastitis. This included standard sickness behaviors (which form part of an adaptive response to conserve energy, minimize heat loss and facilitate recovery following infection and inflammation) and lesser-studied behaviors, that are biologically predicted to change at low-level infection, and therefore particularly relevant for detecting subclinical mastitis (here labelled  luxury behaviors). SAA is a physiological marker of systemic inflammation, with levels in milk and serum already known to increase during subclinical mastitis. We investigated whether the same was true for SAA in cow saliva. Data were collected for n = 34 commercial barn-housed dairy cows, tested in 17 matched-pairs. Each pair comprised a cow with subclinical mastitis (SCM) and a healthy control (CTRL), identified using somatic cell count (SCC) (SCM: SCC >200 x1000 cells/ml; CTRL: SCC <100 x1000 cells/ml). SCM cows were selected for study ad-hoc, at which point they were paired with a CTRL cow, based upon parity and calving date; consequently, the full data set was accrued over several months. Data were collected for each pair over 3 days: SCC (l4:00-l5:00h) Day 1; behavior (24h from 00:00h) Day 2; salivary serum amyloid-A (SAA) Day 3. We report, for the first time, that an increase in salivary SAA occurs during subclinical mastitis; SAA was higher in SCM cows and demonstrated a positive (weak) correlation with SCC. The behavioral comparisons revealed that SCM cows had reductions in activity (behavioral transitions and distance moved), social exploration, social reactivity (here: likelihood to move away/be displaced following receipt of agonism), performance of social grooming and head butts, and the receipt of agonistic non-contact challenges. In addition, SCM cows received more head swipes, and spent a greater proportion of time lying with their head on their flank than CTRL cows. SCM cows also displayed an altered feeding pattern; they spent a greater proportion of feeding time in direct contact with two conspecifics, and a lower proportion of feeding time at self-locking feed barriers, than CTRL cows. Behavioral measures were found to correlate, albeit loosely, with serum SAA in a direction consistent with predictions for sickness behavior. These included positive correlations with both lying duration and the receipt of all agonistic behavior, and negative correlations with feeding, drinking, the performance of all social and all agonistic behavior, and social reactivity. We conclude that changes in salivary SAA, social behavior, and activity offer potential in the detection of subclinical mastitis and recommend further investigation to substantiate and refine our findings. Sensing their plasma membrane curvature allows migrating cells to circumvent obstacles Cell migration is a hallmark out-of-equilibrium process in biology. In addition to persistent self-propelled motion, many cells display remarkable adaptive behaviors when they navigate complex environments within the body. Combining theory and experiments, we identify a curvature-sensing mechanism underlying obstacle avoidance in immune-like cells. The genetic perturbation of this machinery leads to a reduced capacity to evade obstructions combined with faster and more persistent cell migration in obstacle-free environments. We propose that the active polymerization of the actin cytoskeleton at the advancing edge of migrating cells is locally inhibited by the curvature-sensitive BAR protein Snx33 in regions with inward plasma membrane curvature. This coupling between actin and membrane dynamics leads to a mechanochemical instability that generates complex protrusive patterns at the cellular front. Adaptive motility thus arises from two simultaneous curvature-dependent effects, i) the specific reduction of propulsion in regions where external objects deform the plasma membrane and ii) the intrinsic patterning capacity due to the membrane-actin coupling that promotes spontaneous changes in the cells protrusions. Our results show how cells utilize actin- and plasma membrane biophysics to sense their environment, allowing them to adaptively decide if they should move ahead or turn away. On the basis of our findings, we propose that the natural diversity of BAR proteins may allow cells to tune their curvature sensing machinery to match the shape characteristics in their environment. Towards a unified model of naive T cell dynamics across the lifespan Naive CD4 and CD8 T cells are cornerstones of adaptive immunity, but the dynamics of their establishment early in life and how their kinetics change as they mature following release from the thymus are poorly understood. Further, due to the diverse signals implicated in naive T cell survival, it has been a long-held and conceptually attractive view that they are sustained by active homeostatic control as thymic activity wanes. Here we employ multiple experimental systems to identify a unified model of naive CD4 and CD8 T cell population dynamics across the mouse lifespan. We infer that both subsets divide rarely and progressively increase their survival capacity with cell age. Strikingly, this simple model captures naive CD4 T cell dynamics throughout life. In contrast, we find that newly generated naive CD8 T cells are lost more rapidly during the first 3-4 weeks of life, likely due to increased recruitment into memory. We find no evidence for elevated division rates in neonates, or for feedback regulation of naive T cell numbers at any age. We show how confronting mathematical models with diverse datasets can reveal a quantitative and remarkably simple picture of naive T cell dynamics from birth into old age. Structural basis for Nirmatrelvir in vitro efficacy against SARS-CoV-2 variants The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (Mpro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). We evaluated the in vitro catalytic activity and in vitro potency of nirmatrelvir against the main protease (Mpro) of prevalent variants of concern (VOC) or variants of interest (VOI): Alpha (, B.1.1.7), Beta ({beta}, B.1.351), Delta ({delta}, B1.617.2), Gamma ({gamma}, P.1), Lambda ({lambda}, B.1.1.1.37/C37), Omicron (o, B.1.1.529) as well as the original Washington or wildtype strain. These VOC/VOI carry prevalent mutations at varying frequencies in the Mpro specifically for: , {beta},{gamma} (K90R), {lambda} (G15S) and o (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant Mpros demonstrate that they are catalytically comparable to wildtype. Nirmatrelvir has similar potency against each mutant Mpro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the o, {lambda} and {beta} Mpro at 1.63 - 2.09 [A] resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells (1). Suppression of sost/sclerostin and dickkopf-1 promote intervertebral disc structure in mice Intervertebral disc (IVD) degeneration is a leading cause of low back pain and characterized by accelerated extracellular matrix breakdown and IVD height loss but there is no approved pharmacological therapeutic. Deletion of Wnt signaling receptor Lrp5 induces IVD degeneration and suggests that Wnt signaling in the IVD may be responsive to inhibition of Wnt signaling inhibitors sost(gene)/sclerostin(protein) or dickkopf-1 (dkk1). Anti-sclerostin antibody (Scl-Ab) is an FDA-approved bone therapeutic that activates Wnt signaling. We (1) determined if pharmacological neutralization of sclerostin, dkk1 or their combination stimulate Wnt signaling and promote IVD structure and (2) determined the extent of the response of the IVD to global, persistent deletion of sost. Nine-week-old C57Bl/6J female mice (n=6-7/grp) were subcutaneously injected 2x/wk for 5.5 wk with scl-Ab (25 mg/kg), dkk1-Ab (25 mg/kg), 3:1 scl-Ab/dkk1-Ab (18.75:6.25 mg/kg) or vehicle (Veh). Separately, IVD of sost KO and WT (wildtype) mice (n=8, grp) were harvested at 16 weeks of age. First, compared to vehicle, scl-Ab, dkk1-Ab and 3:1 scl-Ab/dkk1-Ab similarly increased lumbar IVD height and {beta}-catenin gene expression. Despite these similarities, scl-Ab decreased cellular stress-related heat shock protein gene expressions while neither dkk1-Ab nor scl-Ab/dkk1-Ab altered the same. Genetically and compared to WT, sost KO increased MRI-determined hydration and proteoglycan staining in the IVD. Notably, persistent deletion of sost was compensated by upregulation of dkk1, which consequently reduced the cell nuclear fraction for Wnt signaling transcription factor {beta}-catenin in whole IVD. Lastly, RNA-sequencing pathway analysis confirmed the compensatory suppression of Wnt signaling and determined a reduction of cellular stress pathways. Together, suppression of sost/sclerostin or dkk1 each promote IVD structure by stimulating Wnt signaling, but sclerostin and dkk1 may differentially regulate cellular stress pathways. Ultimately, postmenopausal women prescribed scl-Ab injections to prevent vertebral fracture may also benefit from a restoration of IVD height and health.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=72 SRC="FIGDIR/small/449486v2_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@db2e17org.highwire.dtl.DTLVardef@1af1bbborg.highwire.dtl.DTLVardef@145cc4corg.highwire.dtl.DTLVardef@167a19e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Suppression of Wnt signaling inhibitors by genetic or pharmacological approaches promotes intervertebral disc structure and hydration by Wnt signaling. However, persistent activation of Wnt signaling induces a compensatory reduction of Wnt signaling that shifts IVD cells toward a chondrocyte-like (CLC) phenotype. AF: annulus fibrosus, NC: notochordal cell, NP: nucleus pulposus, PG: proteoglycan 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. Characterization of regulatory transcriptional mechanisms in hepatocyte lipotoxicity Non-alcoholic fatty liver disease is a continuum of disorders among which non-alcoholic steatohepatitis (NASH) is particularly associated with a negative prognosis. Hepatocyte lipotoxicity is one of the main pathogenic factors of liver fibrosis and NASH. However, the molecular mechanisms regulating this process are poorly understood. The main aim of this study was to dissect transcriptional mechanisms regulated by lipotoxicity in hepatocytes. We achieved this aim by combining transcriptomic, proteomic and chromatin accessibility analyses from human liver and mouse hepatocytes. This integrative approach revealed several transcription factor networks deregulated by NASH and lipotoxicity. To validate these predictions, genetic deletion of the transcription factors MAFK and TCF4 was performed, resulting in hepatocytes that were better protected against saturated fatty acid oversupply. MAFK- and TCF4-regulated gene expression profiles suggest a mitigating effect against cell stress, while promoting cell survival and growth. Moreover, in the context of lipotoxicity, some MAFK and TCF4 target genes were to the corresponding differentially regulated transcripts in human liver fibrosis. Collectively, our findings comprehensively profile the transcriptional response to lipotoxicity in hepatocytes, revealing new molecular insights and providing a valuable resource for future endeavours to tackle the molecular mechanisms of NASH. 2,3-cAMP treatment mimics stress molecular response in Arabidopsis thaliana The role of the RNA degradation product 2,3-cyclic adenosine monophosphate (2,3-cAMP) is poorly understood. Recent studies have identified 2,3-cAMP in plant material and determined its role in stress signaling. The level of 2,3-cAMP increases upon wounding, dark, and heat, and 2,3-cAMP by binding to an RNA-binding protein, Rbp47b, promotes stress granule (SG) assembly. To gain further mechanistic insight into 2,3-cAMP function, we used a multi-omics approach combining transcriptomics, metabolomics, and proteomics to dissect Arabidopsis response to 2,3-cAMP treatment. We demonstrated that 2,3-cAMP is metabolized into adenosine, suggesting that the well-known cyclic nucleotide-adenosine pathway from human cells might also exist in plants. Transcriptomic analysis revealed only minor overlap between 2,3-cAMP-and adenosine-treated plants, suggesting that these molecules act through independent mechanisms. Treatment with 2,3-cAMP changed the levels of hundreds of transcripts, proteins, and metabolites, many previously associated with plant stress responses including protein and RNA degradation products, glucosinolates, chaperones and SG components. Finally, we demonstrated that 2,3-cAMP treatment influences the movement of processing bodies, supporting the role of 2,3-cAMP in the formation and motility of membraneless organelles.