On-person adaptive evolution of Staphylococcus aureus during treatment for atopic dermatitis Genetic variation among bacterial strains can contribute to heterogeneity in the severity of chronic inflammatory diseases 1,2, but the degree of variation created by de novo mutation during colonization is not well understood. The inflamed skin of people with atopic dermatitis (AD) is frequently colonized with Staphylococcus aureus, an opportunistic pathogen associated with both asymptomatic colonization of nasal passages and invasive disease3-6. While genetic risk and barrier disruption are critical to AD initiation 7,8, S. aureus colonization worsens disease severity by promoting skin damage9. Here we show, from tracking 25 children treated for AD over 9 months, that S. aureus adapts via de novo mutations during colonization. Patients S. aureus populations are typically dominated by a single lineage, with infrequent invasion by distant genotypes. Variants emerge within each lineage with mutation rates similar to S. aureus in other contexts. Some of these variants replace their ancestors across the body within months, with signatures of adaptive, rather than neutral, forces. Most strikingly, the capsule synthesis gene capD obtained four parallel mutations within one patient and was involved in mutational sweeps in multiple patients. We confirm that selection for capD negativity is common in AD, but not in other contexts, via reanalysis of public S. aureus genomes from 276 people. Our finding of disease-specific selection raises the possibility that adaptation of pathobionts during colonization prolongs the positive feedback cycle of inflammation.
Atypical integrative element with strand-biased circularization activity assists interspecies antimicrobial resistance gene transfer from Vibrio alfacsensis The exchange of antimicrobial resistance (AMR) genes between aquaculture and terrestrial microbial populations has emerged as a serious public health concern. However, the nature of the mobile genetic elements in marine bacteria is poorly documented. To gain insight into the genetic mechanisms underlying AMR gene transfer from marine bacteria, we mated a multidrug-resistant Vibrio alfacsensis strain with an Escherichia coli strain, and then determined the complete genome sequences of the donor and transconjugant strains. Sequence analysis revealed a multidrug resistance conjugative plasmid in the donor strain, which was integrated into the chromosome of the recipient. The plasmid backbone in the transconjugant chromosome was flanked by two copies of a 7.1 kb unclassifiable integrative element harboring a {beta}-lactamase gene. The 7.1 kb element (named SE-6945) and previously reported Tn6283 share four coding sequences, two of which encode the catalytic R-H-R-Y motif of tyrosine recombinases. Polymerase chain reaction and sequencing experiments revealed that SE-6945 generates a circular copy of its one specific strand without generating an empty donor site, in contrast to the movement of integron gene cassette or ICE/IMEs discovered to date. Copy numbers and location of SE-6945 in the chromosome affected antibiotic resistance levels of the transconjugants. SE-6945/Tn6283-related integrative elements are collectively termed SEs (strand-biased circularizing integrative elements). SEs were identified in the genomes of other Vibrio species. Overall, these results suggest that in addition to conjugative plasmids, SEs are involved in the spread of AMR genes among marine bacteria.
Agrobacterium mediated transformation and deciphering SNPs in TaLr67 gene homeologs for gene editing in wheat. Diseases adversely affect grain yield of crop plants. Leaf rust is a major disease of wheat. As race-specific resistance breaks down, introduction of newer sources of resistance often from older accessions is necessary. Linkage drag from the donor accessions adversely impacts grain yield. As a result, CIMMYT breeding effort has shifted to using durable resistance, also referred to as race non-specific host resistance, which allows slow rusting but maintains grain yield. One of the key genes for durable resistance is Lr67 and the resistant form of this gene is absent in CIMMYT elite lines. Hence, we have initiated efforts to convert the susceptible copy of Lr67 into its resistant form directly in these elite lines using gene editing. This would eliminate backcrossing and thus save time as well as eliminate linkage drag that would accompany the resistant copy of the gene if it were to be introgressed from an older accession. As first steps, we have isolated and sequenced the genomic copies from each of the A, B, and D genome of the Lr67 gene from three elite lines and an experimental line. Identification of more than 50 single nucleotide polymorphisms (SNPs) in the open reading frames (ORF) among these lines would be useful in designing the guide RNA molecules with precision. Further, we have streamlined genetic transformation of these elite lines, a prerequisite step for gene editing. Life history adaptations to fluctuating environments: Combined effects of demographic buffering and lability of demographic parameters Demographic buffering and lability have been identified as adaptive strategies to optimise fitness in a fluctuating environment. These are not mutually exclusive, however we lack efficient methods to measure their relative importance for a given life history. Here, we decompose the stochastic growth rate (fitness) into components arising from nonlinear responses and variance-covariance of demographic parameters to an environmental driver, which allows studying joint effects of buffering and lability. We apply this decomposition for 154 animal matrix population models under different scenarios, to explore how these main fitness components vary across life histories. Faster-living species appear more responsive to environmental fluctuations, either positively or negatively. They have the highest potential for strong adaptive demographic lability, while demographic buffering is a main strategy in slow-living species. Our decomposition provides a comprehensive framework to study how organisms adapt to variability through buffering and lability, and to predict species responses to climate change. Relationship between shear wave velocity and muscle activation is inconsistent across different muscle types There is an increasing use of shear wave ultrasound elastography to quantify mechanical properties of muscles under various conditions such as changes muscle length and levels of activation in healthy and pathological muscle. However, little is known about the variability in shear wave velocity among muscles as most studies investigate one specific muscle. The purpose of this study was to determine if the relationship between SWV and muscle activation is consistent across muscles with different architectures: biceps brachii, tibialis anterior, and medial gastrocnemius, All measures were made at matching levels of activation and approximately at the optimal length for each muscle to control for length-dependent changes in the relationship between activation and force or stiffness. We also conducted a control experiment to determine how the passive force within a muscle alters the relationship between muscle activation and shear wave velocity. The relationship between shear wave velocity-squared and activation above 10% MVC differed across muscles, with biceps brachii and medial gastrocnemius showing a lower slope than tibialis anterior. Shear wave velocity-squared also differed between muscles at the shortest length (p<0.001) and the increase in shear wave velocity-squared with muscle lengthening also differed between muscle types (p = 0.005) Muscle-specific differences could not be explained by the architectural features such as pennation angle, fiber length, and physiological cross-sectional area. Our results demonstrate that there is not a unique relationship between muscle activation and shear wave velocity, highlighting the importance of understanding the many factors contributing to shear wave propagation in muscle before elastography can be used to make quantitative comparisons across muscle types. Polycystins are upregulated in fibrotic kidneys and promote renal fibrosis Autosomal polycystic kidney disease (ADPKD) is caused by mutations in PKD1or PKD2 gene, which encodes polycystin-1 or polycystin-2 respectively. Paradoxically, polycystins are up-regulated in human ADPKD kidneys and overexpression of polycystin-1 even induced renal cyst formation in mice, suggesting that up-regulation of polycystin is pathogenic. In this study, we determined the role of polycystin in renal fibrosis and explored its underlying mechanisms. We found that the expression of polycystin-1 or polycystin-2 was positively correlated with the expression of collagen-I in the kidney of three different fibrotic mouse models. Pharmaceutical inhibition or genetic knockout of polycystin-2 reduced the expression of pro-fibrotic markers in fibrotic kidneys of two animal models. Moreover, conditional knockout of Pkd1 gene attenuated renal fibrosis in two mouse models. Enhancer of Zeste homolog 2 (EZH2) is a methyltransferase which plays important role in renal fibrosis. We found that EZH2 was positively correlated with the expression of polycystins in fibrotic kidneys. Conditional knockout of the Ezh2 attenuated the anti-fibrotic responses induced by Pkd1 deletion in fibrotic kidneys.

In conclusion, polycystins are up-regulated in injured kidneys and polycystins promote renal fibrosis through up-regulation of EZH2 in fibrotic kidneys, suggesting that renal fibrosis in polycystic kidneys is not the direct effect of PKD gene mutations.

SignificanceMutations of polycystins cause autosomal polycystic kidney disease (ADPKD) which is associated with increased renal fibrosis in cystic kidneys. However, the role of polycystin in renal fibrosis remains unclear. Here, we showed that the expression of polycystins is positively correlated with the development of renal fibrosis in injured kidneys. Pharmaceutical inhibition or genetic deletion of polycystin-2 attenuated renal fibrosis in fibrotic kidneys. Moreover, genetic deletion of polycystin-1 inhibited renal fibrosis in two different animal models. Finally, we showed that ablation of the methyltransferase EZH2 mitigated the anti-fibrotic effect of polycystin-1 deletion in fibrotic kidneys. We conclude that polycystin promotes through EZH2 and renal fibrosis in cystic kidney is not the direct effect of PKD gene mutations. Enhancement of CRISPR-mediated homologous recombination-based gene knock-in efficiency by simultaneous HDR overexpression and NHEJ suppression in poplar CRISPR-mediated genome editing has become a powerful tool for genetic modification of biological traits. However, developing an efficient, site-specific, gene knock-in system based on homology-directed DNA repair (HDR) remains a significant challenge in plants, especially in woody species like poplar. Here, we show that simultaneous inhibition of non-homologous end joining (NHEJ) recombination cofactor XRCC4 and overexpression of HDR enhancer factors CtIP and MRE11 can improve the HDR efficiency for gene knock-in. Using this approach, the BleoR gene was integrated onto the 3' end of the MKK2 MAP Kinase gene to generate a BleoR-MKK2 fusion protein. Based on exogenous BleoR expression, the HDR-mediated knock-in efficiency was up to [~]40-fold greater when using a XRCC4 silencing incorporated with a combination of CtIP and MRE11 overexpression compared to no HDR enhancement or NHEJ silencing. Furthermore, this corporation of HDR enhancer overexpression and NHEJ repression also resulted in 7-fold fewer CRISPR-induced Insertions and Deletions (InDels), resulting in no functional effects on MKK2-based salt stress responses in poplar. Therefore, this approach may be useful not only in poplar and plants or crops but also in mammalians for improving CRISPR-mediated gene knock-in efficiency. Mitochondrial pyruvate supports lymphoma proliferation by fueling a non-canonical glutamine metabolism pathway The fate of pyruvate, which is modulated by the activity of the mitochondrial pyruvate carrier (MPC), is a defining metabolic feature in many cancers. Diffuse large B-cell lymphomas (DLBCLs) are a genetically and metabolically heterogeneous cancer. Although MPC expression and activity differed between DLBCL subgroups, mitochondrial pyruvate was uniformly consumed by glutamate pyruvate transaminase 2 (GPT2) to support -ketoglutarate production as part of glutaminolysis. This led us to discover that glutamine exceeds pyruvate as a carbon source for the tricarboxylic acid (TCA) cycle in DLBCLs. Furthermore, we found that MPC inhibition unexpectedly leads to decreased glutaminolysis, which is contrary to previous observations in other cell types. We also discovered that MPC inhibition and depletion only decreased DLBCL proliferation in an extracellular matrix (ECM) environment and in vivo xenografts, but not in the typical DLBCL suspension environment. We also have found that the metabolic profile of DLBCL cells in ECM is markedly different from cells in suspension environment. Thus, we report that besides the canonical glutamate dehydrogenase (GDH)-mediated glutaminolysis, the non-canonical GPT2 mediated consumption and assimilation of glutamine and pyruvate in DLBCLs enables their proliferation in an extracellular environment-dependent manner.

HIGHLIGHTSO_LIGlutamine, but not glucose, is a major carbon source for the tricarboxylic acid cycle in diffuse large B-cell lymphomas.
C_LIO_LIMitochondrial pyruvate supports glutaminolysis in diffuse large B-cell lymphomas by supplying pyruvate for glutamate pyruvate transaminase 2-mediated  -ketoglutarate production.
C_LIO_LIMitochondrial pyruvate carrier inhibition leads to decreased glutaminolysis in diffuse large B-cell lymphomas.
C_LIO_LI -ketoglutarate production is important for diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment.
C_LIO_LIMitochondrial pyruvate carrier activity supports diffuse large B-cell lymphoma proliferation in a solid extracellular matrix environment and in mouse xenografts.
C_LI Repeated subfunctionalization of a modular antimicrobial peptide gene for neural function Antimicrobial peptides (AMPs) are host-encoded antibiotics that combat invading pathogens. These genes commonly encode multiple products as post-translationally cleaved polypeptides. Recent studies have highlighted roles for AMPs in neurological contexts suggesting functions for these defence molecules beyond infection. During our immune study characterizing the antimicrobial peptide gene Baramicin, we recovered multiple Baramicin paralogs in Drosophila melanogaster and other species, united by their N-terminal IM24 domain. Not all paralogs were immune-induced. Here, through careful dissection of the Baramicin familys evolutionary history, we find that these non-immune paralogs result from repeated events of duplication and subsequent truncation of the coding sequence from an immune-inducible ancestor. These truncations leave only the IM24 domain as the prominent gene product. Surprisingly, using mutation and targeted gene silencing we demonstrate that two such genes are adapted for function in neural contexts in D. melanogaster. We also show enrichment in the head for independent Baramicin genes in other species. The Baramicin evolutionary history reveals that the IM24 Baramicin domain is not strictly useful in an immune context. We thus provide a case study for how an AMP-encoding gene might play dual roles in both immune and non-immune processes via its multiple peptide products. We reflect on these findings to highlight a blind spot in the way researchers approach AMP research in in vivo contexts.

Significance statementAntimicrobial peptides are immune proteins recently implicated in neurological roles. To date little attention has been paid to the contributions of different gene products in this function. Here we show that an antimicrobial peptide gene encodes multiple products with either immune-specific or neurological roles. A phenomics approach for in vitro antiviral drug discovery BackgroundThe current COVID-19 pandemic has highlighted the need for new and fast methods to identify novel or repurposed therapeutic drugs. Here we present a method for untargeted phenotypic drug screening of virus-infected cells, combining Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties.

ResultsWe show that the methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with Human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state.

ConclusionsThe method can be used in drug discovery for morphological profiling of novel antiviral compounds on both infected and non-infected cells. Human cells contain myriad excised linear introns with potential functions in gene regulation and as RNA biomarkers By using TGIRT-seq, we identified >8,500 short full-length excised linear intron (FLEXI) RNAs in human cells. Subsets of FLEXIs accumulated in a cell-type specific manner, and [~]200 corresponded to agotrons or mirtrons or encoded snoRNAs. Analysis of CLIP-seq datasets identified potential interactions between FLEXIs and >100 different RNA-binding proteins (RBPs), 53 of which had binding sites in [&ge;]30 different FLEXIs. In addition to proteins that function in RNA splicing, these 53 RBPs included transcription factors, chromatin remodeling proteins, and cellular growth regulators that impacted FLEXI host gene alternative splicing and/or mRNA levels in knockdown datasets. We computationally identified six groups of RBPs whose binding sites were enriched in different subsets of FLEXIs: AGO1-4 and DICER associated with agotrons and mirtrons; AATF, DKC1, NOLCI, and SMNDC1 associated with snoRNA-encoding FLEXIs; two different combinations of alternative splicing factors found in stress granules; and two novel RBP-intron combinations, one including LARP4 and PABC4, which function together in the cytoplasm to regulate ribosomal protein translation. Our results suggest a model in which proteins involved in transcriptional regulation, alternative splicing, or post-splicing secondary functions bind and stabilize cell-type specific subsets of FLEXIs that perform different biological functions and have potential utility as biomarkers. 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. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Resolving when (and where) the Thylacine went extinct Like the Dodo and Passenger Pigeon before it, the predatory marsupial Thylacine (Thylacinus cynocephalus), or  Tasmanian tiger, has become an iconic symbol of human-caused extinction. The last captive animal died in 1936, but even today reports of the Thylacines possible ongoing survival in remote regions of Tasmania are newsworthy and capture the publics imagination. Extirpated from mainland Australia in the mid-Holocene, the large island of Tasmania became the species final stronghold. Following European settlement in the 1800s, the Thylacine was heavily persecuted and pushed to the margins of its range, although many sightings were reported thereafter--even well beyond the 1930s. To gain a new depth of insight into the extinction of the Thylacine, we assembled an exhaustive database of 1,237 observational records from Tasmania (from 1910 onwards), quantified their uncertainty, and charted the patterns these revealed. We also developed a new method to visualize the species 20th-century spatio-temporal dynamics, to map potential post-bounty refugia and pinpoint the most-likely location of the final persisting subpopulation. A direct reading of the high-quality records (confirmed kills and captures, in combination with sightings by past Thylacine hunters and trappers, wildlife professionals and experienced bushmen) implies a most-likely extinction date within four decades following the last capture (i.e., 1940s to 1970s). However, uncertainty modelling of the entire sighting record, where each observation is assigned a probability and the whole dataset is then subject to a sensitivity analysis, suggests that extinction might have been as recent as the late 1980s to early 2000s, with a small chance of persistence in the remote south-western wilderness areas. Beyond the intrinsically fascinating problem of reconstructing the final fate of the Thylacine, the new spatio-temporal mapping of extirpation developed herein would also be useful for conservation prioritization and search efforts for other rare taxa of uncertain status. Placental uptake and metabolism of 25(OH)Vitamin D determines its activity within the fetoplacental unit Pregnancy 25-hydroxyvitamin D (25(OH)D) concentrations are associated with maternal and fetal health outcomes. Using physiological human placental perfusion and villous explants, we investigate the role of the placenta in regulating the relationships between maternal 25(OH)D and fetal physiology. We demonstrate active placental uptake of 25(OH)D3 by endocytosis, placental metabolism of 25(OH)D3 into 24,25-dihydroxyvitamin D3 and active 1,25-dihydroxyvitamin D [1,25(OH)2D3], with subsequent release of these metabolites into both the maternal and fetal circulations. Active placental transport of 25(OH)D3 and synthesis of 1,25(OH)2D3 demonstrate that fetal supply is dependent on placental function rather than simply the availability of maternal 25(OH)D3. We demonstrate that 25(OH)D3 exposure induces rapid effects on the placental transcriptome and proteome. These map to multiple pathways central to placental function and thereby fetal development, independent of vitamin D transfer. Our data suggest that the underlying epigenetic landscape helps dictate the transcriptional response to vitamin D treatment. This is the first quantitative study demonstrating vitamin D transfer and metabolism by the human placenta, with widespread effects on the placenta itself. These data demonstrate a complex interplay between vitamin D and the placenta and will inform future interventions using vitamin D to support fetal development and maternal adaptations to pregnancy. Novel MOG analogues to explore the MCT2 pharmacophore, α-ketoglutarate biology and cellular effects of N-oxalylglycine -ketoglutarate (KG) is a central metabolic node with a broad influence on cellular physiology. The KG analogue N-oxalylglycine (NOG) and its membrane-permeable pro-drug derivative dimethyl-oxalylglycine (DMOG) have been extensively used as tools to study prolyl hydroxylases (PHDs) and other KG-dependent processes. In cell culture media, DMOG is rapidly converted to MOG, which enters cells through monocarboxylate transporter MCT2, leading to intracellular NOG concentrations that are sufficiently high to inhibit glutaminolysis enzymes and cause cytotoxicity. Therefore, the degree of (D)MOG instability together with MCT2 expression levels determine the intracellular targets NOG engages with and, ultimately, its effects on cell viability. Here we designed and characterised a series of MOG analogues with the aims of improving stability and exploring the functional requirements for interaction with MCT2, a relatively understudied member of the SLC16 family. We report MOG analogues that maintain ability to enter cells via MCT2, and identify compounds that do not inhibit glutaminolysis or cause cytotoxicity but can still inhibit PHDs. We use these analogues to show that glutaminolysis-induced activation of mTORC1 can be uncoupled from PHD activity. Therefore, these new compounds can help deconvolute cellular effects that result from the polypharmacological action of NOG. Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants BackgroundGenome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

ResultsBy comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and produce datasets similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jack-knifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

ConclusionsGOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data hold promise for generating novel hypotheses about comparative gene function and traits. Conformational dynamics and allosteric modulation of the SARS-CoV-2 spike Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects cells through binding to angiotensin-converting enzyme 2 (ACE2). This interaction is mediated by the receptor-binding domain (RBD) of the viral spike (S) glycoprotein. Structural and dynamic data have shown that S can adopt multiple conformations, which controls the exposure of the ACE2-binding site in the RBD. Here, using single-molecule Forster resonance energy transfer (smFRET) imaging we report the effects of ACE2 and antibody binding on the conformational dynamics of S from the Wuhan-1 strain and the B.1 variant (D614G). We find that D614G modulates the energetics of the RBD position in a manner similar to ACE2 binding. We also find that antibodies that target diverse epitopes, including those distal to the RBD, stabilize the RBD in a position competent for ACE2 binding. Parallel solution-based binding experiments using fluorescence correlation spectroscopy (FCS) indicate antibody-mediated enhancement of ACE2 binding. These findings inform on novel strategies for therapeutic antibody cocktails. Probabilistic solvers enable a straight-forward exploration of numerical uncertainty in neuroscience models Understanding neural computation on the mechanistic level requires models of neurons and neuronal networks. To analyze such models one typically has to solve coupled ordinary differential equations (ODEs), which describe the dynamics of the underlying neural system. These ODEs are solved numerically with deterministic ODE solvers that yield single solutions with either no, or only a global scalar bound on precision. It can therefore be challenging to estimate the effect of numerical uncertainty on quantities of interest, such as spike-times and the number of spikes. To overcome this problem, we propose to use recently developed sampling-based probabilistic solvers, which are able to quantify such numerical uncertainties. They neither require detailed insights into the kinetics of the models, nor are they difficult to implement. We show that numerical uncertainty can affect the outcome of typical neuroscience simulations, e.g. jittering spikes by milliseconds or even adding or removing individual spikes from simulations altogether, and demonstrate that probabilistic solvers reveal these numerical uncertainties with only moderate computational overhead. A common SEC23B missense mutation in congenital dyserythropoietic anemia type II leads to growth restriction and chronic pancreatitis in mice Human loss-of-function mutations in SEC23B result in congenital dyserythropoietic anemia type II (CDAII). Complete deficiency of SEC23B in mice leads to perinatal death caused by massive degeneration of professional secretory tissues. Functions of SEC23B in postnatal mice are unclear. In this study, we generated knockout (KO) mice with deletion of exons 5 and 6 of Sec23b and knockin (KI) mice with the E109K mutation, the most common CDAII missense mutation. Despite decreased SEC23B protein level, Sec23bki/ki mice showed no obvious abnormalities. Hemizygous (Sec23bki/ko) mice exhibit a partial lethal phenotype, with only half of expected mice survive past weaning. Sec23bki/ko pancreas had exocrine insufficiency and chronic pancreatitis histology, as well as increased ER stress and apoptosis. Moreover, Sec23bki/ko mice exhibited severe growth restriction accompanied by growth hormone (GH) insensitivity, reminiscent of the Laron syndrome. Growth restriction is not associated with hepatocyte-specific Sec23b deletion, suggesting non-liver origin of the phenotype. Inflammation associated with chronic pancreatic deficiency may explain GH insensitivity in Sec23bki/ko mice. Our results indicate that phenotype severities are linked to the residual functions of SEC23B, demonstrating the importance of SEC23B in pancreatic acinar function in adult mice. The Sec23bki/ko mice provide a novel model of chronic pancreatitis and growth retardation. A biochemically-realisable relational model of the self-manufacturing cell As shown by Hofmeyr, the processes in the living cell can be divided into three classes of efficient causes that produce each other, so making the cell closed to efficient causation, the hallmark of an organism. They are the enzyme catalysts of covalent metabolic chemistry, the intracellular milieu that drives the supramolecular processes of chaperone-assisted folding and self-assembly of polypeptides and nucleic acids into functional catalysts and transporters, and the membrane transporters that maintain the intracellular milieu, in particular its electrolyte composition. Each class of efficient cause can be modelled as a relational diagram in the form of a mapping in graph-theoretic form, and a minimal model of a self-manufacturing system that is closed to efficient causation can be constructed from these three mappings using the formalism of relational biology. This Fabrication-Assembly or (F,A)-system serves as an alternative to Robert Rosens replicative Metabolism-Repair or (M,R)-system, which has been notoriously problematic to realise in terms of real biochemical processes. A key feature of the model is the explicit incorporation of formal cause, which arrests the infinite regress that plagues all relational models of the cell. The (F,A)-system is extended into a detailed relational model of the self-manufacturing cell that has a clear biochemical realisation. This (F,A) cell model allows the interpretation and visualisation of concepts such as the metabolism and repair components of Rosens (M,R)-system, John von Neumanns universal constructor, Howard Pattees symbol-function split via the symbol-folding transformation, Marcello Barbieris genotype-ribotype-phenotype ontology, and Tibor Gantis chemoton. Birth Order Differences in Education Originate in Post-Natal Environments Siblings share many environments and much of their genetics. Yet, siblings turn out different. Intelligence and education are influenced by birth order, with earlier-born siblings outperforming later-borns. We investigate whether birth order differences in education are caused by biological differences present at birth, that is, genetic differences or in-utero differences. Using family data that spans two generations, combining registry, survey, and genotype information, this study is based on the Norwegian Mother, Father and Child Cohort Study (MoBa). We show that there are no genetic differences by birth order as captured by polygenic scores (PGSs) for educational attainment. Furthermore, we show that earlier-born have lower birth weight than later-born, indicating worse uterine environments. Educational outcomes are still higher for earlier-born children when we adjust for PGSs and in utero variables, indicating that birth order differences arise post-natally. Finally, we consider potential environmental influences, such as differences according to maternal age, parental educational attainment, and sibling genetic nurture. We show that birth order differences are not biological in origin, but pinning down their specific causes remains elusive. Notch signaling via Hey1 and Id2b regulates Muller glia's regenerative response to retinal injury Zebrafish Muller glia (MG) respond to retinal injury by suppressing Notch signaling and producing progenitors for retinal repair. A certain threshold of injury-derived signal must be exceeded in order to engage MG in a regenerative response (MGs injury-response threshold). Pan-retinal Notch inhibition expands the zone of injury-responsive MG at the site of focal injury, suggesting that Notch signaling regulates MGs injury-response threshold. We found that Notch signaling enhanced chromatin accessibility and gene expression at a subset of regeneration-associated genes in the uninjured retina. Two Notch effector genes, hey1 and id2b, were identified that reflect bifurcation of the Notch signaling pathway, and differentially regulate MGs injury-response threshold and proliferation of MG-derived progenitors. Furthermore, Notch signaling component gene repression in the injured retina suggests a role for Dll4, Dlb, and Notch3 in regulating Notch signaling in MG and epistasis experiments confirm that the Dll4/Dlb-Notch3-Hey1/Id2b signaling pathway regulates MGs injury-response threshold and proliferation. 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. Bulk-assembly of monodisperse coacervates and giant unilamellar vesicles with programmable hierarchical complexity Compartmentalization is an essential step for the emergence of protocellular life on the early earth and the endeavor to bottom-up construct synthetic cells in the lab. Among manifold natural strategies for compartmentalization, assembly of lipid vesicles or coacervate droplets with resemblance to cell envelope or cytoplasm have emerged as two dominant paradigms. However, the spontaneous assembly of lipid vesicles or coacervates generally results in polydisperse microcompartments and lacks the ability to efficiently integrate functional building blocks towards defined complexity, which give rise to the difficulty to achieve the order and complexity towards the scope of that of life via spontaneous assembly. Herein, we show that the interplay of coacervates and colloidal particles provides a bulk-assembly approach to form monodisperse coacervates and giant unilamellar vesicles (GUVs) with programmable hierarchical complexity. The hierarchical complexity of the synthetic cell microcompartments can be delicately engineered from the control of the spatial or temporal organization of components in entity unit to manifold forms of interconnected synthetic cell consortia, achieving a compartment hierarchy analogous to the organelle-cell-tissue structure. The increase in complexity gives rise to emergent properties, for instance, collective morphology evolution of coacervate assemblages, gated permeability of GUVs without sophisticated protein machinery, and remarkable structural and functional stability of GUVs at extreme conditions. This work paves an unprecedented step to form monodisperse, hierarchical coacervates and GUVs via bulk-assembly, which provides a prebiotically plausible approach for the proper integration of inanimate matter towards the order and complexity of life and monodisperse, hierarchical microdroplets for follow-up scientific and technological applications. Proteomic phenotyping of stimulated Müller cells uncovers profound pro-inflammatory signaling and antigen-presenting capacity Muller cells are the main macroglial cells of the retina exerting a wealth of functions to maintain retinal homoeostasis. Upon pathological changes in the retina, they become gliotic with both protective and detrimental consequences. Accumulating data also provide evidence for a pivotal role of Muller cells in the pathogenesis of diabetic retinopathy (DR). While microglial cells, the resident immune cells of the retina are considered as main players in inflammatory processes associated with DR, the implication of activated Muller cells in chronic retinal inflammation remains to be elucidated. In order to assess the signaling capacity of Muller cells and their role in retinal inflammation, we performed in-depth proteomic analysis of Muller cell proteomes and secretomes after stimulation with INF{gamma}, TNF, IL-4, IL-6, IL-10, VEGF, TGF{beta}1, TGF{beta}2 and TGF{beta}3. We used both, primary porcine Muller cells and the human Muller cell line MIO-M1 for our hypothesis generating approach. Our results point towards an intense signaling capacity of Muller cells, which reacted in a highly discriminating manner upon treatment with different cytokines. Stimulation of Muller cells resulted in a primarily pro-inflammatory phenotype with secretion of cytokines and components of the complement system. Furthermore, we observed evidence for mitochondrial dysfunction, implying oxidative stress after treatment with the various cytokines. Finally, both MIO-M1 cells and primary porcine Muller cells showed several characteristics of atypical antigen-presenting cells, as they are capable of inducing MHC class I and MHC class II with co-stimulatory molecules. In line with this, they express proteins associated with formation and maturation of phagosomes. Thus, our findings underline the importance of Muller cell signaling in the inflamed retina, indicating an active role in chronic retinal inflammation underlying the pathogenesis of diabetic retinopathy.