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.
Emergent evolutionary forces in spatial models of microbial growth in the human gut microbiota The genetic composition of the gut microbiota is constantly reshaped by ecological and evolutionary forces. These strain-level dynamics can be challenging to understand because they emerge from complex spatial growth processes that take place within a host. Here we introduce a general population genetic framework to predict how stochastic evolutionary forces emerge from simple models of microbial growth in spatially extended environments like the intestinal lumen. Our framework shows how fluid flow and longitudinal variation in growth rate combine to shape the frequencies of genetic variants in sequenced fecal samples, yielding analytical expressions for the effective generation times, selection coefficients, and rates of genetic drift. We find that the emergent evolutionary dynamics can often be captured by well-mixed models that lack explicit spatial structure, even when there is substantial spatial variation in species-level composition. By applying these results to the human colon, we find that continuous fluid flow and simple forms of wall growth are unlikely to create sufficient bottlenecks to allow large fluctuations in mutant frequencies within a host. We also find that the effective gener-ation times may be significantly shorter than expected from traditional average growth rate estimates. Our results provide a starting point for qua ntifying genetic turnover in spatially extended settings like the gut microbiota, and may be relevant for other microbial ecosystems where unidirectional fluid flow plays an important role.
Embedded 3D printing in self-healing annealable composites for precise patterning of functionally mature human neural constructs Human in vitro models of neural tissue with controllable cellular identity, tunable microenvironment, and defined spatial arrangement are needed to facilitate studies of brain development and disease. Towards this end, embedded printing in jammed microgel supports (i.e., granular gels) holds great promise as it allows precise and programmable patterning of extremely soft and compliant tissue constructs. However, in contrast to the vast material landscape available for bulk hydrogels, granular printing support formulations are restricted to a handful of materials without the ability for facile adjustment of biofunctional properties of the cellular microenvironment. Therefore, there has been a need for novel materials that take advantage of versatile biomimicry of bulk hydrogels while providing high-fidelity support for embedded printing akin to granular gels. To address this need, we present a modular platform for bioengineering of neuronal networks via direct embedded 3D printing of human stem cells inside Self-Healing Annealable Particle-Extracellular matrix (SHAPE) composites. SHAPE composites consist of soft microgels immersed in viscous extracellular-matrix solution to enable precise freeform patterning of human stem cells and consequent generation and long-term maintenance of mature subtype-specific neurons that extend projections within the volume of the annealed support. The developed approach further allows multi-ink deposition, live spatial and temporal monitoring of oxygen levels, as well as creation of vascular channels. Due to its modularity, SHAPE biomanufacturing toolbox not only offers a solution for functional modeling of mechanically sensitive neural constructs, but also has potential to be applied to a wide range of biomaterials with different crosslinking mechanisms to model tissues and diseases where recapitulation of complex architectural features and topological cues is essential. 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. Data-Driven Network Models for Genetic Circuits From Time-Series Data with Incomplete Measurements Synthetic gene networks are frequently conceptualized and visualized as static graphs. This view of biological programming stands in stark contrast to the transient nature of biomolecular interaction, which is frequently enacted by labile molecules that are often unmeasured. Thus, the network topology and dynamics of synthetic gene networks can be difficult to verify in vivo or in vitro, due to the presence of unmeasured biological states. Here we introduce the dynamical structure function as a new mesoscopic, data-driven class of models to describe gene networks with incomplete measurements. We introduce a network reconstruction algorithm and a code base for reconstructing the dynamical structure function from data, to enable discovery and visualization of graphical relationships in a genetic circuit diagram as time-dependent functions rather than static, unknown weights. We prove a theorem, showing that dynamical structure functions can provide a data-driven estimate of the size of crosstalk fluctuations from an idealized model. We illustrate this idea with numerical examples. Finally, we show how data-driven estimation of dynamical structure functions can explain failure modes in two experimentally implemented genetic circuits, a historical genetic circuit and a new E. coli based transcriptional event detector. Fatal COVID-19 ARDS associated with incomplete AEC1 differentiation from the transitional state without senescence or fibrosis ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis. Intrinsic cardiac adrenergic cells contribute to septic cardiomyopathy Occurring independently of cardiac sympathetic nervous system, the intrinsic cardiac adrenergic (ICA) cells have been identified as an important regulator in both of developing and adult cardiac physiological and pathological processes. However, its role in septic cardiomyopathy remains unknown. Herein, we report that lipopolysaccharide (LPS) dose- and time-dependently increased norepinephrine (NE) release from ICA cells, which aggravates myocardial TNF- production and dysfunction. Inhibition of NE synthesis in ICA cells alleviated LPS-elicited cardiac dysfunction as well as TNF- production in Langendorff perfusing hearts. Mechanistically, ICA cell expressed Toll-like receptor 4 (TLR4), activated by LPS, to increase the expression of tyrosine hydroxylase, a key enzyme responsible for NE biosynthesis, via AP-1 binding to its promoter. Surprisingly, LPS-TLR4 signaling triggered no TNF- production in ICA cells due to the elevated Nfkbia and Tnfaip6 expression. In LPS-treated co-culture of ICA cells and cardiomyocytes, the raised NE from ICA cells activated cardiomyocyte {beta}1-adrenergic receptor ({beta}1-AR), driving Ca2+/calmodulin-dependent protein kinase II (CaMKII) to increase the activities of NF-{kappa}B and mitogen-activated protein kinase pathways, which were mimicked by dobutamine. Our findings reveal a cell type-specific TLR4 function triggering NE synthesis, but not TNF- production in inflammatory pathogenesis, and identify ICA cell-derived NE as a paracrine signal in the cross talk among different cardiac cells to enhance myocardial injury during LPS challenge, suggesting that targeting ICA cell-derived NE may be a potential therapeutic strategy for septic cardiomyopathy. The genome of the cereal pest Sitophilus oryzae: a transposable element haven BackgroundThe rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

ResultsWe sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 Kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino-acids required for insect development and cuticle biosynthesis.

ConclusionsHere we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes. A Substrate-triggered μ-Peroxodiiron(III) Intermediate in the 4-Chloro-L-Lysine-Fragmenting Heme-Oxygenase-like Diiron Oxidase (HDO) BesC: Substrate Dissociation from, and C4 Targeting by, the Intermediate The enzyme BesC from the {beta}-ethynyl-L-serine biosynthetic pathway in Streptomyces cattleya fragments 4-chloro-L-lysine (produced from L-Lysine by BesD) to ammonia, formaldehyde, and 4-chloro-L-allylglycine and can analogously fragment L-Lys itself. BesC belongs to the emerging family of O2-activating non-heme-diiron enzymes with the "heme-oxygenase-like" protein fold (HDOs). Here we show that binding of L-Lys or an analog triggers capture of O2 by the proteins diiron(II) cofactor to form a blue {micro}-peroxodiiron(III) intermediate analogous to those previously characterized in two other HDOs, the olefin-installing fatty acid decarboxylase, UndA, and the guanidino-N-oxygenase domain of SznF. The [~] 5- and [~] 30-fold faster decay of the intermediate in reactions with 4-thia-L-Lys and (4RS)-chloro-DL-lysine than in the reaction with L-Lys itself, and the primary deuterium kinetic isotope effects (D-KIEs) on decay of the intermediate and production of L-allylglycine in the reaction with 4,4,5,5-[2H]-L-Lys, imply that the peroxide intermediate or a successor complex with which it reversibly interconverts initiates the oxidative fragmentation by abstracting hydrogen from C4. Surprisingly, the sluggish substrate L-Lys can dissociate after triggering the intermediate to form, thereby allowing one of the better substrates to bind and react. Observed linkage between Fe(II) and substrate binding suggests that the triggering event involves a previously documented (in SznF) ordering of the dynamic HDO architecture that contributes one of the iron sites, a hypothesis consistent with the observation that the diiron(III) product cluster produced upon decay of the intermediate spontaneously degrades, as it has been shown to do in all other HDOs studied to date. Disruption of miR-18a alters proliferation, photoreceptor replacement kinetics, inflammatory signaling and microglia/macrophage numbers during retinal regeneration in zebrafish In mammals, photoreceptor loss causes permanent blindness, but in zebrafish (Danio rerio), photoreceptor loss reprograms Muller glia to function as stem cells, producing progenitors that regenerate photoreceptors. MicroRNAs (miRNAs) regulate CNS neurogenesis, but the roles of miRNAs in injury-induced neuronal regeneration are largely unknown. In the embryonic zebrafish retina, miR-18a regulates photoreceptor differentiation. The purpose of the current study was to determine, in zebrafish, the function of miR-18a during injury-induced photoreceptor regeneration. RT-qPCR, in-situ hybridization and immunohistochemistry showed that miR-18a expression increases throughout the retina between 1 and 5 days post-injury (dpi). To test miR-18a function during photoreceptor regeneration, we used homozygous miR-18a mutants (miR-18ami5012), and knocked down miR-18a with morpholino oligonucleotides. During photoreceptor regeneration, miR-18ami5012 retinas have fewer mature photoreceptors than WT at 7 and 10 dpi, but there is no difference at 14 dpi, indicating that photoreceptor regeneration is delayed. Labeling dividing cells with 5-bromo-2'-deoxyuridine (BrdU) showed that at 7 and 10 dpi, there are excess dividing progenitors in both mutants and morphants, indicating that miR-18a negatively regulates injury-induced proliferation. Tracing 5-ethynyl-2'-deoxyuridine (EdU) and BrdU-labeled cells showed that in miR-18ami5012 retinas excess progenitors migrate to other retinal layers in addition to the photoreceptor layer. Inflammation is critical for photoreceptor regeneration, and RT-qPCR showed that in miR-18ami5012 retinas, inflammatory gene expression and microglia activation are prolonged. Suppressing inflammation with dexamethasone rescues the miR-18ami5012 phenotype. Together, these data show that in the injured zebrafish retina, disruption of miR-18a alters proliferation, inflammation, the microglia/macrophage response, and the timing of photoreceptor regeneration. 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. Colicin E1 opens its hinge to plug TolC The double membrane architecture of Gram-negative bacteria forms a barrier that is effectively impermeable to extracellular threats. Bacteriocin proteins evolved to exploit the accessible, surface-exposed proteins embedded in the outer membrane to deliver cytotoxic cargo. Colicin E1 is a bacteriocin produced by, and lethal to, Escherichia coli that hijacks the outer membrane proteins TolC and BtuB to enter the cell. Here we capture the colicin E1 translocation domain inside its membrane receptor, TolC, by high-resolution cryoEM, the first reported structure of a bacteriocin bound to TolC. Colicin E1 binds stably to TolC as an open hinge through the TolC pore--an architectural rearrangement from colicin E1s unbound conformation. This binding is stable in live E. coli cells as indicated by single-molecule fluorescence microscopy. Finally, colicin E1 fragments binding to TolC plugs the channel, inhibiting its native efflux function as an antibiotic efflux pump and heightening susceptibility to three antibiotic classes. In addition to demonstrating that these protein fragments are useful starting points for developing novel antibiotic potentiators, this method could be expanded to other colicins to inhibit other outer membrane protein functions. PRSS2 stimulates tumor growth by remodeling the TME via repression of Tsp1 In the earliest stages of tumor development, epithelial tumors (carcinomas) are physically confined to the area of the tissue in which they form. These nascent lesions (carcinomas in situ) are sequestered from the tissue parenchyma by the basement membrane. Within the tissue parenchyma lie a myriad of cell types comprised of fibroblasts, immune and inflammatory cells and endothelial cells. Upon invasion across the basement membrane and into the tissue parenchyma, tumors must manipulate the expression of pro- and anti-tumorigenic proteins such that pro-tumorigenic factors are produced in excess to anti-tumorigenic proteins. One such anti-tumorigenic protein is Thrombospondin-1 (Tsp-1). We have previously demonstrated that stimulation of Tsp-1 in the tumor microenvironment (TME) potently inhibits tumor growth and progression and in some cases induces tumor regression. Here, we identify a novel tumor-mediated mechanism to repress the expression of Tsp-1 in the TME via secretion of the serine protease PRSS2. We demonstrate that PRSS2 represses Tsp-1, not via its enzymatic activity, but by binding to low-density lipoprotein receptor-related protein 1 (LRP1). These findings describe a novel activity for PRSS2 through binding to LRP1 and represent a potential therapeutic strategy to treat cancer by blocking the PRSS2-mediated repression of Tsp-1. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Electrocommunication signals indicate motivation to compete during dyadic interactions of an electric fish Animals across species compete for limited resources. While in some species competition behavior is solely based on own abilities, others assess their opponents to facilitate these interactions. Using cues and communication signals, contestants gather information about their opponent, adjust their behavior accordingly, and can thereby avoid high costs of escalating fights. We tracked electrocommunication signals, in particular "rises", and agonistic behaviors of the gymnotiform electric fish Apteronotus leptorhynchus in staged competition experiments. A larger body size relative to the opponent was the sole significant predictor for winners. Sex and the frequency of the continuously emitted electric field were only mildly influencing competition outcome. In males, correlations of body-size and winning were stronger than in females and, especially when losing against females, communication and agonistic interactions were enhanced, hinting towards males being more motivated to compete. Fish that lost competitions emitted the majority of rises in dependence on the competitors relative size and sex. The emission of rises was costly since it provoked ritualized biting or chasing behaviors by the other fish. Despite winners being accurately predictable based on rise numbers already after the initial 25 minutes, losers continued to emit rises. The number of rises emitted by losers and the duration of chasing behaviors depended in similar ways on physical attributes of contestants. The detailed evaluation of these correlations hint towards A. leptorhynchus adjusting their competition behavior according to mutual assessment, where rises could signal a losers motivation to continue assessment through ritualized fighting. Mouse oocytes do not contain a Balbiani body Oocytes spend the majority of their lifetime in a primordial, dormant state. Unlike many somatic cell types and mature oocytes, the cellular and molecular biology of primordial oocytes is largely unexplored. Yet, studying their cellular biology is necessary to understand the mechanisms through which oocytes maintain cellular fitness for decades, and why they eventually fail with age.

A hallmark of primordial oocytes in many species is the Balbiani body, a non-membrane bound compartment that contains the majority of mitochondria in the oocyte cytoplasm. The Balbiani body has been proposed to be essential for maintaining mitochondria in a healthy state during long-lasting dormancy, however, the architecture and function of the mammalian Balbiani body remains unknown.

Here, we develop enabling methods for live-imaging based comparative characterization of Xenopus, mouse and human primordial oocytes. We show that primordial oocytes in all three vertebrate species contain active mitochondria, Golgi apparatus and lysosomes. We further demonstrate that human and Xenopus oocytes have a Balbiani body characterized by a dense accumulation of mitochondria in their cytoplasm. However, despite previous reports, we did not find a Balbiani body in mouse oocytes. Instead, we demonstrate what was previously used as a marker for the Balbiani body in mouse primordial oocytes is in fact a ring-shaped Golgi apparatus that is not functionally associated with oocyte dormancy. Our work provides the first insights into the organisation of the cytoplasm in mammalian primordial oocytes, and clarifies relative advantages and limitations of choosing different model organisms for studying oocyte dormancy. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. spatialGE: Quantification and visualization of the tumor microenvironment heterogeneity using spatial transcriptomics SummarySpatially-resolved transcriptomics promises to increase our understanding of the tumor microenvironment and improve cancer prognosis and therapies. Nonetheless, analytical methods to explore associations between the spatial heterogeneity of the tumor and clinical data are not available. Hence, we have developed spatialGE, a software that provides visualizations and quantification of the tumor microenvironment heterogeneity through gene expression surfaces, spatial heterogeneity statistics (SThet) that can be compared against clinical information, spot-level cell deconvolution, and spatially-informed clustering (STclust), all using a new data object to store data and resulting analyses simultaneously.

Availability and implementationThe R package and tutorial/vignette are available at https://github.com/FridleyLab/spatialGE. A script to reproduce the analyses in this manuscript is available in Supplementary information.

ContactFridley.Lab@Moffitt.org or Brooke.Fridley@Moffitt.org

Supplementary informationAvailable at Bioinformatics online.



O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/454023v2_figG1.gif" ALT="Figure 1">
View larger version (73K):
org.highwire.dtl.DTLVardef@191cbe7org.highwire.dtl.DTLVardef@e1b6c7org.highwire.dtl.DTLVardef@a6f68forg.highwire.dtl.DTLVardef@1853601_HPS_FORMAT_FIGEXP  M_FIG O_FLOATNOGraphical abstractC_FLOATNO Overview of spatialGE features. A. The STList data object from spatialGE can be creared from several sources, including comma- or tab-separated files containing gene counts and spatial coordinates. The object can also be created directly from Visium outputs, Seurat objects, or GeoMx outputs. B. Users can optionally provide a metadata file, containing information associated with each sample (one row per sample, or per ROI if GeoMx data). C. Methods for quality control of data are provided by spatialGE, including visualizations of counts and genes per spot, as well as filtering of spots or genes within user-determined thresholds. D. A novel method (STclust) performs spatially informed clustering of spots and tissue domain identification. E. spatialGE provides different types of data visualization, including gene expression at each spot ("quilt plots"), as well as adaptation of spatial interpolation ("kriging") to spatial transcriptomics data (transcriptomic surface). F. spatialGE also leverages spatial statistics (Morans I, Gearys C, Getis-Ord Gi) to quantitatively describe heterogeneity within the tumor microenvironment and to explore associations between spatial heterogeneity and clinical oucomes. G. Gene expression deconvolution can also be applied to each spot to detect immune cell types (xCell) and classification of spots as tumor or stroma (ESTIMATE).

C_FIG Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans. The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP3 produced by PI3Ks. Excess PIP3 promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP2 and PIP3. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP3. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP3 allosterically unfolding the N-terminal PIP2 binding domain, positioning it favorably for the polybasic motif interaction with PIP2. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer. Regulatory perturbations of ribosome allocation optimize the growth proteome with a trade-off in adaptation capacity Bacteria regulate their cellular resource allocation to enable fast growth-adaptation to a variety of environmental niches. We studied the ribosomal allocation, growth and expression profiles of two sets of fast-growing mutants of Escherichia coli K-12 MG1655 in glucose minimal medium. Mutants with only 3 of the seven copies of ribosomal RNA operons grew faster than the wild-type strain in minimal media and show similar phenotype to previously studied rpoB mutants. Higher growth rates due to increased ribosome content affected resource allocation. Expression profiles of fast-growing mutants shared downregulation of hedging functions and upregulated growth functions. Mutants showed longer diauxic shifts and reduced activity of gluconeogenic promoters during glucose-acetate shifts, suggesting reduced availability of the RNA Polymerase for expressing hedging proteome. These results show that the regulation of ribosomal allocation underlies the growth/hedging phenotypes obtained from laboratory evolution experiments. We show how two different regulatory perturbations (rRNA promoters or rpoB mutations) reshape the proteome for growth with a concomitant fitness cost

HighlightsMutants with only 3 ribosomal operons grow faster than wild-type in minimal medium

 {Delta}4 rrn and rpoB mutants share phenotypic traits

Faster growth of mutants is achieved by increased ribosome content

Fast-growing mutants display reduced hedging expression and adaptation trade-offs

Despite similar ribosomal content in rich medium the mutants present growth defects Glypiated FGF directs contact-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. In vitro maturation of Toxoplasma gondii bradyzoites in human myotubes and their metabolomic characterization The apicomplexan parasite Toxoplasma gondii forms bradyzoite-containing tissue cysts that cause chronic and drug-tolerant infections. However, current in vitro models do not allow long-term culture of these cysts to maturity. Here, we developed a human myotube-based in vitro culture model of functionally mature tissue cysts that are orally infectious to mice and tolerate exposure to a range of antibiotics and temperature stresses. Metabolomic characterization of purified cysts reveals global changes that comprise increased levels of amino acids and decreased abundance of nucleobase- and tricarboxylic acid cycle-associated metabolites. In contrast to fast replicating tachyzoite forms of T. gondii these tissue cysts tolerate exposure to the aconitase inhibitor sodium fluoroacetate. Direct access to persistent stages of T. gondii under defined cell culture conditions will be essential for the dissection of functionally important host-parasite interactions and drug evasion mechanisms. It will also facilitate the identification of new strategies for therapeutic intervention. Methylome dynamic upon proteasome inhibition by the Pseudomonas syringae virulence factor Syringolin A DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, where modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effectors molecule, several of which act as proteasome inhibitors. Here we investigated the effect of proteasome inhibition by the bacterial virulence factor Syringolin A on genome-wide DNA methylation. We show that Syringolin A treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH DMRs that are enriched in the proximity of transcriptional start sites. Syringolin A treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway: AGO3, AGO9 and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of a epi-genomic arms race against pathogens. Are microbes growing on flowers evil? Effects of old flower microbes on fruit set in a wild ginger with one-day flowers, Alpinia japonica (Zingiberaceae) Flowers are colonized and inhabited by diverse microbes. Plants rapidly replace flowers with short lifespan, and old flowers senesce. This may contribute to avoiding adverse effects of the microbes. In this study, we investigate if the flower microbial community on old flowers impedes fruit and seed production in a wild ginger with one-day flowers. We inoculated newly opened flowers with old flower microbes, and monitored the effects on fruit and seed set. We also assessed prokaryotic communities on the flowers using amplicon sequencing. We found six bacterial amplicon sequence variants (ASVs) whose proportions were increased on the inoculated flowers. These ASVs were also found on flower buds and flowers that were bagged by net or paper during anthesis, suggesting that they had been present in small numbers prior to flowering. Fruit set was negatively associated with the proportions of these ASVs, while seed set was not. The results suggest that old flowers harbor microbial communities different from those at anthesis, and that the microbes abundant on old flowers negatively affect plant reproduction. Though the short lifespan of flowers has gotten little attention, it might be an essential defense mechanism to cope with antagonistic microbes that rapidly proliferate on the flowers.