The elevational ascent and spread of exotic annual grass dominance in the Great Basin, USA AimIn the western US, sagebrush (Artemisia spp.) and salt desert shrublands are rapidly transitioning to communities dominated by exotic annual grasses, a novel and often self-reinforcing state that threatens the economic sustainability and conservation value of rangelands. Climate change is predicted to directly and indirectly favor annual grasses, potentially pushing transitions to annual grass dominance into higher elevations and north-facing aspects. We sought to quantify the expansion of annual grass-dominated vegetation communities along topographic gradients over the past several decades.

LocationOur analysis focused on rangelands among three ecoregions in the Great Basin of the western US, where several species of exotic annual grasses are widespread among shrub and perennial grass-dominated vegetation communities.

MethodsWe used recently developed remote sensing-based rangeland vegetation data to produce yearly maps of annual grass-dominated vegetation communities spanning the period 1990-2020. With these maps, we quantified the rate of spread and characterized changes in the topographic distribution (i.e., elevation and aspect) of areas transitioning to annual grass dominance.

ResultsWe documented more than an eight-fold increase in annual grass-dominated area (to >77,000 km2) occurring at an average rate of >2,300 km2 yr-1. In 2020, annual grasses dominated one fifth (19.8%) of Great Basin rangelands. This rapid expansion is associated with a broadening of the topographic niche, with widespread movement into higher elevations and north-facing aspects.

Main conclusionsAccelerated, strategic intervention is critically needed to conserve the fragile band of rangelands being compressed between annual grassland transitions at lower elevations and woodland expansion at higher elevations.
Multilevel multinomial models as a tool to extract population spatial networks from capture-recapture data O_LIEstimating the impacts of anthropogenic disturbances requires an understanding of the habitat use patterns of individuals within a population. This is especially the case when disturbances are localized within a populations spatial range, as variation in habitat-use within a population can drastically alter the distribution of impacts.
C_LIO_LIHere, we illustrate the potential for multilevel multinomial models to generate spatial networks from capture-recapture data, a common data source use in wildlife studies to monitor population dynamics and habitat use. These spatial networks capture which regions of a populations spatial distribution share similar/dissimilar individual usage patterns, and can be especially useful for detecting structured habitat use within the populations spatial range.
C_LIO_LIUsing simulations and 18 years of capture-recapture data from St. Lawrence Estuary (SLE) beluga, we show that this approach can successfully estimate the magnitude of similarities/dissimilarities in individual usage patterns across sectors, and identify sectors that share similar individual usage patterns that differ from other sectors, i.e., structured habitat use. In the case of SLE beluga, this method identified multiple clusters of individuals, each preferentially using restricted areas within their summer range of the SLE.
C_LIO_LISynthesis and applications. Multilevel multinomial models can be effective at estimating spatial structure in habitat use within wildlife populations sampled by capture-recapture of individuals. Our finding of a structured habitat use within the SLE beluga summer range has direct implications for estimating individual exposures to localized stressors, such as underwater noise from shipping or other activities.
C_LI
Mutations in coral soma and sperm imply lifelong stem cell renewal and cell lineage selection In many animals, the germline differentiates early in embryogenesis, so only mutations that accumulate in germ cells are inherited by offspring1. Exceptions to this developmental process may indicate that other mechanisms have evolved to limit the effects of deleterious mutation accumulation2. Stony corals are animals that can live for hundreds of years3 and have long been thought to produce gametes from somatic tissue4. To clarify conflicting evidence about germline-soma distinction in corals, we sequenced high coverage, full genomes with technical replicates for parent coral branches and their sperm pools. We identified post-embryonic single nucleotide variants (SNVs) unique to each parent branch, then checked if each SNV was shared by the respective sperm pool: 26% of post-embryonic SNVs were shared by the sperm but 74% were not. We also identified germline SNVs, those that were present in the sperm but not in the parent. These data suggest that self-renewing stem cells in corals differentiate into germ and soma throughout the adult life of the colony, with SNV rates and patterns differing markedly in stem, soma, and germ lineages. In addition to informing the important place in the evolutionary spectrum from non-Weismannian to Weismmanian animals that corals occupy, these insights inform how corals may generate adaptive diversity necessary in the face of global climate change. Creation of a single cell RNA seq meta-atlas to define human liver immune homeostasis The liver is unique in both its ability to maintain immune homeostasis and in its potential for immune tolerance following solid organ transplantation. Single-cell RNA sequencing (scRNA seq) is a powerful approach to generate highly dimensional transcriptome data to understand cellular phenotypes. However, when scRNA data is produced by different groups, with different data models, different standards, and samples processed in different ways, it can be challenging to draw meaningful conclusions from the aggregated data. The goal of this study was to establish a method to combine  human liver scRNA seq datasets by 1) characterizing the heterogeneity between studies and 2) using the meta-atlas to define the dominant phenotypes across immune cell subpopulations in healthy human liver. Publicly available scRNA seq data generated from liver samples obtained from a combined total of 17 patients and [~]32,000 cells were analyzed. Liver-specific immune cells (CD45+) were extracted from each dataset, and immune cell subpopulations (myeloid cells, NK and T cells, plasma cells, and B cells) were examined using dimensionality reduction (UMAP), differential gene expression, and ingenuity pathway analysis. All datasets co-clustered, but cell proportions differed between studies. Gene expression correlation demonstrated similarity across all studies, and canonical pathways that differed between datasets were related to cell stress and oxidative phosphorylation rather than immune-related function. Next, a meta-atlas was generated via data integration and compared against PBMC data to define gene signatures for each hepatic immune subpopulation. This analysis defined key features of hepatic immune homeostasis, with decreased expression across immunologic pathways and enhancement of pathways involved with cell death. This method for meta-analysis of scRNA seq data provides a novel approach to broadly define the features of human liver immune homeostasis. Specific pathways and cellular phenotypes described in this human liver immune meta-atlas provide a critical reference point for further study of immune mediated disease processes within the liver. 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 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. Feasibility of Administering Human Pancreatic Cancer Chemotherapy in a Spontaneous Pancreatic Cancer Mouse Model BackgroundBoth modified FOLFIRINOX (mFFX) and gemcitabine/nab-paclitaxel chemotherapy regimens have been shown to improve clinical outcomes in patients with pancreatic cancer, and are often used interchangeably as the standard of care. Preclinical studies often do not use these regimens, since administering these multiagent approaches can be difficult. In this study, we assessed the feasibility of administering these two chemotherapy regimens in spontaneous pancreatic tumors using KPC mice with the ultimate goal of advancing preclinical studies.

MethodsKPC mice were created by breeding KrasLSL-G12D/+ to Trp53fl/fl;Ptf1Cre/+, resulting in KrasLSL-G12D/+;p53fl/+;Ptf1Cre/+ mice. At 14 weeks of age, mice were palpated for spontaneous tumor growth that was verified using ultrasounds. Mice with tumors under 15 mm in diameter were used. The mice were assigned to one of seven treatment regimens: 1 cycle of mFOLFIRINOX (FFX X1), 2 cycles of mFOLFIRINOX (FFX X2), 1 cycle of mFOLFIRINOX with 40Gy SBRT (FFX SBRT), 1 cycle of gemcitabine/nab-paclitaxel (GEM/AB X1), 2 cycles of gemcitabine/nab-paclitaxel (GEM/AB X2), 2 cycles of gemcitabine/nab-paclitaxel with 40Gy SBRT (GEM/AB SBRT), or no chemotherapy treatment (control).

ResultsIn total, 92 mice were included. The median OS in the FFX X2 group was slightly longer that the median OS in the FFX X1 group (15 days vs 11 days, P=0.003). Mice in the GEM/AB X2 group had longer OS when compared to mice in the GEM/AB X1 group (33.5 vs 13 days, P=0.001). Mice treated with chemotherapy survived longer than untreated control animals (median OS: 6.5 days, P <0.001). Moreover, in mice treated with chemotherapy, mice that received 2 cycles of GEM/AB X2 had the longest survival, while the FFX X1 group had the poorest OS (P <0.001). Lastly, chemotherapy followed by consolidative SBRT trended towards increased local control and survival.

ConclusionsWe demonstrate the utility and feasibility of clinicall relevant mFOLFIRINOX and gemcitabine/nab-paclitaxel in preclinical models of pancreatic cancer. Feature-reweighted representational similarity analysis: A method for improving the fit between computational models, brains, and behavior Representational Similarity Analysis (RSA) has emerged as a popular method for relating representational spaces from human brain activity, behavioral data, and computational models. RSA is based on the comparison of representational (dis-)similarity matrices (RDM or RSM), which characterize the pairwise (dis-)similarities of all conditions across all features (e.g. fMRI voxels or units of a model). However, classical RSA treats each feature as equally important. This  equal weights assumption contrasts with the flexibility of multivariate decoding, which reweights individual features for predicting a target variable. As a consequence, classical RSA may lead researchers to underestimate the correspondence between a model and a brain region and, in case of model comparison, may lead them to select an inferior model. The aim of this work is twofold: First, we sought to broadly test feature-reweighted RSA (FR-RSA) applied to computational models and reveal the extent to which reweighting model features improves RSM correspondence and affects model selection. Previous work suggested that reweighting can improve model selection in RSA but it has remained unclear to what extent these results generalize across datasets and data modalities. To draw more general conclusions, we utilized a range of publicly available datasets and three popular deep neural networks (DNNs). Second, we propose voxel-reweighted RSA, a novel use case of FR-RSA that reweights fMRI voxels, mirroring the rationale of multivariate decoding of optimally combining voxel activity patterns. We found that reweighting individual model units markedly improved the fit between model RSMs and target RSMs derived from several fMRI and behavioral datasets and affected model selection, highlighting the importance of considering FR-RSA. For voxel-reweighted RSA, improvements in RSM correspondence were even more pronounced, demonstrating the utility of this novel approach. We additionally show that classical noise ceilings can be exceeded when FR-RSA is applied and propose an updated approach for their computation. Taken together, our results broadly validate the use of FR-RSA for improving the fit between computational models, brain, and behavioral data, possibly allowing us to better adjudicate between competing computational models. Further, our results suggest that FR-RSA applied to brain measurement channels could become an important new method to assess the correspondence between representational spaces. TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines. Identification of a conserved neutralizing epitope present on spike proteins from all highly pathogenic coronaviruses Three pathogenic human coronaviruses have emerged within the last 20 years, with SARS-CoV-2 causing a global pandemic. Although therapeutic antibodies targeting the SARS-CoV-2 spike currently focus on the poorly conserved receptor-binding domain, targeting essential neutralizing epitopes on the more conserved S2 domain may provide broader protection. We report an antibody binding an epitope conserved in the pre-fusion core of MERS-CoV, SARS-CoV and SARS-CoV-2 spike S2 domains. Antibody 3A3 binds a conformational epitope with ~2.5 nM affinity and neutralizes spike from SARS-CoV, SARS-CoV-2 and variants of concern in in vitro pseudovirus assays. Hydrogen-deuterium exchange mass spectrometry identified residues 980-1006 in the flexible hinge region at the S2 apex as the 3A3 epitope, suggesting 3A3 prevents the S2 conformational rearrangements required for conversion to the spike post-fusion state and virus-host cell fusion. This work defines a conserved vulnerable site on the SARS-CoV-2 S2 domain and guides the design of pan-protective spike immunogens. Dynamics of SARS-CoV-2 mutations reveals regional-specificity and similar trends of N501 and High-Frequency mutation N501Y in different levels of control measures Coronavirus disease 2019 (COVID-19) is a contagious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This disease has spread globally, causing more than 161.5 million cases and 3.3 million deaths to date. Surveillance and monitoring of new mutations in the virus genome are crucial to our understanding of the adaptation of SARS-CoV-2. Moreover, how the temporal dynamics of these mutations is influenced by control measures and non-pharmaceutical interventions (NPIs) is poorly understood. Using 1 058 020 SARS-CoV-2 from sequenced COVID-19 cases from 98 countries (totaling 714 country-month combinations), we perform a normalization by COVID-19 cases to calculate the relative frequency of SARS-CoV-2 mutations and explore their dynamics over time. We found 115 mutations estimated to be present in more than 3 % of global COVID-19 cases and determined three types of mutation dynamics: High-Frequency, Medium-Frequency, and Low-Frequency. Classification of mutations based on temporal dynamics enable us to examine viral adaptation and evaluate the effects of implemented control measures in virus evolution during the pandemic. We showed that Medium-Frequency mutations are characterized by high prevalence in specific regions and/or in constant competition with other mutations in several regions. Finally, taking N501Y mutation as representative of High-Frequency mutations, we showed that level of control measure stringency negatively correlates with the effective reproduction number of SARS-CoV-2 with High-Frequency or not-High-Frequency and both follows similar trends in different levels of stringency. TFAP2 paralogs facilitate chromatin access for MITF at pigmentation genes but inhibit expression of cell-cell adhesion genes independently of MITF In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that Tfap2 facilitates MITFs ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa- expressing cells that also express Tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the latter necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A coupled mechano-biochemical framework for root meristem morphogenesis Plants grow roots to adjust their bodies to dynamic changes in the surrounding environment. How do plant roots emerge from the stem cell reservoir during embryogenesis is however poorly understood. Here, we present a bottom-up strategy to address this challenge by combining empirical observations with advanced computer modeling techniques. We demonstrate that the anisotropy of root growth results from differential growth rates of adjacent tissues, whereas the root meristem development incorporates a multi-level feedback loop between complex transport network of phytohormone auxin, auxin-dependent cell growth and cytoskeleton rearrangements. In silico model predictions are in close agreement with in vivo patterns of anisotropic growth, auxin distribution, and cell polarity, as well as several root phenotypes caused by chemical, mechanical, or genetic perturbations. Our findings reveal a minimal set of design principles connecting tissue mechanics, cell anisotropy, and directional transport that are sufficient for self-organization of the root meristem shape. A mobile auxin signal transported through immobile cells orchestrates polarity and growth mechanics to instruct the morphogenesis of an independent organ. Unevolved proteins from modern and prebiotic amino acids manifest distinct structural profiles Natural proteins represent numerous but tiny structure/function islands in a vast ocean of possible protein sequences, most of which has not been explored by either biological evolution or research. Recent studies have suggested this uncharted sequence space possesses surprisingly high structural propensity, but development of an understanding of this phenomenon has been awaiting a systematic high-throughput approach.

Here, we designed, prepared, and characterized two combinatorial protein libraries consisting of randomized proteins, each 105 residues in length. The first library constructed proteins from the entire canonical alphabet of 20 amino acids. The second library used a subset of only 10 residues (A,S,D,G,L,I,P,T,E,V) that represent a consensus view of plausibly available amino acids through prebiotic chemistry. Our study shows that compact conformations resistant to proteolysis are (i) abundant (up to 40%) in random sequence space, (ii) independent of general Hsp70 chaperone system activity, and (iii) not granted solely by "late" and complex amino acid additions. The Hsp70 chaperone system effectively increases solubility and refoldability of the canonical alphabet but has only a minor impact on the "early" library. The early alphabet proteins are inherently more soluble and refoldable, possibly assisted by the cell-like environment in which these assays were performed.

Our work indicates that both early and modern amino acids are predisposed to supporting protein structure (either in forms of oligomers or globular/molten globule structures) and that protein structure may not be a unique outcome of evolution. AutoGGN: A Gene Graph Network AutoML Tool for Multi-Omics Research Omics data identifies biological characteristics from genetic to phenotypic levels during the life span. Molecular interaction networks have a fundamental impact on life activities. Integrating omics data and molecular interaction networks will help researchers delve into comprehensive information underlying the data. Here, we proposed a new multimodal method called AutoGGN to aggregate multi-omics data and molecular interaction networks based on graph convolutional neural networks. We evaluated AutoGGN using two different tasks: cancer type classification and single-cell stage classification. On both tasks, AutoGGN showed better performance compared to other methods, the trend is relevant to the ability of utilizing much more information from biological data. The phenomenon indicated AutoGGN has the potential to incorporate valuable information from molecular interaction networks and multi-omics data effectively. Furthermore, in order to provide a better understanding of the mechanism of prediction results, we assessed the explanation using SHAP module and identified the key genes contributing to the prediction of classification, which will provide insights for the downstream design of biological experiments. Phylogenetic analysis of molecular pathways in Tregs Suppressor Function Treg-mediated suppression of conventional T cells is a fundamental step in regulating the adaptive immune response. It is known that Treg first appeared in vertebrates. However, little is known about the of major suppression pathways mediated by Tregs. We employed artificial intelligence text mining system to highlight the suppression pathways currently known to be utilized by Tregs. Our system identified various pathways mediated by CTLA4, calcium signaling, NfkB and NFAT. We simultaneously employed detailed phylogenetic analysis including multiple sequence alignment, phylogenetic tree building, ancestral sequence reconstruction, neutrality tests and positive selection test to investigate the evolution of Treg mediated pathways. We found that CTLA4 first appeared in vertebrates possibly arising from an IGV containing protein in cartilaginous fish. Conversely, we found that Tregs repurposed ancient pathways such as Calcineurin and CAMP Response Element Modulator. Interestingly, these two pathways were highly conserved between vertebrates and lower invertebrates indicating conservation of function. Taken together, our research indicate that Tregs developed its regulatory systems through evolution in vertebrates as well as reusing conserved ancient regulatory systems that are related to the innate immune system. Maternal microbiota Bifidobacterium promotes placental morphogenesis, nutrient transport and fetal growth in mice The gut microbiota plays a central role in regulating host metabolism. However, while substantial progress has been made in discerning how the microbiota influences host functions post birth and beyond, little is known about how key members of the maternal gut microbiota can influence feto-placental growth. Here, using germ-free and specific-pathogen-free mice, we demonstrate that the bacterium Bifidobacterium breve UCC2003 modulates maternal body adaptations, placental structure and nutrient transporter capacity, with implications for fetal metabolism and growth. The effects of B. breve UCC2003 on feto-placental growth are mediated, in part, by changes in the maternal and placental metabolome (i.e. acetate and carnitine). Histological analysis of the placenta confirmed that Bifidobacterium modifies placental structure via changes in Igf2P0, Dlk1, Mapk1 and Mapk14 expression. Additionally, B. breve UCC2003, acting through Slc2a1 and Fatp1-4 transporters, was shown to restore fetal glycaemia and fetal growth in association with changes in the fetal hepatic transcriptome. This study provides knowledge towards a novel and safe therapeutic strategy for treating pregnancy disorders via modulation of the maternal gut microbiota. Single-molecule imaging of cytoplasmic dynein in cellulo reveals the mechanism of motor activation and cargo capture Cytoplasmic dynein 1 (dynein) is the primary minus end-directed motor protein in most eukaryotic cells. Dynein remains in an inactive conformation until the formation of a tripartite complex comprising dynein, its regulator dynactin and a cargo adaptor. How this process of dynein activation occurs is un-clear, since it entails the formation of a three-protein complex inside the crowded environs of a cell. Here, we employed live-cell, single-molecule imaging to visualise and track fluorescently tagged dynein. First, we observed that only [~]30% of dynein molecules that bound to the microtubule (MT) engaged in minus end-directed movement, and that too for a short duration of [~]0.6 s. Next, using high-resolution imaging in live and fixed cells, and using correlative light and electron microscopy, we discovered that dynactin remained persistently attached to MTs, and endosomal cargo remained in proximity to the MTs and dynactin. Finally, we employed two-colour imaging to visualise cargo movement effected by single motor binding. We then used these discoveries as the basis for a stochastic model incorporating dynamic motors binding to cargo located along MTs, and also developed a coarse-grained 3-state run-and-tumble particle (RTP) model for the cargo that quantitatively recapitulates the emergent statistics of cargo movement. Taken together, we discovered a search mechanism that is facilitated by dyneins frequent MT binding-unbinding kinetics: (1) in a futile event when dynein does not encounter cargo anchored in proximity to the MT, dynein dissociates and diffuses into the cytoplasm, (2) when dynein encounters cargo and dynactin upon MT-binding, it moves cargo in a short run. Several of these short runs are undertaken in succession for long-range directed movement. In conclusion, we demonstrate that dynein activation and cargo capture are coupled in a step that relies on the reduction of dimensionality to enable minus end-directed transport in cellulo, and that complex cargo behaviour emerges from stochastic motor-cargo interactions. Comparative toxicogenomics of glyphosate and Roundup herbicides by mammalian stem cell-based genotoxicity assays and molecular profiling in Sprague-Dawley rats Whether glyphosate-based herbicides (GBHs) are more potent than glyphosate alone at activating cellular mechanisms, which drive carcinogenesis remains controversial. As GBHs are more cytotoxic that glyphosate, we reasoned they may also be more capable of activating carcinogenic pathways. We tested this hypothesis by comparing the effects of glyphosate with Roundup GBHs both in vitro and in vivo. First, glyphosate was compared with representative GBHs namely MON 52276 (EU), MON 76473 (UK) and MON 76207 (USA) using the mammalian stem cell-based ToxTracker system. Here, MON 52276 and MON 76473, but not glyphosate and MON 76207, activated oxidative stress and unfolded protein responses. Second, molecular profiling of liver was performed in female Sprague-Dawley rats exposed to glyphosate or MON 52276 (both at 0.5, 50, 175 mg/kg bw/day glyphosate) for 90 days. MON 52276 but not glyphosate increased hepatic steatosis and necrosis. MON 52276 and glyphosate altered the expression of genes in liver reflecting TP53 activation by DNA damage and circadian rhythm regulation. Genes most affected in liver were similarly altered in kidneys. Small RNA profiling in liver showed decreased amounts of miR-22 and miR-17 from MON 52276 ingestion. Glyphosate decreased mir-30 while miR-10 levels were increased. DNA methylation profiling of liver revealed 5,727 and 4,496 differentially methylated CpG sites between the control and glyphosate and MON 52276 exposed animals respectively. Apurinic/apyrimidinic DNA damage formation in liver was increased with glyphosate exposure. Altogether, our results show that Roundup formulations cause more biological changes linked with carcinogenesis than glyphosate. Non-Newtonian Blood Rheology Impacts Left Atrial Stasis in Patient-Specific Simulations The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects can be subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered in calculating patient-specific blood stasis indices by computational fluid dynamics. Different contribution of the monkey prefrontal and premotor dorsal cortex in decision making during a transitive inference task Several studies have reported similar neural modulations between brain areas of the frontal cortex, such as the dorsolateral prefrontal (DLPFC) and the premotor dorsal (PMd) cortex, in tasks requiring encoding of the abstract rules for selecting the proper action. Here we compared the neuronal modulation of the DLPFC and PMd of monkeys trained to choose the higher rank from a pair of abstract images (target item), selected from an arbitrarily rank-ordered set (A>B>C>D>E>F) in the context of a transitive inference task. Once acquired by trial-and-error, the ordinal relationship between pairs of adjacent images (i.e., A>B; B>C; C>D; D>E; E>F), monkeys were tested in indicating the ordinal relation between items of the list not paired during learning. During these decisions, we observed that the choice accuracy increased and the reaction time decreased as the rank difference between the compared items enhanced. This result is in line with the hypothesis that after learning, the monkeys built an abstract mental representation of the ranked items, where rank comparisons correspond to the items position comparison on this representation. In both brain areas, we observed higher neuronal activity when the target item appeared in a specific location on the screen with respect to the opposite position and that this difference was particularly enhanced at lower degrees of difficulty. By comparing the time evolution of the activity of the two areas, we observed that the neural encoding of target item spatial position occurred earlier in the DLPFC than in the PMd. Increased early sodium current causes familial atrial fibrillation and dampens effect of flecainide Atrial fibrillation (AF) is the most common cardiac arrhythmia. Pathogenic variants in genes encoding ion channels are associated with familial AF. The point mutation M1875T in the SCN5A gene, which encodes the -subunit of the cardiac sodium channel Nav1.5, has been associated with increased atrial excitability and familial AF. Designing a new murine model carrying the Scn5a-M1875T mutation enabled us to study the effects of the Nav1.5 mutation in detail.

Atrial cardiomyocytes from newly generated Scn5a-M1875T+/- mice showed a selective increase in the early (peak) cardiac sodium current, measured by patch clamp. Atrial microelectrode recordings revealed larger action potential amplitudes and a faster peak upstroke velocity. Conduction was studied using optical mapping. Conduction slowing by the sodium channel blocker flecainide was less pronounced in Scn5a-M1875T+/- compared to wildtype atria. In vivo analysis using electrocardiograms and echocardiography, as well as cardiac histology, excluded overt hypertrophy or heart failure in young adult Scn5a-M1875T+/- mice.

The Scn5a-M1875T point mutation causes a cardiac sodium channel gain-of-function. Our results suggest increased atrial peak sodium current as a potential trigger for increased atrial excitability and atrial fibrillation.

Statement of SignificanceThe observed changes induced by this point mutation suggest a selective increase in peak sodium current as a cause of familial AF. Our findings provide a possible explanation for the variable effectiveness of sodium channel blockers in patients with AF. Carriers of such sodium channel gain-of-function mutations may benefit more from tailored treatments. Deletion of Rap-Phr systems in Bacillus subtilis influences in vitro biofilm formation and plant root colonization Natural isolates of the soil-dwelling bacterium Bacillus subtilis form robust biofilms under laboratory conditions and colonize plant roots. B. subtilis biofilm gene expression displays phenotypic heterogeneity that is influenced by a family of Rap-Phr regulatory systems. Most Rap-Phr systems in B. subtilis have been studied independently, in different genetic backgrounds and under distinct conditions, hampering true comparison of the Rap-Phr systems impact on bacterial differentiation. Here, we investigated each of the 12 Rap-Phr systems of B. subtilis NCIB 3610 for their effect on biofilm formation. By studying single {Delta}rap-phr mutants, we show that despite redundancy between the cell-cell communication systems, deletion of each of the 12 Rap-Phr systems influences matrix gene expression, which could possibly enable fine-tuning of the timing and level of matrix production in response to specific conditions. Furthermore, some of the {Delta}rap-phr mutants demonstrated altered biofilm formation in vitro and colonization of Arabidopsis thaliana roots, but not necessarily similarly in both processes, indicating that the pathways regulating matrix gene expression and other factors important for biofilm formation may be differently regulated under these distinct conditions.

Significance StatementNatural isolates of Bacillus subtilis form robust biofilms in vitro and on plant roots. The formation of these heterogeneous populations is regulated by diverse Rap-Phr systems. However, most Rap-Phr systems in B. subtilis have been studied independently and in different genetic backgrounds. Here, we report that all 12 Rap-Phr systems affect matrix gene expression, while some of them affect development of in vitro biofilms and plant root colonization. Our study highlights the importance of the Rap-Phr systems in environmental adaptation of B. subtilis, specifically during biofilm formation in the rhizosphere. 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. Kinetics and Mechanism of Fentanyl Dissociation from the mu-Opioid Receptor Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid over-dose is resuscitated by the use of naloxone, which competitively binds and antagonizes the {micro}-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H2976.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the N{delta}-protonated H2976.52, fentanyl can descend as much as 10 [A] below the level of the conserved D1473.32 before escaping the receptor, and has a calculated residence time{tau} of 38 s. In contrast, with the N - and doubly protonated H2976.52, the calculated{tau} are 2.6 s and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental{tau} of [~]4 min for fentanyl and the role of TM6 in mOR activation, we suggest that the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform the regulatory decision. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics.