Scaling effects of temperature on parasitism from individuals to host-parasite systems Parasitism is expected to change in a warmer future, but whether warming leads to substantial increases in parasitism remains unclear. Understanding how warming effects on parasitism in individual hosts (e.g., parasite load) translate to effects on population-level parasitism (e.g., prevalence, R0) remains a major knowledge gap. We analyzed the temperature dependence of parasitism at both host and population levels in thirteen empirical vector-borne host-parasite systems and found a strong, significant positive correlation between the thermal optima of individual- and population-level parasitism. We also found a significant, positive correlation in eleven environmentally-transmitted parasite systems, though several of these systems exhibited thermal optima >5{o}C apart between individual and population levels. Similarly, parasitism thermal optima were close to host performance thermal optima in vector-borne systems but not in environmentally-transmitted systems. We then adapted and simulated simple models for both transmission modes and found a similar pattern to the empirical systems: thermal optima in vector-borne systems were more strongly correlated across scales compared to environmentally-transmitted systems. Generally, our results suggest that information on the temperature-dependence, and specifically the thermal optimum, at either the individual- or population-level should provide a useful--though not quantitatively exact--baseline for predicting temperature dependence at the other level, especially in vector-borne parasite systems. Environmentally-transmitted parasitism may operate by a different set of rules, in which temperature-dependence is decoupled in some systems, requiring the need for trait-based studies of temperature-dependence at individual and population levels.
Defining coarse-grainability in a model of structured microbial ecosystems Despite their complexity, microbial ecosystems appear to be at least partially "coarse-grainable" in that some properties of interest can be adequately described by effective models of dimension much smaller than the number of interacting lineages (frequently in the dozens or hundreds). This is especially puzzling since recent studies demonstrate that a surprising amount of functionally relevant diversity is present at all levels of resolution, down to strains differing by 100 nucleotides or fewer. Rigorously defining coarse-grainability and understanding the conditions for its emergence is of critical importance for understanding microbial ecosystems. To begin addressing these questions, we propose a minimal model for investigating hierarchically structured ecosystems within the framework of resource competition. We use our model to operationally define coarse-graining quality based on reproducibility of the outcomes of a specified experiment and show that a coarse-graining can be operationally valid despite grouping together functionally diverse strains. We further show that, at least within our model, a high diversity of strains (while nominally more complex) may in fact facilitate coarse-grainability. However, this only applies if the pool of interacting strains is sampled from the appropriate "native" environment, as we discuss.
Bioinformatic and mechanistic analysis of the palmerolide PKS-NRPS biosynthetic pathway from the microbiome of an Antarctic ascidian Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the [~]75 kb trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicohabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression. Within-host evolutionary dynamics of seasonal and pandemic human influenza A viruses in young children The evolution of influenza viruses is fundamentally shaped by within-host processes. However, the within-host evolutionary dynamics of influenza viruses remain incompletely understood, in part because most studies have focused on within-host virus diversity of infections in otherwise healthy adults based on single timepoint data. Here, we analysed the within-host evolution of 82 longitudinally-sampled individuals, mostly young children, infected with A/H3N2 or A/H1N1pdm09 viruses between 2007 and 2009. For A/H1N1pdm09 infections during the 2009 pandemic, nonsynonymous changes were common early in infection but decreased or remained constant throughout infection. For A/H3N2 viruses, early infection was dominated by purifying selection. However, as infections progressed, nonsynonymous variants increased in frequencies even though within-host virus titres decreased, leading to the maintenance of virus diversity via mutation-selection balance. Our findings suggest that this maintenance of genetic diversity in these children combined with their longer duration of infection may provide important opportunities for within-host virus evolution. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. Effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI Intragenomic rearrangements of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses Variation of the betacoronavirus SARS-CoV-2 has been the bane of COVID-19 control. Documented variation includes point mutations, deletions, insertions, and recombination among closely or distantly related coronaviruses. Here, we describe yet another aspect of genome variation by beta- and alphacoronaviruses. Specifically, we report numerous genomic insertions of 5-untranslated region sequences into coding regions of SARS-CoV-2, other betacoronaviruses, and alphacoronaviruses. To our knowledge this is the first systematic description of such insertions. In many cases, these insertions change viral protein sequences and further foster genomic flexibility and viral adaptability through insertion of transcription regulatory sequences in novel positions within the genome. Among human Embecorivus betacoronaviruses, for instance, from 65% to all of the surveyed sequences in publicly available databases contain 5-UTR-derived inserted sequences. In limited instances, there is mounting evidence that these insertions alter the fundamental biological properties of mutant viruses. Intragenomic rearrangements add to our appreciation of how variants of SARS-CoV-2 and other beta- and alphacoronaviruses may arise.

SignificanceUnderstanding mechanisms of variation in coronaviruses is vital to control of their associated diseases. Beyond point mutations, insertions, deletions and recombination, we here describe for the first time intragenomic rearrangements and their relevance to changes in transmissibility, immune escape and/or virulence documented during the SARS-CoV-2 pandemic. A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Single cell profiling of functionally cured Chronic Hepatitis B patients reveals the emergence of activated innate and an altered adaptive immune response in the intra-hepatic environment Hepatitis B surface antigen (HBsAg) loss or functional cure (FC), is considered the desirable therapeutic outcome for chronic hepatitis B (CHB) patients. However, the immuno-pathological biomarkers and underlying mechanisms remain unclear. Here we present a comprehensive single cell-transcriptomic atlas together with immune-phenotyping of disease-associated cell states (DACS) isolated from intra-hepatic tissue and matched PBMCs of either CHB or FC patients. We find that the intra-hepatic environment displays specific cell identities and molecular signatures that are distinct from PBMCs. FC is associated with emergence of an altered adaptive immune response marked by CD4 cytotoxic T lymphocytes (CD4-CTLs), and an activated innate response represented by liver-resident natural killer (LR-NK) cells. Overall, these findings provide novel insights into immuno-pathological cell states associated with FC that could serve as prognostic biomarkers. Phylogenomic analyses of Alismatales shed light into adaptations to aquatic environments Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptation to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptation, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the backbone of the three main clades that could be explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and whole-genome duplications using phylogenomics. Drug screening at single-organoid resolution via bioprinting and interferometry High-throughput drug screening is an established approach to investigate tumor biology and identify therapeutic leads. Traditional platforms for high-throughput screening use two-dimensional cultures of immortalized cell lines which do not accurately reflect the biology of human tumors. More clinically relevant model systems, such as three-dimensional tumor organoids, can be difficult to screen and scale. For example, manually seeded organoids coupled to destructive endpoint assays allow for the characterization of response to treatment, but do not capture the transitory changes and intra-sample heterogeneity underlying clinically observed resistance to therapy. We therefore developed a pipeline to generate bioprinted tumor organoids linked to label-free, real-time imaging via high-speed live cell interferometry (HSLCI) and machine learning-based quantitation of individual organoids. Bioprinting cells gives rise to 3D organoid structures that preserve tumor histology and gene expression. HSLCI imaging in tandem with machine learning-based image segmentation and organoid classification tools enables accurate, label-free parallel mass measurements for thousands of bioprinted organoids. We demonstrate that our method quantitatively identifies individual organoids as insensitive, transiently sensitive, or persistently sensitive to specific treatments. This opens new avenues for rapid, actionable therapeutic selection using automated tumor organoid screening. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Nuclear alpha-synuclein is present in the human brain and is modified in dementia with Lewy bodies Dementia with Lewy bodies is pathologically defined by the cytoplasmic accumulation of alpha-synuclein within neuronal cells in the brain. Alpha-synuclein is predominately pre-synaptic, but has been reported present in various subcellular compartments in cell and animal models. In particular, nuclear alpha-synuclein is evident in-vitro and in disease models and has been associated with altered DNA integrity, gene transcription, nuclear homeostasis. However, owing to various factors, the presence of alpha-synuclein in the nuclei of human brain cells remains controversial, as does its role in synucleinopathies.

Here, we close this gap and provide a unique demonstration confirming the presence of nuclear alpha-synuclein in post-mortem brain tissue obtained from cases of dementia with Lewy bodies as well as from controls via immunohistochemistry, immunoblot, and label-free mass-spectrometry.

Discrete intra-nuclear alpha-synuclein puncta reactive against phosphorylated serine 129-alpha-synuclein and pan-alpha-synuclein antibodies were observed in cortical neurons and non-neuronal cells in fixed brain sections and in isolated nuclear preparations from Dementia with Lewy bodies cases and matched controls. Subsequent biochemical analysis of subcellular fractionated tissue confirmed alpha-synuclein as present in a nuclear fraction at levels ~ 10-fold lower than in the cytoplasm. Critically, however, an increase in monomeric nuclear alpha-synuclein phosphorylated as serine 129 was observed in cases of dementia with Lewy bodies alongside higher molecular weight pan- and phosphorylation reactive alpha-synuclein species, consistent with the formation of intranuclear phosphorylated alpha-synuclein oligomers. Furthermore, the presence of nuclear alpha-synuclein was confirmed via label free mass spectrometry, as 6 unique alpha-synuclein derived peptide sequences were identified in nuclear fractions (71.4% sequence coverage).

Collectively, our data confirm the presence of nuclear alpha-synuclein in human brain tissue and describe nuclear pathology associated with dementia with Lewy bodies. These findings address a major controversy in the synucleinopathy field by confirming the presence of nuclear alpha-synuclein in autoptic human brain tissue and, for the first time, identify that alpha-synuclein is aggregated into novel and potentially pathological assemblies in the nucleus as part of the disease process associated with dementia with Lewy bodies and thus may contribute to the disease phenotype. FGFR2 fusion protein-driven mouse models of intrahepatic cholangiocarcinoma unveil a necessary role for Erk signaling Background and aimsAbout 15% of intrahepatic cholangiocarcinoma (iCCA) express fibroblast growth factor receptor 2 (FGFR2) fusion proteins (FFs), most often in concert with mutationally inactivated TP53, CDKN2A or BAP1. FFs span residues 1-768 of FGFR2 fused to sequences encoded by any of a long list (>60) of partner genes, a configuration sufficient to ignite oncogenic FF activation. In line, FGFR-specific tyrosine kinase inhibitors (F-TKI) were shown to provide clinical benefit in FF+ iCCA, although responses were partial and/or limited by resistance mechanisms, including the FF V565F gatekeeper mutation. Herein we present an FF-driven murine iCCA model and exploit its potential for pre-clinical studies on FF therapeutic targeting.

MethodsFour iCCA FFs carrying different fusion sequences were expressed in Tp53-/- mouse liver organoids. Tumorigenic properties of genetically modified liver organoids were assessed by intrahepatic/subcutaneous transplantation in immuno-deficient mice. Cellular models derived from neoplastic lesions were exploited for pre-clinical studies.

ResultsTumors diagnosed as CCA were obtained upon transplantation of FF-expressing liver organoids. The penetrance of this tumorigenic phenotype was influenced by FF identity. Tumor organoids and 2D cell lines derived from CCA lesions were addicted to FF signaling via Ras-Erk, regardless of FF identity or presence of V565F mutation. Double blockade of FF-Ras-Erk pathway by concomitant pharmacological inhibition of FFs and Mek1/2 provided greater therapeutic efficacy than single agent F-TKI in vitro and in vivo.

ConclusionsFF-driven iCCA pathogenesis was successfully modelled in murine Tp53-/- background. This model revealed biological heterogeneity among structurally different FFs. Double blockade of FF-Erk signaling deserves consideration for improving precision-based approaches against human FF+ iCCA.

Abbreviations used in this paper: ANOVA, analysis of variance; Bap1, BRCA1-Associated-Protein 1; Cdkn2a, cyclin-dependent kinase inhibitor 2A; Cftr, cystic fibrosis transmembrane conductance regulator; Ck19, cytokeratin 19; Cyp3A, cytochrome P450, family 3, subfamily A; EGF, Epidermal growth factor; EGFR, Epidermal growth factor receptor; EpCAM, epithelial cell adhesion molecule; Erk, extracellular signal-regulated kinase; FGFR2, fibroblast growth factor receptor 2; FRS2, fibroblast growth factor receptor substrate 2; GRB2, growth factor receptor-bound 2; GSEA, gene set enrichment analysis; GSVA, gene set variation analysis; H&E, hematoxylin and eosin; HepPar1, hepatocyte Paraffin 1; Hnf4-7, hepatocyte nuclear factor 4 alpha; Hprt, hypoxanthine-guanine phosphoribosyl transferase; LGR5, leucine-rich repeat-containing G-protein coupled receptor 5; NTRK, neurotrophic Tyrosine Kinase; Parp, Poly (ADP-ribose) polymerase; RECIST, response evaluation criteria in solid tumors; SHP2, Src homology phosphatase 2; Ttr, transthyretin. The relative speed of the glucorticoid stress response varies independently of scope and is predicted by environmental variability and longevity across birds. The acute glucocorticoid response is a key mediator of the coordinated vertebrate response to unpredictable challenges. Rapid glucocorticoid increases initiate changes that allow animals to cope with stressors. The scope of the glucocorticoid response -- defined here as the absolute increase in glucocorticoids -- is associated with individual differences in performance and varies across species with environment and life history. In addition to varying in scope, responses can differ enormously in speed; however, relatively little is known about whether speed and absolute glucocorticoid levels covary, how selection shapes speed, or what aspects of speed are important. We used corticosterone samples collected at 5 time points from 1,750 individuals of 60 species of birds to ask i) how the speed and scope of the glucocorticoid response covary and ii) whether variation in absolute or relative speed is predicted by environmental context or life history. As predicted by a recent model, faster absolute glucocorticoid responses were strongly associated with a larger scope. Despite this covariation, the relative speed of the glucocorticoid response (standardized within species) varied independently of absolute scope, suggesting that selection could operate on both features independently. Species with faster relative glucocorticoid responses lived in locations with more variable temperature and had shorter lifespans. Our results suggest that rapid changes associated with the speed of the glucocorticoid response, such as those occurring through non-genomic receptors, might be an important determinant of coping ability and we emphasize the need for studies designed to measure speed independently of absolute glucocorticoid levels. Sablefish (Anoplopoma fimbria Pallas, 1814) plasma biochemistry and hematology reference intervals including blood cell morphology Plasma biochemistry and hematology reference intervals are integral health assessment tools in all medical fields, including aquatic animal health. As sablefish (Anoplopoma fimbria) are becoming aquaculturally and economically more important, this manuscript provides essential reference intervals (RI) for their plasma biochemistry and hematology along with reference photomicrographs of blood cells in healthy, fasted sablefish. Blood cell morphology can differ between fish species. In addition, blood cell counts and blood chemistry can vary between fish species, demographics, water conditions, seasons, diets, and culture systems, which precludes the use of RIs from other fish species. For this study, blood was collected for plasma biochemistry and hematology analysis between June 20 and July 18, 2019, from healthy, yearling sablefish, hatched and reared in captivity on a commercial diet. Overnight fast of 16-18 hours did not sufficiently reduce lipids in the blood, which led to visible lipemia and frequent rupture of blood cells during analysis. Therefore, sablefish should be fasted for 24 to 36 hours before blood is collected to reduce hematology artifacts or possible reagent interference in plasma biochemistry analysis. Lymphocytes were the most dominant leukocytes (98%), while eosinophils were rare, and basophils were not detected in sablefish. Neutrophils were very large cells with Dohle bodies. In mammals and avian species, Dohle bodies are usually signs of toxic change from inflammation, but no such association was found in these fish. In conclusion, lipemia can interfere with sablefish blood analysis, and available removal methods should be evaluated as fasting for up to 36 h might not always be feasible. Also, more studies are required to establish RI for different developmental stages and rearing conditions. Localization of signaling receptors maximizes cellular informationacquisition in spatially-structured natural environments Cells in natural environments like tissue or soil sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization maximizes information acquisition in simulated natural contexts, including tissue and soil. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold. Broadly, our work shows how cells can maximize the fidelity of information transfer by adapting the spatial organization of signaling molecules to the spatial structure of the environment. Sex differences and risk factors for bleeding in Alagille syndrome Spontaneous lethal bleeds are major cause of death in the pediatric liver disease Alagille syndrome (ALGS), yet risk factors and screening methods have not been established. We performed a systematic review and identified significantly more female than male patients with idiopathic intracranial hemorrhage (10:1). We investigated bleeding and vasculature in patients and a mouse model for ALGS (Jag1Ndr/Ndr mice) and asked whether phenotypes identified in mice could be detected in patients non-invasively. Jag1Ndr/Ndr mice bled spontaneously, exhibiting a thin skull and vascular defects including artery-vein crossings, tortuous vessels, capillary breakdown and CADASIL-like sparse vascular smooth muscle cell coverage which was aggravated by hypertension. Retinographs from patients confirmed tortuous blood vessels and artery-vein crossings in ALGS. In conclusion, Jag1Ndr/Ndr mice could be used to develop interventions for vascular defects in ALGS, and retinography could provide a non-invasive method for vascular analysis in these pediatric patients. Epithelial and neutrophil interactions and coordinated response to Shigella in a human intestinal enteroid-neutrophil co-culture model Polymorphonuclear neutrophils (PMN) are recruited to the gastrointestinal mucosa in response to inflammation, injury, and infection. Herein, we report the development and the characterization of an ex vivo tissue co-culture model consisting of human primary intestinal enteroid monolayers and PMN, and a mechanistic interrogation of PMN-epithelial cell interaction and response to Shigella, a primary cause of childhood dysentery. Cellular adaptation and tissue integration, barrier function, PMN phenotypic and functional attributes, and innate immune responses were examined. PMN within the enteroid monolayers acquired a distinct activated/migratory phenotype that was influenced by direct epithelial cell contact as well as by molecular signals. Seeded on the basal side of the intestinal monolayer, PMN intercalated within the epithelial cells and moved paracellularly toward the apical side. Co-cultured PMN also increased basal secretion of IL-8. Shigella added to the apical surface of the monolayers evoked additional PMN phenotypic adaptations, including increased expression of cell surface markers associated with chemotaxis and cell degranulation (CD47, CD66b, and CD88). Apical Shigella infection triggered rapid transmigration of PMN to the luminal side, NET formation as well as bacterial phagocytosis and killing. Shigella infection modulated cytokine production in the co-culture; apical MCP-1, TNF-, and basolateral IL-8 production were downregulated, while basolateral IL-6 secretion was increased. We demonstrated, for the first time, PMN phenotypic adaptation, mobilization, and coordinated epithelial cell-PMN innate response upon Shigella infection in the human intestinal environment. The enteroid monolayer-PMN co-culture represents a technical innovation for mechanistic interrogation of gastrointestinal physiology, host-microbe interaction, innate immunity, and evaluation of preventive/therapeutic tools.

ImportanceStudies of mucosal immunity and microbial host cell interaction have traditionally relied on animal models and in vitro tissue culture using immortalized cancer cell lines, which render non-physiological and often unreliable results. Herein we report the development and characterization of an ex vivo enteroid-PMN co-culture consisting of normal human intestinal epithelium and a mechanistic interrogation of PMN and epithelial cell interaction and function in the context of Shigella infection. We demonstrated tissue-driven phenotypic and functional adaptation of PMN and a coordinated epithelial cell and PMN response to Shigella, a primary cause of dysentery in young children in the developing world. Thermodynamic coupling between neighboring binding sites in homo-oligomeric ligand sensing proteins from mass resolved ligand dependent population distributions Homo-oligomeric ligand-activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand binding sites. Understanding the basis for this regulation requires both quantifying the free energy {Delta}G transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring-shaped homo-oligomeric trp RNA binding attenuation protein, TRAP. First, we developed nearest-neighbor statistical thermodynamic binding models comprising microscopic free energies for ligand binding to isolated sites {Delta}GN0, and for coupling between one or both adjacent sites, {Delta}GN1 and {Delta}GN2. Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2N possible liganded states. We then experimentally monitored liganddependent population shifts using conventional spectroscopic and calorimetric methods, and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand-dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2N liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure-thermodynamic linkage that drives cooperativity.

TOC Figure (draft)

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=74 SRC="FIGDIR/small/484990v3_ufig1.gif" ALT="Figure 1">
View larger version (30K):
org.highwire.dtl.DTLVardef@2ce1f3org.highwire.dtl.DTLVardef@a6541borg.highwire.dtl.DTLVardef@1dbd08aorg.highwire.dtl.DTLVardef@58297e_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOTOC Figure.C_FLOATNO Ligand (Trp) binding to multiple sites on homo-oligomeric ring-shaped proteins like TRAP alters their functional states. Homotropic cooperativity is expected to alter the activation pathway in response to cellular ligand concentration. In the presence of positive nearest-neighbor cooperativity, ligand binding is favored at adjacent sites, whereas in the absence of cooperativity, a random "Normal" distribution is expected.

C_FIG Movement decoding using spatio-spectral features of cortical and subcortical local field potentials The application of machine learning to intracranial signal analysis has the potential to revolutionize deep brain stimulation (DBS) by personalizing therapy to dynamic brain states, specific to symptoms and behaviors. Most decoding pipelines for movement decoding in the context of adaptive DBS are based on single channel frequency domain features, neglecting spatial information available in multichannel recordings. Such features are extracted either from DBS lead recordings in the subcortical target and/or from electrocorticography (ECoG). To optimize the simultaneous use of both types of signals, we developed a supervised online-compatible movement decoding pipeline based on multichannel and multiple site recordings. We found that adding spatial information to the model has the potential to improve decoding. In addition, we demonstrate movement decoding from spatio-spectral features derived from cortical and subcortical oscillations. We demonstrate between-patients variability of the spatial neural maps and its relationship to feature decoding performance. This application of spatial filters to decode movement from combined cortical and subcortical recordings is an important step in developing machine learning approaches for intelligent DBS systems. Accurate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion Single cell RNA-seq (scRNA-seq) has been widely used to uncover cellular heterogeneity, however, the constraints of cost make it impractical as a routine on large patient cohorts. Here we present ENIGMA, a method that accurately deconvolute bulk tissue RNA-seq into single cell-type resolution given the knowledge gained from scRNA-seq. ENIGMA applies a matrix completion strategy to minimize the distance between mixture transcriptome and weighted combination of cell type-specific expression, allowing quantification of cell type proportions and reconstruction of cell type-specific transcriptome. The superior performance of ENIGMA was validated in simulated and realistic datasets, including disease-related tissues, demonstrating its ability in novel biological findings. Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.) Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Significance StatementUsing traditional and automated/high-throughput (using a novel deep learning approach) phenotyping methods we studied leaf anatomical variation in sunflower. Genome-wide association (GWA) analyses identified numerous genomic regions underlying individual trait variation and regions underlying major multivariate axes of phenotypic variation. These results illustrate the value of employing a multivariate approach to GWA analyses and shed light on the extent to which leaf trait (co-)variation can be genetically decoupled to explore novel phenotypic space. Liberation from fibrogenesis or tumorigenesis via cellular senescence by a novel small molecule BackgroundUncontrolled proliferative diseases such as fibrosis or cancer can be fatal. We previously found that a compound containing a chromone scaffold, ONG41008, had potent anti-fibrotic effects in diseased human lung myofibroblasts but not normal human lung fibroblasts.

MethodsWe investigated the effects of ONG41008 on tumor cells, and compared these effects with those in pathogenic myofibrotic cells and normal fibroblasts cells.

FindingsStimulation of A549 lung carcinoma epithelial cells with ONG41008 resulted in cellular senescence, indicating that dysregulated cell proliferation is common to fibrotic cells and tumor cells. Replicative senescence of A549 cells resulted in multinucleation, which was followed by oncogene-induced senescence. There was significant upregulation of expression and nuclear translocation of p-TP53 and p16 in ONG41008-treated A549 cells, and all cells died after 72 hr. Similar effects occurred after ONG41008 treatment in several human aggressive cancer cell lines such as PANC1, MCF7, PC3, or primary non-small cell lung carcinoma cells. Unlike cisplatin, ONG41008 was not toxic to normal human lung fibroblasts or primary prostate epithelial cells, suggesting ONG41008 can distinguish the intracellular microenvironment between normal cells and aged or diseased cells. This effect might occur as a result of the increased NAD/NADH ratio or increased lactate dehydrogenase levels in aged or diseased cells.

InterpretationTo our best knowledge, this is the first study to show that a small molecule can arrest uncontrolled proliferation during fibrogenesis or tumorigenesis in vitro. ONG41008 could be a potential drug for a broad range of fibrotic or tumorigenic diseases.

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSThe notion that aging is a disease and that diseases occur as a consequence of aging was first put forward by David Sinclair and colleagues (Aging Cell;14(4):497-510). There are biological systems that provide evidence for this notion; for example, somatic cells can revert to embryonic cells, producing younger somatic cells. This phenomenon underlies induced pluripotent stem cells. Another example is that some types of jellyfish can live forever. These examples suggest that a counter-aging program exists in animals. Human diseases are the manifestations of cell aging generated by the accumulation of somatic mutations. Aged and pathogenic cells are senescent, so a drug that specifically targeted senescent cells might initiate a cellular program that could ameliorate age-associated disease. Indeed, the kinase inhibitor dasatinib induces cellular senescence (Clin Ther 2007 29:2289-2308). In 2017, two drugs that target senescent cells were identified: quercetin and fisetin. These drugs selectively kill senescent cells, and are referred to as senolytic drugs or senotherapeutics (Aging 2017 8;9(3):955-963). Although it is established that senescent cells accumulate in cancer and idiopathic pulmonary fibrosis (Nat Commun 2017 23;8 :14532), the effect of senolytic drugs in these diseases is largely unknown.

Added value of this studyThis study characterized a novel drug, termed ONG41008, which was found to have both senogenic and senolytic effects in cell-based assays. ONG41008 induced senescence in myofibroblasts and several cancer cell lines representative of aggressive human cancers, which was followed by cell death. Importantly, ONG41008 exhibited essentially no toxicity on normal human lung fibroblasts or primary prostate epithelial cells.

Implication of all the evidenceBased on our results, we believe that ONG41008 is a potent inducer of cellular senescence (replicative senescence and oncogene-induced senescence) and causes arrest of uncontrolled, pathogenic proliferation of myofibroblasts or cancer cells. 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. RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub.