Agonist efficiency from concentration-response curves: structural implications and applications Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC50) that indicates potency, a high-concentration asymptote that indicates efficacy and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency ({eta}), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that {eta} can be calculated from EC50 and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of {eta} values is bimodal with population means at 51% (including ACh, nornicotine and DMPP) and 40% (including epibatidine, varenicline and cytisine). The value of {eta} is related inversely to the size of the agonists head-group, with high-versus low-efficiency ligands having an average volume of 70 [A]3 versus 102 [A]3. Most binding site mutations have only a small effect on ACh efficiency except for Y190A (35%), W149A (60%) and those at G153 (42%). If {eta} is known, the midpoint and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC50 of a CRC that has been normalized to 1. Given {eta}, the level of constitutive activity can be estimated from a single CRC.

Statement of significanceReceptors are molecular machines that convert chemical energy from agonist binding into mechanical energy of a global conformational change that generates a cell response. Agonists are distinguished by their potency (proportional to affinity) and efficacy, but also by the efficiency at which their binding energy is applied to receptor activation. Here we show that agonist efficiency can be estimated from a single concentrationresponse curve, and estimate efficiencies of 20 nicotinic receptor agonists. These have a bimodal distribution with larger agonists belonging to the lower-efficiency population. We further show that mutations of some binding site residues alter efficiency, and that knowledge of agonist efficiency simplifies and extends dose-response curve analysis.
ATP binding facilitates target search of SWR1 chromatin remodeler by promoting one-dimensional diffusion on DNA One-dimensional (1D) target search is a well characterized phenomenon for many DNA binding proteins but is poorly understood for chromatin remodelers. Herein, we characterize the 1D scanning properties of SWR1, a yeast chromatin remodeler that performs histone exchange on +1 nucleosomes adjacent to a nucleosome depleted region (NDR) at promoters. We demonstrate that SWR1 has a kinetic binding preference for DNA of NDR length as opposed to gene-body linker length DNA. Using single and dual color single particle tracking on DNA stretched with optical tweezers, we directly observe SWR1 diffusion on DNA. We found that various factors impact SWR1 scanning, including ATP which promotes diffusion through nucleotide binding rather than ATP hydrolysis. A DNA binding subunit, Swc2, plays an important role in the overall diffusive behavior of the complex, as the subunit in isolation retains similar, although faster, scanning properties as the whole remodeler. ATP-bound SWR1 slides until it encounters a protein roadblock, of which we tested dCas9 and nucleosomes. The median diffusion coefficient, 0.024 m2/sec, in the regime of helical sliding, would mediate rapid encounter of NDR-flanking nucleosomes at length scales found in cells.
The limits of evolutionary convergence in sympatry: reproductive interference and historical constraints leading to local diversity in warning traits Mutualistic interactions between defended species represent a striking case of evolutionary convergence in sympatry, driven by the increased protection against predators brought by mimicry in warning traits. However, such convergence is often limited: sympatric defended species frequently display different or imperfectly similar warning traits. The phylogenetic distance between sympatric species may indeed prevent evolution towards the exact same signal. Moreover, warning traits are also involved in mate recognition, so that trait convergence might result in heterospecific courtship and mating. Here, we develop a mathematical model to investigate the strength and direction of evolution of warning trait in defended species with different ancestral traits. Specifically, we determine the effect of phenotypic distances between ancestral trait states of sympatric defended species and of costs of heterospecific sexual interactions on imperfect mimicry and trait divergence. Our analytical results confirm that reproductive interference and historical constraints limit the convergence of warning traits, leading to either complete divergence, or imperfect mimicry. Our model reveals that imperfect mimicry evolves only when ancestral trait values differ between species, due to historical constraints and highlights the importance of female and predator discriminations in the evolution of such imperfect mimicry. Our study thus provides new predictions on how reproductive interference interacts with historical constraints and may promote the emergence of novel warning traits, enhancing mimetic diversity. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Greater spear nosed bats in Panama do not use social proximity to improve foraging efficiency Group foraging animals may change their spatial relationships in response to the resource landscape. The strength of social network ties across foraging groups could then reflect the energy balance of individual group members through their effect on foraging efficiency. Greater spear-nosed bats (Phyllostomus hastatus) live in social groups where unrelated females develop group-specific social calls, defend each others offspring, and are hypothesized to coordinate social foraging to feed on balsa flowers. To test if social groups forage together, and if this increases foraging efficiency, we tracked 18 individuals from three social groups of P. hastatus on Isla Colon, Panama using high frequency GPS for 1 - 3 days each. We found that individuals commuted alone but then joined group members in small foraging patches with high densities of flowering balsas, 25 km from their roost on the mainland. Proximity-based social network measures did not predict foraging efficiency or energy expenditure, but individuals with higher social associations spent more time outside of the cave. Individuals did not coordinate movement, but close proximity among individuals increased resting bout duration and synchrony, particularly when in the same patch. These results show that stronger social proximity associations were more closely tied to resting behavior, and indicate that factors other than increased feeding efficiency may structure network relationships of group members while foraging.

Lay SummarySocial foraging can increase foraging efficiency on hard-to-find food sources. Miniaturized high-frequency GPS allowed us to test this on groups of greater spear-nosed bats. We estimated food intake rates, energy expenditure, and the strength of bat social relationships. In contrast to our expectations, we did not find a relationship between bat social relationships and energy expenditure or foraging efficiency. However, bats took longer resting breaks when group members were close. Why these bats commute 25 km individually across the ocean to then forage near each other thus requires further investigation. A partially nested cortical hierarchy of neural states underlies event segmentation in the human brain A fundamental aspect of human experience is that it is segmented into discrete events. This may be underpinned by transitions between distinct neural states. Using an innovative data-driven state segmentation method, we investigate how neural states are organized across the cortical hierarchy and where in the cortex neural state boundaries and perceived event boundaries overlap. Our results show that neural state boundaries are organized in a temporal cortical hierarchy, with short states in primary sensory regions, and long states in lateral and medial prefrontal cortex. State boundaries are shared within and between groups of brain regions that resemble well-known functional networks. Perceived event boundaries overlap with neural state boundaries across large parts of the cortical hierarchy, particularly when those state boundaries demarcate a strong transition or are shared between brain regions. Taken together, these findings suggest that a partially nested cortical hierarchy of neural states forms the basis of event segmentation. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. Limited selection towards lighter seeds in Danish grasslands over an eight-year period Time-series analysis of cover data from 236 Danish grassland sites demonstrated community selection on seed mass. Across four grassland habitat types and during an eight-year period there was significant community selection towards lighter seeds. In the analysis, the continuous seed size variable was used to group plant species into functional types. This method was chosen in order to account for the sampling process of the abundance data, and is in contrast to most other analyses of trait selection, where the community weighted mean of the traits is used as the dependent variable. Generally, little is known on the underlying causes of possible plant trait community selection in semi-natural grasslands that are not subjected to successional processes. Localization of signaling receptors maximizes cellular informationacquisition in spatially-structured natural environments Cells in natural environments like tissue or soil sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization maximizes information acquisition in simulated natural contexts, including tissue and soil. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold. Broadly, our work shows how cells can maximize the fidelity of information transfer by adapting the spatial organization of signaling molecules to the spatial structure of the environment. 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. Weighted transitivity scores that account for triadic edge similarity in undirected graphs The graph transitivity measures the probability that adjacent vertices in a network are interconnected, thus revealing the existence of tightly connected neighborhoods playing a role in information and pathogen circulation. When the connections vary in strength, focusing on whether connections exist or not can be reductive. I score the weighted transitivity according to the similarity between the weights of the three possible links in each triad. I illustrate the biological relevance of that information with two reanalyses of animal contact networks. In the rhesus macaque Macaca mulatta, a species in which kin relationships strongly predict social relationships, the new metrics revealed striking similarities in the configuration of grooming networks in captive and free-ranging groups, but only as long as the matrilines were preserved. In the barnacle goose Branta leucopsis, in an experiment designed to test the long-term effect of the goslings social environment, the new metrics uncovered an excess of weak triplets closed by strong links in males compared to females, and consistent with the triadic process underlying goose dominance relationships. A network approach to genetic circuit designs As genetic circuits become more sophisticated, the size and complexity of data about their designs increases. This data captured goes beyond monolithic genetic sequences and towards circuit modularity and functional details, which are beneficial for analyzing circuit performance and establishing design automation techniques. However, the accessibility, visualisation and usability of design data (and metadata) have received relatively little attention to date. Here, we present a method to turn circuit designs into networks and showcase its potential to enhance the utility of design data. Since networks are dynamic structures, initial graphs can be interactively shaped into sub-networks of relevant information based on requirements such as abstraction, hierarchy and protein interactions. Additionally, several visual changes can be applied, such as colouring or clustering nodes based on types (e.g., genes or promoters), resulting in easier comprehension from a user perspective. This approach allows circuit designs to be coupled to other networks, such as metabolic pathways or implementation protocols captured in graph-like formats. Therefore, we advocate using networks to structure, access and improve synthetic biology information. 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. A chemical screen based on an interruption of zebrafish gastrulation identifies the HTR2C inhibitor Pizotifen as a suppressor of EMT-mediated metastasis Metastasis is responsible for approximately 90% of cancer-associated mortality but few models exist that allow for rapid and effective screening of anti-metastasis drugs. Current mouse models of metastasis are too expensive and time consuming to use for rapid and high-throughput screening. Therefore, we created a unique screening concept utilizing conserved mechanisms between zebrafish gastrulation and cancer metastasis for identification of potential anti-metastatic drugs. We hypothesized that small chemicals that interrupt zebrafish gastrulation might also suppress metastatic progression of cancer cells and developed a phenotype-based chemical screen to test the hypothesis. The screen used epiboly, the first morphogenetic movement in gastrulation, as a marker and enabled 100 chemicals to be tested in five hours. The screen tested 1280 FDA-approved drugs and identified Pizotifen, an antagonist for serotonin receptor 2C (HTR2C) as an epiboly-interrupting drug. Pharmacologic and genetic inhibition of HTR2C suppressed metastatic progression in a mouse model. Blocking HTR2C with Pizotifen restored epithelial properties to metastatic cells through inhibition of Wnt-signaling. In contrast, HTR2C induced epithelial to mesenchymal transition (EMT) through activation of Wnt-signaling and promoted metastatic dissemination of human cancer cells in a zebrafish xenotransplantation model. Taken together, our concept offers a novel platform for discovery of anti-metastasis drugs. Experimental evolution reveals post-transcriptional regulation as a novel driver of Leishmania fitness gain The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network driving parasite fitness, with genome instability causing highly reproducible, gene dosage-dependent changes in protein abundance linked to post-transcriptional regulation. These in turn were associated with a gene dosage-independent reduction in abundance of flagellar transcripts and a coordinated increase in abundance of coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of fitness-adapted ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings. Deep Heterogeneous Dilation of LSTM for Transient-phase Gesture Prediction through High-density Electromyography: Towards Application in Neurorobotics Deep networks have been recently proposed to estimate motor intention using conventional bipolar surface electromyography (sEMG) signals for myoelectric control of neurorobots. In this regard, Deepnets are generally challenged by long training times (affecting practicality and calibration), complex model architectures (affecting the predictability of the outcomes), and a large number of trainable parameters (increasing the need for big data). Capitalizing on our recent work on homogeneous temporal dilation in a Recurrent Neural Network (RNN) model, this paper proposes, for the first time, heterogeneous temporal dilation in an LSTM model and applies that to high-density surface electromyography (HD-sEMG), allowing for the decoding of dynamic temporal dependencies with tunable temporal foci. In this paper, a 128-channel HD-sEMG signal space is considered due to the potential for enhancing the spatiotemporal resolution of human-robot interfaces. Accordingly, this paper addresses a challenging motor intention decoding problem of neurorobots, namely, transient intention identification. Our approach uses only the dynamic and transient phase of gesture movements when the signals are not stabilized or plateaued, which can significantly enhance the temporal resolution of human-robot interfaces. This would eventually enhance seamless real-time implementations. Additionally, this paper introduces the concept of "dilation foci" to modulate the modeling of temporal variation in transient phases. In this work a high number (e.g., 65) of gestures is included, which adds to the complexity and significance of the understudied problem. Our results show state-of-the-art performance for gesture prediction in terms of accuracy, training time, and model convergence. Polycomb-mediated repression of paternal chromosomes maintains haploid dosage in diploid embryos of Marchantia Complex mechanisms regulate gene dosage throughout eukaryotic life cycles. Mechanisms controlling gene dosage have been extensively studied in animals, however it is unknown how generalizable these mechanisms are to diverse eukaryotes. Here, we use the haploid plant Marchantia polymorpha to assess gene dosage control in its short-lived diploid embryo. We show that throughout embryogenesis, paternal chromosomes are repressed resulting in functional haploidy. The paternal genome is targeted for genomic imprinting by the Polycomb mark H3K27me3 starting at fertilization, rendering the maternal genome in control of embryogenesis. Maintaining haploid gene dosage by this new form of imprinting is essential for embryonic development. Our findings illustrate how haploid-dominant species can regulate gene dosage through paternal chromosome inactivation and initiates the exploration of the link between life cycle history and gene dosage in a broader range of organisms. Common coupling map advances GPCR-G protein selectivity Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein  couplome. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on GPCR-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling towards the exploitation of G protein signaling bias in design of safer drugs. Human 14-3-3 proteins site-selectively bind the mutational hotspot region of SARS-CoV-2 nucleoprotein modulating its phosphoregulation Phosphorylated within its Ser/Arg-rich region, the SARS-CoV-2 nucleoprotein (N) recruits the phosphopeptide-binding human 14-3-3 proteins that play a well-recognized role in replication of many viruses. Here we use genetic code expansion to demonstrate that phosphorylation of SARS-CoV-2 N at either of two pseudo-repeats centered at Ser197 and Thr205 is sufficient for 14-3-3 binding. According to fluorescence anisotropy, the pT205-motif, present in SARS-CoV-2 but not in SARS-CoV, is preferred over the pS197-motif by all seven human 14-3-3 isoforms, which display unforeseen pT205/pS197 binding selectivity hierarchy. Crystal structures demonstrate that pS197 and pT205 are mutually exclusive 14-3-3-binding sites, whereas SAXS and biochemical data indicate 14-3-3 binding occludes the Ser/Arg-rich region, inhibiting its dephosphorylation. This Ser/Arg-rich region of N is highly prone to mutations, as exemplified by the Omicron and Delta variants, with our data suggesting how the strength of its 14-3-3 binding can be linked with the replicative fitness of the virus. Multiple origins of nematode-Wolbachia symbiosis in supergroup F and convergent loss of bacterioferritin in filarial Wolbachia The intracellular endosymbiotic proteobacteria Wolbachia have evolved across the phyla nematoda and arthropoda. In Wolbachia phylogeny, supergroup F is the only clade with members from both arthropod and filarial nematode hosts and therefore can provide unique insights into their evolution and biology. In this study, 4 new supergroup F Wolbachia genomes have been assembled using a metagenomic assembly and binning approach, wMoz and wMpe from the human filarial parasites Mansonella ozzardi and Mansonella perstans, and wOcae and wMoviF from the blue mason bee Osmia caerulescens and the sheep ked Melophagus ovinus respectively. A comprehensive phylogenomic analysis revealed two independent origins of filarial Wolbachia in supergroup F, likely from ancestral arthropod hosts, as well as a convergent pseudogenization and loss of the bacterioferritin gene across all filarial Wolbachia. The new genomes provide a valuable resource for further studies on symbiosis, evolution, and the discovery of new antibiotics to treat mansonellosis.

Data depositionRaw Illumina reads for the M. perstans isolate Mpe1 and the M. ozzardi isolate Moz2 have been deposited in the NCBI SRA database under BioProject accession numbers PRJNA666672 and PRJNA666671 respectively. The assembled genome and PGAP annotations for the Wolbachia wMpe and wMoz are available from the NCBI GenBank database with accession numbers JACZHU000000000 and JADAKK000000000 respectively. The annotated genomes for wOcae and wMoviF are available as GenBank accessions JAHDJT000000000 and JAHRBA000000000 respectively, deposited under BioProject accessions PRJNA244005 and PRJNA739984 respectively. RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub. Quantifying ERK-activity in response to inhibition of the BRAFV600E-MEK-ERK cascade using mathematical modelling BackgroundSimultaneous inhibition of multiple components of the BRAF-MEK-ERK cascade (vertical inhibition) has become a standard of care for treating BRAF-mutant melanoma. However, the molecular mechanism of how vertical inhibition synergistically suppresses intracellular ERK-activity, and consequently cell proliferation, are yet to be fully elucidated.

MethodsWe develop a mechanistic mathematical model that describes how the mutant BRAF-inhibitor, dabrafenib, and the MEK-inhibitor, trametinib, affect BRAFV600E-MEK-ERK signalling. The model is based upon a system of chemical reactions that describes cascade signalling dynamics. Using mass action kinetics, the chemical reactions are re-expressed as ordinary differential equations which are parameterised by in vitro data and solved numerically to obtain the temporal evolution of cascade component concentrations.

ResultsThe model provides a quantitative method to compute how dabrafenib and trametinib can be used in combination to synergistically inhibit ERK-activity in BRAFV600E-mutant melanoma cells. The model elucidates molecular mechanisms of vertical inhibition of the BRAFV600E-MEK-ERK cascade and delineates how elevated BRAF concentrations generate drug resistance to dabrafenib and trametinib. The computational simulations further suggest that elevated ATP levels could be a factor in drug resistance to dabrafenib.

ConclusionsThe model can be used to systematically motivate which dabrafenib-trametinib dose-combinations, for treating BRAFV600E-mutated melanoma, warrant experimental investigation. Classification of daily crop phenology in PhenoCams using deep learning and hidden markov models Near surface cameras, such as those in the PhenoCam network, are a common source of ground truth data in modelling and remote sensing studies. Despite having locations across numerous agricultural sites, few studies have used near surface cameras to track the unique phenology of croplands. Due to management activities, crops do not have a natural vegetation cycle which many phenological extraction methods are based on. For example, a field may experience abrupt changes due to harvesting and tillage throughout the year. A single camera can also record several different plants due to crop rotations, fallow fields, and cover crops. Current methods to estimate phenology metrics from image time series compress all image information into a relative greenness metric, which discards a large amount of contextual information. This can include the type of crop present, whether snow or water is present on the field, the crop phenology, or whether a field lacking green plants consists of bare soil, fully senesced plants, or plant residue. Here we developed a modelling workflow to create a daily time series of crop type and phenology, while also accounting for other factors such as obstructed images and snow covered fields. We used a mainstream deep learning image classification model, VGG16. Deep learning classification models do not have a temporal component, so to account for temporal correlation among images our workflow incorporates a hidden markov model in the post-processing. The initial image classification model had out of sample F1 scores of 0.83-0.85, which improved to 0.86-0.91 after all post-processing steps. The resulting time series show the progression of crops from emergence to harvest, and can serve as a daily, local scale dataset of field states and phenological stages for agricultural research. Discovery and functional assessment of a novel adipocyte population driven by intracellular Wnt/β-catenin signaling in mammals Wnt/{beta}-catenin signaling has been well established as a potent inhibitor of adipogenesis. Here, we identified a population of adipocytes that exhibit persistent activity of Wnt/{beta}-catenin signaling, as revealed by the Tcf/Lef-Gfp reporter allele, in embryonic and adult mouse fat depots, named as Wnt+ adipocytes. We showed that the {beta}-catenin-mediated signaling activation in these cells is Wnt ligand- and receptor-independent but relies on AKT/mTOR pathway and is essential for cell survival. Such adipocytes are distinct from classical ones in transcriptomic and genomic signatures and can be induced from various sources of mesenchymal stromal cells including human cells. Genetic lineage-tracing and targeted cell ablation studies revealed that these adipocytes convert into beige adipocytes directly and are also required for beige fat recruitment under thermal challenge, demonstrating both cell autonomous and non-cell autonomous roles in adaptive thermogenesis. Furthermore, mice bearing targeted ablation of these adipocytes exhibited glucose intolerance, while mice receiving exogenously supplied such cells manifested enhanced glucose utilization. Our studies uncover a unique adipocyte population in regulating beiging in adipose tissues and systemic glucose homeostasis. Liberation from fibrogenesis or tumorigenesis via cellular senescence by a novel small molecule BackgroundUncontrolled proliferative diseases such as fibrosis or cancer can be fatal. We previously found that a compound containing a chromone scaffold, ONG41008, had potent anti-fibrotic effects in diseased human lung myofibroblasts but not normal human lung fibroblasts.

MethodsWe investigated the effects of ONG41008 on tumor cells, and compared these effects with those in pathogenic myofibrotic cells and normal fibroblasts cells.

FindingsStimulation of A549 lung carcinoma epithelial cells with ONG41008 resulted in cellular senescence, indicating that dysregulated cell proliferation is common to fibrotic cells and tumor cells. Replicative senescence of A549 cells resulted in multinucleation, which was followed by oncogene-induced senescence. There was significant upregulation of expression and nuclear translocation of p-TP53 and p16 in ONG41008-treated A549 cells, and all cells died after 72 hr. Similar effects occurred after ONG41008 treatment in several human aggressive cancer cell lines such as PANC1, MCF7, PC3, or primary non-small cell lung carcinoma cells. Unlike cisplatin, ONG41008 was not toxic to normal human lung fibroblasts or primary prostate epithelial cells, suggesting ONG41008 can distinguish the intracellular microenvironment between normal cells and aged or diseased cells. This effect might occur as a result of the increased NAD/NADH ratio or increased lactate dehydrogenase levels in aged or diseased cells.

InterpretationTo our best knowledge, this is the first study to show that a small molecule can arrest uncontrolled proliferation during fibrogenesis or tumorigenesis in vitro. ONG41008 could be a potential drug for a broad range of fibrotic or tumorigenic diseases.

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSThe notion that aging is a disease and that diseases occur as a consequence of aging was first put forward by David Sinclair and colleagues (Aging Cell;14(4):497-510). There are biological systems that provide evidence for this notion; for example, somatic cells can revert to embryonic cells, producing younger somatic cells. This phenomenon underlies induced pluripotent stem cells. Another example is that some types of jellyfish can live forever. These examples suggest that a counter-aging program exists in animals. Human diseases are the manifestations of cell aging generated by the accumulation of somatic mutations. Aged and pathogenic cells are senescent, so a drug that specifically targeted senescent cells might initiate a cellular program that could ameliorate age-associated disease. Indeed, the kinase inhibitor dasatinib induces cellular senescence (Clin Ther 2007 29:2289-2308). In 2017, two drugs that target senescent cells were identified: quercetin and fisetin. These drugs selectively kill senescent cells, and are referred to as senolytic drugs or senotherapeutics (Aging 2017 8;9(3):955-963). Although it is established that senescent cells accumulate in cancer and idiopathic pulmonary fibrosis (Nat Commun 2017 23;8 :14532), the effect of senolytic drugs in these diseases is largely unknown.

Added value of this studyThis study characterized a novel drug, termed ONG41008, which was found to have both senogenic and senolytic effects in cell-based assays. ONG41008 induced senescence in myofibroblasts and several cancer cell lines representative of aggressive human cancers, which was followed by cell death. Importantly, ONG41008 exhibited essentially no toxicity on normal human lung fibroblasts or primary prostate epithelial cells.

Implication of all the evidenceBased on our results, we believe that ONG41008 is a potent inducer of cellular senescence (replicative senescence and oncogene-induced senescence) and causes arrest of uncontrolled, pathogenic proliferation of myofibroblasts or cancer cells. Identification of age and ethnicity specific gene expression biomarkers for immune aging Human immune system functions over an entire lifetime, yet how and why the immune system becomes less effective with age are not well understood. Here, we characterize peripheral blood mononuclear cells transcriptome from 172 healthy adults with 21~90 years of age using RNA-seq and the weighted gene correlation network analyses (WGCNA). These data reveal a set of insightful gene expression modules and representative gene biomarkers for human immune system aging from Asian and Caucasian ancestry, respectively. Among them, the aging-specific modules show an age-related gene expression variation spike around early-seventies. In addition, it is not known whether Asian and Caucasian immune systems go through similar gene expression changes throughout their lifespan, and to what extent these aging-associated changes are shared among ethnicities. We find the top hub genes including NUDT7, CLPB, OXNAD1 and MLLT3 are shared between Asian and Caucasian aging related modules and further validated in human PBMCs from different age groups. Overall, the impact of age and race on transcriptional variation elucidated from this study provide insights into the transcriptional driver of immune aging.