Essential Regression - a generalizable framework for inferring causal latent factors from multi-omic human datasets High-dimensional cellular and molecular profiling of human samples highlights the need for analytical approaches that can integrate multi-omic datasets to generate predictive biomarkers and prioritized causal inferences. Current methods are limited by high dimensionality of the combined datasets, the differences in their data distributions and their integration to infer causal relationships. Here we present Essential Regression (ER), an interpretable machine learning approach for high-dimensional multi-omic datasets, that addresses these problems by identifying latent factors and their likely cause-effect relationships with the system-wide outcome/properties of interest. ER is a novel data-distribution-free latent-factor regression model that integrates multi-omic datasets and identifies latent factors significantly associated with an outcome. ER outperforms a range of state-of-the-art methods in terms of prediction performance on simulated datasets. ER can be coupled with probabilistic graphical modeling thereby strengthening the causal inferences. ER generates novel cellular and molecular predictions, using multi-omic human systems immunology datasets, pertaining to immunosenescence and immune dysregulation.
Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments.
Implications of the (co)evolution of agriculture and resource foraging for the maintenance of species diversity and community structure Agriculture is found in numerous taxa such as humans, ants, beetles, fishes and even bacteria. This type of niche construction has evolved independently from hunting, though many species remain primarily predators. When a consumer has a positive effect on its resource, we can expect an allocative cost of agriculture, as the agricultural care diverts time and energy from other activities. Defending the resource against predators may divert time from its consumption (exploitation cost). The cost may also occur on the foraging of alternative resources, for instance if the consumer spends more time nearby the farmed resource and underexploiting resources elsewhere (opportunity cost). We here investigate transitions from predation to agriculture in a simple three-species model of a farmer that consumes two resources and has a positive effect on one. We study the conditions for the (co)evolution of the investment into agriculture and specialization on the two resources, and its consequences on the ecological dynamics of the community. Eco-evolutionary dynamics generate a feedback between the evolution of agriculture and specialization on the helped resource, that can lead to varying selected intensity of agriculture, from generalist strategies with no agriculture, to specialist farmers, with possible coexistence between these two extreme strategies. 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. "Neighbourhood watch" model: embryonic epiblast cells assess positional information in relation to their neighbours In many developing and regenerating systems, tissue pattern is established through gradients of informative morphogens, but we know little about how cells interpret these. Using experimental manipulation of early chick embryos including misexpression of an inducer (VG1 or ACTIVIN) and an inhibitor (BMP4), we test two alternative models for their ability to explain how the site of primitive streak formation is positioned relative to the rest of the embryo. In one model, cells read morphogen concentrations cell-autonomously. In the other, cells sense changes in morphogen status relative to their neighbourhood. We find that only the latter model can account for the experimental results, including some counter-intuitive predictions. This mechanism (which we name "neighbourhood watch" model) illuminates the classic "French Flag Problem" and how positional information is interpreted by a sheet of cells in a large developing system.

Summary statementIn a large developing system, the chick embryo before gastrulation, cells interpret gradients of positional signals relative to their neighbours to position the primitive streak, establishing bilateral symmetry. Early epigenetic influence on EGR1, KLF2 and KLF4 transcriptional network associates with alteration of HSCs quiescence Excessive fetal growth is associated with DNA methylation alterations in human hematopoietic stem and progenitor cells (HSPC), but their functional impact remains elusive. We implemented an integrative analysis combining single-cell epigenomics, single-cell transcriptomics, and in vitro analyses to functionally link DNA methylation changes to putative alterations of HSPC functions. We showed in hematopoietic stem cells (HSC) from large for gestational age neonates that both DNA hypermethylation and chromatin rearrangement target a specific network of transcription factors known to sustain stem cell quiescence. In parallel, we found a decrease expression of key genes regulating HSC differentiation including EGR1, KLF2, SOCS3, and JUNB. Our functional analyses showed that this epigenetic programming was associated with a decreased ability for HSCs to stay quiescent. Taken together, our multimodal approach using single-cell (epi)genomics showed that human fetal overgrowth affects hematopoietic stem cells quiescence maintenance via epigenetic programming. Time- and cell-dependent atypia and cell death are caused by progressive deficiency in DNA replication Pleiotropy caused by single-gene mutations is common and poorly understood. A zebrafish null mutant of DNA polymerase  subunit B, huli hutu (hht), evolves a complex pleiotropy associated with DNA damage and S phase arrest across multiple organ systems over 5-7 days, including nuclear atypia, a common cellular feature in human cancers and pre-cancers, in gastrointestinal organs, and nuclear fragmentation in the eye and brain. The pleiotropic pattern of hht phenotypes is explained by progressive loss of wild-type maternal pola2 function in homozygous mutant embryos whose pola2 mRNA becomes undetectable by 24 hours post-fertilization (hpf). Inhibition of DNA synthesis by aphidicolin or hydroxyurea in wild-type embryos from 24 hpf phenocopied the pleiotropic pattern of hht. These results are consistent with a model in which time-sensitive, reduced capacity for DNA synthesis results in cell death in fast-replicating cells, and nuclear atypia in tissues with fewer and larger cells. Analysis of science journalism reveals gender and regional disparities in coverage Science journalism is a critical way in which the public can remain informed and benefit from new scientific findings. Such journalism also shapes the publics view of the current state of scientific findings and legitimizes experts. Those covering science can only cite and quote a limited number of sources. Sources may be identified by the journalists research or by recommendations by other scientists. In both cases, biases may influence who is identified and ultimately included as an expert. To analyze possible biases in science journalism, we analyzed 22,001 non-research articles published by Nature. We chose to analyze Nature non-research articles since its research articles provide a natural comparator. Our analysis considered two possible sources of disparity: gender and name origin. To explore these sources of disparity, we extracted cited authors names as well as extracted names of quoted speakers. While citations and quotations within a piece do not reflect the entire information-gathering process, they can provide insight into the demographics of visible sources. We then used the extracted names to predict gender and name origin of the cited authors and speakers.

In order to appropriately quantify the level of difference, we must identify a suitable reference set for comparison. We chose first and last authors within primary research articles in Nature and a subset of Springer Nature articles in the same time period as our comparator. In our analysis, we found a skew towards male quotation in Nature science journalism-related articles. However, quotation is trending toward equal representation at a faster rate than first and last authorship in academic publishing. Interestingly, we found that the gender disparity in Nature quotes was column-dependent, with the "Career Features" column reaching gender parity. Our name origin analysis found a significant over-representation of names with predicted Celtic/English origin and under-representation of names with a predicted East Asian origin. This finding was observed both in extracted quotes and journal citations, but dampened in citations. The combined treatment of Molnupiravir and Favipiravir results in a marked potentiation of efficacy in a SARS-CoV2 hamster infection model through an increased frequency of mutations in the viral genome. Favipiravir and Molnupiravir, orally available antivirals, have been reported to exert antiviral activity against SARS-CoV2. In recent days preliminary efficacy data have been reported in COVID-19 patients. We here studied the combined antiviral effect of the drugs in the SARS-CoV2 hamster infection model. We first demonstrate that Molnupiravir can reduce infectious virus titers in lungs of infected animals in a dose-dependent manner by up to 3.5 log10 which is associated with a marked improvement of virus-induced lung pathology. When animals are treated with a combination of suboptimal doses of Molnupiravir and Favipiravir (that each alone result in respectively a 1.3 log10 and 1.1 log10 reduction of infectious virus titers in the lungs), a marked combined potency is observed. Infectious virus titers in the lungs of animals treated with the combo are on average reduced by 4.5 log10 and infectious virus are no longer detected in the lungs of 60% of treated infected animals. Both drugs result in an increased mutation frequency of the remaining viral RNA recovered from the lungs. In the combo-treated hamsters an increased frequency of C-to-T and G-to-A mutations in the viral RNA is observed as compared to the single treatment groups which may explain the pronounced antiviral potency of the combination. Our findings may lay the basis for the design of clinical studies to test the efficacy of the combination of Molnupiravir and Favipiravir in the treatment of COVID-19. 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. Getting more out of FLAG-Tag co-immunoprecipitation mass spectrometry experiments using FAIMS. Co-immunoprecipitation of proteins coupled to mass spectrometry is critical for the understanding of protein interaction networks. In instances where a suitable antibody is not available, it is common to graft synthetic tags onto target protein sequences and allowing the use of commercially available antibodies for affinity purification. A common approach is through FLAG-Tag co-immunoprecipitation. To allow the selective elution of protein complexes, competitive displacement using a large molar excess of the tag peptides is often carried out. Yet, this creates downstream challenges for the mass spectrometry analysis due to the presence of large quantities of these peptides. Here, we demonstrate that Field Asymmetric Ion Mobility Spectrometry (FAIMS), a gas phase ion separation device prior to mass spectrometry analysis can be applied to FLAG-Tag co-immunoprecipitation experiment to increase the depth of protein coverage. By excluding these abundant tag peptides, we were able to observe deeper coverage of interacting proteins and as a result, deeper biological insights, without the need for additional sample handling or altering sample preparation protocols. Loss of Nuclear DNA Ligase III Reverts PARP Inhibitor Resistance in BRCA1/53BP1 Double-deficient Cells by Exposing ssDNA Gaps Inhibitors of poly(ADP-ribose) (PAR) polymerase (PARPi) have entered the clinic for the treatment of homologous recombination (HR)-deficient cancers. Despite the success of this approach, preclinical and clinical research with PARPi has revealed multiple resistance mechanisms, highlighting the need for identification of novel functional biomarkers and combination treatment strategies. Functional genetic screens performed in cells and organoids that acquired resistance to PARPi by loss of 53BP1 identified loss of LIG3 as an enhancer of PARPi toxicity in BRCA1-deficient cells. Enhancement of PARPi toxicity by LIG3 depletion is dependent on BRCA1 deficiency but independent of the loss of 53BP1 pathway. Mechanistically, we show that LIG3 loss promotes formation of MRE11-mediated post-replicative ssDNA gaps in BRCA1-deficient and BRCA1/53BP1 double-deficient cells exposed to PARPi, leading to an accumulation of chromosomal abnormalities. LIG3 depletion also enhances efficacy of PARPi against BRCA1-deficient mammary tumors in mice, suggesting LIG3 as a potential therapeutic target. Collective polymerase dynamics emerge from DNA supercoiling during transcription All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts and transcription stops. DNA supercoiling can also mediate non-local interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton) whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.

PACS numbers: mTORC2 contributes to murine systemic autoimmunity The development of many systemic autoimmune diseases, including systemic lupus erythematosus, is associated with overactivation of the type I interferon (IFN) pathway, lymphopenia, and increased follicular helper T (Tfh) cell differentiation. However, the cellular and molecular mechanisms underlying these immunological perturbations remain incompletely understood. Here we show that the mechanistic target of rapamycin complex 2 (mTORC2) promotes Tfh differentiation and disrupts Treg homeostasis. Inactivation of mTORC2 in total T cells, but not in Tregs, greatly ameliorated the immunopathology in a systemic autoimmunity mouse model. This was associated with reduced Tfh differentiation, B cell activation, and reduced T cell glucose metabolism. Finally, we show that type I IFN can synergize with TCR ligation to activate mTORC2 in T cells, which partially contributes to T cell lymphopenia. These data indicate that mTORC2 may act as downstream of type I IFN, TCR, and costimulatory receptor ICOS, to promote glucose metabolism, Tfh differentiation, and T cell lymphopenia, but not to suppress Treg function in systemic autoimmunity. Our results suggest that mTORC2 might be a rational target for systemic autoimmunity treatment. A Novel Proposal for an Index for Regional Cerebral Perfusion Pressure - A Theoretical Approach Using Fluid Dynamics Cerebral blood flow (CBF) / cerebral blood volume (CBV) ratio derived by [15O] H2O/ CO2 and CO positron emission tomography (PET) examination has been used empirically as an index for cerebral perfusion pressure (CPP). However, this equation lacks a theoretical basis and furthermore, has not been confirmed to be proportionate to CPP, as direct measurement of local CPP is not practical.

We have developed a new index for CPP using the Poiseuille equation based on a simple model. Our model suggests that CBF/CBV2 is proportionate to CPP and that it is mathematically a more accurate index than CBF/CBV. 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. Non-specific effects of a CINNAMATE-4-HYDROXYLASE inhibitor on auxin homeostasis Chemical inhibitors are often implemented for the functional characterization of genes to overcome the limitations associated with genetic approaches. Although being a powerful tool, off-target effects of these inhibitors are easily overlooked in a complex biological setting. Here we illustrate the implications of such secondary effects by focusing on piperonylic acid (PA), an inhibitor of CINNAMATE-4-HYDROXYLASE (C4H) that is often used to investigate the involvement of lignin during plant growth and development. When supplied to plants, we found that PA is recognized as a substrate by GRETCHEN HAGEN 3.6 (GH3.6), an amido synthetase involved in the formation of the auxin catabolite indole-3-acetic acid (IAA)-Asp. By competing for the same enzyme, PA interferes with auxin conjugation, resulting in an increase in cellular auxin concentrations. These increased auxin levels likely further contribute to an increase in adventitious rooting previously observed upon PA-treatment. Despite the focus on GH3.6 in this report, PA is conjugated by an array of enzymes and their subsequent reduced activity on native substrates could potentially affect a whole set of physiological processes in the plant. We conclude that surrogate occupation of the endogenous conjugation machinery in the plant by exogenous compounds is likely a more general phenomenon that is rarely considered in pharmacological studies. Our results hereby provide an important basis for future reference in studies using chemical inhibitors. Distributions of extinction times from fossil ages and tree topologies: the example of mid-Permian synapsid extinctions Given a phylogenetic tree that includes only extinct, or a mix of extinct and extant taxa, where at least some fossil data are available, we present a method to compute the distribution of the extinction time of a given set of taxa under the Fossilized-Birth-Death model.

Our approach differs from the previous ones in that it takes into account (i) the possibility that the taxa or the clade considered may diversify before going extinct and (ii) the whole phylogenetic tree to estimate extinction times, whilst previous methods do not consider the diversification process and deal with each branch independently. Because of this, our method can estimate extinction times of lineages represented by a single fossil, provided that they belong to a clade that includes other fossil occurrences.

We assess and compare our new approach with a standard previous one using simulated data. Results show that our method provides more accurate confidence intervals.

This new approach is applied to the study of the extinction time of three Permo-Carboniferous synapsid taxa (Ophiacodontidae, Edaphosauridae, and Sphenacodontidae) that are thought to have disappeared toward the end of the Cisuralian (early Permian), or possibly shortly thereafter. The timing of extinctions of these three taxa and of their component lineages supports the idea that the biological crisis in the late Kungurian/early Roadian consisted of a progressive decline in biodiversity throughout the Kungurian. Structural effects driven by rare point mutations in amylin hormone, the type II diabetes-associated peptide BackgroundAmylin is a 37-amino-acid peptide hormone co-secreted with insulin, which participates in glucose homeostasis. This hormone is able to aggregate in a {beta}-sheet conformation and deposit in islet amyloids, a hallmark in type II diabetes. Since amylin is a gene-encoded hormone, this peptide has variants caused by point mutations that can impact its functions.

MethodsHere, we analyzed the structural effects caused by S20G and G33R point mutations which, according to the 1000 Genomes Project, have frequency in East Asian and European populations, respectively. The analyses were performed by means of aggrescan server, SNP functional effect predictors, and molecular dynamics.

ResultsWe found that both mutations have aggregation potential and cause changes in the monomeric forms when compared with wild-type amylin. Furthermore, comparative analyses with pramlintide, an amylin drug analogue, allowed us to infer that second -helix maintenance may be related to the aggregation potential.

ConclusionsThe S20G mutation has been described as pathologically related, which is in agreement with our findings. In addition, our data suggest that the G33R mutation might have a deleterious effect. The data presented here also provide new therapy opportunities, whether for creating more effective drugs for diabetes or implementing specific treatment for patients with these mutations.

General SignificanceOur data could help to better understand the impact of mutations on the wild-type amylin sequence, as a starting point for the evaluation and characterization of other variations. Moreover, these findings could improve the health of patients with type II diabetes. A liquid +TIP-network regulates microtubule dynamics through tubulin condensation Regulation of microtubule dynamics is essential for diverse cellular functions, and proteins that bind to dynamic microtubule ends can regulate network dynamics. Here we show, that two conserved microtubule end-binding proteins, CLIP-170 and EB3, undergo phase separation and form dense liquid-networks. When CLIP-170 and EB3 act together the multivalency of the network increases, which synergistically increases the amount of protein in the dense phase. In vitro and in cells these liquid networks can condense tubulin. In vitro in the presence of microtubules, EB3/CLIP-170 phase separation can co-condense tubulin all along the microtubule. At this condition microtubule growth speed increases up to two-fold and depolymerization events are strongly reduced, compared to conditions with phase separation deficient networks. Our data show that phase separated EB3/CLIP-170 networks impact microtubule growth dynamics beyond direct protein-microtubule interactions. RNF213 loss of function reshapes vascular transcriptome and spliceosome leading to disrupted angiogenesis and aggravated vascular inflammatory responses. RationaleMoyamoya disease (MMD) is a rare cerebrovascular occlusive disease that affects Asian population more often. The pathogenesis of MMD is related to mutation in RNF213 gene. However, why, and how RNF213 mutation leads to MMD is still not fully understood.

ObjectiveAnalyze the impact of RNF213 loss of function on vascular cells and the observed changes correlate with MMD.

Methods and resultsRNF213 KD was conducted in HUVEC and vascular smooth muscle cells (vSMCs). First, HUVEC cells showed alteration of angiogenesis, migration under LPS stimulation, Leukocyte endothelial transmigration, and endothelial-to-vSMCs communication.

Transcriptome analysis revealed downregulation of genes regulating cell division and mitosis. Interestingly, Alternative splicing (AS) analysis revealed hundreds of AS events to occur after RNF213 KD in various types of AS. The alternatively spliced genes showed minimal overlap with transcriptome profiling. Pathway analysis revealed many processes and pathways to be regulated by AS events observed. Transcriptome profiling was also performed after LPS treatment and revealed the basis for increased sensitivity to LPS observed in our analysis. AS changes were also observed after LPS treatment in RNF213 KD HUVEC. Different types of AS showed different patterns of convergence or divergence in terms of the regulated pathways after LPS treatment. Finally, transcriptome and AS analysis was performed in vSMCs and showed various processes that impact vSMCs function and phenotype in RNF213 loss of function.

ConclusionOur data provide a wealth of information on RNF213 gene function in vascular tissues and shed important light on different processes that contribute to MMD pathogenesis. Our results signify an immune-driven and hemodynamically linked process of MMD initiation and progression. Royal knifefish generate powerful suction feeding through large neurocranial elevation and high epaxial muscle power Suction feeding in ray-finned fishes involves powerful buccal cavity expansion to accelerate water and food into the mouth. Previous XROMM studies in largemouth bass (Micropterus salmoides), bluegill sunfish (Lepomis macrochirus), and channel catfish (Ictalurus punctatus) have shown that more than 90% of suction power in high performance strikes comes from the axial musculature. Thus, the shape of the axial muscles and skeleton may impact suction feeding mechanics. Royal knifefish (Chitala blanci) have an unusual postcranial morphology, with a ventrally flexed vertebral column and relatively large mass of epaxial muscle. Based on their body shape, we hypothesized that royal knifefish would generate high power strikes by utilizing large neurocranial elevation, vertebral column extension, and epaxial shortening. As predicted, C. blanci generated high suction expansion power compared to the other three species studied to date (up to 160 W), which was achieved by increasing both the rate of volume change and the intraoral subambient pressure. The large epaxial muscle (25% of body mass) shortened at high velocities to produce large neurocranial elevation and vertebral extension (up to 41 deg, combined), as well as high muscle mass-specific power (up to 800 W kg-1). For the highest power strikes, axial muscles generated 95% of the power, and 64% of the axial muscle mass consisted of the epaxial muscles. The epaxial-dominated suction expansion of royal knifefish supports our hypothesis that postcranial morphology may be a strong predictor of suction feeding biomechanics.

SUMMARY STATEMENTRoyal knifefish rely on their distinct postcranial morphology--with a curved vertebral column and large dorsal body muscles--to produce large neurocranial elevation and powerful suction feeding. Critical re-evaluation of the slope factor of the operational model of agonism Background and PurposeAlthough being a relative term, agonist efficacy is a cornerstone in the proper assessment of agonist selectivity and signalling bias. The operational model of agonism (OMA) has become successful in the determination of agonist efficacies and ranking them. In 1985, Black et al. introduced the slope factor to the OMA to make it more flexible and allow for fitting steep as well as flat concentration-response curves.

Experimental ApproachFunctional analysis of OMA demonstrates that the slope factor implemented by Black et al. affects relationships among parameters of the OMA, leading to wrong estimates of operational efficacy and affinity of an agonist. In contrast, functional analysis of OMA with Hill coefficient as a slope factor has not this flaw.

Key ResultsFitting of the OMA with Black et al. slope factor to concentration-response curves of several models (homotropic allosteric modulation, substrate inhibition and non-competitive auto-inhibition) gives wrong estimates of operational efficacy and affinity. In contrast, fitting of the OMA modified by the Hill coefficient to concentration-response curves of these several models gives correct estimates of operational efficacy and affinity.

Conclusion and ImplicationsThe OMA modified by the Hill coefficient can be used for correct ranking of agonism and subsequent analysis, like quantification of signalling bias. Movement planning and feedback are sufficient for sensory prediction error computation and motor adaptation Prediction errors guide many forms of learning, providing teaching signals that help us improve our performance. Implicit motor adaptation, for instance, is driven by sensory prediction errors (SPEs), which occur when the expected and observed consequences of a movement differ. Traditionally, SPE computation is thought to require movement execution. However, recent work suggesting that the brain generates and accounts for sensory predictions based on motor imagery or planning alone calls this assumption into question. Here, by measuring implicit adaptation during a visuomotor task, we tested whether motor planning and well-timed sensory feedback are sufficient for SPE computation. Human participants were cued to reach to a target and were, on a subset of trials, rapidly cued to withhold these movements. Errors displayed both on trials with and without movements induced single-trial implicit learning. Learning following trials without movements persisted even when movement trials had never been paired with errors, and when the direction of movement and sensory feedback trajectories were decoupled. These observations demonstrate that the brain can compute SPEs without generating overt movements, leading to the adaptation of planned movements even when they are not performed. Differences in gait stability and acceleration characteristics between healthy young and older females Our aim was to evaluate differences in gait acceleration intensity, variability and stability of feet and trunk between older females and young females using inertial sensors. Twenty older females (OF; mean age 68.4, SD 4.1 years) and eighteen young females (YF; mean age 22.3, SD 1.7 years) were asked to walk straight for 100 meters at their preferred speed, while wearing inertial sensors on heels and lower back. We calculated spatiotemporal measures, foot and trunk acceleration characteristics and their variability, as well as trunk stability using the local divergence exponent (LDE). Two-way analysis of variance (including the factors foot and age), Students t-test, and Mann-Whitney U test were used to compare statistical differences of measures between groups. Cohens d effects were calculated for each variable. Foot maximum vertical acceleration and amplitude, trunk-foot vertical acceleration attenuation, as well as their variability were significantly smaller in OF than in YF. In contrast, trunk mediolateral acceleration amplitude, maximum vertical acceleration, and amplitude, as well as their variability were significantly larger in OF than in YF. Moreover, OF showed lower stability (i.e., higher LDE values) in mediolateral acceleration, mediolateral and vertical angular velocity of the trunk. Even though we measured healthy older females, these participants showed lower vertical foot accelerations with higher vertical trunk acceleration, lower trunk-foot vertical acceleration attenuation, less gait stability, and more variability of the trunk, and hence, were more likely to fall. These findings suggest that instrumented gait measurements may help for early detection of changes or impairments in gait performance, even before this can be observed by clinical eye or gait speed.