Pulsatile contractions and pattern formation in excitable actomyosin cortex The actin cortex is an active adaptive material, embedded with complex regulatory networks that can sense, generate, and transmit mechanical forces. The cortex exhibits a wide range of dynamic behaviours, from generating pulsatory contractions and travelling waves to forming organised structures. Despite the progress in characterising the biochemical and mechanical components of the actin cortex, the emergent dynamics of this mechanochemical system is poorly understood. Here we develop a reaction-diffusion model for the RhoA signalling network, the upstream regulator for actomyosin assembly and contractility, coupled to an active actomyosin gel, to investigate how the interplay between chemical signalling and mechanical forces regulate stresses and patterns in the cortex. We demonstrate that mechanochemical feedback in the cortex acts to destabilise homogeneous states and robustly generate pulsatile contractions. By tuning active stress in the system, we show that the cortex can generate propagating contraction pulses, form network structures, or exhibit topological turbulence.
Cryptic pathogen-sugar interactions revealed by universal saturation transfer analysis Many host pathogen interactions such as human viruses (including non-SARS-coronaviruses) rely on attachment to host cell-surface glycans. There are conflicting reports about whether the Spike protein of SARS-CoV-2 binds to sialic acid commonly found on host cell-surface N-linked glycans. In the absence of a biochemical assay, the ability to analyze the binding of glycans to heavily- modified proteins and resolve this issue is limited. Classical Saturation Transfer Difference (STD) NMR can be confounded by overlapping sugar resonances that compound with known experimental constraints. Here we present  universal saturation transfer analysis (uSTA), an NMR method that builds on existing approaches to provide a general and automated workflow for studying protein-ligand interactions. uSTA reveals that B-origin-lineage-SARS-CoV-2 spike trimer binds sialoside sugars in an  end on manner and modelling guided by uSTA localises binding to the spike N-terminal domain (NTD). The sialylated-polylactosamine motif is found on tetraantennary human N-linked-glycoproteins in deeper lung and may have played a role in zoonosis. Provocatively, sialic acid binding is abolished by mutations in some subsequent SARS- CoV-2 variants-of-concern. A very high resolution cryo-EM structure confirms the NTD location and  end on mode; it rationalises the effect of NTD mutations and the structure-activity relationship of sialic acid analogues. uSTA is demonstrated to be a robust, rapid and quantitative tool for analysis of binding, even in the most demanding systems.

Extended AbstractThe surface proteins found on both pathogens and host cells mediate entry (and exit) and influence disease progression and transmission. Both types can bear host-generated post- translational modifications such as glycosylation that are essential for function but can confound biophysical methods used for dissecting key interactions. Several human viruses (including non- SARS-coronaviruses) attach to host cell-surface N-linked glycans that include forms of sialic acid (sialosides). There remains, however, conflicting evidence as to if or how SARS-associated coronaviruses might use such a mechanism. Here, we demonstrate quantitative extension of  saturation transfer protein NMR methods to a complete mathematical model of the magnetization transfer caused by interactions between protein and ligand. The method couples objective resonance-identification via a deconvolution algorithm with Bloch-McConnell analysis to enable a structural, kinetic and thermodynamic analysis of ligand binding beyond previously-perceived limits of exchange rates, concentration or system. Using an automated and openly available workflow this  universal saturation transfer analysis (uSTA) can be readily-applied in a range of even heavily-modified systems in a general manner to now obtain quantitative binding interaction parameters (KD, kEx). uSTA proved critical in mapping direct interactions between natural sialoside sugar ligands and relevant virus-surface attachment glycoproteins - SARS-CoV-2-spike and influenza-H1N1-haemagglutinin variants - by quantitating ligand signal in spectral regions otherwise occluded by resonances from mobile protein glycans (that also include sialosides). In B- origin-lineage-SARS-CoV-2 spike trimer  end on-binding to sialoside sugars was revealed contrasting with  extended surface-binding for heparin sugar ligands; uSTA-derived constraints used in structural modelling suggested sialoside-glycan binding sites in a beta-sheet-rich region of spike N-terminal domain (NTD). Consistent with this, uSTA-glycan binding was minimally- perturbed by antibodies that neutralize the ACE2-binding domain (RBD) but strongly disrupted in spike from the B1.1.7/alpha and B1.351/beta variants-of-concern, which possess hotspot mutations in the NTD. Sialoside binding in B-origin-lineage-NTD was unequivocally pinpointed by cryo-EM to a site that is created from residues that are notably deleted in variants (e.g. H69,V70,Y145 in alpha). An analysis of beneficial genetic variances in cohorts of patients from early 2020 suggests a model in which this site in the NTD of B-origin-lineage-SARS-CoV-2 (but not in alpha/beta-variants) may have exploited a specific sialylated-polylactosamine motif found on tetraantennary human N-linked-glycoproteins in deeper lung. Together these confirm a novel binding mode mediated by the unusual NTD of SARS-CoV-2 and suggest how it may drive virulence and/or zoonosis via modulation of glycan attachment. Since cell-surface glycans are widely relevant to biology and pathology, uSTA can now provide ready, quantitative, widespread analysis of complex, host-derived and post-translationally modified proteins with putative ligands relevant to disease even in previously confounding complex systems.
Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Citation needed? Wikipedia and the COVID-19 pandemic BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Risk modelling of single cell populations revealed the heterogeneity of immune infiltration in hepatocellular carcinoma Hepatocellular carcinoma (HCC) is one of the most common primary hepatic malignancies. E2F transcription factors play an important role in the tumorigenesis and progression of HCC, mainly through the RB/E2F pathway. Prognostic models for HCC based on gene signatures have been developed rapidly in recent years; however, their discriminating ability at the single cell level remains elusive, which could reflect the underlying mechanisms driving the sample bifurcation. In this study, we constructed and validated a predictive model based on E2F expression, successfully stratifying HCC patients into two groups with different survival risks. Then we used a single-cell dataset to test the discriminating ability of the predictive model on infiltrating T cells, demonstrating remarkable cellular heterogeneity as well as altered cell fates. We identified distinct cell subpopulations with diverse molecular characteristics. We also found that the distribution of cell subpopulations varied considerably across onset stages among patients, providing fundamental basis for patient-oriented precision evaluation. Moreover, single-sample gene set enrichment analysis revealed that subsets of CD8+ T cells with significantly different cell adhesion level could be associated with different patterns of tumor cell dissemination. Therefore, our findings linked the conventional prognostic gene signature to the immune microenvironment and cellular heterogeneity at the single-cell level, thereby providing deeper insights into the understanding of HCC tumorigenesis. wnt16 exerts pleiotropic effects on bone and lean mass in zebrafish Osteoporosis, a disease of bone fragility, and sarcopenia, a condition of reduced muscle mass and strength, commonly occur in the same individual--a condition termed osteosarcopenia. Genetic variants at the CPED1-WNT16 locus have pleiotropic effects on human bone mineral density (BMD) and lean tissue mass. However, our current understanding of Wnt16 and its influence on bone fails to explain effects on lean mass at this locus. In this study, we sought to better understand the biological mechanism underlying pleiotropy at the CPED1-WNT16 locus on BMD and lean mass. Using single cell analysis, microCT imaging, and genetic approaches in zebrafish, we show that wnt16 exerts dual influence on muscle and bone, and provide evidence that this occurs through a function in somite development, a period in which precursors for both of these tissues are specified. We show that mutations in wnt16 alter vertebral size and lean mass during axial skeletogenesis. We show that wnt16+ cells transcriptionally resemble cells within the amniote dermomyotome, and provide evidence they function as myogenic precursors during embryonic development. We also provide evidence that WNT16 is a gene of major effect on lean mass at the CPED-WNT16 locus. Our studies support the potential for wnt16 to exert pleiotropic effects on bone and lean mass through a function in myogenic precursors. Lipid membranes modulate the activity of RNA through sequence-dependent interactions RNA is a ubiquitous biomolecule that can serve as both catalyst and information carrier. Understanding how RNA bioactivity is controlled is crucial for elucidating its physiological roles and potential applications in synthetic biology. Here we show that lipid membranes can act as RNA organization platforms, introducing a novel mechanism for ribo-regulation. The activity of R3C ribozyme can be modified by the presence of lipid membranes, with direct RNA-lipid interactions dependent on RNA nucleotide content, base pairing and length. In particular, the presence of guanine in short RNAs is crucial for RNA-lipid interactions, and G-quadruplex formation further promotes lipid binding. Lastly, by artificially modifying the R3C substrate sequence to enhance membrane binding we generated a lipid-sensitive ribozyme reaction with riboswitch-like behavior. These findings introduce RNA-lipid interactions as a tool for developing synthetic riboswitches and novel RNA-based lipid biosensors, and bear significant implications for RNA World scenarios for the origin of life. POPULATION GENOMIC ANALYSIS OF CRYPTOCOCCUS BRAZILIAN ISOLATES REVEALS AN AFRICAN TYPE SUBCLADE DISTRIBUTION The genomes of a large number of Cryptococcus neoformans isolates have been sequenced and analyzed in recent years. These genomes have been used to understand the global population structure of this opportunistic pathogen. However, only a small number of South American isolates have been considered in these studies, and the population structure of C. neoformans in this part of the world remains elusive. Here, we analyzed the genomic sequences of 53 Brazilian Cryptococcus isolates and deciphered the C. neoformans population structure in this country. Our data reveal an African-like structure that suggested repeated intercontinental transports from Africa to South America. We also identified a mutator phenotype in one VNBII Brazilian isolate, exemplifying how fast-evolving isolates can shape the Cryptococcus population structure. Finally, phenotypic analyses revealed wide diversity but not lineage specificity in the expression of classical virulence traits within the set of isolates. Accounting for errors in data improves timing in single-cell cancer evolution Single-cell sequencing provides a new way to explore the evolutionary history of cells. Compared to traditional bulk sequencing, where a population of heterogeneous cells is pooled to form a single observation, single-cell sequencing isolates and amplifies genetic material from individual cells, thereby preserving the information about the origin of the sequences. However, single-cell data is more error-prone than bulk sequencing data due to the limited genomic material available per cell. Here, we present error and mutation models for evolutionary inference of single-cell data within a mature and extensible Bayesian framework, BEAST2. Our framework enables integration with biologically informative models such as relaxed molecular clocks and population dynamic models. Our simulations show that modeling errors increase the accuracy of relative divergence times and substitution parameters. We reconstruct the phylogenetic history for a colorectal cancer patient and a healthy patient from single-cell DNA sequencing data. We find that the estimated times of terminal splitting events are shifted forward in time compared to models which ignore errors. We observed that not accounting for errors can overestimate the phylogenetic diversity in single-cell DNA sequencing data. We estimate that 30-50% of the apparent diversity can be attributed to error. Our work enables a full Bayesian approach capable of accounting for errors in the data within the integrative Bayesian software framework BEAST2. Diverse classes of constraints enable broader applicability of a linear programming-based dynamic metabolic modeling framework Current metabolic modeling tools suffer from a variety of limitations, from scalability to simplifying assumptions, that preclude their use in many applications. We recently created a modeling framework, LK-DFBA, that addresses a key gap: capturing metabolite dynamics and regulation while retaining a potentially scalable linear programming structure. Key to this frameworks success are the linear kinetics and regulatory constraints imposed on the system. However, while the linearity of these constraints reduces computational complexity, it may not accurately capture the behavior of many biochemical systems. Here, we developed three new classes of LK-DFBA constraints to better model interactions between metabolites and the reactions they regulate. We tested these new approaches on several synthetic and biological systems, and also performed the first-ever comparison of LK-DFBA predictions to experimental data. We found that no single constraint approach was optimal across all systems examined, and systems with the same topological structure but different parameters were often best modeled by different types of constraints. However, we did find that the optimal constraint approach was generally robust to local perturbations of the system, indicating that just a single wild-type dataset could allow identification of the ideal constraint for a given system. These results suggest that the availability of multiple constraint approaches will allow LK-DFBA to model a wider range of metabolic systems. GLOBAL LAND-USE AND LAND-COVER DATA: HISTORICAL, CURRENT AND FUTURE SCENARIOS Land-use land-cover (LULC) data are important predictors of species occurrence and biodiversity threat. Although there are LULC datasets available under current conditions, there is a lack of such data under historical and future climatic conditions. This hinders, for example, projecting niche and distribution models under global change scenarios at different time scenarios. The Land Use Harmonization Project (LUH2) is a global terrestrial dataset at 0.25{degrees} spatial resolution that provides LULC data from 850 to 2300 for 12 LULC state classes. The dataset, however, is compressed in a file format (NetCDF) that is incompatible for many analyses and requires layer extractions and transformations that are intractable for most researchers. Here we selected and transformed the LUH2 in a standard GIS format data to make it more user-friendly. We provide LULC for every year from 850 to 2100, and from 2015 on, the LULC dataset is provided under two Shared Socioeconomic Pathways (SSP2-4.5 and SSP5-8.5). We provide two types of files for each year: separate files with continuous values for each of the 12 LULC state classes, and a single categorical file with all state classes combined. To create the categorical layer, we assigned the state with the highest value in a given pixel among the 12 continuous data. LUH2 predicts a pronounced decrease in primary forest, particularly noticeable in the Amazon, the Brazilian Atlantic Forest, the Congo Basin and the boreal forests, an equally pronounced increase in secondary forest and non-forest lands, and in croplands in the Brazilian Atlantic Forest and sub-Saharan Africa. The final dataset provides LULC data for 1251 years that will be of interest for macroecology, ecological niche modeling, global change analysis, and other applications in ecology and conservation. Hepatic MIR20B promotes nonalcoholic fatty liver disease by suppressing PPARA BackgroundNon-alcoholic fatty liver disease (NAFLD) is characterized by excessive lipid accumulation and imbalances in lipid metabolism in the liver. Although nuclear receptors (NRs) play a crucial role in hepatic lipid metabolism, the underlying mechanisms of NR regulation in NAFLD remain largely unclear.

MethodsUsing network analysis and RNA-seq to determine the correlation between NRs and microRNA in human NAFLD patients, we revealed that MIR20B specifically targets PPARA. MIR20B mimic and anti-MIR20B were administered to human HepG2 and Huh-7 cells and mouse primary hepatocytes as well as high fat diet (HFD)- or methionine-deficient diet (MCD)-fed mice to verify the specific function of MIR20B in NAFLD. We tested the inhibition of the therapeutic effect of a PPAR agonist, fenofibrate, by MIR20B and the synergic effect of combination of fenofibrate with anti-MIR20B.

ResultsWe revealed that MIR20B specifically targets PPARA through miRNA regulatory network analysis of nuclear receptor genes in NAFLD. The expression of MIR20B was upregulated in free fatty acid (FA)- treated hepatocytes and the livers of both obesity-induced mice and NAFLD patients. Overexpression of MIR20B significantly increased hepatic lipid accumulation and triglyceride levels. Furthermore, MIR20B significantly reduced FA oxidation and mitochondrial biogenesis by targeting PPARA. In MIR20B-introduced mice, the effect of fenofibrate to ameliorate hepatic steatosis was significantly suppressed. Finally, inhibition of MIR20B significantly increased FA oxidation and uptake, resulting in improved insulin sensitivity and a decrease in NAFLD progression. Moreover, combination of fenofibrate and anti-MIR20B exhibited the synergic effect on improvement of NAFLD in MCD-fed mice.

ConclusionsTaken together, our results demonstrate that the novel MIR20B targets PPARA, plays a significant role in hepatic lipid metabolism, and present an opportunity for the development of novel therapeutics for NAFLD.

FundingThis research was funded by Korea Mouse Phenotyping Project (2016M3A9D5A01952411), the National Research Foundation of Korea (NRF) grant funded by the Korea government (2020R1F1A1061267, 2018R1A5A1024340, NRF-2021R1I1A2041463, 2020R1I1A1A01074940), and the Future-leading Project Research Fund (1.210034.01) of UNIST. Reuse of a wing venation gene-regulatory network in patterning the eyespot rings of butterflies Novel organismal traits might reuse ancestral gene-regulatory networks (GRNs) in their development, but data supporting this mechanism are still sparse. Here we propose that a conserved venation gene network, used in patterning insect wings, has been recruited to pattern the sharp and distinct rings of color in butterfly eyespot development. This network involves Decapentaplegic (Dpp), acting as a central morphogen, activating optix and spalt at different concentration thresholds. It also involves Spalt cross-regulating optix resulting in the formation of a sharp boundary between the two genes. First, using in situ hybridizations, we show that dpp is expressed in the eyespot centers during the early pupal stage and that CRISPR-Cas9 embryonic knockouts of dpp, optix, and spalt result in the loss of eyespots, loss of the orange ring, and loss of the black scales, respectively. Furthermore, using RNAi during the pupal stage, we show that dpp down-regulation results in the complete or partial loss of eyespots. Using CRISPR disruptions of optix or spalt, followed by immunostainings, we show that Spalt represses optix in the central region of the eyespot, limiting optix expression to a more peripheral ring. We further illustrate that optix disruptions alter scale pigmentation and scale ultrastructure as they transform orange to brown scales. Altogether, the regulatory interactions between dpp, spalt and optix mimic those found in earlier larval developmental stages, in the anterior compartment of an insects wing, where these genes play a role in venation patterning. The network similarities suggest that the venation GRN was co-opted to aid in the differentiation of the eyespot rings. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. 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. Oxford Nanopore R10.4 long-read sequencing enables near-perfect bacterial genomes from pure cultures and metagenomes without short-read or reference polishing Long-read Oxford Nanopore sequencing has democratized microbial genome sequencing and enables the recovery of highly contiguous microbial genomes from isolates or metagenomes. However, to obtain near-perfect genomes it has been necessary to include short-read polishing to correct insertions and deletions derived from homopolymer regions. Here, we show that Oxford Nanopore R10.4 can be used to generate near-perfect microbial genomes from isolates or metagenomes without shortread or reference polishing. 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. Differential fate of acellular vascular scaffolds in vivo involves macrophages and T-cell subsets Biological scaffold or implant is a popular choice for the preparation of tissue-engineered organs and has the potential to address donor shortage in the clinics. However, biological scaffolds prepared by physical or chemical agents cause damage to the extracellular matrix by potentially inducing immune responses after implantation. The current study explores an alternative route for the preparation of acellular scaffolds and explores the fate of the prepared scaffolds in a milieu of immune cells following implantation without using immunosuppressant. Using the syngeneic (Lewis male-Lewis female) and allogeneic (Brown Norway male-Lewis female) models and different tissue routes (subcutaneous vs omentum) for transplantation, normal blood vascular scaffolds were implanted which was converted to acellular vascular scaffolds by in vivo natural decellularization at the end of 2 months of observation. We also prepared chemically decellularized acellular scaffolds from normal untreated blood vascular scaffolds using a cocktail of chemicals which was also similarly placed in subcutaneous and omentum sites. Here, we applied in-depth quantitative proteomics along with histology and image analysis to comprehensively describe and compare the proteome of the natural and chemically decellularized scaffold. Our data confirm that site-specific advantages exist in modulating the ECM and regulating the immune responses (macrophage and T cells) following implantation, which possibly led to the production of an acellular scaffold (natural decellularization) under in vivo conditions. The current approach opens up the possibility to create tailor-made acellular scaffolds to build functional blood vessels. In addition, the identification of different tissue sites facilitates differential immune response against the scaffolds. This study provides a rich resource aimed toward an enhanced mechanistic understanding to study immune responses under similar settings in the field of transplantation and regenerative medicine.

Impact statementThe development of a scaffold helps in the preparation of a functional organ in the clinics. In the current study, we prepared an acellular vascular scaffold by utilizing site specific tissue changes and vis-a-vis compared with a conventionally chemically prepared biological scaffold at genomic and protein level, which helped us to identify immunological trigger following implantation. The current study which was carried out without any immunosuppressive agents could help to establish (a) alternative strategies for preparing biological scaffolds as well as (b) implantable sites as potential bioreactors to circumvent any adverse immune reactions for acceptance of the scaffold/implant post implantation. A terrain treadmill to study animal locomotion through large obstacles A major challenge to understanding locomotion in complex 3-D terrain with large obstacles is to create tools for controlled, systematic lab experiments. Existing terrain arenas only allow observations at small spatiotemporal scales ([~]10 body length, [~]10 stride cycles). Here, we create a terrain treadmill to enable high-resolution observations of animal locomotion through large obstacles over large spatiotemporal scales. An animal moves through modular obstacles on an inner sphere, while a rigidly-attached, concentric, transparent outer sphere rotated with the opposite velocity via closed-loop feedback to keep the animal on top. During sustained locomotion, a discoid cockroach moved through pillar obstacles for 25 minutes ({approx} 2500 strides) over 67 m ({approx} 1500 body lengths), and was contained within a radius of 4 cm (0.9 body length) for 83% of the duration, even at speeds of up to 10 body length/s. The treadmill enabled observation of diverse locomotor behaviors and quantification of animal-obstacle interaction.

Summary StatementA device keeps a small animal on top of a sphere while it traverses obstacles, creating a "terrain treadmill" to study locomotion over a long time and distance at high-resolution. Chromosome breaks in breast cancers occur near herpes tumor virus sequences and explain why cancer drugs stop working One explanation for why breast cancer has a relentless tendency to come back after treatment is that breast cancer phenomes are fluid, so it is difficult for therapy to target new phenotypes that arise. In breast cancer, viral interactions cause chromosome breakage and fusion, producing anaphase bridges containing two centromeres. Cell division then leads to more chromosome breaks, more mutations, and new cancer phenotypes. Analyses of public data from about 2100 breast cancers implicate Epstein-Barr virus (EBV or Human Herpes Virus 4) in the cycles that break chromosomes and produce anaphase bridges. EBV generates breakage cycles with highly fragmented chromosomes (chromothripsis) as an inherent feature. These episodes continually produce new cancer driver mutations. DNA at chromothripsis boundaries in breast cancers consistently matches EBV tumor variants. Breakpoints are often near known EBV binding sites, and breast cancer shares an EBV methylation signature with other known EBV cancers. Breast cancer breakpoints on all chromosomes cluster around the same positions as in cancers definitively associated with EBV infection (nasopharyngeal cancer and Burkitts lymphoma). The status of EBV in the microbiome and the presence of anaphase bridges are phenotypic markers of breast cancer that are not explicitly targeted by chemotherapy. The breast cancer phenome continually supplies new cancer driver genes that foster breast cancer recurrence.

SummaryBreaks in breast cancer chromosomes and in fragments of anaphase bridges occur mainly near regions that resemble Epstein-Barr virus (EBV) tumor variants

Extensive breaks are an inherent outcome of anaphase bridges and the reassembly of highly fragmented chromosomes continually generates new cancer driver mutations and new cancer phenotypes.

Treatment does not explicitly target the underlying breakage and reassembly cycles so breast cancer often comes back.

The breast cancer phenome is not fixed and can continuously generate new phenotypes

Graphical Abstract

O_FIG O_LINKSMALLFIG WIDTH=154 HEIGHT=200 SRC="FIGDIR/small/467751v4_ufig1.gif" ALT="Figure 1">
View larger version (38K):
org.highwire.dtl.DTLVardef@6537cforg.highwire.dtl.DTLVardef@132fa50org.highwire.dtl.DTLVardef@1cfcb04org.highwire.dtl.DTLVardef@11979aa_HPS_FORMAT_FIGEXP  M_FIG C_FIG A TIR-1/SARM1 phase transition underlies p38 immune pathway activation in the C. elegans intestine Intracellular signaling regulators can be concentrated into membrane-free, higher-ordered protein assemblies to initiate protective responses during stress -- a process known as phase transition. Here, we show that a phase transition of the Caenorhabditis elegans Toll/interleukin-1 receptor domain protein (TIR-1), an NAD+ glycohydrolase homologous to mammalian sterile alpha and TIR motif-containing 1 (SARM1), underlies p38 PMK-1 immune pathway activation in C. elegans intestinal epithelial cells. Through visualization of fluorescently labeled TIR-1/SARM1 protein, we demonstrate for the first time that physiologic stresses, both pathogen and non-pathogen, induce multimerization of TIR-1/SARM1 into visible puncta within intestinal epithelial cells. In vitro enzyme kinetic analyses revealed that, like mammalian SARM1, the NAD+ glycohydrolase activity of C. elegans TIR-1 is dramatically potentiated by protein oligomerization and a phase transition. Accordingly, C. elegans with genetic mutations that specifically block either multimerization or the NAD+ glycohydrolase activity of TIR-1/SARM1 fail to induce p38 PMK phosphorylation, are unable to increase immune effector expression, and are dramatically susceptible to bacterial infection. Finally, we demonstrate that the TIR-1/SARM1 phase transition is modified by dietary cholesterol, revealing a new adaptive response that allows a metazoan host to anticipate pathogen threats during micronutrient deprivation, a time of relative susceptibility to infection. When cholesterol is limited, TIR-1/SARM1 oligomerizes into puncta in intestinal epithelial cells and engages its NAD+ glycohydrolase activity, which increases p38 PMK-1 phosphorylation, and primes immune effector induction in a manner that promotes pathogen clearance from the intestine during a subsequent infection. Thus, a phase transition of TIR-1/SARM1 as a prerequisite for its NAD+ glycohydrolase activity is strongly conserved across millions of years of evolution and is essential for diverse physiological processes in multiple cell types. A novel Fiji/ImageJ plugin for the rapid analysis of blebbing cells When confined, cells have recently been shown to undergo a phenotypic switch to what has been termed, fast amoeboid (leader bleb-based) migration. However, as this is a nascent area of research, few tools are available for the rapid analysis of cell behavior. Here, we demonstrate that a novel Fiji/ImageJ-based plugin, Analyze_Blebs, can be used to quickly obtain cell migration parameters and morphometrics from time lapse images. As validation, we show that Analyze_Blebs can detect significant differences in cell migration and morphometrics, such as the largest bleb size, upon introducing different live markers of F-actin, including F-tractin and LifeAct tagged with green and red fluorescent proteins. We also demonstrate, using flow cytometry, that live markers increase total levels of F-actin. Furthermore, that F-tractin increases cell stiffness, which was found to correlate with a decrease in migration, thus reaffirming the importance of cell mechanics as a determinant of Leader Bleb-Based Migration (LBBM).

Highlight summaryA new plugin, Analyze_Blebs, enables the rapid analysis of cell migration and morphometrics of fast amoeboid cells. Morphometrics combined with cell stiffness are found to be predictive of confined, fast amoeboid (leader bleb-based), migration. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Controllability boosts neural and cognitive signatures of information seeking In uncertain environments, seeking information about alternative choice options is essential for adaptive learning and decision-making. However, information seeking is usually confounded with changes-of-mind about the reliability of the preferred option. Here, we exploited the fact that information seeking requires control over which option to sample to isolate its behavioral and neurophysiological signatures. We found that changes-of-mind occurring with control require more evidence against the current option, are associated with reduced confidence, but are nevertheless more likely to be confirmed on the next decision. Multimodal neurophysiological recordings showed that these changes-of-mind are preceded by stronger activation of the dorsal attention network in magnetoencephalography (MEG), and followed by increased pupil-linked arousal during the presentation of decision outcomes. Together, these findings indicate that information seeking increases the saliency of evidence perceived as the direct consequence of ones own actions.

HighlightsO_LIExploration-exploitation dilemmas confound information seeking with changes-of-mind.
C_LIO_LIUnlike changes-of-mind, information seeking requires control over information sampling.
C_LIO_LIInformation seeking is associated with reduced confidence and active hypothesis testing.
C_LIO_LIStronger neurophysiological correlates of attention and arousal during information seeking.
C_LI A revised view on the evolution of glutamine synthetase isoenzymes in plants Glutamine synthetase (GS) is a key enzyme responsible for the incorporation of inorganic nitrogen in the form of ammonium into the amino acid glutamine. The genes encoding GS are among the oldest existing genes in living organisms. In plants, two groups of functional GS enzymes are found: eubacterial GSIIb (GLN2) and eukaryotic GSIIe (GLN1/GS). Phylogenetic analyses have shown that the GLN2 group originated from bacteria following horizontal gene transfer. Only GLN1/GS genes are found in vascular plants, which suggests that they are involved in the final adaptation of plants to terrestrial life. The present phylogenetic study reclassifies the different GS of seed plants into three clusters: GS1a, GS1b and GS2. The presence of genes encoding GS2 has been expanded to Cycadopsida gymnosperms, which suggests the origin of this gene in a common ancestor of Cycadopsida, Ginkgoopsida and angiosperms. GS1a genes have been identified in all gymnosperms, basal angiosperms and some Magnoliidae species. Previous studies in conifers and the gene expression profiles obtained in ginkgo and magnolia in the present work could explain the absence of GS1a in more recent angiosperm species (e.g., monocots and eudicots) due to the redundant roles of GS1a and GS2 in photosynthetic cells. Altogether, the results provide a better understanding of the evolution of plant GS isoenzymes and their physiological roles, which is valuable for improving crop nitrogen use efficiency and productivity. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain.