Species-specific developmental timing dictates expansion of the avian wing skeletal pattern A fundamental question in biology is how species-specific rates of embryonic development are controlled in differently sized species. To address this problem, we compared wing development in the quail and the larger chick. We reveal that pattern formation (humerus to digits) is faster in the quail wing as determined by the earlier activation of 5 Hox genes, termination of developmental organisers (Shh and Fgf8) and the laying down of the skeleton (Sox9). Using interspecies tissue grafts, we show that developmental timing can be reset during a critical window in which retinoic acid signalling is active. Accordingly, retinoic acid can switch developmental timing between the quail and chick, and also between the chick and the larger turkey. We show that the incremental growth rate is equivalent between all three species, suggesting that differences in the expansion of the skeletal pattern are governed primarily by the tempo of development. Our findings have implications for how development is timed both within and between different species.
Loss of Resf1 reduces the efficiency of embryonic stem cell self-renewal and germline entry Retroelement silencing factor 1 (RESF1) interacts with the key regulators of mouse embryonic stem cells (ESCs) OCT4 and NANOG, and its absence results in sterility of mice. However, the function of RESF1 in ESCs and germ line specification is poorly understood. In this study, we used Resf1 knockout cell lines to determine the requirements of RESF1 for ESCs self-renewal and for in vitro specification of ESCs into primordial germ cell-like cells (PGCLCs). We found that deletion of Resf1 in ESCs cultured in serum and LIF reduces self-renewal potential whereas episomal expression of RESF1 has a modest positive effect on ESC self-renewal. In addition, RESF1 is not required for the capacity of NANOG and its downstream target ESRRB to drive self-renewal in the absence of LIF. However, Resf1 deletion reduces efficiency of PGCLC differentiation in vitro. These results identify Resf1 as a novel player in the regulation of pluripotent stem cells and germ cell specification.
Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. NOMA: a high-throughput microchip for robust, sequential measurements of secretions from the same single-cells Despite recent advances in single-cell analysis technologies, lacking simple methods to keep the live single-cells traceable for longitudinal detection reliably poses a significant obstacle in single-cell secretion analysis. Here we reported the high-density NOMA (narrow-opening microwell array) microchip that realized the retention of[&ge;] 97% of trapped single cells in dedicated spatial locations during repetitive detection procedures, verified with both adherent and suspension cells by two researchers independently. We applied it in monitoring single-cell protein secretions sequentially from the same single cells, and we found the digital protein secretion patterns dominate the protein secretion. We also demonstrated the microchip for longitudinally tracking IL-8 and the CD81+EV secretions from the same single-cells over days, which revealed the presence of "super secretors" within the cell population be more persistent to secrete protein or extracellular vesicle for an extended period. The NOMA platform reported here is simple, robust, and easy to operate for tracking sequential measurements from the same single cells, representing a novel and informative tool to inspire new observations in biomedical research.

Table of Contents Graphic

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/427049v4_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@aa69c7org.highwire.dtl.DTLVardef@6bc79org.highwire.dtl.DTLVardef@3fcbbforg.highwire.dtl.DTLVardef@597f4b_HPS_FORMAT_FIGEXP  M_FIG C_FIG Antibody-induced glomerulonephritis is amplified by RTEC-intrinsic IL-17 signaling and restrained by IL-17-mediated induction of the endoribonuclease Regnase-1 (Zc3h12a) Antibody-mediated glomerulonephritis (AGN) is a clinical manifestation of many autoimmune kidney diseases for which few effective treatments exist. Chronic inflammatory circuits in renal glomerular and tubular cells lead to tissue damage in AGN. These cells are targeted by the cytokine IL-17, which has recently been shown to be a central driver of the pathogenesis of AGN. However, surprisingly little is known about the regulation of pathogenic IL-17 signaling in the kidney. Here, using a well characterized mouse model of AGN, we show that IL-17 signaling in renal tubular epithelial cells (RTECs) is necessary for AGN development. We also show that Regnase-1, an RNA binding protein with endoribonuclease activity, is a negative regulator of IL-17 signaling in RTECs. Accordingly, mice with a selective Regnase-1 deficiency in RTECs exhibited exacerbated kidney dysfunction in AGN. Mechanistically, Regnase-1 inhibits IL-17-driven expression of the transcription factor I{kappa}B{xi} and consequently its downstream gene targets including Il6 and Lcn2. Moreover, deletion of Regnase-1 in human RTECs reduced inflammatory gene expression in an I{kappa}B{xi}-dependent manner. Overall, these data identify an IL-17-driven inflammatory circuit in RTECs during AGN that is constrained by Regnase-1. NPEPPS regulates intracellular import and sensitivity to cisplatin by interaction with volume regulated anion channels Despite routine use of platinum-based chemotherapeutics in cancer treatment, there remains a need to improve efficacy and patient selection. Multi-omic assessment of human bladder cancer cell lines and their cisplatin resistant derivatives and whole-genome CRISPR screens identified NPEPPS, the puromycin-sensitive aminopeptidase as a novel driver of cisplatin resistance. NPEPPS depletion increased cisplatin import and sensitization of resistant cells in vitro and in vivo. Pharmacologic inhibition of NPEPPS in cells and chemoresistant, patient-derived tumor organoids improved response to cisplatin. NPEPPS was found in complex with volume regulated anion channel (VRAC) subunits LRRC8A and LRRC8D, whose loss is known to enhance resistance to cisplatin. Depletion of LRRC8A the only obligate subunit for normal VRAC function abrogated the effect of NPEPPS-mediated cisplatin import. Our findings describe the first mechanism by which VRACs can be targeted for therapeutic benefit.

HIGHLIGHTSO_LIMulti-omic screening found NPEPPS is a driver of cisplatin resistance in cancer
C_LIO_LINPEPPS is in protein complex with LRRC8A and D, volume regulated anion channel (VRAC) subunits
C_LIO_LILRRC8A and D loss increases cisplatin resistance and reduce its intracellular levels
C_LIO_LINPEPPS regulates cisplatin import and sensitivity through VRACs
C_LIO_LINPEPPS function is inhibited by small molecules offering potential for clinical translation
C_LI 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. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. Ineffectual AEC1 differentiation from KRT8hi transitional state without fibrosis is associated with fatal COVID-19 ARDS 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. DetecDiv, a generalist deep-learning platform for automated cell division tracking and survival analysis Automating the extraction of meaningful temporal information from sequences of microscopy images represents a major challenge to characterize dynamical biological processes. So far, strong limitations in the ability to quantitatively analyze single-cell trajectories have prevented large-scale investigations to assess the dynamics of entry into replicative senescence in yeast. Here, we have developed DetecDiv, a microfluidic-based image acquisition platform combined with deep learning-based software for high-throughput single-cell division tracking. We show that DetecDiv can automatically reconstruct cellular replicative lifespans with high accuracy and performs similarly with various imaging platforms and geometries of microfluidic traps. In addition, this methodology provides comprehensive temporal cellular metrics using time-series classification and image semantic segmentation. Last, we show that this method can be further applied to automatically quantify the dynamics of cellular adaptation and the real-time cell survival upon exposure to environmental stress. Hence, this methodology provides an all-in-one toolbox for high-throughput phenotyping for cell cycle, stress response, and replicative lifespan assays. XRN2-mediated regulation of let-7 microRNAs is of critical pathophysiological significance in Humans MicroRNAs are critical regulators of diverse developmental and physiological processes, and they themselves get regulated both at the level of biogenesis and turnover. We demonstrate that the ribonuclease XRN2 can degrade the mature forms of certain let-7 family members in multiple human cancer cell lines, without affecting their precursors. XRN2 depletion results in a reduction in the expression of a number of oncogenes, and diminishes the proliferative and metastatic potential of cancer cells. The clinical relevance of these observations is also verified in tumour transcriptomics data from public RNA-sequencing datasets, where XRN2 mRNA expression is inversely correlated with the levels of a large number of miRNAs, including let-7 members, and high XRN2 mRNA levels are associated with poor survival in hepatocellular carcinoma, lung adenocarcinoma, and glioblastoma. We demonstrate that the miRNA is released by an as-yet unidentified proteinaceous  miRNA release factor from the grasp of Argonaute before its degradation. Our analyses of the patient-derived transcriptomics data also show that XRN2, via its regulation of let-7, affects multiple pathways in a consistent manner across epithelial and glial cell lineages, and thus, is of critical pathophysiological significance. Stimulated oxygen transport in tissue by magnetic needle acupuncture AimTo understand physiology of magnetic needle acupuncture by investigating O2 transport in tissue during needle intervention.

MethodsO2 transport in tissue is modeled by utilizing COMSOL with magnetic needle inserted into muscle tissue in a 2D porous media. The damaged tissue has been mimicked by an extracted tissue block with 1st order O2 consumption rate. The convection-diffusion O2 transport in the damaged tissue has been further evaluated by varying magnetic flux density B0 of the needle (0-1 T), myoglobin concentration (0-1 mM), O2 tension (5-100 Torr), O2 consumption rates, tissue permeability (10-12-10-6 m2) and porosity (0.1-0.9).

Results1) Active O2 transport carried by interstitial flow is enhanced with the intervention of a magnetic needle by generating a high gradient magnetic field around the tip, which exerts a strong force (104 N/m3) on the diamagnetic interstitial fluid to accelerate the flow. 2) This interstitial flow can reach 30 m/s at B0 = 1T and strongly correlates to the needle B0 and tissue permeability. 3) The needle stirs the interstitial flow can pump O2 flux by 1-2 orders of magnitude compared to that without magnetic field. 4) The enhancement of active O2 transport by magnetic needle is site-specific to the tissue in the vicinity of the tip. This enhancement is more effective in edema condition with a high tissue permeability (>10-9 m2).

ConclusionsThe dramatic enhancement of O2 transport to restore the O2 mitochondria metabolism for dysfunctional muscle tissues is the fundamental physiological mechanism of magnetic needle acupuncture. Phenotypic plasticity is aligned with phenological adaptation on both micro- and macroevolutionary timescales In seasonally-variable environments, phenotypic plasticity in phenology may be critical for adaptation to fluctuating environmental (temperature) conditions. Using an 18-generation longitudinal dataset from natural insect (damselfly) populations, we show that phenology has strongly advanced. Individual fitness data suggest this is likely an adaptive response towards a thermally-dependent fitness optimum. A laboratory experiment revealed that developmental plasticity qualitatively matches environmental-dependence of selection, partially explaining observed phenological advance. Expanding our analysis to the macroevolutionary level, we use a database of over 1-million occurrence records and spatiotemporally-matched temperature data from 49 Swedish Odonate species to infer macroevolutionary dynamics of phenology. Phenological plasticity was more closely aligned with adaptation for species that have recently colonized northern latitudes, with more phenological mismatch at lower latitudes. Our results show that phenological plasticity plays a key role in microevolutionary dynamics within a single species, and such plasticity may have facilitated post-Pleistocene range expansion in this insect clade. Membrane-assisted assembly and selective autophagy of enteroviruses Enteroviruses are non-enveloped positive-sense RNA viruses that cause diverse diseases in humans. Their rapid multiplication depends on remodeling of cytoplasmic membranes for viral genome replication. It is unknown how virions assemble around these newly synthesized genomes and how they are then loaded into autophagic membranes for release through secretory autophagy. Here, we use cryo-electron tomography of infected cells to show that poliovirus assembles directly on replication membranes. Pharmacological untethering of capsids from membranes abrogates RNA encapsidation. Our data directly visualize a membrane-bound half-capsid as a prominent virion assembly intermediate. Assembly progression past this intermediate depends on the class III phosphatidylinositol 3-kinase VPS34, a key host-cell autophagy factor. On the other hand, the canonical autophagy initiator ULK1 is shown to restrict virion production since its inhibition leads to increased accumulation of virions in vast intracellular arrays, followed by an increased vesicular release at later time points. Finally, we identify multiple layers of selectivity in virus-induced autophagy, with a strong selection for RNA-loaded virions over empty capsids and the segregation of virions from other types of autophagosome contents. These findings provide an integrated structural framework for multiple stages of the poliovirus life cycle. The visual cortex produces gamma band echo in response to broadband visual flicker The aim of this study is to uncover the network dynamics of the human visual cortex by driving it with a broadband random visual flicker. We here applied a broadband flicker (1-720 Hz) while measuring the MEG and then estimated the temporal response function (TRF) between the visual input and the MEG response. This TRF revealed an early response in the 40-60 Hz gamma range as well as in the 8-12 Hz alpha band. While the gamma band response is novel, the latter has been termed the alpha band perceptual echo. The gamma echo preceded the alpha perceptual echo. The dominant frequency of the gamma echo was subject-specific thereby reflecting the individual dynamical properties of the early visual cortex. To understand the neuronal mechanisms generating the gamma echo, we implemented a pyramidal-interneuron gamma (PING) model that produces gamma oscillations in the presence of constant input currents. Applying a broadband input current mimicking the visual stimulation allowed us to estimate TRF between the input current and the population response (akin to the local field potentials). The TRF revealed a gamma echo that was similar to the one we observed in the MEG data. Our results suggest that the visual gamma echo can be explained by the dynamics of the PING model even in the absence of sustained gamma oscillations.

Author SummaryThe properties of the neuronal dynamics governing the visual system are highly debated. While some emphasize the neuronal firing rate and evoked activity in response to visual stimuli, others emphasize the oscillatory neuronal dynamics. To investigate the dynamical properties of the visual system, we recorded the magnetoencephalography while stimulating the visual system using a broadband (1-720 Hz) visual flicker. By estimating the temporal response function (similar to cross-correlation) between the visual input and neuronal activity, we demonstrated a clear response in the gamma band that we term the gamma echo. We then constructed a physiologically realistic network model that could generate gamma-band oscillations by a pyramidal-interneuron gamma (PING) mechanism. This model allowed us to account for empirically observed response in the gamma band, and to provide novel insight on the neuronal dynamics governing the early visual system. The stage is now set for further investigating how the gamma echo is modulated by tasks such as spatial attention as well as uncovering how the echo might propagate in the visual hierarchy. Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Nej1XLF interacts with Sae2CTIP to inhibit Dna2 mediated resection at DNA doublestrand break The two major pathways of DNA double strand break (DSB) repair, non-homologous end-joining (NHEJ) and homologous recombination (HR), are highly conserved from yeast to mammals. The regulation of 5 DNA resection controls repair pathway choice and influences repair outcomes. Nej1 was first identified as a canonical NHEJ factor involved in stimulating the ligation of broken DNA ends and more recently, it was shown to be important for DNA end-bridging and inhibiting 5 resection mediated by Dna2-Sgs1. Nej1 interacts with Sae2 and this impacts DSB repair in three ways. First, Nej1 inhibits MRX-Sae2 interactions and Sae2 localization to a DSB. Second, Nej1 inhibits Sae2-dependent recruitment of Dna2 in the absence of Sgs1. Third, NEJ1 and SAE2 showed an epistatic relationship for DNA end-bridging, an event that restrains the broken ends and reduces the frequency of genomic deletions from developing at the DSB. Deletion of NEJ1 suppressed the synthetic lethality of sae2{Delta} sgs1{Delta} and was dependent on the nuclease activity of Dna2. These Nej1 functions promote end-joining DSB repair, but could also be relevant for controlling resection initiation during HR repair.

HighlightsO_LINej1 physically interacts with Sae2 and inhibits end-resection at a DSB.
C_LIO_LINej1 inhibits Sae2 interactions with the MRX complex.
C_LIO_LINej1 inhibits Sae2-dependent recruitment of Dna2 to a DSB.
C_LIO_LINEJ1 and SAE2 are epistatic for DNA end-bridging.
C_LIO_LIDeletion of NEJ1 suppresses the synthetic lethality of sae2{Delta} sgs1{Delta}, which is dependent on Dna2 nuclease activity.
C_LI Temperature Increase Drives Critical Slowing Down of Fish Ecosystems Fish ecosystems perform ecological functions that are critically important for the sustainability of marine ecosystems, such as global food security and carbon stock. During the 21st century, significant global warming caused by climate change has created pressing challenges for fish ecosystems that threaten species existence and global ecosystem health. Here, we study a coastal fish community in Maizuru Bay, Japan, and investigate the relationships between fluctuations of ST, abundance-based species interactions and salient fish biodiversity. Observations show that a local 20% increase in temperature from 2002 to 2014 underpins a long-term reduction in fish diversity ([~]25%) played out by some native and invasive species (e.g. Chinese wrasse) becoming exceedingly abundant; this causes a large decay in commercially valuable species (e.g. Japanese anchovy) coupled to an increase in ecological productivity. The fish community is analyzed considering five temperature ranges to understand its atemporal seasonal sensitivity to ST changes, and long-term trends. An optimal information flow model is used to reconstruct species interaction networks that emerge as topologically different for distinct temperature ranges and species dynamics. Networks for low temperatures are more scale-free compared to ones for intermediate (15-20{degrees}C) temperatures in which the fish ecosystem experiences a first-order phase transition in interactions from locally stable to metastable and globally unstable for high temperatures states as suggested by abundance-spectrum transitions. The dynamic dominant eigenvalue of species interactions shows increasing instability for competitive species (spiking in summer due to intermediate-season critical transitions) leading to enhanced community variability and critical slowing down despite higher time-point resilience. Native competitive species whose abundance is distributed more exponentially have the highest total directed interactions and are keystone species (e.g. Wrasse and Horse mackerel) for the most salient links with cooperative decaying species. Competitive species, with higher eco-climatic memory and synchronization, are the most affected by temperature and play an important role in maintaining fish ecosystem stability via multitrophic cascades (via cooperative-competitive species imbalance), and as bioindicators of change. More climate-fitted species follow temperature increase causing larger divergence divergence between competitive and cooperative species. Decreasing dominant eigenvalues and lower relative network optimality for warmer oceans indicate fishery more attracted toward persistent oscillatory states, yet unpredictable, with lower cooperation, diversity and fish stock despite the increase in community abundance due to non-commercial and venomous species. We emphasize how changes in species interaction organization, primarily affected by temperature fluctuations, are the backbone of biodiversity dynamics and yet for functional diversity in contrast to taxonomic richness. Abundance and richness manifest gradual shifts while interactions show sudden shift. The work provides data-driven tools for analyzing and monitoring fish ecosystems under the pressure of global warming or other stressors. Abundance and interaction patterns derived by network-based analyses proved useful to assess ecosystem susceptibility and effective change, and formulate predictive dynamic information for science-based fishery policy aimed to maintain marine ecosystems stable and sustainable. In vitro culture of freshly isolated Trypanosoma brucei brucei bloodstream forms results in gene copy-number changes. Most researchers who study unicellular eukaryotes work with an extremely limited number of laboratory-adapted isolates that were obtained from the field decades ago, but the effects of passage in laboratory rodents, and adaptation to in vitro culture, have been little studied. For example, the vast majority of studies of Trypanosoma brucei biology have concentrated on just two strains, Lister 427 and EATRO1125, which were taken from the field over half a century ago and have since have undergone innumerable passages in rodents and culture. We here describe two new Trypanosoma brucei brucei strains. MAK65 and MAK98, which have undergone only 3 rodent passages since isolation from Ugandan cattle. High-coverage sequencing revealed that adaptation of the parasites to culture was accompanied by changes in gene copy numbers. T. brucei has so far been considered to be uniformly diploid, but we also found trisomy of chromosome 5 not only in one Lister 427 culture, but also in the MAK98 field isolate. Trisomy of chromosome 6, and increased copies of other chromosome segments, were also seen in established cultured lines. The two new T. brucei strains should be useful to researchers interested in trypanosome differentiation and pathogenicity. Initial results suggested that the two strains have differing infection patterns in rodents. MAK65 is uniformly diploid and grew more reproducibly in bloodstream-form culture than MAK98. Curation, inference, and assessment of a globally reconstructed gene regulatory network for Streptomyces coelicolor Streptomyces coelicolor A3(2) is a model microorganism for the study of Streptomycetes, antibiotic production, and secondary metabolism in general. Even though S. coelicolor has an outstanding variety of regulators among bacteria, little effort to globally study its transcription has been made. We manually curated 29 years of literature and databases to assemble a meta-curated experimentally-validated gene regulatory network (GRN) with 5386 genes and 9707 regulatory interactions ([~]41% of the total expected interactions). This provides the most extensive and up-to-date reconstruction available for the regulatory circuitry of this organism. Only [~]6% (534/9707) are supported by experiments confirming the binding of the transcription factor to the upstream region of the target gene, the so-called "strong" evidence. While for the remaining interactions there is no confirmation of direct binding. To tackle network incompleteness, we performed network inference using several methods (including two proposed here) for motif identification in DNA sequences and GRN inference from transcriptomics. Further, we contrasted the structural properties and functional architecture of the networks to assess the reliability of the predictions, finding the inference from DNA sequence data to be the most trustworthy approach. Finally, we show two applications of the inferred and the curated networks. The inference allowed us to propose novel transcription factors for the key Streptomyces antibiotic regulatory proteins (SARPs). The curated network allowed us to study the conservation of the system-level components between S. coelicolor and Corynebacterium glutamicum. There we identified the basal machinery as the common signature between the two organisms. The curated networks were deposited in Abasy Atlas (https://abasy.ccg.unam.mx/) while the inferences are available as Supplementary Material. Crowdsourced identification of multi-target kinase inhibitors for RET- and TAU- based disease: the Multi-Targeting Drug DREAM Challenge A continuing challenge in modern medicine is the identification of safer and more efficacious drugs. Precision therapeutics, which have one molecular target, have been long promised to be safer and more effective than traditional therapies. This approach has proven to be challenging for multiple reasons including lack of efficacy, rapidly acquired drug resistance, and narrow patient eligibility criteria. An alternative approach is the development of drugs that address the overall disease network by targeting multiple biological targets ( polypharmacology). Rational development of these molecules will require improved methods for predicting single chemical structures that target multiple drug targets. To address this need, we developed the Multi-Targeting Drug DREAM Challenge, in which we challenged participants to predict single chemical entities that target pro-targets but avoid anti-targets for two unrelated diseases: RET-based tumors and a common form of inherited Tauopathy. Here, we report the results of this DREAM Challenge and the development of two neural network-based machine learning approaches that were applied to the challenge of rational polypharmacology. Together, these platforms provide a potentially useful first step towards developing lead therapeutic compounds that address disease complexity through rational polypharmacology.

Author SummaryMany modern drugs are developed with the goal of modulating a single cellular pathway or target. However, many drugs are, in fact,  dirty; they target multiple cellular pathways or targets. This phenomenon is known as multi-targeting or polypharmacology. While some strive to develop  cleaner therapeutics that eliminate secondary targets, recent work has shown that multi-targeting therapeutics have key advantages for a variety of diseases. However, while multi-targeting drugs that affect a precisely-defined profile of targets may be more effective, it is difficult to computationally predict which molecules have desirable target profiles. Here, we report the results of a competitive crowdsourcing project (the Multi-Targeting Drug DREAM Challenge), where we challenged participants to predict chemicals that have desired target profiles for cancer and neurodegenerative disease. Cell-specific Bioorthogonal Tagging of Glycoproteins Altered glycosylation is an undisputed corollary of cancer development. Understanding these alterations is paramount but hampered by limitations underlying cellular model systems. For instance, the intricate interactions between tumour and host cannot be adequately recapitulated in monoculture of tumour-derived cell lines. More complex co-culture models usually rely on sorting procedures for proteome analyses and rarely capture the details of protein glycosylation. Here, we report a strategy termed Bio-Orthogonal Cell line-specific Tagging of Glycoproteins (BOCTAG). Cells are equipped by transfection with an artificial biosynthetic pathway that transforms bioorthogonally tagged sugars into the corresponding nucleotide-sugars. Only transfected cells incorporate bioorthogonal tags into glycoproteins in the presence of non-transfected cells. We employ BOCTAG as an imaging technique and to annotate cell-specific glycosylation sites in mass spectrometry-glycoproteomics. We demonstrate application in co-culture and mouse models, allowing for profiling of the glycoproteome as an important modulator of cellular function. Defining An Expanded RAS Conformational Landscape Based on Over 700 Experimentally Determined Structures of KRAS, NRAS, and HRAS For many human cancers and tumor-associated diseases, mutations in the RAS isoforms (KRAS, NRAS, and HRAS) are the most common oncogenic alterations, making these proteins high-priority therapeutic targets. Effectively targeting the RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structure-guided catalogue of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. We present an expanded classification of RAS conformations based on analyzing their catalytic switch 1 (SW1) and switch 2 (SW2) loops. From all 721 available human KRAS, NRAS, and HRAS structures in the Protein Data Bank (PDB) (206 RAS-protein complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of residue Y32 in SW1 and residue Y71 in SW2. Subsequently, we defined additional conformational subsets (some previously undescribed) by clustering all well modeled SW1 and SW2 loops using a density-based machine learning algorithm with a backbone dihedral-based distance metric. In all, we identified three SW1 conformations and nine SW2 conformations, each which are associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation can be further subdivided based on the hydrogen (H)-bond type made between residue Y32 and the GTP {gamma}-phosphate: water-mediated, direct, or no H-bond. Further analyzing these structures clarified the catalytic impact of the G12D and G12V RAS mutations, and the inhibitor chemistries that bind to each druggable RAS conformation. To facilitate future RAS structural analyses, we have created a web database, called Rascore, presenting an updated and searchable dataset of human KRAS, NRAS, and HRAS structures in the PDB, and which includes a page for analyzing user uploaded RAS structures by our algorithm (http://dunbrack.fccc.edu/rascore/).

SignificanceAnalyzing >700 experimentally determined RAS structures helped define an expanded landscape of active, inactive and druggable RAS conformations, the structural impact of common RAS mutations, and previously uncharacterized RAS-inhibitor binding modes. Plant-environment microscopy tracks interactions of Bacillus subtilis with plant roots across the entire rhizosphere Our understanding of plant-microbe interactions in soil is limited by the difficulty of observing processes at the microscopic scale throughout plants large volume of influence. Here, we present the development of 3D live microscopy for resolving plant-microbe interactions across the environment of an entire seedling growing in a transparent soil in tailor-made mesocosms, maintaining physical conditions for the culture of both plants and microorganisms. A tailor made dual-illumination light-sheet system acquired scattering signals from the plant whilst fluorescence signals were captured from transparent soil particles and labelled microorganisms, allowing the generation of quantitative data on samples approximately 3600 mm3 in size with as good as 5 m resolution at a rate of up to one scan every 30 minutes. The system tracked the movement of Bacillus subtilis populations in the rhizosphere of lettuce plants in real time, revealing previously unseen patterns of activity. Motile bacteria favoured small pore spaces over the surface of soil particles, colonising the root in a pulsatile manner. Migrations appeared to be directed towards the root cap, the point "first contact", before subsequent colonisation of mature epidermis cells. Our findings show that microscopes dedicated to live environmental studies present an invaluable tool to understand plant-microbe interactions. Assessing the Role of Calmodulin's Linker Flexibility in Target Binding 1Calmodulin (CaM) is a universal Ca2+ binding protein known to bind at least 300 targets. The selectivity and specificity towards these targets are partially attributed to the proteins flexible alpha-helical linker that connects its N- and C-domains. How this flexible linker mediates the driving forces guiding CaMs binding to regulatory targets is not well-established. Therefore, we utilized the Martini coarse-grained (CG) molecular dynamics simulations to probe interrelationships between CaM/target assembly and the role of its linker region. As a model system, we simulated the binding of CaM to the CaM binding region (CaMBR) of calcineurin (CaN). The simulations were conducted assuming a  wild-type calmodulin with normal flexibility of its linker, as well as a labile, highly flexible linker variant. For the wild-type model, 98% of the 600 simulations across three ionic strengths adopted a bound complex within 2 {micro}s of simulation time; of these, 1.7% sampled the fully-bound state observed in experimentally-determined crystallographic structure. By calculating the mean-first-passage-time for these simulations, we estimated the association rate to be ka = 5.9 x 108 M-1 s-1, which is similar to the experimentally-determined rate of 2.2 x 108 M-1 s-1 [1]. Further, our simulations recapitulated the inverse relationship between the association rate and solution ionic strength reported in the literature. In contrast, although over 97% of the labile linker simulations formed tightly-bound complexes, only 0.3% achieved the fully-bound configuration. This effect appears to stem from a difference in the ensembles of extended and collapsed states controlled by the linker properties. Specifically, the labile linker variant samples fewer extended states compatible with target peptide binding. Therefore, our simulations suggest that variations in the CaM linkers propensity for alpha-helical secondary structure can modulate the kinetics of target binding. This finding is important, as the linker region houses several CaM variants sites for post-translational modifications, that may alter the proteins normal regulatory functions.