LSD and psilocybin flatten the brain's energy landscape: insights from receptor-informed network control theory Psychedelics like lysergic acid diethylamide (LSD) and psilocybin offer a powerful window into the function of the human brain and mind, by temporarily altering subjective experience through their neurochemical effects. A recent model postulates that serotonin 2a (5-HT2a) receptor agonism allows the brain to explore its dynamic landscape more readily, as reflected by more diverse (entropic) brain activity. We postulate that this increase in entropy may arise in part from a flattening of the brains control energy landscape, which can be observed using network control theory to quantify the energy required to transition between recurrent brain states measured using functional magnetic resonance imaging (fMRI) in individuals under LSD, psilocybin, and placebo conditions. We show that LSD and psilocybin reduce the amount of control energy required for brain state transitions, and, furthermore, that, across individuals, LSDs reduction in control energy correlates with more frequent state transitions and increased entropy of brain state dynamics. Through network control analysis that incorporates the spatial distribution of 5-HT2a receptors from publicly available (non-drug) positron emission tomography (PET) maps, we demonstrate the specific role of this receptor in reducing control energy. Our findings provide evidence that 5-HT2a receptor agonist compounds allow for more facile state transitions and more temporally diverse brain activity. More broadly, by combining receptor-informed network control theory with pharmacological modulation, our work highlights the potential of this approach in studying the impacts of targeted neuropharmacological manipulation on brain activity dynamics.

Significance StatementWe present a multi-modal framework for quantifying the effects of two psychedelic drugs (LSD and psilocybin) on brain dynamics by combining functional magnetic resonance imaging (fMRI), diffusion MRI (dMRI), positron emission tomography (PET) and network control theory. Our findings provide evidence that psychedelics flatten the brains control energy landscape, allowing for more facile state transitions and more temporally diverse brain activity. We also demonstrate that the spatial distribution of serotonin 2a receptors - the main target of LSD and psilocybin - is optimized for generating these effects. This approach could be used to understand how drugs act on different receptors in the brain to influence brain activity dynamics.
Synergistic effect of mutations in the Parkinson's disease-linked endocytic proteins synaptojanin 1 (PARK20) and auxilin (PARK19) Parkinsons disease (PD) is a neurodegenerative disorder characterized by defective dopaminergic (DAergic) input to the striatum. Mutations in two genes encoding synaptically-enriched clathrin-uncoating factors, synaptojanin 1 (SJ1) and auxilin, have been implicated in PD. SJ1 knock-in (SJ1-KIRQ) mice carrying PD mutation display neurological manifestation reminiscent of PD. Here we report that auxilin knockout (Aux-KO) mice display dystrophic changes of a subset of nigrostriatal DAergic terminals similar to those of SJ1-KIRQ mice. Furthermore, Aux-KO/SJ1-KIRQ double mutant mice have shorter lifespan and more severe synaptic defects than single mutant mice. These include increase in dystrophic striatal nerve terminals positive for DAergic markers and for the PD risk protein SV2C, as well as adaptive changes in striatal interneurons. The synergistic effect of the two mutations demonstrates a special lability of DAergic neurons to defects in clathrin uncoating, with implications for PD pathogenesis in at least some forms of this condition.
Climate change and chill accumulation: implications for tree fruit production in cold-winter regions Winter chill accumulation is critical for the productivity and profitability of perennial tree fruit systems. Several studies have quantified the impacts of global warming on chill accumulation in the warmer production regions of the world, where insufficient chill events occur and their frequency is increasing. In contrast, we focus on a region with relatively cold winters-the Pacific Northwest United States (PNW)-where insufficient chill events are currently absent, and quantify the potential for introduction of these risks under climate change. Results identified spatial variation within the PNW, with chill accumulation projected to increase in some areas but decrease in others. There was also spatiotemporal variation in the driving factors of changes to chill accumulation. Even with decreases in chill accumulations, there are likely minimal issues with insufficient chill accumulation. However, delayed chill accumulation in combination with advances in the onset of heat accumulation can potentially shift the region from one where spring phenology is primarily forcing-driven to one where the dynamic interplay between chilling and forcing processes become important. These interactions might create production risks for varieties with high chill requirements, post mid-21st-century under high emissions scenarios. Future work should focus on understanding, modeling, and projecting responses across these overlapping chilling and forcing processes. Additionally, given significant spatial differences across a relatively small geographic range, it is also critical to understand and model these dynamics at a local landscape resolution for regions such as the PNW. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Variation in female songbird state determines signal strength needed to evoke copulation It is the female response to male signals that determines courtship success. In most songbirds, females control reproduction via the copulation solicitation display (CSD), an innate, stereotyped posture produced in direct response to male displays. Because CSD can be elicited in the absence of males by the presentation of recorded song, CSD production enables investigations into the effects of underlying signal features and behavioral state on female mating preferences. Using computer vision to quantify CSD trajectory in female brown-headed cowbirds (Molothrus ater), we show that both song quality and a females internal state predict CSD production, as well as the onset latency and duration of the display. We also show that CSD can be produced in a graded fashion based on both signal strength and internal state. These results emphasize the importance of underlying receiver state in determining behavioral responses and suggest that female responsiveness acts in conjunction with male signal strength to determine the efficacy of male courtship. 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. Bone marrow haematopoietic stem cells influence liver homeostatic networks and cancer development In fatty liver disease, systemic homeostasis is perturbed. While pre-clinical models are used to understand its pathogenesis, translating this knowledge to patients is difficult. However, by focusing on the most preserved homeostasis systems between species and models, novel disease dimensions can be unearthed. We interrogated core liver gene co-expression networks in a mouse model of liver cancer following dietary challenge. Immunometabolic modules showed temporal changes under the influence of diet duration and aging. The behaviour of immune network in tumours mirrored their counterparts in non-tumour liver. A high immune response network was associated with a lower tumour burden in mice and humans. This module in mice was enriched for genes related to haematopoietic cell differentiation. Consistently, the bone marrow haematopoietic stem and progenitor cells response was reflective of the liver immune response. Linking haematopoiesis to hepatic homeostasis uncovers a hitherto unexplored dimension of tissue crosstalk that can inform pathogenesis. Evaluation of methods for estimating coalescence times using ancestral recombination graphs The ancestral recombination graph (ARG) is a structure that describes the joint genealogies of sampled DNA sequences along the genome. Recent computational methods have made impressive progress towards scalably estimating whole-genome genealogies. In addition to inferring the ARG, some of these methods can also provide ARGs sampled from a defined posterior distribution. Obtaining good samples of ARGs is crucial for quantifying statistical uncertainty and for estimating population genetic parameters such as effective population size, mutation rate, and allele age. Here, we use standard neutral coalescent simulations to benchmark the estimates of pairwise coalescence times from three popular ARG inference programs: ARGweaver, Relate, and tsinfer+tsdate. We compare 1) the true coalescence times to the inferred times at each locus; 2) the distribution of coalescence times across all loci to the expected exponential distribution; 3) whether the sampled coalescence times have the properties expected of a valid posterior distribution. We find that inferred coalescence times at each locus are most accurate in ARGweaver, and often more accurate in Relate than in tsinfer+tsdate. However, all three methods tend to overestimate small coalescence times and underestimate large ones. Lastly, the posterior distribution of ARGweaver is closer to the expected posterior distribution than Relates, but this higher accuracy comes at a substantial trade-off in scalability. The best choice of method will depend on the number and length of input sequences and on the goal of downstream analyses, and we provide guidelines for the best practices. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. The isoforms of pyruvate kinase act as nutrient sensors for the β-cell KATP channel Insulin release from pancreatic {beta}-cells is triggered by the metabolism-dependent closure of ATP-sensitive K+ (KATP) channels. Here, we used isoform-specific deletion of pyruvate kinase (PK) to identify the unique roles of PKm1 and PKm2 in the {beta}-cell nutrient response. In the absence of PK activity, cytosolic ATP/ADP was insufficient to close KATP channels. Both PK isoforms were found on the plasma membrane, where they endow KATP channels with nutrient sensitivity and specificity. Although PKm1 accounts for >90% of cellular activity, PKm2 controls glucosedependent KATP closure. Amino acids, which generate mitochondrial phosphoenolpyruvate (PEP) without allosterically activating PKm2, depend on PKm1 to raise ATP/ADP and close KATP channels. Shifting the {beta}-cell from PKm2 to PKm1 correspondingly increased amino acid sensitivity and boosted secretory output. In addition to these membrane depolarizing "on-switches," {beta}-cell deletion of PCK2 exposed a mitochondrial PEP-dependent "off-switch" that facilitates Ca2+ extrusion. Restricting PCK2 prolonged Ca2+ influx and increased insulin secretion. In summary, it is not the KATP channel alone, but the PK-KATP microcompartment that confers nutrient sensitivity to the {beta}-cell plasma membrane. The differential response of the PK isoforms to the glycolytic and mitochondrial sources of PEP influences {beta}-cell nutrient preference, and controls the oscillatory cycle regulating insulin secretion. Opposing mechanisms of competition determine species invasions and functional diversity Understanding how species phenotypic differences affect competition is key to explaining community assembly and predicting biodiversity responses. Many studies overlook the variable effects that species trait differences can have on the direction of competitive exclusion, which reverses depending on the specific mechanism at play. We performed a comprehensive trait-based study of an ant invasion integrating morphological, dietary, physiological and behavioral analyses. We found that trait differences between invasive and resident species were not only associated with niche differences which promoted the coexistence of dissimilar species, but also competitive ability differences which acted in the opposite fashion. Furthermore, competition along separate trait axes led to complex and contrasting patterns in the invaded assemblages, where species were at once similar (clustered) in some traits but also dissimilar (overdispersed) in others. Our results reveal that different aspects of phenotype may distinctly modulate the effect of competition in structuring ecological communities and functional diversity. The extracellular loops of OmpA control the slow rate of in vitro folding Outer membrane proteins are all beta barrels and these barrels have a variety of well-documented loop conformations. Here we test the effect of three different loop types on outer membrane protein A (OmpA) folding. We designed twelve 5-residue loops and experimentally tested the effect of replacing the long loops of outer membrane protein OmpA with the designed loops. Our studies succeeded in creating the smallest known outer membrane barrel. We find that significant changes in OmpA loops do not have a strong overall effect on OmpA folding. However, when decomposing folding into a fast rate and a slow rate we find that changes in loops strongly affect the slow rate of OmpA folding. Extracellular loop types with higher levels of hydrogen bonds had more instances of increasing the slow folding rate and extracellular loop types with low levels of hydrogen bonds had more instances of decreasing the slow folding rate. Having the slow rate affected by loop composition is consistent with the slow rate being associated with the insertion step of outer membrane protein folding. Feed efficiency of lactating Holstein cows was not as repeatable across diets as within diet over subsequent lactation stages BackgroundImproving feed efficiency has become a common target for dairy farmers to meet the requirement of producing more milk with fewer resources. To improve feed efficiency, a prerequisite is to ensure that the cows identified as most or least efficient will remain as such, independently of diet composition. Therefore, the current research analysed the ability of lactating dairy cows to maintain their feed efficiency while changing the energy density of the diet by changing its concentration in starch and fibre. A total of 60 lactating Holstein cows, including 33 primiparous cows, were first fed a high starch diet-low fibre (diet S+F-), then switched over to a low starch diet-high fibre (diet S-F+). To know if diet affect feed efficiency, we compared the ability of feed efficiency to be maintained within a diet over subsequent lactation stages, known as repeatability of feed efficiency, with its ability to be maintained across diets, known as reproducibility of feed efficiency. To do so we used two indicators: the estimation of the error of repeatability/reproducibility, which is commonly used in metrology, and the coefficient of correlation of concordance (CCC), which is used in biology. The effect of diet change could also lead to a change in cows sorting behaviour which could potentially affect feed efficiency if for example the most efficient cows select more concentrate than the least efficient. We therefore analysed the relationship between the differences in individual feed refusals composition and the differences in feed efficiency. To do so, the composition of each feed refusal was described with its near infra-red (NIR) spectroscopy and was performed on each individual feed ingredient, diet and refusals and used as composition variable. The variability of the NIR spectra of the refusals was described with its principal components thanks to a principal component analysis (PCA). The Pearson correlation was estimated to check the relationship between feed efficiency and refusals composition, i.e. sorting behaviour.

ResultsThe error of reproducibility of feed efficiency across diets was 2.95 MJ/d. This error was significantly larger than the errors of repeatability estimated within diet, which were 2.01 MJ/d within diet S-F+ and 2.40 MJ/d within diet S+F-. The CCC was 0.64 between feed efficiency estimated within diet S+F- and feed efficiency estimated within diet S-F+. This CCC was smaller than the one observed for feed efficiency estimated within diet between two subsequent lactation stages (CCC = 0.72 within diet S+F- and 0.85 within diet S-F+). Feed efficiency was poorly correlated to the first two principal components, which explained 90% of the total variability of the NIR spectra of the individual refusals. This suggests that feed sorting behaviour did not explain differences in feed efficiency.

ConclusionsFeed efficiency was significantly less reproducible across diets than repeatable within the same diet over subsequent lactation stages, but cows ranking for feed efficiency was not significantly affected by diet change. This loss in repeatability across diets could be due to a more pronounced feed sorting subsequent to the change in diet composition. However, the differences in sorting behaviour between cows were not associated to feed efficiency differences in this trial neither with the S+F- diet nor with the S-F+ diet. Those results have to be confirmed on diets having different forage to concentrate ratios to ensure that the least and most efficient cows will not change. Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific 4s subunit replaces the 4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82-92%) during differentiation, presumably due to the shift from 4 to 4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S, and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between 4 and 4s. We propose that these differences induced by 4s incorporation result in significant changes in the way the s20S interacts with its partners, and dictate its role in germ cell differentiation. A cell-based screen identifies HDAC inhibitors as activators of RIG-I signaling Enhancing the immune microenvironment in cancer by targeting the nucleic acid sensors is becoming a potent therapeutic strategy. Among the nucleic acid sensors, activation of the RNA sensor Retinoic Acid-inducible Gene (RIG-I) using small hairpin RNAs has been shown to elicit powerful innate and adaptive immune responses. Given the challenges inherent in pharmacokinetics and delivery of RNA based agonists, we set out to discover small molecule agonists of RIG-I using a cell-based assay. To this end, we established and validated a robust high throughput screening assay based on a commercially available HEK293 reporter cell line with a luciferase reporter downstream of tandem interferon stimulated gene 54 (ISG54) promoter elements. We first confirmed that the luminescence in this cell line is dependent on RIG-I and the interferon receptor using a hairpin RNA RIG-I agonist. We established a 96-well and a 384-well format HTS based on this cell line and performed a proof-of-concept screen using an FDA approved drug library of 1200 compounds. Surprisingly, we found two HDAC inhibitors Entinostat, Mocetinostat and the PLK1 inhibitor Volasertib significantly enhanced ISG-luciferase activity. This luminescence was substantially diminished in the null reporter cell line indicating the increase in signaling was dependent on RIG-I expression. Combination treatment of tumor cell lines with Entinostat increased RIG-I induced cell death in a mammary carcinoma cell line that is resistant to either Entinostat or RIG-I agonist alone. Taken together, our data indicates an unexpected role for HDAC1,-3 inhibitors in enhancing RIG-I signaling and highlight potential opportunities for therapeutic combinations. Strength of interactions in the Notch gene regulatory network determines patterning and fate in the notochord Development of multicellular organisms requires the generation of gene expression patterns that determines cell fate and organ shape. Groups of genetic interactions known as Gene Regulatory Networks (GRNs) play a key role in the generation of such patterns. However, how the topology and parameters of GRNs determine patterning in vivo remains unclear due to the complexity of most experimental systems. To address this, we use the zebrafish notochord, an organ where coin-shaped precursor cells are initially arranged in a simple unidimensional geometry. These cells then differentiate into vacuolated and sheath cells. Using newly developed transgenic tools together with in vivo imaging, we identify jag1a and her6/her9 as the main components of a Notch GRN that generates a lateral inhibition pattern and determines cell fate. Making use of this experimental system and mathematical modeling we show that lateral inhibition patterning requires that ligand-receptor interactions are stronger within the same cell than in neighboring cells. Altogether, we establish the zebrafish notochord as an experimental system to study pattern generation, and identify and characterize how the properties of GRNs determine self-organization of gene patterning and cell fate. Small extracellular vesicles from young mice prevent frailty, improve healthspan and decrease epigenetic age in old mice. Aging is associated with an increased risk of frailty, disability, comorbidities, institutionalization, falls, fractures, hospitalization, and mortality. Searching for strategies to delay the degenerative changes associated with aging and frailty is interesting. We treated old animals intravenously with small extracellular vesicles (sEVs) derived from adipose mesenchymal stem cells (ADSCs) of young animals, and we found an improvement of several functional parameters usually altered with aging, such as motor coordination, grip strength, fatigue resistance, fur regeneration, and renal function. Frailty index analysis showed that 40% of old control mice were frail, whereas none of the old ADSCs-sEVs treated mice were. Molecular and structural benefits in muscle and kidney accompanied this functional improvement. ADSCs-sEVs induced pro-regenerative effects and a decrease in oxidative stress, inflammation, and senescence markers. Moreover, predicted epigenetic age was lower in tissues of old mice treated with ADSCs-sEVs and their metabolome changed to a youth-like pattern. Finally, we gained some insight into the miRNAs contained in sEVs that might be, at least in part, responsible for the effects observed. We propose that young sEVs treatment can be beneficial against frailty and therefore can promote healthy aging. Mesoscale organization in the cell envelope of Deinococcus radiodurans S-layers are highly ordered coats of proteins localized on the cell surface of many bacterial species. In these structures, one or more proteins form elementary units that self-assemble into a crystalline monolayer tiling the entire cell surface. Here, the cell envelope of the radiation-resistant bacterium Deinococcus radiodurans was studied by high-resolution cryo-electron microscopy finding the crystalline regularity of the S-layer extended into the layers below. The cell envelope appears to be highly packed and resulting from a three-dimensional crystalline distribution of protein complexes organized in close continuity but allowing different degrees of voidness in the entire thickness. These insights grade S-layers to mesoscale hubs behaving as structural and functional architraves essential for the entire cell body. Recurring adaptive introgression of a supergene variant that determines social organization Introgression has been proposed as an essential source of adaptive genetic variation. However, a key barrier to adaptive introgression is that recombination can break down combinations of alleles that underpin many traits. This barrier might be overcome in supergene regions, where suppressed recombination leads to joint inheritance across many loci. Here, we study the evolution of a large supergene region that determines a major social and ecological trait in Solenopsis fire ants: whether colonies have one queen or multiple queens. Using coalescent-based phylogenies built from the genomes of 365 haploid fire ant males, we show that the supergene variant responsible for multiple-queen colonies evolved in one species and repeatedly spread to other species through introgressive hybridization. This finding highlights how supergene architecture can enable a complex adaptive phenotype to recurrently permeate species boundaries. Deficiency of gluconeogenic enzyme PCK1 promotes murine NASH progression and fibrosis by activation PI3K/AKT/PDGF axis Background and AimsNonalcoholic steatohepatitis (NASH) is a chronic liver disease characterized by hepatic lipid accumulation, inflammation, and progressive fibrosis. However, the pathomechanisms underlying NASH are incompletely explored. Phosphoenolpyruvate carboxykinase 1 (PCK1) catalyzes the first rate-limiting step of gluconeogenesis. This study was designed to determine the role of PCK1 in regulating NASH progression.

MethodsLiver metabolism, hepatic steatosis, and fibrosis were evaluated at 24 weeks in liver-specific Pck1-knockout mice fed with NASH diet or chow diet. Gain- or loss-of-function approaches were used to explore the underlying mechanism in vitro. AKT and RhoA inhibitors were evaluated for NASH treatment in vivo.

ResultsHepatic PCK1 was downregulated in patients with NASH and mouse models of NASH. Mice with liver Pck1 deficiency displayed hepatic lipid disorder and liver injury fed with normal diet, while showed aggravated fibrosis and inflammation when fed NASH diet. Mechanistically, PCK1 deficiency upregulated genes involved in fatty acid transport and lipid droplet formation. Moreover, metabolomics analysis showed the accumulation of glycerol 3-phosphate, the substrate of triglyceride synthesis. Furthermore, the loss of hepatic PCK1 could activate the RhoA/PI3K/AKT pathway, which leads to increased secretion of PDGF-AA and promotes the activation of hepatic stellate cells. Accordingly, treatment with RhoA and AKT inhibitors alleviated NASH progression in the presence of Pck1 deletion in vivo.

ConclusionsPCK1 deficiency plays a key role in the development of hepatic steatosis and fibrosis by facilitating the RhoA/PI3K/AKT/PDGF-AA axis. These findings provide a novel insight into therapeutic approaches for the treatment of NASH. Tackling polypharmacology of kinase inhibitors and selecting chemical kinase probes by transcription factor activity profiling. Protein kinase inhibitors (PKI) evolved as promising drugs for multiple diseases. However, PKIs promiscuity causes polypharmacological effects challenging therapeutic development. We show that the polypharmacology of PKIs can be inferred from their impact on transcription factors (TF) lying at the apexes of signaling pathways. Using a high-content reporter assay, we evaluated TF activity profiles (TFAP) in cells treated with Akt, CDK, Aurora, Raf, MEK, and ERK inhibitors. For each kinase family, we found kinase-specific consensus TFAP signatures signifying the on-target PKI activity. Remarkably, proximal kinases had the highest similarity of their PKIs consensus signatures. Furthermore, we show that individual PKIs exhibited the consensus, "on-target signatures" within certain "specificity windows" and distinct "off-target signatures" at other concentrations. We also show that the TFAP signatures permit pinpointing the off-target PKI activity at irrelevant kinases and non-kinase effectors. Therefore, the TFAP approach provides clear-cut quantitative metrics for assessing PKI activity and concentration ranges where the on-target activity dominates cell response. By illuminating PKI biology invisible to target-based techniques, this effect-based approach streamlines the selection of kinase chemical probes and PKI drug prioritization. Native American Ancestry and Pigmentation Allele Contributions to Skin Color in a Caribbean Population Interest in the genetic basis of variation in skin pigmentation in Native American populations led us to seek indigenous populations of the Western Hemisphere with African and minimal European admixture to study the effect of Native American ancestry on skin color. Admixture analysis from DNA collected from 458 individuals in the Kalinago territory of the Commonwealth of Dominica showed shared ancestry with East Asians at K=3 and 55% Native American, 32% African, and 11% European ancestry at K=6, the highest Native American ancestry of Caribbean populations. Skin pigmentation was 20 to 80 melanin units, averaging 46. Three albino individuals were homozygous for multi-nucleotide polymorphism OCA2NW273KV of African origin, whose population allele frequency was 0.03 and single allele effect size was -8 melanin units. Hypopigmenting allele frequencies for SLC24A5A111T and SLC45A2L374F were 0.14 and 0.05, whose single allele effect sizes were -6 and -3, respectively. Skin color plots of individuals lacking known hypopigmenting alleles suggests that Native American Ancestry reduced pigmentation by more than 20 melanin units (low and high estimates 21.8 and 28.5). Shared ancestry with East Asians at K=3 suggests potential sharing of one or more pigmentation alleles. Re-convolving the compositional landscape of primary and recurrent glioblastoma reveals prognostic andtargetable tissue states Glioblastoma is an aggressive diffusely infiltrating neoplasm that spreads beyond surgical resection margins, where it intermingles with non-neoplastic brain cells. This complex microenvironment harbouring infiltrating glioma and non-neoplastic brain cells is the origin of tumor recurrence. Thus, understanding the cellular and molecular features of the glioma microenvironment is therapeutically and prognostically important. Here, we used single-nucleus RNA sequencing (snRNAseq) of primary and recurrent glioma to define three compositional  tissue-states rooted in cohabitation of cell-types and transcriptional states. These comprise a state featuring A) abundant normal brain microenvironment cells, B) reactive/inflamed brain microenvironment, and C) cellular/proliferative tumor. All these states exhibited variable degrees of infiltration by glioma cells. Spatial transcriptomics confirmed that the cell-types and transcriptomics states which compositionally cohabitate tissue states also colocalize in space. Tissue states are clinically significant because they correlate with radiographic, histopathologic, and prognostic features. We found that enrichment of tissue state B signature correlated with shorter survival. Importantly, we found that our compositionally defined tissue states are enriched in distinct metabolic pathways. One such pathway is fatty acid biosynthesis, which was enriched in tissue state B - a state enriched in recurrent glioblastoma and composed of astrocyte-like/mesenchymal glioma cells, reactive astrocytes resembling those seen in neurodegeneration, and monocyte-like myeloid cells. We showed that treating acute slices of GBM with a fatty acid synthesis inhibitor is sufficient to deplete the transcriptional signature of tissue state B. Our findings define a novel compositional approach to the glioma-infiltrated tissue which allows us to discover prognostic and targetable features, paving the way to new mechanistic and therapeutic discoveries. A protein microarray-based in vitro transglutaminase assay platform for epitope mapping and peptide discovery: implication for transglutaminase-mediated immunodominant determination Transglutaminases (TGs) are a family of crosslinking enzymes catalyzing the formation of intra- and inter-molecular glutamine-lysine isopeptide bonds in a calcium dependent manner. Activation of transglutaminases is pathogenically associated with severe human diseases including neurodegenerations, cardiovascular diseases, and autoimmune diseases. Although continuous efforts determining the enzymes substrate preference have been witnessed, a high-throughput assay platform with the "omic" efficiency is still missing for the global identification of substrate-specific TG modification sites. In this study we report a protein microarray-based in vitro TG assay platform for rapid and large-scale (up to 30000 reactions per chip) determination of the glutamine (Q)-bearing TG modification motifs. With this platform we identified the Q16 in superoxide dismutase 1 and Q109 in alpha-synuclein as the primary modification sites for tissue transglutaminase (TG2), the most ubiquitous member of the enzyme family. Of particular interest, we found a close match between the modification motifs and published vaccine epitope sequences in alpha-synuclein, implying an essential and intrinsic role transglutaminase might play in the determination of immunodominant epitopes. Our data collectively suggest the glutamine and its follow-up five residues on the C terminal of a protein compose a minimal determinant motif for TG2 modification and the TG2 modification motifs determined by our platform could finally correspond to the substrates epitope sequences in antigen processing. To screen for site-specific interfering peptides and assist gene editing-based protein engineering, we optimized through onchip amino-acid scanning the TG2 modification motif QQIV in the extracellular matrix protein fibronectin and obtained 14 variants with significantly higher TG2 reactivity that might serve as the competitive inhibitor peptides and 1 with lower reactivity. We further confirmed this optimization approach with 12-mer peptides, the longest ones that could be synthesized on the chip. Taken together, our synthetic transglutaminase assay platform might be able to deliver a precise epitope blueprint for immunotherapeutic targeting and provide proof-of-concept and directional studies for TG-based peptide discovery and protein design. LinearTurboFold: Linear-Time Global Prediction of Conserved Structures for RNA Homologs with Applications to SARS-CoV-2 The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in SARS-CoV-2 genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length, and are thus infeasible for coronaviruses, which possess the longest genomes ([~]30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurbo-Folds purely in silico prediction not only is close to experimentally-guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5 and 3 UTRs ([~]29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies novel conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, siRNAs, CRISPR-Cas13 guide RNAs and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies, and will be a useful tool in fighting the current and future pandemics.

Significance StatementConserved RNA structures are critical for designing diagnostic and therapeutic tools for many diseases including COVID-19. However, existing algorithms are much too slow to model the global structures of full-length RNA viral genomes. We present LinearTurboFold, a linear-time algorithm that is orders of magnitude faster, making it the first method to simultaneously fold and align whole genomes of SARS-CoV-2 variants, the longest known RNA virus ([~]30 kilobases). Our work enables unprecedented global structural analysis and captures long-range interactions that are out of reach for existing algorithms but crucial for RNA functions. LinearTurboFold is a general technique for full-length genome studies and can help fight the current and future pandemics.