bluuebunny/crossref_metadata_2025 huggingface dataset:
A Random Forest Classifier for Protein-Protein Docking Models
Authors: Barradas-Bautista Didier , Cao Zhen , Vangone Anna , Oliva Romina , Cavallo Luigi
Published: June 2021

Herein, we present the results of a machine learning approach we developed to single out correct 3D docking models of protein-protein complexes obtained by popular docking software. To this aim, we generated a set of ≈7×106 docking models with three different docking programs (HADDOCK, FTDock and ZDOCK) for the 230 complexes in the protein-protein interaction benchmark, version 5 (BM5). Three different machine-learning approaches (Random Forest, Supported Vector Machine and Perceptron) were used to train classifiers with 158 different scoring functions (features). The Random Forest algorithm outperformed the other two algorithms and was selected for further optimization. Using a features selection algorithm, and optimizing the random forest hyperparameters, allowed us to train and validate a random forest classifier, named CoDES (COnservation Driven Expert System). Testing of CoDES on independent datasets, as well as results of its comparative performance with machine-learning methods recently developed in the field for the scoring of docking decoys, confirm its state-of-the-art ability to discriminate correct from incorrect decoys both in terms of global parameters and in terms of decoys ranked at the top positions.
Genomic history and forensic characteristics of Sherpa highlanders on the Tibetan Plateau inferred from high-resolution genome-wide InDels and SNPs
Authors: Wang Mengge , Du Weian , Tang Renkuan , Liu Yan , Zou Xing , Yuan Didi , Wang Zheng , Liu Jing , Guo Jianxin , Yang Xiaomin , Chen Jing , Yang Meiqing , Zhang Xianpeng , Wei Lan-Hai , Yuan Haibing , Yeh Hui-Yuan , Wang Chuan-Chao , Liu Chao , He Guanglin
Published: June 2021

Sherpa people, one of the high-altitude hypoxic adaptive populations, mainly reside in Nepal and the southern Tibet Autonomous Region. The genetic origin and detailed evolutionary profiles of Sherpas remain to be further explored and comprehensively characterized. Here we analyzed the newly-generated InDel genotype data from 628 Dingjie Sherpa people by merging with 4222 worldwide InDel profiles and collected genome-wide SNP data (approximately 600K SNPs) from 3324 individuals in 382 modern and ancient populations to explore and reconstruct the fine-scale genetic structure of Sherpas and their relationships with nearby modern and ancient East Asians based on the shared alleles and haplotypes. The forensic parameters of 57 autosomal InDels (A-InDels) included in our used new-generation InDel amplification system showed that this updated InDel panel is informative and polymorphic in Sherpas, suggesting that it can be used as the supplementary tool for forensic personal identification and parentage testing in the highland East Asians. Descriptive findings from the PCA, ADMIXTURE and TreeMix-based phylogeny suggested that Sherpas showed excess allele sharing with neighboring Tibeto-Burman Tibetans. Furthermore, patterns of allele sharing in f-statistics demonstrated that Sherpa people had a different evolutionary history compared with their neighbors from Nepal (Newar and Gurung) but showed genetic similarity with 2700-year-old Chokhopani and modern Tibet Tibetans. QpAdm/qpGraph-based admixture sources and models further showed that Sherpa, core Tibetans and Chokhopani formed one clade which could be fitted as having the main ancestry from late Neolithic Qijia millet farmers and other deep ancestries from early Asians. Chromosome painting profiles and shared IBD fragments inferred from FineStructure and ChromoPainter not only confirmed the abovementioned genomic affinity patterns but also revealed the fine-scale microstructures among Sino-Tibetan speakers. Finally, natural-selection signals revealed via iHS, nSL, and iHH12 showed signatures associated with disease susceptibility in Sherpa people. Generally, we provided the comprehensive landscape of admixture and evolutionary history of Sherpa people based on the shared alleles and haplotypes from the low-density forensic markers and high-density genome-wide SNP data. The more detailed genetic landscape of Sherpa people should be further confirmed and characterized via ancient genomes or single-molecule real-time sequencing technology.
Sex-specific effects of polygenic risk for schizophrenia on lifespan cognitive functioning in healthy individuals
Authors: Koch Elise , Nyberg Lars , Lundquist Anders , Pudas Sara , Adolfsson Rolf , Kauppi Karolina
Published: June 2021

Polygenic risk for schizophrenia has been associated with lower cognitive ability and age-related cognitive change in healthy individuals. Despite well-established neuropsychological sex differences in schizophrenia patients, genetic studies on sex differences in schizophrenia in relation to cognitive phenotypes are scarce. Here, we investigated whether the effect of a polygenic risk score (PRS) for schizophrenia on childhood, midlife and late life cognitive function in healthy individuals is modified by sex, and if PRS is linked to accelerated cognitive decline. Using a longitudinal data set from healthy individuals aged 25-100 years (N = 1,459) spanning a 25-year period, we found that PRS was associated with lower cognitive ability (episodic memory, semantic memory, visuospatial ability), but not with accelerated cognitive decline. A significant interaction effect between sex and PRS was seen on cognitive task performance, and sex-stratified analyses showed that the effect of PRS was male-specific. In a sub-sample, we observed a male-specific effect of the PRS on school performance at age 12 (N = 496). Our findings of sex-specific effects of schizophrenia genetics on cognitive functioning across the life-span indicate that the effects of underlying disease genetics on cognitive functioning is dependent on biological processes that differ between the sexes.
Local temporal Rac1-GTP nadirs and peaks restrict cell protrusions and retractions
Authors: Hu Jianjiang , Gong Xiaowei , Strömblad Staffan
Published: June 2021

Spatiotemporal coordination of the GTP-binding activity of Rac1 and RhoA initiates and reinforces cell membrane protrusions and retractions during cell migration1–7. However, while protrusions and retractions form cycles that cells use to efficiently probe their microenvironment8–10, the control of their finite lifetime remains unclear. To examine if Rac1 or RhoA may also control protrusion and retraction lifetimes, we here define the relation of their spatiotemporal GTP-binding levels to key protrusion and retraction events, as well as to cell-ECM mechanical forces in fibrosarcoma cells grown on collagen of physiologically relevant stiffness. We identified temporal Rac1-GTP nadirs and peaks at the maximal edge velocity of local membrane protrusions and retractions, respectively, followed by declined edge velocity. Moreover, increased local Rac1-GTP consistently preceded increased cell-ECM traction force. This suggests that Rac1-GTP nadirs and peaks may restrain the lifetime of protrusions and retractions, possibly involving the regulation of local traction forces. Functional testing by optogenetics validated this notion, since local Rac1-GTP elevation applied early in the process prolonged protrusions and restrained retractions, while local Rac1-GTP inhibition acted in reverse. Optogenetics also defined Rac1-GTP as a promotor of local traction force. Together, we show that Rac1 plays a fundamental role in restricting the size and durability of protrusions and retractions, plausibly in part through controlling traction forces.
LSD induces increased signalling entropy in rats’ prefrontal cortex
Authors: Savino Aurora , Nichols Charles D.
Published: June 2021

Psychedelic drugs are gaining attention from the scientific community as potential new compounds for the treatment of psychiatric diseases such as mood and substance use disorders. The 5-HT2A receptor has been identified as the main molecular target, and early studies pointed to an effect on the expression of neuroplasticity genes. Analysing RNA-seq data from the prefrontal cortex of rats chronically treated with lysergic acid diethylamide (LSD), we describe the psychedelic-induced rewiring of gene co-expression networks, which become less centralized but more complex, with an overall increase in signalling entropy, typical of highly plastic systems. Intriguingly, signalling entropy mirrors, at the molecular level, the increased brain entropy reported through neuroimaging studies in human, suggesting the underlying mechanisms of higher-order phenomena. Moreover, from the analysis of network topology we identify potential transcriptional regulators and imply different cell types in psychedelics’ activity.
Probing the rules of cell coordination in live tissues by interpretable machine learning based on graph neural networks
Authors: Yamamoto Takaki , Cockburn Katie , Greco Valentina , Kawaguchi Kyogo
Published: June 2021

Robustness in developing and homeostatic tissues is supported by various types of spatiotemporal cell-to-cell interactions. Although live imaging and cell tracking are powerful in providing direct evidence of cell coordination rules, extracting and comparing these rules across many tissues with potentially different length and timescales of coordination requires a versatile framework of analysis. Here we demonstrate that graph neural network (GNN) models are suited for this purpose, by showing how they can be applied to predict cell fate in tissues and utilized to infer the cell interactions governing the multicellular dynamics. Analyzing the live mammalian epidermis data, where spatiotemporal graphs constructed from cell tracks and cell contacts are given as inputs, GNN discovers distinct neighbor cell fate coordination rules that depend on the region of the body. This approach demonstrates how the GNN framework is powerful in inferring general cell interaction rules from live data without prior knowledge of the signaling involved.
Panel Informativity Optimizer (PIO): an R package to improve cancer NGS panel informativity
Authors: Alcazer Vincent , Sujobert Pierre
Published: June 2021

Mutation detection by next generation sequencing (NGS) is routinely used for cancer diagnosis. Selecting an optimal set of genes for a given cancer is not trivial as it has to optimize informativity (i.e. the number of patients with at least one mutation in the panel), while minimizing panel length in order to reduce sequencing costs and increase sensitivity. We propose herein Panel Informativity Optimizer (PIO), an open-source software developed as an R package with a user-friendly graphical interface to help optimize cancer NGS panel informativity. Using patient-level mutational data from either private datasets or preloaded dataset of 91 independent cohort from 31 different cancer type, PIO selects an optimal set of genomic intervals to maximize informativity and panel size in a given cancer type. Different options are offered such as the definition of genomic intervals at the gene or exon level, and the use of optimization strategy at the patient or patient per kilobase level. PIO can also propose an optimal set of genomic intervals to increase informativity of custom panels. A panel tester function is also available for panel benchmarking. Using public databases, as well as data from real-life settings, we demonstrate that PIO allows panel size reduction of up to 1000kb, and accurately predicts the performance of custom or commercial panels. PIO is available online at https://vincentalcazer.shinyapps.io/Panel_informativity_optimizer/ or can be set on a locale machine from https://github.com/VincentAlcazer/PIO.
Entamoeba histolytica adaption to auranofin: a phenotypic and multi-omics characterization
Authors: Shaulov Yana , Sarid Lotem , Trebicz-Geffen Meirav , Ankri Serge
Published: June 2021

Auranofin (AF), an antirheumatic agent, targets mammalian thioredoxin reductase (TrxR), an important enzyme controlling redox homeostasisis, AF is also very effective against a diversity of pathogenic bacteria and protozoan parasites. Here, we report about the resistance of the parasite Entamoeba histolytica to 2 μM of AF that has been acquired by gradual exposure of the parasite to increasing amount of the drug. AF adapted E.histolytica trophozoites (AFAT) has an impaired growth, cytopathic activity and they are more sensitive to oxidative stress (OS), nitrosative stress (NS) and metronidazole (MTZ) than wild type (WT) trophozoites. Integrated transcriptomics and redoxomics analyses showed that many upregulated genes in AFAT, including genes encoding for dehydrogenase and cytoskeletal proteins, have their product oxidized in wild type trophozoites exposed to AF (acute AF trophozoites) but not in AFAT. We also showed that the level of reactive oxygen species (ROS) and oxidized proteins (OXs) in AFAT is lower than that of acute AF trophozoites. Overexpression of E.histolytica TrxR (EhTrxR) did not protect the parasite against AF which suggests that EhTrxR is not central is the mechanism of adaptation to AF.
A live-cell marker to visualize the dynamics of stable microtubules
Authors: Jansen Klara I. , Burute Mithila , Kapitein Lukas C.
Published: June 2021

The microtubule (MT) cytoskeleton underlies processes such as intracellular transport and cell division. Immunolabeling for post-translational modifications of tubulin has revealed the presence of different MT subsets, which are believed to differ in stability and function. Whereas dynamic MTs can readily be studied using live-cell plus-end markers, the dynamics of stable MTs have remained obscure due to a lack of tools to directly visualize these MTs in living cells. Here, we present a live-cell marker to visualize stable MTs and explore their dynamics. We demonstrate that a rigor mutant of kinesin-1 binds selectively to acetylated MTs without affecting MT organization and organelle transport. These MTs are long-lived, do not depolymerize upon nocadozale-treatment or laser-based severing, and display rich dynamics, including undulation, looping and sliding. This marker will help to explore how different MT subsets contribute to cellular organization and transport.
SARS-CoV-2 neutralising antibodies in Dogs and Cats in the United Kingdom
Authors: Smith Shirley L. , Anderson Enyia. R. , Cansado-Utrilla Cintia , Prince Tessa , Farrell Sean , Brant Bethaney , Smyth Stephen , Noble Peter-John M. , Pinchbeck Gina L. , Marshall Nikki , Roberts Larry , Hughes Grant L. , Radford Alan D. , Patterson Edward I.
Published: June 2021

Companion animals are susceptible to SARS-CoV-2 infection and sporadic cases of pet infections have occurred in the United Kingdom. Here we present the first large-scale serological survey of SARS-CoV-2 neutralising antibodies in dogs and cats in the UK. Results are reported for 688 sera (454 canine, 234 feline) collected by a large veterinary diagnostic laboratory for routine haematology during three time periods; pre-COVID-19 (January 2020), during the first wave of UK human infections (April-May 2020) and during the second wave of UK human infections (September 2020-February 2021). Both pre-COVID-19 sera and those from the first wave tested negative. However, in sera collected during the second wave, 1.4% (n=4) of dogs and 2.2% (n=2) cats tested positive for neutralising antibodies. The low numbers of animals testing positive suggests pet animals are unlikely to be a major reservoir for human infection in the UK. However, continued surveillance of in-contact susceptible animals should be performed as part of ongoing population health surveillance initiatives.
How Communication Pathways Bridge Local and Global Conformations in an IgG4 Antibody: a Molecular Dynamics Study
Authors: Tarenzi Thomas , Rigoli Marta , Potestio Raffaello
Published: June 2021

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recently, however, evidence has accumulated that points toward a nontrivial relation between the CDR and distal sites on the antibody structure: variations in the binding strengths have been observed upon mutating amino acids separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4μsof MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR; a pivotal role in this relationship is played by the relatively short hinge, which mechanically couples Fab and Fc domains. These results support the hypothesis that pembrolizumab behaves as a complex machinery, with a multi-scale hierarchy of global and local conformational changes that communicate with one another. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

Synopsis

Antigen binding restricts the conformational variability of the therapeutic antibody pembrolizumab in an interplay between the paratope and hinge region, mediated by a full-scale interaction network.

Graphical TOC Entry
A rare natural lipid induces neuroglobin expression to prevent amyloid oligomers toxicity and retinal neurodegeneration
Authors: Oamen Henry Patrick , Romero Nathaly Romero , Knuckles Philip , Saarikangas Juha , Dong Yuhong , Caudron Fabrice
Published: June 2021

Most neurodegenerative diseases such as Alzheimer’s disease are proteinopathies linked to the toxicity of amyloid oligomers. Treatments to delay or cure these diseases are lacking. Using budding yeast, we report that the natural lipid tripentadecanoin induces expression of the nitric oxide oxidoreductase Yhb1 to prevent the formation of protein aggregates during aging and extends replicative lifespan. In mammals, tripentadecanoin induces expression of the Yhb1 orthologue, neuroglobin, to protect neurons against amyloid toxicity. Tripentadecanoin also rescues photoreceptors in a mouse model of retinal degeneration and retinal ganglion cells in a Rhesus monkey model of optic atrophy. Together, we propose that tripentadecanoin affects p-bodies to induce neuroglobin expression and offers a potential treatment for proteinopathies and retinal neurodegeneration.

One Sentence Summary

The natural lipid tripentadecanoin is cytoprotective against amyloid oligomer toxicity and retinal neurodegeneration by inducing YHBI/neuroglobin expression in yeast and mammals.
Designing BH3-mimetic Peptide Inhibitors for the Viral Bcl-2 Homologs A179L and BHRF1: Importance of long-range electrostatic interactions
Authors: Reddy C. Narendra , Sankararamakrishnan Ramasubbu
Published: June 2021

Viruses have evolved strategies to prevent apoptosis of infected cells at early stages of infection. The viral proteins (vBcl-2s) from specific viral genes adopt a helical fold that is structurally similar to that of mammalian anti-apoptotic Bcl-2 proteins and exhibit little sequence similarity. Hence vBcl-2 homologs are attractive targets to prevent viral infection. However, very few studies have focused on developing inhibitors for vBcl-2 homologs. In this study, we have considered two vBcl-2 homologs, A179L from African swine fever virus and BHRF1 from Epstein-Barr virus. We generated two sets of 8000 randomized BH3-like sequences from eight wild-type pro-apoptotic BH3 peptides. During this process, the four conserved hydrophobic residues and an Asp residue were retained at their respective positions and all other positions were substituted randomly without any bias. We constructed 8000 structures each for A179L and BHRF1 in complex with BH3-like sequences. Histograms of interaction energies calculated between the peptide and the protein resulted in negatively skewed distributions. The BH3-like peptides with high helical propensities selected from the negative tail of respective interaction energy distributions exhibited more favorable interactions with A179L and BHRF1 and they are rich in basic residues. Molecular dynamics studies and electrostatic potential maps further revealed that both acidic and basic residues favorably interact with A179L while only basic residues have the most favorable interactions with BHRF1. As in mammalian homologs, the role of long range interactions and non-hotspot residues have to be taken into account while designing specific BH3-mimetic inhibitors for vBcl-2 homologs.
A Bayesian Multivariate Mixture Model for Spatial Transcriptomics Data
Authors: Allen Carter , Chang Yuzhou , Neelon Brian , Chang Won , Kim Hang J. , Li Zihai , Ma Qin , Chung Dongjun
Published: June 2021

High throughput spatial transcriptomics (HST) is a rapidly emerging class of experimental technologies that allow for profiling gene expression in tissue samples at or near single-cell resolution while retaining the spatial location of each sequencing unit within the tissue sample. Through analyzing HST data, we seek to identify sub-populations within a tissue sample that reflect distinct cell types or states. Existing methods either ignore the spatial heterogeneity in gene expression profiles, fail to account for important statistical features such as skewness, or are heuristic network-based clustering methods that lack the inferential benefits of statistical modeling. To address this gap, we develop SPRUCE: a Bayesian spatial multivariate finite mixture model based on multivariate skew-normal distributions, which is capable of identifying distinct cellular sub-populations in HST data. We further implement a novel combination of Pólya–Gamma data augmentation and spatial random effects to infer spatially correlated mixture component membership probabilities without relying on approximate inference techniques. Via a simulation study, we demonstrate the detrimental inferential effects of ignoring skewness or spatial correlation in HST data. Using publicly available human brain HST data, SPRUCE outperforms existing methods in recovering expertly annotated brain layers. Finally, our application of SPRUCE to human breast cancer HST data indicates that SPRUCE can distinguish distinct cell populations within the tumor microenvironment.
Mice are not automatons; subjective experience in premotor circuits guides behavior
Authors: Schreiner Drew C. , Cazares Christian , Renteria Rafael , Gremel Christina M.
Published: June 2021

Subjective experience is a powerful driver of decision-making and continuously accrues. However, most neurobiological studies constrain analyses to task-related variables and ignore how continuously and individually experienced internal, temporal, and contextual factors influence adaptive behavior during decision-making and the associated neural mechanisms. We show mice rely on learned information about recent and longer-term subjective experience of variables above and beyond prior actions and reward, including checking behavior and the passage of time, to guide self-initiated, self-paced, and self-generated actions. These experiential variables were represented in secondary motor cortex (M2) activity and its projections into dorsal medial striatum (DMS). M2 integrated this information to bias strategy-level decision-making, and DMS projections used specific aspects of this recent experience to plan upcoming actions. This suggests diverse aspects of experience drive decision-making and its neural representation, and shows premotor corticostriatal circuits are crucial for using selective aspects of experiential information to guide adaptive behavior.
Automated classification of bat echolocation call recordings with artificial intelligence
Authors: Tabak Michael A. , Murray Kevin L. , Lombardi John A. , Bay Kimberly J.
Published: June 2021

Acoustic recorders are commonly used to remotely monitor and collect data on bats (Order Chiroptera). These efforts result in many acoustic recordings that must be classified by a bat biologist with expertise in call classification in order to obtain useful information. The rarity of this expertise and time constraints have prompted efforts to automatically classify bat species in acoustic recordings using a variety of learning methods. There are several software programs available for this purpose, but they are imperfect and the United States Fish and Wildlife Service often recommends that a qualified acoustic analyst review bat call identifications even if using these software programs. We sought to build a model to classify bat species using modern computer vision techniques. We used images of bat echolocation calls (i.e., plots of the pulses) to train deep learning computer vision models that automatically classify bat calls to species. Our model classifies 10 species, five of which are protected under the Endangered Species Act. We evaluated our models using standard model validation procedures, but we also performed two out-of-distribution tests. For the out-of-distribution tests, an entire dataset was withheld from the procedure before splitting the data into training and validation sets. We found that our validation accuracy (93%) and out-of-distribution accuracy (90%) were higher than when we used Kaleidoscope Pro and BCID software (65% and 61% accuracy, respectively) to evaluate the same calls. Our results suggest that our approach is effective at classifying bat species from acoustic recordings, and our trained model will be incorporated into new bat call identification software: WEST-EchoVision.
Hippocampal place cells use vector computations to navigate
Authors: Ormond Jake , O’Keefe John
Published: June 2021

One function of the Hippocampal Cognitive Map is to provide information about salient locations in familiar environments such as those containing reward or danger, and to support navigation towards or away from those locations1. Although much is known about how the hippocampus encodes location in world-centred coordinates, how it supports flexible navigation is less well understood. We recorded from CA1 place cells while rats navigated to a goal or freely foraged on the honeycomb maze2. The maze tests the animal’s ability to navigate using indirect as well as direct paths to the goal and allows the directionality of place cells to be assessed at each choice point during traversal to the goal. Place fields showed strong directional polarization in the navigation task, and to a lesser extent during random foraging. This polarization was characterized by vector fields which converged to sinks distributed throughout the environment. The distribution of these convergence sinks was centred near the goal location, and the population vector field converged on the goal, providing a strong navigational signal. Changing the goal location led to the movement of ConSinks and vector fields towards the new goal and within-days, the ConSink distance to the goal decreased with continued training. The honeycomb maze allows the independent assessment of spatial representation and spatial action in place cell activity and shows how the latter depends on the former. The results suggest a vector-based model of how the hippocampus supports flexible navigation, allowing animals to select optimal paths to destinations from any location in the environment.
Modular, Cascade-like Transcriptional Program of Regeneration inStentor
Authors: Sood Pranidhi , Lin Athena , Yan Connie , McGillivary Rebecca , Diaz Ulises , Makushok Tatyana , Marshall Wallace F.
Published: June 2021

The giant ciliateStentor coeruleusis a classical model system for studying regeneration and morphogenesis at the level of a single cell. The anterior of the cell is marked by an array of cilia, known as the oral apparatus, which can be induced to shed and regenerate in a series of reproducible morphological steps, previously shown to require transcription. If a cell is cut in half, each half will regenerate an intact cell, including a new oral apparatus in the posterior half. We used RNAseq to assay the dynamic changes inStentor’s transcriptome during regeneration, after both oral apparatus shedding and bisection, allowing us to identify distinct temporal waves of gene expression including kinases, RNA binding proteins, centriole biogenesis factors, and orthologs of human ciliopathy genes implicated in Meckel and Joubert syndromes. By comparing transcriptional profiles of different regeneration events in the same species, we were able to identify distinct modules of gene expression corresponding to oral apparatus regeneration, posterior holdfast regeneration, and recovery after wounding. By measuring gene expression in cells in which translation is blocked, we show that the sequential waves of gene expression involve a cascade mechanism in which later waves of expression are triggered by translation products of early-expressed genes. Among the early-expressed genes, we identified an E2F transcription factor and the conserved RNA binding protein Pumilio as potential regulators of regeneration based on the expression pattern of their predicted target genes. RNAi mediated knockdown experiments indicate that Pumilio is required for regenerating oral structures of the correct size. We show that E2F is involved in the completion of regeneration but is dispensable for earlier steps. This work allows us to classify regeneration genes into groups based on their potential role for regeneration in distinct cell regeneration paradigms, and provides insight into how a single cell can coordinate complex morphogenetic pathways to regenerate missing structures.
E484K and N501Y SARS-CoV 2 Spike Mutants Increase ACE2 Recognition but Reduce Affinity for Neutralizing Antibody
Authors: Chakraborty Sandipan
Published: June 2021

SARS-CoV2 mutants emerge as variants of concern (VOC) due to altered selection pressure and rapid replication kinetics. Among them, lineages B.1.1.7, B.1.351, and P.1 contain a key mutation N501Y. B.1.135 and P.1 lineages have another mutation, E484K. Here, we decode the effect of these two mutations on the host receptor, ACE2, and neutralizing antibody (B38) recognition. The gain in binding affinity for the N501Y RBD mutant to the ACE2 is attributed to improved π-π stacking and π-cation interactions. The enhanced receptor affinity of the E484K mutant is caused due to the formation of a specific hydrogen bond and salt-bridge interaction with Glu75 of ACE2. Notably, both the mutations reduce the binding affinity for B38 due to the loss of several hydrogen-bonding interactions. The insights obtained from the study are crucial to interpret the increased transmissibility and reduction in the neutralization efficacy of rapidly emerging SARS-CoV2 VOCs.
A Multi-omics Method Enabled by Sequential Metabolomics and Proteomics for Human Pluripotent Stem Cell-derived Cardiomyocytes
Authors: Bayne Elizabeth F. , Simmons Aaron D. , Roberts David S. , Zhu Yanlong , Aballo Timothy J. , Wancewicz Benjamin , Palecek Sean P. , Ge Ying
Published: June 2021

Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) show immense promise for patient-specific disease modeling, cardiotoxicity screening, and regenerative therapy development. However, hPSC-CMs in culture have not recapitulated the structural or functional properties of adult CMs in vivo thus far. To gain global insight into hPSC-CM biology, we established a multi-omics method for analyzing the hPSC-CM metabolome and proteome from the same cell culture, creating multi-dimensional profiles of hPSC-CMs. Specifically, we developed a sequential extraction to capture metabolites and proteins from the same hPSC-CM monolayer cultures, and analyzed these extracts using high-resolution mass spectrometry (MS). Using this method, we annotated 205 metabolites/lipids and 4,008 proteins from 106 cells with high reproducibility. We further integrated the proteome and metabolome measurements to create network profiles of molecular phenotypes for hPSC-CMs. Out of 310 pathways identified using metabolomics and proteomics, 40 pathways were considered significantly overrepresented (FDR-corrected p ≤ 0.05). Highly populated pathways included those involved in protein synthesis (ribosome, spliceosome), ATP generation (oxidative phosphorylation), and cardiac muscle contraction. This multi-omics method achieves deep coverage of metabolites and proteins, creating a multidimensional view of the hPSC-CM phenotype, which provides a strong technological foundation to advance the understanding of hPSC-CM biology.
Bioaffinity-based surface-immobilization of antibodies to capture endothelial colony-forming cells
Authors: Boulanger Mariève D. , Elkhodiry Mohamed A. , Bashth Omar S. , Laroche Gaétan , Hoesli Corinne A.
Published: June 2021

Maximizing the re-endothelialization of vascular implants such as prostheses or stents has the potential to significantly improve their long-term performance. Endothelial progenitor cell capture stents with surface-immobilized antibodies show significantly improved endothelialization in the clinic. However, most current antibody-based stent surface modification strategies rely on antibody adsorption or direct conjugation via amino or carboxyl groups which leads to poor control over antibody surface concentration and/or molecular orientation, and ultimately bioavailability for cell capture. Here, we assess the utility of a bioaffinity-based surface modification strategy consisting of a surface-conjugated cysteine-tagged protein G molecules that immobilize Immunoglobulin G (IgG) antibodies via the Fc domain to capture circulating endothelial colony-forming cells (ECFCs). The cysteine-tagged protein G was grafted onto aminated substrates at different concentrations as detected by an enzyme-linked immunosorbent assay and fluorescence imaging. Different IgG antibodies were successfully immobilized on the protein G-modified surfaces and higher antibody surface concentrations were achieved compared to passive adsorption methods. Surfaces with immobilized antibodies targeting endothelial surface proteins, such as CD144, significantly enhanced the capture of circulating ECFCsin vitrocompared to surfaces with non-endothelial specific antibodies such as anti-CD14. This work presents a potential avenue for enhancing the clinical performance of vascular implants by using covalent grafting of protein G to immobilize IgG antibodies more effectively.

Table of Contents

Antibody immobilization via surface-conjugated recombinant cysteine-protein G provides an effective approach to capture circulating therapeutic cells.
Genetic Loci and Metabolic States Associated With Murine Epigenetic Aging
Authors: Mozhui Khyobeni , Lu Ake T. , Li Caesar Z , Haghani Amin , Sandoval-Sierra Jose Vladimir , Wu Yibo , Williams Robert W. , Horvath Steve
Published: June 2021

Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a ‘pan-mammalian’ microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs, and show that high fat diet augments the age-associated changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate rates of epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and a second on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed associations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging.
Aberrant FGF signaling promotes granule neuron precursor expansion in SHH subgroup infantile medulloblastoma
Authors: Yabut Odessa R. , Gomez Hector , Arela Jessica , Castillo Jesse Garcia , Ngo Thomas , Pleasure Samuel J.
Published: June 2021

Mutations in Sonic Hedgehog (SHH) signaling pathway genes, e.g.,Suppressor of Fused(SUFU), drive granule neuron precursors (GNP) to form medulloblastomas (MBSHH). However, how different molecular lesions in the Shh pathway drive transformation is frequently unclear, andSUFUmutations in the cerebellum seem distinct. In this study, we show that fibroblast growth factor 5 (FGF5) signaling is integral for many infantile MBSHHcases and thatFGF5expression is uniquely upregulated in infantile MBSHHtumors. Similarly, mice lacking SUFU (Sufu-cKO) ectopically expressFgf5specifically along the secondary fissure where GNPs harbor preneoplastic lesions and show that FGFR signaling is also ectopically activated in this region. Treatment with an FGFR antagonist rescues the severe GNP hyperplasia and restores cerebellar architecture. Thus, direct inhibition of FGF signaling may be a promising and novel therapeutic candidate for infantile MBSHH.
Parallel Evaluation of Nucleophilic and Electrophilic Chemical Probes for Sulfenic Acid: Reactivity, Selectivity and Biocompatibility
Authors: Shi Yunlong , Carroll Kate S.
Published: June 2021

S-sulfenylation of cysteine thiols (Cys-SOH) is a regulatory posttranslational modification in redox signaling and an important intermediate to other cysteine chemotypes. Owing to the dual chemical nature of the sulfur in sulfenic acid, both nucleophilic and electrophilic chemical probes have been developed to react with and detect Cys-SOH; however, the efficiency of existing probes has not been evaluated in a side-by-side comparison. Here, we employ small-molecule and protein models of Cys-SOH and compare the chemical probe reactivity. These data clearly show that 1,3-diketone-based nucleophilic probes react more efficiently with sulfenic acid as compared to strained alkene/alkyne electrophilic probes. Kinetic experiments that rigorously address the selectivity of the 1,3-diketone-based probes are also reported. Consideration of these data alongside relative cellular abundance, indicates that biological electrophiles, including cyclic sulfenamides, aldehydes, disulfides and hydrogen peroxide, are not meaningful targets of 1,3-diketone-based nucleophilic probes, which still remain the most viable tools for the bioorthogonal detection of Cys-SOH.
Repertoire analyses reveal TCR sequence features that influence T cell fate
Authors: Lagattuta Kaitlyn A. , Kang Joyce B. , Nathan Aparna , Pauken Kristen E. , Jonsson Anna Helena , Rao Deepak A. , Ishigaki Kazuyoshi , Raychaudhuri Soumya
Published: June 2021

T cells acquire a regulatory phenotype when their T cell receptors (TCRs) experience an intermediate-high affinity interaction with a self-peptide presented on MHC. Using TCR sequences from FACS-sorted human cells, we identified TCR features that shape affinity to these self-peptide-MHC complexes, finding that 1) CDR3β hydrophobicity and 2) certain TRBV genes promote Treg fate. We developed a scoring system for TCR-intrinsic regulatory potential (TiRP) and found that within the tumor microenvironment clones exhibiting Treg-Tconv plasticity had higher TiRP than expanded clones maintaining the Tconv phenotype. To elucidate drivers of these predictive TCR features, we examined the two elements of the Treg TCR ligand separately: the self-peptide via murine Tregs, and the human MHC II molecule via human memory Tconvs. These analyses revealed that CDR3β hydrophobicity promotes reactivity to self-peptides, while TRBV gene usage shapes the TCR’s general propensity for MHC II-restricted activation.
Why are SNAREpins rod-shaped?
Authors: Zeng Jin , McDargh Zachary , An Dong , O’Shaughnessy Ben
Published: June 2021

SNARE proteins are the core components of the cellular machineries that fuse membranes for neurotransmitter or hormone release and other fundamental processes. Fusion is accomplished when SNARE proteins hosted by apposing membranes form SNARE complexes called SNAREpins, but the mechanism of fusion remains unclear. Computational simulations of SNARE-mediated membrane fusion are challenging due to the millisecond timescales of physiological membrane fusion. Here we used ultra-coarse-grained (UCG) simulations to investigate the minimal requirements for a molecular intracellular fusogen, and to elucidate the mechanisms of SNARE-mediated fusion. We find fusion by simple body forces that push vesicles together is highly inefficient. Inter-vesicle fusogens with different aspect ratios can fuse vesicles only if they are rodlike, of sufficient length to clear the fusogens from the fusion site by entropic forces. Simulations with rod-shaped SNAREpin-like fusogens fused 50-nm vesicles on ms timescales, driven by entropic forces along a reproducible fusion pathway. SNARE-SNARE and SNARE-membrane entropic forces cleared the fusion site and pressed the vesicles into an extended contact zone (ECZ), drove stalk nucleation at the high curvature ECZ boundary, and expanded the stalk into a long-lived hemifusion diaphragm in which a simple pore completed fusion. Our results provide strong support for the entropic hypothesis of SNARE-mediated membrane fusion, and implicate the rodlike structure of the SNAREpin complex as a necessity for entropic force production and fusion.
Cardiovascular disease causes proinflammatory microvascular changes in the human right atrium
Authors: Linna-Kuosmanen Suvi , Schmauch Eloi , Galani Kyriakitsa , Ojanen Johannes , Boix Carles A. , Örd Tiit , Toropainen Anu , Singha Prosanta K. , Moreau Pierre R. , Harju Kristiina , Blazeski Adriana , Segerstolpe Åsa , Lahtinen Veikko , Hou Lei , Kang Kai , Meibalan Elamaran , Agudelo Leandro Z. , Kokki Hannu , Halonen Jari , Jalkanen Juho , Gunn Jarmo , MacRae Calum A. , Hollmén Maija , Hartikainen Juha , Kaikkonen Minna U. , García-Cardeña Guillermo , Tavi Pasi , Kiviniemi Tuomas , Kellis Manolis
Published: June 2021

Ischemic heart disease is globally the leading cause of death. It plays a central role in the electrical and structural remodeling of the right atrium, predisposing to arrhythmias, heart failure, and sudden death. Here, we provide the first dissection of the gene expression changes in the live right atrial tissue, using single-nuclei RNA-seq and spatial transcriptomics. We investigate matched samples of the tissue and pericardial fluid and reveal substantial differences in disease- associated gene expression in all cell types, leading to inflammatory microvascular dysfunction and changes in the tissue composition. Our study demonstrates the importance of creating high- resolution cellular maps and partitioning disease signals beyond epicardial coronary arteries and ischemic left ventricle to identify candidate mechanisms leading to more severe types of human cardiovascular disease.

One-Sentence Summary

Single-cell dissection ofex vivoheart biopsies and pericardial fluid in ischemic heart disease and heart failure
Dynamic aspects of plant water potential revealed by a microtensiometer
Authors: Pagay Vinay
Published: June 2021

Water potential is a fundamental thermodynamic parameter that describes the activity of water. In this paper, we describe the continuous measurement of plant water potential, a reliable indicator of its water status, using a novelin situsensor known as a ‘microtensiometer’ in mature grapevines under field conditions. The microtensiometer operates on the principle of equilibration of water potentials of internal liquid water with an external vapour or liquid phase. We characterised the seasonal and diurnal dynamics of trunk water potentials (Ψtrunk) obtained from microtensiometers installed in two grapevine cultivars, Shiraz and Cabernet Sauvignon, and compared these values to pressure chamber-derived stem (Ψstem) and leaf (Ψleaf) water potentials as well as leaf stomatal conductance. Diurnal patterns of Ψtrunkmatched those of Ψstemand Ψleafunder low vapour pressure deficit (VPD) conditions, but diverged under high VPD conditions. The highest diurnal values of Ψtrunkwere observed shortly after dawn, while the lowest values were typically observed in the late afternoon. Differential responses of Ψtrunkto VPD were observed between cultivars, with Shiraz more sensitive than Cabernet to increasing VPD over long time scales, and both cultivars had a stronger VPD response than soil moisture response. On a diurnal basis, however, time cross correlation analysis revealed that Shiraz Ψtrunklagged Cabernet Ψtrunkin response to changing VPD. Microtensiometers were shown to operate reliably under field conditions over several months. To be useful for irrigation scheduling of woody crops, new thresholds of Ψtrunkneed to be developed.
Multiple mechanisms impact fluconazole resistance of mutant Erg11 proteins in Candida glabrata
Authors: Vu Bao Gia , Moye-Rowley W. Scott
Published: June 2021

Azoles, the most commonly used antifungal drugs, specifically inhibit the fungal lanosterol α-14 demethylase enzyme, which is referred to as Erg11. Inhibition of Erg11 ultimately leads to a reduction in ergosterol production, an essential fungal membrane sterol. Many Candida species, such as Candida albicans, develop mutations in this enzyme which reduces the azole binding affinity and results in increased resistance. Candida glabrata is also a pathogenic yeast that has low intrinsic susceptibility to azole drugs and easily develops elevated resistance. These azole resistant mutations typically cause hyperactivity of the Pdr1 transcription factor and rarely lie within the ERG11 gene. Here, we generated C. glabrata ERG11 mutations that were analogous to azole resistance alleles from C. albicans ERG11. Three different Erg11 forms (Y141H, S410F, and the corresponding double mutant (DM)) conferred azole resistance in C. glabrata with the DM Erg11 form causing the strongest phenotype. The DM Erg11 also induced cross-resistance to amphotericin B and caspofungin. Resistance caused by the DM allele of ERG11 imposed a fitness cost that was not observed with hyperactive PDR1 alleles. Crucially, the presence of the DM ERG11 allele was sufficient to activate the Pdr1 transcription factor in the absence of azole drugs. Our data indicate that azole resistance linked to changes in ERG11 activity can involve cellular effects beyond an alteration in this key azole target enzyme. Understanding the physiology linking ergosterol biosynthesis with Pdr1-mediated regulation of azole resistance is crucial for ensuring the continued efficacy of azole drugs against C. glabrata.

Importance

Azole drugs target the Erg11 enzyme and lead to a reduction in fungal ergosterol, a vital sterol in yeast. Mutations in Erg11 are common among azole resistant Candida albicans clinical isolates, but not in C. glabrata, a major human pathogen. In this study, we showed that ERG11 mutations were tolerated in C. glabrata, and these mutations could confer azole resistance. We found that the strongest azole-resistant allele of ERG11 led to induction of the Pdr1 transcription factor and Cdr1 ATP-binding cassette transporter protein in the absence of drug. ERG11 mutations can cause azole resistance via altered enzymatic properties but also by triggering induction of other resistance systems owing to impacts on ergosterol biosynthesis. These data illustrate the deep connections between ergosterol biosynthesis and regulation of membrane transporter proteins via Pdr1 and the ergosterol-responsive transcription factor Upc2A.