Sound localization in world and head-centered space in ferrets The location of sounds can be described in multiple coordinate systems that are defined relative to ourselves, or the world around us. Evidence from neural recordings in animals point towards the existence of both head-centered and world-centered representations of sound location in the brain; however, it is unclear whether such neural representations have perceptual correlates in the sound localization abilities of non-human listeners. Here, we establish novel behavioral tests to determine the coordinate systems in which ferrets can localize sounds. We found that ferrets could learn to discriminate between sound locations that were fixed in either world-centered or head-centered space, across wide variations in sound location in the alternative coordinate system. Using probe sounds to assess broader generalization of spatial hearing, we demonstrated that in both head and world-centered tasks, animals used continuous maps of auditory space to guide behavior. Single trial responses of individual animals were sufficiently informative that we could then model sound localization using speaker position in specific coordinate systems and accurately predict ferrets actions in held-out data. Our results indicate that an animal model in which neurons are known to be tuned to sound location in egocentric and allocentric reference frames can also localize sounds in multiple head and world-centered spaces.

Significance StatementHumans can describe the location of sounds either relative to themselves, or in the world, independent of their momentary position. These different spaces are also represented in the activity of neurons in animals, but its not clear whether non-human listeners also perceive both head and world-centered sound location. Here, we designed two behavioral tasks in which ferrets had to discriminate between two sounds using their position in the world, or relative to the head. Subjects learnt to solve both problems and showed the ability to generalize sound location in each space when presented with infrequent probe sounds. These findings reveal a perceptual correlate of neural sensitivity previously observed in the ferret brain and establish that, like humans, ferrets can access an auditory map of their local environment.
Flexible control of Pavlovian-instrumental transfer based on expected reward value The Pavlovian-instrumental transfer (PIT) paradigm is widely used to assay the motivational influence of reward-predictive cues, reflected by their ability to invigorate instrumental behavior. Leading theories of incentive learning assume that cues are assigned motivational properties based on the value of the rewards they predict. We outline an alternative account of PIT which recognizes that reward-predictive cues may suppress rather than motivate instrumental behavior under certain conditions. We posit that cues capable of triggering a strong reward expectancy will tend to inhibit instrumental behavior, which is exploratory by nature, in order to facilitate retrieval of the expected reward before it is lost or stolen. According to this view, the motivation to engage in instrumental behavior during a cue should be inversely related to the value of the predicted reward, since there is more to lose by failing to secure a high-value reward than a low-value reward. We tested this hypothesis in rats using a high-conflict PIT protocol which elicits competing instrumental and food-port entry responses. In Experiment 1, cues signaling different reward magnitudes (1, 3, or 9 pellets) differentially influenced these behaviors. Whereas the 1-pellet cue increased instrumental behavior, cues signaling 3 or 9 pellets suppressed instrumental behavior and elicited high levels of food-port activity. Experiment 2 found that reward-predictive cues suppressed instrumental behavior and increased food-port activity in a flexible manner that was disrupted by post-training reward devaluation. We discuss how the PIT task provides a useful tool for studying cognitive control over cue-motivated behavior in rodents.
Sharpened and mechanically durable carbon fiber electrode arrays for neural interfacing Bioelectric medicine treatments target disorders of the nervous system unresponsive to pharmacological methods. While current stimulation paradigms effectively treat many disorders, the underlying mechanisms are relatively unknown, and current neuroscience recording electrodes are often limited in their specificity to gross averages across many neurons or axons. Here, we develop a novel, durable carbon fiber electrode array adaptable to many neural structures for precise neural recording. Carbon fibers were sharpened using a reproducible blowtorch method that uses the reflection of fibers against the surface of a water bath. The arrays were developed by partially embedding carbon fibers in medical-grade silicone to improve durability. We recorded acute spontaneous electrophysiology from the rat cervical vagus nerve (CVN), feline dorsal root ganglia (DRG), and rat brain. Blowtorching resulted in fibers of 72.3 {+/-} 33.5 degree tip angle with 146.8 {+/-} 17.7 m exposed carbon. Silicone-embedded carbon fiber arrays were durable with 87.5% of fibers remaining after 50,000 passes. Observable neural clusters were recorded using sharpened carbon fiber electrodes from rat CVN (41.8 Vpp), feline DRG (101.1 Vpp), and rat brain (80.7 Vpp). Recordings from the feline DRG included physiologically-relevant signals from increased bladder pressure and cutaneous brushing. These results suggest that this carbon fiber array is a uniquely durable and adaptable neural recording device. In the future, this device may be useful as a bioelectric medicine tool for diagnosis and closed-loop neural control of therapeutic treatments and monitoring systems. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. An Essential Role of UBXN3B in B Lymphopoiesis Hematopoiesis is finely regulated to enable timely production of the right numbers and types of mature immune cells to maintain tissue homeostasis. Dysregulated hematopoiesis may compromise antiviral immunity and/or exacerbate immunopathogenesis. Herein, we report an essential role of UBXN3B in maintenance of hematopoietic homeostasis and restriction of immunopathogenesis during respiratory viral infection. Ubxn3b deficient (Ubxn3b-/-) mice are highly vulnerable to SARS-CoV-2 and influenza A infection, characterized by more severe lung immunopathology, lower virus-specific IgG, significantly fewer B cells, but more myeloid cells than Ubxn3b+/+ littermates. This aberrant immune compartmentalization is recapitulated in uninfected Ubxn3b-/- mice. Mechanistically, UBXN3B controls precursor B-I (pre-BI) transition to pre-BII and subsequent proliferation in a cell-intrinsic manner, by maintaining BLNK protein stability and pre-BCR signaling. These results reveal an essential role of UBXN3B for the early stage of B cell development. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. 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. Postmitotic accumulation of histone variant H3.3 in new cortical neurons establishes neuronal chromatin, transcriptome, and identity Histone variants, which can be expressed outside of S-phase and deposited DNA synthesis-independently, provide long-term histone replacement in postmitotic cells, including neurons. Beyond replenishment, histone variants also play active roles in gene regulation by modulating chromatin states or enabling nucleosome turnover. Here, we uncover crucial roles for the histone H3 variant H3.3 in neuronal development. We find that newborn cortical excitatory neurons, which have only just completed replication-coupled deposition of canonical H3.1 and H3.2, substantially accumulate H3.3 immediately post mitosis. Co-deletion of H3.3-encoding genes H3f3a and H3f3b from newly postmitotic neurons abrogates H3.3 accumulation, markedly alters the histone posttranslational modification (PTM) landscape, and causes widespread disruptions to the establishment of the neuronal transcriptome. These changes coincide with developmental phenotypes in neuronal identities and axon projections. Thus, preexisting, replication-dependent histones are insufficient for establishing neuronal chromatin and transcriptome; de novo H3.3 is required. Stage-dependent deletion of H3f3a and H3f3b from (1) cycling neural progenitor cells, (2) neurons immediately post mitosis, or (3) several days later, reveals the first postmitotic days to be a critical window for de novo H3.3. After H3.3 accumulation within this developmental window, co-deletion of H3f3a and H3f3b does not lead to immediate H3.3 loss, but causes progressive H3.3 depletion over several months without widespread transcriptional disruptions or cellular phenotypes. Our study thus uncovers key developmental roles for de novo H3.3 in establishing neuronal chromatin, transcriptome, identity, and connectivity immediately post mitosis that are distinct from its role in maintaining total histone H3 levels over the neuronal lifespan.

SignificanceDNA is packaged around histones into chromatin, which compacts the genome, but also restricts access to DNA. Gene transcription thus requires chromatin reorganization that is precisely regulated, including via variant forms of histones. Here, we find that during a critical developmental window for establishing postmitotic neuronal identity, newly generated cortical excitatory neurons substantially accumulate the histone H3 variant H3.3. Conditional deletion of H3.3-encoding genes from new neurons abrogates de novo H3.3 accumulation, and broadly disrupts neuronal histone modifications, gene expression, subtype identity, and axon projections. Thus, preexisting H3 histones are insufficient for establishing neuronal chromatin and transcriptome; de novo H3.3 is essential. This developmental requirement for H3.3 is distinct from H3.3 contribution to long-term maintenance of histones in mature neurons. Community Evaluation of Glycoproteomics Informatics Solutions Reveals High-Performance Search Strategies of Serum N- and O-Glycopeptide Data Glycoproteome profiling (glycoproteomics) is a powerful yet analytically challenging research tool. The complex tandem mass spectra generated from glycopeptide mixtures require sophisticated analysis pipelines for structural determination. Diverse software aiding the process have appeared, but their relative performance remains untested. Conducted through the HUPO Human Proteome Project - Human Glycoproteomics Initiative, this community study, comprising both developers and users of glycoproteomics software, evaluates the performance of informatics solutions for system-wide glycopeptide analysis. Mass spectrometry-based glycoproteomics datasets from human serum were shared with all teams. The relative team performance for N- and O-glycopeptide data analysis was comprehensively established and validated through orthogonal performance tests. Excitingly, several high-performance glycoproteomics informatics solutions were identified. While the study illustrated that significant informatics challenges remain, as indicated by a high discordance between annotated glycopeptides, lists of high-confidence (consensus) glycopeptides were compiled from the standardised team reports. Deep analysis of the performance data revealed key performance-associated search variables and led to recommendations for improved "high coverage" and "high accuracy" glycoproteomics search strategies. This study concludes that diverse software for comprehensive glycopeptide data analysis exist, points to several high-performance search strategies, and specifies key variables that may guide future software developments and assist informatics decision-making in glycoproteomics. 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. First evidence of in vitro cytotoxic effects of marine microlitter on Merluccius merluccius and Mullus barbatus, two Mediterranean commercial fish species Marine litter is composed mainly of plastics and is recognized as a serious threats to marine ecosystems. Ecotoxicological approaches have started elucidating the potential severity of microplastics (MPs) in controlled laboratory studies with pristine materials but no information exist on marine environmental microlitter as a whole. Here, we characterized the litter in the coastal Northern Tyrrhenian sea and in the stomach of two fish species of socio-economic importance, and exposed primary cell cultures of mucosal and lymphoid organs to marine microlitter for evaluating possible cytotoxic effects. An average of 0.30 {+/-} 0.02 microlitter items m-3 was found in water samples. FT-IR analysis revealed that plastic particles, namely HDPE, polyamide and polypropylene were present in 100% and 83.3% of Merluccius merluccius and Mullus barbatus analyzed, which overall ingested 14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual, respectively. Moreover, microlitter was confirmed as a vector of microorganisms. Lastly, the apical end-point of viability was found to be significantly reduced in splenic cells exposed in vitro to two microlitter conditions. Considering the role of the spleen in the mounting of adaptive immune responses, our results warrant more in-depth investigations for clarifying the actual susceptibility of these two species to anthropogenic microlitter.

HighlightsO_LI0.30 {+/-} 0.02 microlitter items m-3 were found at the surface of coastal Northern Tyrrhenian sea
C_LIO_LI14.67 {+/-} 4.10 and 5.50 {+/-} 1.97 items/individual were retrieved from the stomach of hakes and mullets
C_LIO_LIThe ingested microlitter contained plastic items
C_LIO_LIMicrolitter was validated as a carrier of bacteria, fungi and flagellates
C_LIO_LISplenic cells exposed to two microlitter conditions for 72 hours suffered cytotoxicity
C_LI Symbiodiniaceae conduct under natural bleaching stress during advanced gametogenesis stages of a mesophotic coral The mesophotic coral Alveopora allingi from the northern Gulf of Eilat/Aqaba, Red Sea, is affected by year-round partial coral-bleaching events. During these events, the migration of Symbiodiniaceae takes place from the coral-host mesoglea to the developed oocytes in bleached parts of colonies of A. allingi but not in the non-bleached parts. Additionally, these oocytes are abnormal, missing part of the structural material of the peripheral areas and are also significantly larger in the bleached areas of the colonies. Hence, we suggest a parasitic behavior of the symbionts or a commensalism relationship which enhance symbionts needs during bleaching periods and may boost the gametogenesis development in these corals. We propose that evolutionarily, this behavior may greatly contribute to the symbiont community survival throughout the bleaching period, and it can also be beneficial for the hosts persistence and adaptation to bleaching through the acquisition of a specific symbiont community following the bleaching event. Partial exogastrulation due to apical-basal polarity of F-actin distribution disruption in sea urchin embryo by omeprazole Gastrulation is a universal process in the morphogenesis of many animal embryos. Although morphological and molecular events in gastrulation have been well studied, the mechanical driving forces and underlying regulatory mechanisms are not fully understood. Here, we investigated the gastrulation of embryos of a sea urchin, Hemicentrotus pulcherrimus, which involves the invagination of a single-layered vegetal plate into the blastocoel. We observed that omeprazole, a proton pump inhibitor capable of perturbing the left-right asymmetry of sea urchin embryo, induced "partial exogastrulation" where the secondary invagination proceeds outward. During early gastrulation, intracellular apical-basal polarity of F-actin distribution in vegetal half were higher than those in animal half, while omeprazole treatment disturbed the apical-basal polarity of F-actin distribution in vegetal half. Furthermore, gastrulation stopped and even partial exogastrulation did not occur when F-actin polymerization or degradation in whole embryo was partially inhibited via RhoA or YAP1 knockout. A mathematical model of the early gastrulation reproduced the shapes of both normal and exogastrulating embryos using cell-dependent cytoskeletal features based on F-actin. Additionally, such cell position-dependent intracellular F-actin distributions might be regulated by intracellular pH distributions. Therefore, apical-basal polarity of F-actin distribution disrupted by omeprazole may induce the partial exogastrulation via anomalous secondary invagination. Extracellular matrix and spores impart structural hierarchy and compositional heterogeneity in Bacillus subtilis biofilms Biofilms are multicellular microbial communities that encase themselves in an extracellular matrix (ECM) of secreted biopolymers and attach to surfaces and interfaces. Bacterial biofilms are detrimental in hospital and industrial settings, but they can be beneficial in agricultural contexts. An essential property of biofilms that grants them with increased survival relative to planktonic cells is phenotypic heterogeneity; the division of the biofilm population into functionally distinct subgroups of cells. Phenotypic heterogeneity in biofilms can be traced to the cellular level, however, the molecular structures and elemental distribution across whole biofilms as well as possible linkages between them remain unexplored. Mapping X-ray diffraction (XRD) across intact biofilms in time and space, we revealed the dominant structural features in Bacillus subtilis biofilms, stemming from matrix components, spores and water. By simultaneously following the X-ray fluorescence (XRF) signal of biofilms and isolated matrix components, we discovered that the ECM preferentially binds calcium ions over other metal ions, specifically, zinc, manganese and iron. These ions, remaining free to flow below macroscopic wrinkles that act as water channels, eventually accumulate and lead to sporulation. The possible link between ECM properties, regulation of metal ion distribution and sporulation across whole intact biofilms unravels the importance of molecular-level heterogeneity in shaping biofilm physiology and development.

Significance StatementBiofilms are multicellular soft microbial communities that are able to colonize synthetic surfaces as well as living organisms. To survive sudden environmental changes and efficiently share their common resources, cells in a biofilm divide into subgroups with distinct functions, leading to phenotypic heterogeneity. Here, by studying intact biofilms by synchrotron X-ray diffraction and fluorescence, we revealed correlations between biofilm macroscopic architectural heterogeneity and the spatio-temporal distribution of extracellular matrix, spores, water and metal ions. Our findings demonstrate that biofilm heterogeneity is not only affected by local genetic expression and cellular differentiation, but also by passive effects resulting from the physicochemical properties of the molecules secreted by the cells, leading to differential distribution of nutrients that propagates through macroscopic length scales. The multi-level phenotypic impact of synonymous substitutions: heterologous gene expression in human cells. Redundancy in the genetic code allows for differences in transcription and/or translation efficiency between sequences carrying synonymous polymorphisms, potentially leading to phenotypic changes. It is commonly admitted that the evolution codon usage bias ("CUB", the over-representation of certain codons in a genome, a gene or in positions along a gene) are driven by a combination of neutral and selective processes, but their relative contribution is a matter of debate, especially in mammals. Particularly, integrative studies quantifying the phenotypic impact of CUB at different molecular and cellular levels are lacking. Here we report a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells.

Six synonymous versions of the shble antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was assessed by means of: i) mRNA-to-DNA and protein-to-mRNA ratios for each shble version; ii) cellular fluorescence, using flow cytometry, as a proxy for single-cell level expression; and iii) real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness.

We show that differences in CUB strongly impact the molecular and the cellular phenotype: i) they result in large differences in mRNA and in protein levels, as well in mRNA-to-protein ratio; ii) they introduce splicing events not predicted by current algorithms; iii) they lead to reproducible phenotypic heterogeneity; iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression.

We interpret that CUB modulate mRNA availability and suitability for translation in human cells, leading to differences in protein levels and eventually eliciting phenotypic differences.

AUTHOR SUMMARYThe genetic code is redundant, with several codons encoding for the same amino acid. These synonymous codons are not used with equal frequencies. Instead, codon usage bias (CUB) varies between species, genes, and even positions along a gene. At each of these levels, CUB are shaped by the overall balance between mutational biases and selection forces. To shed light on the molecular mechanisms underlying molecular and organism phenotypes, integrative studies quantifying the phenotypic impact of CUB at different levels of biological integration are necessary. Here, we monitored the multiscale changes induced by six synonymous versions of an antibiotic resistance gene independently expressed in a human cell line. We show that: 1. both mRNA levels, protein levels, and their ratios are affected by CUB; 2. potential effects on mRNA (splicing and availability for translation) seem to constitute the main level of action; 3. cell fitness is severely impacted by a trade-off between the burden of heterologous expression and the benefits of antibiotic resistance. This integrative study provides new insights on translation regulation and the associated phenotypic impact in human cells, associated to CU. Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the -subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitutions effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases.

Significance StatementIndividual amino acid residues within a protein work in concert to produce a functionally coherent structure that must be maintained even as orthologous proteins in different species diverge over time. Given this dependence, we expect identical mutations to have more similar effects on protein function in more closely related species. We tested this hypothesis by performing protein-engineering experiments on ATP1A, an enzyme mediating target-site insensitivity to cardiotonic steroids (CTS) in diverse animals. These experiments reveal that the phenotypic effects of substitutions can sometimes be background-dependent, but also that the magnitude of these phenotypic effects does not correlate with overall levels of ATP1A sequence divergence. Our results suggest that epistatic constraints are determined by states at a small number of sites, potentially explaining the frequent convergent CTS resistance substitutions among Na,K-ATPases of highly divergent taxa. FACT regulates pluripotency through distal regulation of gene expression in murine embryonic stem cells The FACT complex is a conserved histone chaperone with essential roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. To determine this mechanism, we used inducible depletion of FACT subunit SPT16 in murine embryonic stem cells paired with genomic factor localization, nascent transcription, and chromatin accessibility analyses. Over a timecourse of SPT16 depletion, nucleosomes invade loci bound by master pluripotency factors and gene-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (OCT4), Sox2, Nanog, and enhancer RNAs produced at the genes' associated enhancers are downregulated, suggesting that FACT regulates expression of the pluripotency factors themselves. We find that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes IntroductionIncreasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero.

MethodsWe propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F -penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

ResultsWe evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F -penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

Discussion/ConclusionWe have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach. A dynamical systems treatment of transcriptomic trajectories in hematopoiesis Inspired by Waddingtons illustration of an epigenetic landscape, cell-fate transitions have been envisioned as bifurcating dynamical systems, wherein the dynamics of an exogenous signal couples to a cells enormously complex signaling and transcriptional machinery, eliciting a qualitative transition in the collective state of a cell - its fate. It remains unclear, however, whether the dynamical systems framework can go beyond a word-based caricature of the system and provide sharp quantitative insights that further our understanding of differentiation. Single-cell RNA sequencing (scRNA-seq), which measures the distributions of possible transcriptional states in large populations of differentiating cells, provides an alternate view, in which development is marked by the individual concentration variations of a myriad of genes. Here, starting from formal mathematical derivations, we challenge these transcriptomic trajectories to a rigorous statistical evaluation of whether they display signatures consistent with bifurcations. After pinpointing bifurcations along transcriptomic trajectories of the neutrophil branch of hematopoeitic differentiation we are able to further leverage the primitive features of a linear instability to identify the single-direction in gene expression space along which the bifurcation unfolds and identify possible gene contributors. This scheme identifies transcription factors long viewed to play a crucial role in the process of neutrophil differentiation in addition to identifying a host of other novel genetic players. Most broadly speaking, we provide evidence that, though very high-dimensional, a bifurcating dynamical systems formalism might be appropriate for the process of cellular differentiation and that it can be leveraged to provide insights. Ambitiously, our work attempts to take a step beyond data-analysis and towards the construction of falsifiable mathematical models that describe the dynamics of the entire transcriptome. Living with high potassium: how plants balance nutrient acquisition during K-induced salt stress O_LIHigh potassium (K) in the growth medium is more toxic to plants than Na at similar concentrations. However, the molecular mechanisms underlying plant responses to K-induced salt stress are virtually unknown.
C_LIO_LIWe examined Arabidopsis thaliana and its extremophyte relative Schrenkiella parvula, using a comparative multi-omics approach to identify cellular processes affected by excess K and understand which deterministic regulatory pathways are active to avoid tissue damage while sustaining growth.
C_LIO_LIA. thaliana showed limited capacity to curb excess K accumulation and prevent nutrient depletion contrasting to S. parvula which could limit excess K accumulation without restricting nutrient uptake. Facilitated by a targeted transcriptomic response, promoting nitrogen uptake along with other key nutrients and uninterrupted N assimilation into primary metabolites during excess K-stress allowed S. parvula to boost its antioxidant and osmolyte pools concurrently leading to sustained growth. Antithetically, A. thaliana showed transcriptional responses indicative of a poor balance between stress signaling, increased ROS levels, and reduced photosynthesis, subsequently leading to inhibited growth.
C_LIO_LIThe ability to regulate independent nutrient uptake and a coordinated transcriptomic response to avoid non-specific stress signaling are two main deterministic steps towards building stress resilience to excess K+-induced salt stress.
C_LI SARS-CoV-2 triggers Golgi fragmentation via down-regulation of GRASP55 to facilitate viral trafficking The ongoing COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), an enveloped RNA virus. Despite the high economic and life losses caused by SARS-CoV-2, the detailed viral cycle, especially how it assembles and traffics in the secretory pathway, remains largely unknown. Here, we show that SARS-CoV-2 infection induces global alterations of the host endomembrane system, including dramatic Golgi fragmentation. Disrupting Golgi function with small molecules strongly inhibits viral infection. Furthermore, expression of several SARS-CoV-2 proteins individually is sufficient to trigger Golgi fragmentation. Significantly, SARS-CoV-2 infection down-regulates GRASP55 but up-regulates TGN46 expression, while expression of GRASP55 or knockdown of TGN46 reduces the infection rate of both USA-WA1 and Delta variants of SARS-CoV-2. Our study reveals that SARS-CoV-2 modulates Golgi structure and function via altering GRASP55 and TGN46 expression to facilitate viral trafficking, indicating the Golgi as a novel therapeutic target to block SARS-CoV-2 infection. Gene expression and DNA methylation consistently predict cancer mutation states In studies of cellular function in cancer, researchers are increasingly able to choose from many -omics assays as functional readouts. Choosing the correct readout for a given study can be difficult, and which layer of cellular function is most suitable to capture the relevant signal may be unclear. In this study, we consider prediction of cancer mutation status (presence or absence) from functional -omics data as a representative problem. Since functional signatures of cancer mutation have been identified across many data types, this problem presents an opportunity to quantify and compare the ability of different -omics readouts to capture signals of dysregulation in cancer. The TCGA Pan-Cancer Atlas contains genetic alteration data including somatic mutations and copy number variants (CNVs), as well as several -omics data types. From TCGA, we focus on RNA sequencing, DNA methylation arrays, reverse phase protein arrays (RPPA), microRNA, and somatic mutational signatures as -omics readouts.

Across a collection of genes recurrently mutated in cancer, RNA sequencing tends to be the most effective predictor of mutation state. Surprisingly, we found that for many of the genes we considered, one or more other data types are approximately equally effective predictors. Performance was more variable between mutations than it was between data types for the same mutation, and there was often little difference between the top data types. We also found that combining data types into a single multi-omics model provided little or no improvement in predictive ability over the best individual data type. Based on our results, for the design of studies focused on the functional outcomes of cancer mutations, there are often multiple -omics types that can serve as effective readouts, although gene expression seems to be a reasonable default option. 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. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Inter-tissue convergence of gene expression during ageing suggests age-related loss of tissue and cellular identity Developmental trajectories of gene expression may reverse in their direction during ageing, a phenomenon previously linked to cellular identity loss. Our analysis of cerebral cortex, lung, liver and muscle transcriptomes of 16 mice, covering development and ageing intervals, revealed widespread but tissue-specific ageing-associated expression reversals. Cumulatively, these reversals create a unique phenomenon: mammalian tissue transcriptomes diverge from each other during postnatal development, but during ageing, they tend to converge towards similar expression levels, a process we term Divergence followed by Convergence, or DiCo. We found that DiCo was most prevalent among tissue-specific genes and associated with loss of tissue identity, which is confirmed using data from independent mouse and human datasets. Further, using publicly available single-cell transcriptome data, we showed that DiCo could be driven both by alterations in tissue cell type composition and also by cell-autonomous expression changes within particular cell types. CHRM1 is a Druggable Melanoma Target Whose Endogenous Activity is Determined by Inherited Genetic Variation in DOPA Production Melanoma risk is 30 times higher in people with lightly pigmented skin compared to those with darkly pigmented skin. Here we show that this difference results from more than melanin pigment and its ultraviolet radiation (UVR) shielding effect. Using primary human melanocytes representing the full human skin pigment continuum and several preclinical melanoma models, we show that cell-intrinsic differences between dark and light melanocytes regulate melanocyte proliferative capacity, overall cellular differentiation state, and susceptibility to malignant transformation, independently of melanin and UV exposure. We determined that these differences result from dihydroxyphenylalanine (DOPA), a melanin precursor synthesized at higher levels in melanocytes from dark skin. Although DOPA was not previously known to have specific signaling activity, we used both high throughput pharmacologic and genetic in vivo CRISPR screens to determine that DOPA limits melanocyte and melanoma cell proliferation by directly inhibiting the muscarinic acetylcholine receptor M1 (CHRM1), a G Protein-Coupled Receptor (GPCR) not previously known to bind DOPA, nor to affect melanoma pathobiology. Pharmacologic CHRM1 antagonism in melanoma leads to depletion of c-Myc and FOXM1, both of which are proliferation drivers associated with aggressive melanoma. In preclinical mouse melanoma models using both immune deficient and syngeneic immune competent mice, pharmacologic inhibition of CHRM1 or FOXM1 inhibited tumor growth. CHRM1 and FOXM1 may be new therapeutic targets for melanoma.