Value, confidence, deliberation: a functional partition of the medial prefrontal cortex demonstrated across rating and choice tasks Deciding about courses of action involves minimizing costs and maximizing benefits. Decision neuroscience studies have implicated both the ventral and dorsal medial prefrontal cortex (vmPFC and dmPFC) in signaling goal value and action cost, but the precise functional role of these regions is still a matter of debate. Here, we suggest a more general functional partition that applies not only to decisions but also to judgments about goal value (expected reward) and action cost (expected effort). In this conceptual framework, cognitive representations related to options (reward value and effort cost) are dissociated from metacognitive representations (confidence and deliberation) related to solving the task (providing a judgment or making a choice). Thus, we used an original approach with the goal of identifying consistencies across several preference tasks, from likeability ratings to binary decisions involving both attribute integration and option comparison. FMRI results confirmed the vmPFC as a generic valuation system, its activity increasing with reward value and decreasing with effort cost. In contrast, more dorsal regions were not concerned with the valuation of options but with metacognitive variables, confidence level being reflected in mPFC activity and deliberation time in dmPFC activity. Thus, there was a dissociation between the effort attached to choice options (represented in the vmPFC) and the effort invested in deliberation (represented in the dmPFC), the latter being expressed in pupil dilation. More generally, assessing commonalities across preference tasks might help reaching a unified view of the neural mechanisms underlying the cost/benefit tradeoffs that drive human behavior.
Delineation of insula dynorphin and kappa opioid receptor circuit elements in alcohol consumption Alcohol use disorder is complex and multi-faceted, involving the engagement of multiple signaling systems across numerous brain regions to drive pathological behavior. Previous work has indicated that both the insular cortex and dynorphin (DYN)/Kappa opioid receptor (KOR) systems contribute to excessive alcohol use. More recently, we identified a microcircuit in the medial aspect of the insular cortex that signals through DYN/KOR. Here, we explored the role of insula DYN/KOR circuit elements on alcohol intake in a long-term intermittent access (IA) procedure. Using a combination of conditional knockout strategies and site-directed pharmacology, we discovered distinct and sex-specific roles for insula DYN and KOR systems in alcohol drinking and related behavior. Our findings show that insula DYN deletion blocked escalated consumption and decreased overall intake of and preference for alcohol in male and female mice. This effect was specific to alcohol in male mice, as DYN deletion did not impact their sucrose intake. In male mice, insula KOR antagonism reduced alcohol intake and preference during the early phase of IA only. Alcohol consumption as not affected by insula KOR knockout. In addition, we found that long-term IA decreased the intrinsic excitability of DYN and deep layer pyramidal neurons (DLPN) in the insula of male mice. Excitatory synaptic transmission was also impacted by IA, as it drove an increase in excitatory synaptic drive in both DYN neurons and DLPN. Combined, our findings suggest that there is a dynamic interplay between excessive alcohol consumption and insula DYN/KOR circuit elements.
Acetyltransferases GCN5 and PCAF are required for B lymphocyte maturation in mice B lymphocyte development has two DNA recombination processes: V(D)J recombination of the immunoglobulin (Igh) gene variable region and class switching of the Igh constant regions from IgM to IgG, IgA, or IgE. V(D)J recombination is required for successful maturation of B cells from pro-B to pre-B to immature-B and then to mature B cells in the bone marrow. CSR occurs outside of the bone marrow when mature B cells migrate to peripheral lymphoid organs, such as spleen and lymph nodes. Both V(D)J recombination and CSR depend on an open chromatin state that makes DNA accessible to specific enzymes, recombination activating gene (RAG), and activation-induced cytidine deaminase (AID). Acetyltransferases GCN5 and PCAF possess redundant functions acetylating histone H3 lysine 9 (H3K9). Here, we generated a mouse model that lacks both GCN5 and PCAF in B cells. We found that double-deficient mice possess low levels of mature B cells in the bone marrow and peripheral organs, an accumulation of pro-B cells in bone marrow, and reduced CSR levels. We conclude that both GCN5 and PCAF are required for B cell development in vivo. RNA-guided Cas12a and Cas9 knockout and knockin of the omega-1 locus of Schistosoma mansoni We compared the efficiency of gene knockout (KO) and precision of insertion (knock-in, KI) of the RNA-guided AsCas12a nuclease of Acidaminococcus sp. with that of SpCas9 from Streptococcus pyogenes, aiming to enhance the functional genomics toolkit for Schistosoma mansoni. Programmed DNA cleavages catalyzed by Cas12a and Cas9 result in staggered and blunt ended strand breaks, respectively. TTTV, the optimal protospacer adjacent motif for AsCas12a would occur frequently within the AT-rich genome of this platyhelminth. We deployed optimized conditions for the ratio of guide RNAs to the nuclease, donor templates, and electroporation parameters, to target a key enzyme termed omega-1 that is secreted by the schistosome egg. AsCas12a was more efficient than SpCas9 for gene knockout of omega-1 as determined by tracking of indels by decomposition (P < 0.001). Resulting from CRISPREsso2 analysis, most mutations were deletions; SpCas9 induced short deletions of 3 nt in length whereas AsCas12a induced deletions of 2 to 26 nt. Knockout efficiency of both nucleases markedly increased in the presence of short, single stranded oligodeoxynucleotide (ssODN) donor templates. With AsCas12a, ssODNs representative of both the non-CRISPR target (NT) and target (T) strands of the targeted gene were tested, resulting in KO efficiencies of 15.67, 28.71 and 21.43% in the SpCas9 plus donor ssODN, AsCas12a plus NT-ssODN, and AsCas12a plus T-ssODN groups, respectively. Trans cleavage activity against the ssODNs by activated AsCas12a was not apparent in vitro. Programmed SpCas9 editing led to more precise transgene insertion than AsCas12a, with KI efficiencies of 17.07% for the KI_SpCas9 group, 14.58% for KI_AsCas12a-NT-ssODN and 12.37% for KI_AsCas12a-T-ssODN. Although AsCas12a induced fewer mutations per genome than SpCas9, the phenotypic impact on transcription and expression of omega-1 was similar for both nucleases. These findings revealed that AsCas12a and SpCas9 both provide tractable routes for RNA-guided programmed mutation of the genome of the schistosome egg. A family of unusual immunoglobulin superfamily genes in an invertebrate histocompatibility complex Most colonial marine invertebrates are capable of allorecognition, the ability to distinguish between themselves and conspecifics. One long-standing question is whether invertebrate allorecognition genes are homologous to vertebrate histocompatibility genes. In the cnidarian Hydractinia symbiolongicarpus, allorecognition is controlled by at least two genes, Allorecognition 1 (Alr1) and Allorecognition 2 (Alr2), which encode highly polymorphic cell surface proteins that serve as markers of self. Here, we show that Alr1 and Alr2 are part of a family of 41 Alr genes, all of which reside a single genomic interval called the Allorecognition Complex (ARC). Using sensitive homology searches and highly accurate structural predictions, we demonstrate that the Alr proteins are members of the immunoglobulin superfamily (IgSF) with V-set and I-set Ig domains unlike any previously identified in animals. Specifically, their primary amino acid sequences lack many of the motifs considered diagnostic for V-set and I-set domains, yet they adopt secondary and tertiary structures nearly identical to canonical Ig domains. Thus, the V-set domain, which played a central role in the evolution of vertebrate adaptive immunity, was present in the last common ancestor of cnidarians and bilaterians. Unexpectedly, several Alr proteins also have immunoreceptor tyrosine-based activation motifs (ITAMs) and immunoreceptor tyrosine-based inhibitory motifs (ITIMs) in their cytoplasmic tails, suggesting they could participate in pathways homologous to those that regulate immunity in humans and flies. This work expands our definition of the IgSF with the addition of a family of unusual members, several of which play a role in invertebrate histocompatibility.

Significance StatementThe immunoglobulin superfamily (IgSF) is one of the largest and most functionally versatile domain families in animal genomes. Although their amino acid sequences can vary considerably, IgSF domains have been traditionally defined by conserved residues at several key positions in their fold. Here, we sequenced an invertebrate histocompatibility complex and discovered a family of IgSF genes with amino acid sequences that lack most of these residues yet are predicted to adopt folds virtually identical to canonical V-set and I-set IgSF domains. This work broadens the definition of the IgSF and shows that the V-set domain was present earlier in animal evolution than previously appreciated. Interaction among bacterial cells triggers exit from lag phase The relationship between the number of cells colonizing a new environment and time for resumption of growth is a subject of long-standing interest. In microbiology this is known as the "inoculum effect". Its mechanistic basis is unclear with possible explanations ranging from the independent actions of individual cells, to collective actions of populations of cells. Progress requires precise measurement of lag time distributions while at the same time, experimentally controlling inoculum size. Here we use a millifluidic droplet device in which the growth dynamics of hundreds of populations founded by different numbers of Pseudomonas fluorescens cells, ranging from a single cell, to one thousand cells, were followed in real time. Our data show that lag phase decreases with inoculum size. The average decrease, variance across droplets, and distribution shapes, follow predictions of extreme value theory, where the inoculum lag time is determined by the minimum value sampled from the single-cell distribution. Our experimental results show that exit from lag phase depends on strong interactions among cells, consistent with a "leader-cell" triggering end of lag phase for the entire population. First detection of herpesvirus and frequency of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Nanocrown electrodes for reliable and robust intracellular recording of cardiomyocytes and cardiotoxicity screening Drug-induced cardiotoxicity arises primarily when a compound alters the electrophysiological properties of cardiomyocytes. Features of intracellular action potentials (iAPs) are powerful biomarkers that predict proarrhythmic risks. However, the conventional patch clamp techniques for measuring iAPs are either laborious and low throughput or not suitable for measuring electrically connected cardiomyocytes. In the last decade, a number of vertical nanoelectrodes have been demonstrated to achieve parallel and minimally-invasive iAP recordings. Nanoelectrodes show great promise, but the large variability in success rate, signal strength, and the low throughput of device fabrication have hindered them from being broadly adopted for proarrhythmia drug assessment. In this work, we developed vertically-aligned and semi-hollow nanocrown electrodes that are mechanically robust and made through a scalable fabrication process. Nanocrown electrodes achieve >99% success rates in obtaining intracellular access through electroporation, allowing reliable and simultaneous iAP recordings from up to 57 human pluripotent stem-cell-derived cardiomyocytes (hPSC-CMs). The accuracy of nanocrown electrode recordings is validated by simultaneous patch clamp recording from the same cell. Nanocrown electrodes enable prolonged iAP recording for continual monitoring of the same cells upon the sequential addition of four to five incremental drug doses. In this way, the dose-response data is self-referencing, which avoids the cell-to-cell variations inherent to hPSC-CMs. We are hopeful that this technology development is a step towards establishing an iAP screening assay for preclinical evaluation of drug-induced arrhythmogenicity. Ecological interactions in breast cancer: Cell facilitation promotes growth and survival under drug pressure The interplay of positive and negative interactions between drug sensitive and resistant cells determines the efficacy of treatment in heterogeneous cancer cell populations. We developed isogenic estrogen receptor positive (ER+) breast cancer cell lineages sensitive and resistant to ribociclib-induced CDK4/6 inhibition to investigate cell interactions under selective pressure. During treatment with ribociclib alone or in combination with endocrine therapies, sensitive cells survive and proliferate better when grown together with resistant cells than when grown in monoculture. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) assays detected estrone, a less active metabolite, in medium taken from each lineage but detected highly active estradiol only from resistant cells. Because sensitive cells are more dependent on estrogen signaling than resistant cells, they benefit from the increase in local estradiol from resistant cells, leading to facilitated growth when grown under treatment in coculture. We developed a suite of mathematical models to quantify the strength of competition and facilitation between cell lines and evaluate the effects of treatments that modify these interactions. These analyses support the conclusion that resistant cell lineages facilitate growth of sensitive lineages during therapy through the increased production and sharing of local estradiol and themselves suffer from increased competition and suppression as a result. These models show that blocking facilitation has the potential to control both resistant and sensitive cell populations and inhibit the emergence of a refractory population. This combination of experiments and models shows that interactions within heterogeneous cancer cell populations can combine positive and negative interactions to shape the dynamics of resistance. In vivo analysis of the contribution of proprotein convertases to the processing of FGF23 Fibroblast growth factor 23 (FGF23) is a hormone secreted from fully differentiated osteoblasts and osteocytes that inhibits phosphate reabsorption by kidney proximal tubules. The full-length (i.e., intact) protein mediates FGF23 endocrine functions, while endoproteolytic cleavage at a consensus cleavage sequence for the proprotein convertases (PCs) inactivates FGF23. Two PCs, furin and PC5, were shown to cleave FGF23 in vitro at RHTR179{downarrow}, but whether they are fulfilling this function in vivo is currently unknown. To address this question we used here mice lacking either or both furin and PC5 in cell-specific manners and mice lacking the paired basic amino acid-cleaving enzyme 4 (PACE4) in all cells. Our analysis shows that furin inactivation in osteoblasts and osteocytes results in a 25% increase in circulating intact FGF23, without any significant impact on serum phosphate levels, whether mice are maintained on a normal or a low phosphate diet. Under conditions of iron deficiency, FGF23 is normally processed in control mice, but its processing is impaired in mice lacking furin in osteoblasts and osteocytes. In contrast, FGF23 is normally cleaved following erythropoietin or IL-1{beta} injections in mice lacking furin or both furin and PC5, and in PACE4-deficient mice. Altogether, these studies suggest that furin is only partially responsible for FGF23 cleavage under certain conditions in vivo. The processing of FGF23 may therefore involve the redundant action of multiple PCs or of other peptidases in osteoblasts, osteocytes and hematopoietic cells. Non-invasive imaging of gene expression and protein secretion dynamics in living mice: identification of ectopic prothrombin expression as driver of thrombosis in cancer The topology of gene expression and protein localization is a crucial characteristic of life, where the spatiotemporal dynamic of secretory proteins instruct higher order organization, including the orchestration of developmental and adaptive programs. However tools to non-invasively interrogate the fate of secretory proteins in vivo are scarce. Here we introduce a genetic tagging strategy for in vivo imaging of the secretion and expression dynamics of secretory proteins in living animals. Applying this to a prototypical liver-derived secretory protein, we demonstrate that this approach, combined with optical in-vivo imaging, uncovers extrahepatic prothrombin expression in multiple novel anatomical sites (including testes, placenta, brain, kidney, heart and lymphatic system) and in emerging tumors, resulting in significant amounts of tumor-derived prothrombin in the blood with procoagulant properties. Syngeneic cell lines from this mouse model enable unravelling regulatory mechanisms in high resolution, and in a scalable format ex vivo. Beyond discovering new functions of proteins in a targeted manner, this model allows identifying rheostats in the cross-talk between gene expression and availability of a secretory protein. It is also a valuable resource for uncovering novel (tissue-specific) therapeutic vulnerabilities. Correlational networking guides the discovery of cryptic natural product biosynthetic enzymes Bacterial natural product biosynthetic genes, canonically clustered, have been increasingly found to rely on cryptic enzymes encoded elsewhere in the genome for completion of biosynthesis. The study and application of lanthipeptides are frequently hindered by cryptic proteases required for final maturation. Here, we establish a global correlation network bridging the gap between lanthipeptide precursors and cryptic proteases. Applying our analysis to 161,954 bacterial genomes, we establish 5,209 correlations between precursors and cryptic proteases, with 91 prioritized. We use network predictions and co-expression analysis to reveal a previously missing protease for the maturation of class I lanthipeptide paenilan. We further discover widely distributed bacterial M16B metallopeptidases of previously unclear biological function as a new family of lanthipeptide proteases. We show the involvement of a pair of bifunctional M16B proteases in the production of novel class III lanthipeptides with high substrate specificity. Together, these results demonstrate the strength of our correlational networking approach to the discovery of cryptic lanthipeptide proteases and potentially other missing enzymes for natural products biosynthesis.

Graphical abstract O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY 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. A Bayesian method to cluster single-cell RNA sequencing data using Copy Number Alterations MotivationCancers are composed by several heterogeneous subpopulations, each one harbouring different genetic and epigenetic somatic alterations that contribute to disease onset and therapy response. In recent years, copy number alterations leading to tumour aneuploidy have been identified as potential key drivers of such populations, but the definition of the precise makeup of cancer subclones from sequencing assays remains challenging. In the end, little is known about the mapping between complex copy number alterations and their effect on cancer phenotypes.

ResultsWe introduce CONGAS, a Bayesian probabilistic method to phase bulk DNA and single-cell RNA measurements from independent assays. CONGAS jointly identifies clusters of single cells with subclonal copy number alterations, and differences in RNA expression. The model builds statistical priors leveraging bulk DNA sequencing data, does not require a normal reference and scales fast thanks to a GPU backend and variational inference. We test CONGAS on both simulated and real data, and find that it can determine the tumour subclonal composition at the single-cell level together with clone-specific RNA phenotypes in tumour data generated from both 10x and Smart-Seq assays.

AvailabilityCONGAS is available as 2 packages: CONGAS (https://github.com/caravagnalab/congas), which implements the model in Python, and RCONGAS (https://caravagnalab.github.io/rcongas/), which provides R functions to process inputs, outputs, and run CONGAS fits. The analysis of real data and scripts to generate figures of this paper are available via RCONGAS; code associated to simulations is available at https://github.com/caravagnalab/rcongas_test.

Contactgcaravagna@units.it

Supplementary informationSupplementary data are available at Bioinformatics online. 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. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. The Hippocampus uses cryptographic principles rather than plasticity to enable animals to secretly cache and retrieve their food For animals, the ability to hide and retrieve valuable information, such as the location of food, can mean the difference between life and death. Here, we propose that to achieve this, their brain uses spatial cells similarly to how we utilize encryption for data security. Some animals are able to cache hundreds of thousands of food items annually by each individual and later retrieve most of what they themselves stashed. Rather than memorizing their cache locations as previously suggested, we propose that they use a single cryptographic-like mechanism during both caching and retrieval. The model we developed is based on hippocampal spatial cells, which respond to an animals positional attention, such as when the animal enters a specific region (place-cells) or gazes at a particular location (spatial-view-cells). We know that the region that activates each spatial cell remains consistent across subsequent visits to the same area but not between areas. This remapping, combined with the uniqueness of cognitive maps, produces a persistent crypto-hash function for both food caching and retrieval. We also show that the model stores temporal information that helps animals in food caching order preference as previously observed. This behavior, which we refer to as crypto-taxis, might also explain consistent differences in decision-making when animals are faced with a large number of alternatives such as in foraging. 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. Functional binding dynamics relevant to the evolution of zoonotic spillovers in endemic and emergent Betacoronavirus strains Comparative functional analysis of the dynamic interactions between various Betacoronavirus mutant strains and broadly utilized target proteins such as ACE2 and CD26, is crucial for a more complete understanding of zoonotic spillovers of viruses that cause diseases such as COVID-19. Here, we employ machine learning to replicated sets of nanosecond scale GPU accelerated molecular dynamics simulations to statistically compare and classify atom motions of these target proteins in both the presence and absence of different endemic and emergent strains of the viral receptor binding domain (RBD) of the S spike glycoprotein. Machine learning was used to identify functional binding dynamics that are evolutionarily conserved from bat CoV-HKU4 to human endemic/emergent strains. Conserved dynamics regions of ACE2 involve both the N-terminal helices, as well as a region of more transient dynamics encompassing K353, Q325 and a novel motif AAQPFLL 386-92 that appears to coordinate their dynamic interactions with the viral RBD at N501. We also demonstrate that the functional evolution of Betacoronavirus zoonotic spillovers involving ACE2 interaction dynamics are likely pre-adapted from two precise and stable binding sites involving the viral bat progenitor strains interaction with CD26 at SAMLI 291-5 and SS 333-334. Our analyses further indicate that the human endemic strains hCoV-HKU1 and hCoV-OC43 have evolved more stable N-terminal helix interactions through enhancement of an interfacing loop region on the viral RBD, whereas the highly transmissible SARS-CoV-2 variants (B.1.1.7, B.1.351 and P.1) have evolved more stable viral binding via more focused interactions between the viral N501 and ACE2 K353 alone. A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Exploring the natural origins of SARS-CoV-2 in the light of recombination The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (Sarbecoviruses) have been identified, is fuelling speculation on the natural origins of SARS-CoV-2. We have comprehensively analysed phylogenetic relations between SARS-CoV-2, and the related bat and pangolin Sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Clustering of the inferred recombination events is non-random with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift in the ancestry of bat Sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm horseshoe bats, Rhinolophus, are the likely SARS-CoV-2 progenitor reservoir species. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and co-circulation of these viruses ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years or so. We confirm that a direct proximal ancestor to SARS-CoV-2 is yet to be sampled, since the closest relative shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for more wildlife sampling to (i) pinpoint the exact origins of SARS-CoV-2s animal progenitor, and (ii) survey the extent of the diversity in the related Sarbecoviruses phylogeny that present high risk for future spillover.

HighlightsO_LIThe origin of SARS-CoV-2 can be traced to horseshoe bats, genus Rhinolophus, with ranges in both China and Southeast Asia.
C_LIO_LIThe closest known relatives of SARS-CoV-2 exhibit frequent transmission among their Rhinolophus host species.
C_LIO_LISarbecoviruses have undergone extensive recombination throughout their evolutionary history.
C_LIO_LIAccounting for the mosaic patterns of these recombinants is important when inferring relatedness to SARS-CoV-2.
C_LIO_LIBreakpoint patterns are consistent with recombination hotspots in the coronavirus genome, particularly upstream of the pike open reading frame with a coldspot in S1.
C_LI Differed cell division angle, position of cell proliferative area, and localization of ANGUSTIFOLIA3 in lateral organs Leaf meristem is a cell proliferative zone present in the lateral organ primordia. In this study, we investigated how the proliferative zone affects the final morphology of the lateral organs. We examined how cell proliferative zones differ in the primordia of planar floral organs and polar auxin transport inhibitor (PATI)-treated leaves from normal foliage leaf primordia of Arabidopsis thaliana with a focus on the spatial accumulation pattern of ANGUSTIFOLIA3 (AN3), a key element for leaf meristem positioning. We found that organ shape changes by PATI treatment were correlated to the angle of the cell division plane relative to the leaf primordia axis in the leaf meristem (cell division angle), but not with leaf-meristem positioning, size of the leaf meristem, or the localization pattern of AN3 protein. In contrast, different shapes between sepals and petals compared with foliage leaves were associated with both altered meristem position associated with altered AN3 expression patterns and different distributions of cell division angles. These results suggest that lateral organ shapes are regulated via two aspects: position of meristem and cell division angles

Summary statementDifferent lateral organs with different morphology possess different properties of meristems; cell division angles, position of cell proliferative area and AN3 localization patterns. Sperm Histone H3 Lysine 4 tri-methylation serves as a metabolic sensor of paternal obesity and is associated with the inheritance of metabolic dysfunction ObjectiveParental environmental exposures can strongly influence descendant risks for adult disease. How paternal obesity changes the sperm chromatin leading to the acquisition of metabolic disease in offspring remains controversial and ill-defined. The objective of this study was to assess: (1) whether obesity induced by a high-fat diet alters sperm histone methylation; (2) whether paternal obesity can induce metabolic disturbances across generations; (3) whether there could be cumulative damage to the sperm epigenome leading to enhanced metabolic dysfunction in descendants; and (4) whether obesity-sensitive regions associate with embryonic epigenetic and transcriptomic profiles. Using a genetic mouse model of epigenetic inheritance, we investigated the role of histone H3 lysine 4 methylation (H3K4me3) in the paternal transmission of metabolic dysfunction. This transgenic mouse overexpresses the histone demethylase enzyme KDM1A in the developing germline and has an altered sperm epigenome at the level of histone H3K4 methylation. We hypothesized that challenging transgenic sires with a high-fat diet would further erode the sperm epigenome and lead to enhanced metabolic disturbances in the next generations.

MethodsTo assess whether paternal obesity can have inter- or transgenerational impacts, and if so, to identify potential mechanisms of this non-genetic inheritance, we used wildtype C57BL/6NCrl and transgenic males with a pre-existing altered sperm epigenome. To induce obesity, sires were fed either a control or high-fat diet (10% or 60% kcal fat, respectively) for 10-12 weeks, then bred to wildtype C57BL/6NCrl female fed a regular diet. F1 and F2 descendants were characterized for metabolic phenotypes by examining the effects of paternal obesity by sex, on body weight, fat mass distribution, the liver transcriptome, intraperitoneal glucose and insulin tolerance tests. To determine whether obesity altered the F0 sperm chromatin, native chromatin immunoprecipitation-sequencing targeting H3K4me3 was performed. To gain insight into mechanisms of paternal transmission, we compared our sperm H3K4me3 profiles with embryonic and placental chromatin states, histone modification and gene expression profiles.

ResultsObesity-induced alterations in H3K4me3 occurred at genes implicated in metabolic, inflammatory, and developmental processes. These processes were associated with offspring metabolic dysfunction and corresponded to genes enriched for H3K4me3 in embryos, and overlapped embryonic and placenta gene expression profiles. Transgenerational susceptibility to metabolic disease was only observed when obese F0 had a pre-existing modified sperm epigenome. This coincided with increased H3K4me3 alterations in sperm and more severe phenotypes affecting their offspring.

ConclusionsOur data suggest sperm H3K4me3 might serve as a metabolic sensor that connects paternal diet with offspring phenotypes via the placenta. This non-DNA based knowledge of inheritance has the potential to improve our understanding of how environment shapes heritability and may lead to novel routes for the prevention of disease. This study highlights the need to further study the connection between the sperm epigenome, placental development and childrens health. PRSS2 stimulates tumor growth by remodeling the TME via repression of Tsp1 In the earliest stages of tumor development, epithelial tumors (carcinomas) are physically confined to the area of the tissue in which they form. These nascent lesions (carcinomas in situ) are sequestered from the tissue parenchyma by the basement membrane. Within the tissue parenchyma lie a myriad of cell types comprised of fibroblasts, immune and inflammatory cells and endothelial cells. Upon invasion across the basement membrane and into the tissue parenchyma, tumors must manipulate the expression of pro- and anti-tumorigenic proteins such that pro-tumorigenic factors are produced in excess to anti-tumorigenic proteins. One such anti-tumorigenic protein is Thrombospondin-1 (Tsp-1). We have previously demonstrated that stimulation of Tsp-1 in the tumor microenvironment (TME) potently inhibits tumor growth and progression and in some cases induces tumor regression. Here, we identify a novel tumor-mediated mechanism to repress the expression of Tsp-1 in the TME via secretion of the serine protease PRSS2. We demonstrate that PRSS2 represses Tsp-1, not via its enzymatic activity, but by binding to low-density lipoprotein receptor-related protein 1 (LRP1). These findings describe a novel activity for PRSS2 through binding to LRP1 and represent a potential therapeutic strategy to treat cancer by blocking the PRSS2-mediated repression of Tsp-1. 'rbioacc': an R-package to analyse toxicokinetic data O_LI rbioacc is an  R-package dedicated to the analysis of experimental data collected from bioaccumulation tests during which organisms are exposed to a chemical (exposure phase) and then put into a clean media (depuration phase). Internal concentrations are regularly measured over time all along the experiment.
C_LIO_LI rbioacc provides ready-to-use functions to visualize and fully analyze such data. Under a Bayesian framework, this package fits a generic one-compartment toxicokinetic (TK) model automatically built from the data. It provides TK parameter estimates (appropriate uptake and elimination rates) and bioaccumulation metrics (e.g., BCF, BSAF, BMF). All parameter estimates, bioaccumulation metrics as well as predictions of internal concentrations into organisms are delivered with their uncertainty.
C_LIO_LIBioaccumulation metrics are primarily provided in support of environmental risk assessment, in full compliance with regulatory requirements required to approve marketing applications of chemical substances.
C_LIO_LIThis paper gives brief worked examples of the use of  rbioacc from data collected through standard bioaccumulation tests, and publicly available within the scientific literature. These examples constitute step-by-step user-guides to analyze any new data set, uploaded in the right format.
C_LI