title: Whole blood metabolomics of dementia patients reveal classes of disease-linked metabolites

abstract:
Dementia is caused by factors that damage neurons. We quantified small molecular markers in whole blood of dementia patients, using non-targeted liquid chromatography-mass spectroscopy (LC-MS). Thirty-three metabolites, classified into 5 groups (A-E), differed significantly in dementia patients, compared with healthy elderly subjects. Seven Group A metabolites present in plasma, including quinolinic acid, kynurenine, and indoxyl-sulfate, increased. Possibly they act as neurotoxins in the central nervous system (CNS). The remaining 26 compounds (Groups B-E) decreased, possibly causing a loss of support or protection of the brain in dementia. Six Group B metabolites, normally enriched in red blood cells (RBCs) of healthy subjects, all contain trimethylated ammonium moieties. These metabolites include ergothioneine and structurally related compounds have scarcely been investigated as dementia markers, validating the examination of RBC metabolites. Ergothioneine, a potent anti-oxidant, is significantly decreased in various cognition-related disorders, such as mild cognitive impairment and frailty. Group C compounds, also include some oxidoreductants and are normally abundant in RBCs (NADP+, glutathione, ATP, pantothenate, S-adenosyl-methionine, and gluconate). Their decreased levels in dementia patients may also contribute to depressed brain function. Groups D (12) contains plasma compounds, such as amino acids, glycerophosphocholine, dodecanoyl-carnitine, 2-hydroxybutyrate, which normally protect the brain, but their diminution in dementia may reduce that protection. Seven Group D compounds have been identified previously as dementia markers. Group B-E compounds may be critical to maintain the CNS by acting directly or indirectly. How RBC metabolites act in the CNS and why they diminish so significantly in dementia remain to be determined.

Significance StatementDementia is a slowly progressing, chronic, and usually irreversible decline in cognitive function. Mechanistic causes and definitive treatments remain elusive. Using comprehensive metabolomics, we identified 5 groups of metabolites (A-E), 21 of which are novel, possibly useful for diagnosis and therapy of forms of dementia, such as Alzheimers disease. Seven Group A compounds may act as neurotoxins, whereas Group B-E compounds may protect the CNS against oxidative stress, maintain energy reserves, supply nutrients and neuroprotective factors. Five metabolites, ergothioneine, S-methyl-ergothioneine, trimethyl-histidine, methionine, and tryptophan identified in this study overlap with those reported for frailty. Interventions for cognitive diseases involving these dementia metabolomic markers may be accomplished either by inhibiting Group A compounds or by supplementing Group B-E compounds in patients.


title: Spatial transcriptomics using combinatorial fluorescence spectral and lifetime encoding, imaging and analysis

abstract: Multiplexed mRNA profiling in the spatial context provides important new information enabling basic research and clinical applications. Unfortunately, most existing spatial transcriptomics methods are limited due to either low multiplexing or assay complexity. Here, we introduce a new spatialomics technology, termed Multi Omic Single-scan Assay with Integrated Combinatorial Analysis (MOSAICA), that integrates in situ labeling of mRNA and protein markers in cells or tissues with combinatorial fluorescence spectral and lifetime encoded probes, spectral and time-resolved fluorescence imaging, and machine learning-based target decoding. This technology is the first application combining the biophotonic techniques, Spectral and Fluorescence Lifetime Imaging and Microscopy (FLIM), to the field of spatial transcriptomics. By integrating the time dimension with conventional spectrum-based measurements, MOSAICA enables direct and highly-multiplexed in situ spatial biomarker profiling in a single round of staining and imaging while providing error correction removal of background autofluorescence. We demonstrate mRNA encoding using combinatorial spectral and lifetime labeling and target decoding and quantification using a phasor-based image segmentation and machine learning clustering technique. We then showcase MOSAICAs multiplexing scalability in detecting 10-plex targets in fixed colorectal cancer cells using combinatorial labeling of only five fluorophores with facile error-correction and removal of autofluorescent moieties. MOSAICAs analysis is strongly correlated with sequencing data (Pearsons r = 0.9) and was further benchmarked using RNAscopeand LGC Stellaris. We further apply MOSAICA for multiplexed analysis of clinical melanoma Formalin-Fixed Paraffin-Embedded (FFPE) tissues that have a high degree of tissue scattering and intrinsic autofluorescence, demonstrating the robustness of the approach. We then demonstrate simultaneous co-detection of protein and mRNA in colorectal cancer cells. MOSAICA represents a simple, versatile, and scalable tool for targeted spatial transcriptomics analysis that can find broad utility in constructing human cell atlases, elucidating biological and disease processes in the spatial context, and serving as companion diagnostics for stratified patient care.

category: bioengineering


title: metaMIC: reference-free Misassembly Identification and Correction of de novo metagenomic assemblies

abstract: Evaluating the quality of metagenomic assemblies is important for constructing reliable metagenome-assembled genomes and downstream analyses. Here, we present metaMIC (https://github.com/ZhaoXM-Lab/metaMIC), a machine-learning based tool for identifying and correcting misassemblies in metagenomic assemblies. Benchmarking results on both simulated and real datasets demonstrate that metaMIC outperforms existing tools when identifying misassembled contigs. Furthermore, metaMIC is able to localize the misassembly breakpoints, and the correction of misassemblies by splitting at misassembly breakpoints can improve downstream scaffolding and binning results.

category: bioinformatics


title: SNAREs and Synaptotagmin cooperatively determine the Ca2+ sensitivity of neurotransmitter release in fixed stoichiometry modules

abstract: Neurotransmitter release is accomplished by a multi-component machinery including the membrane-fusing SNARE proteins and Ca2+-sensing Synaptotagmin molecules. However, the Ca2+ sensitivity of release was found to increase or decrease with more or fewer SNARE complexes at the release site, respectively, while the cooperativity is unaffected (Acuna et al., 2014; Arancillo et al., 2013), suggesting that there is no simple division of labor between these two components. To examine the mechanisms underlying these findings, we developed molecular dynamics simulations of the neurotransmitter release machinery, with variable numbers of Synaptotagmin molecules and assembled SNARE complexes at the release site. Ca2+ uncaging simulations showed that increasing the number of SNARE complexes at fixed stoichiometric ratio of Synaptotagmin to SNAREs increased the Ca2+ sensitivity without affecting the cooperativity. The physiological cooperativity of ~4-5 was reproduced with 2-3 Synaptotagmin molecules per SNARE complex, suggesting that Synaptotagmin and SNAREs cooperate in fixed stoichiometry modules. In simulations of action potential-evoked release, increased numbers of Synaptotagmin-SNARE modules increased release probability, consistent with experiment. Our simulations suggest that the final membrane fusion step is driven by SNARE complex-mediated entropic forces, and by vesicle-tethering forces mediated by the long Synaptotagmin linker domains. In consequence, release rates are increased when more SNARE complexes and Synaptotagmin monomers are present at the fusion site.

category: neuroscience


title: Vigilance and behavioral state-dependent modulation of cortical neuronal activity throughout the sleep/wake cycle

abstract: GABAergic inhibitory neurons, through their molecular, anatomic and physiological diversity, provide a substrate for the modulation of ongoing cortical circuit activity throughout the sleep-wake cycle. Here, we investigated neuronal activity dynamics of parvalbumin (PV), vasoactive intestinal polypeptide (VIP) and somatostatin (SST) neurons in naturally-sleeping head-restrained mice at the level of layer 2/3 of the primary somatosensory barrel cortex of mice. Through calcium-imaging and targeted single-unit loose-patch or whole-cell recordings, we found that PV action potential (AP) firing activity was largest during both NREM (non-rapid eye movement) and REM sleep stages, that VIP neurons were activated during REM sleep and that the overall activity of SST neurons remained stable throughout the sleep/wake cycle. Analysis of neuronal activity dynamics uncovered rapid decreases in PV cell firing at wake onset followed by a progressive recovery during wake. Simultaneous local field potential (LFP) recordings further revealed that, except for SST neurons, a large proportion of neurons were modulated by ongoing delta and theta waves. During NREM sleep spindles, PV and SST activity increased and decreased, respectively. Finally, we uncovered the presence of whisking behavior in mice during REM sleep and show that the activity of VIP and SST is differentially modulated during awake and sleeping whisking bouts, which may provide a neuronal substrate for internal brain representations occurring during sleep.


title: A double-hit in vivo model of GBA1 viral microRNA-mediated downregulation and human alpha-synuclein overexpression demonstrates nigrostriatal degeneration

abstract: Preclinical and clinical studies support a strong association between mutations in the GBA1 gene that encodes {beta}-glucocerebrosidase (GCase) (EC 3.2.1.45; glucosylceramidase beta) and Parkinsons disease (PD). Alpha-synuclein (AS), a key player in PD pathogenesis, and GBA1 mutations may independently and synergistically cause lysosomal dysfunction and thus, embody clinically well-validated targets of the neurodegenerative disease process in PD. However, double-hit in vivo models, recapitulating pathological features of PD that can be used to dissect the nature of the complex relationship between GCase and AS on the nigrostriatal axis, the region particularly vulnerable in PD, are direly needed. To address this, we implemented a bidirectional approach in mice to examine the effects of: 1) GCase overexpression (wildtype and mutant N370S) on endogenous AS levels and 2) downregulation of endogenous GCase combined with AS overexpression. Striatal delivery of viral-mediated GCase overexpression revealed minimal effects on cortical and nigrostriatal AS tissue levels and no significant effect on dopaminergic system integrity. On the other hand, microRNA (miR)-mediated GBA1 downregulation (miR GBA), combined with virus-mediated human AS overexpression (+AS), yields decreased GCase activity in the cortex, mimicking levels seen in GBA1 heterozygous carriers (30-40%), increased astrogliosis and microgliosis, decreased striatal dopamine levels (50% compared to controls) and loss of nigral dopaminergic neurons (~33%)-effects that were all reversible with miR rescue. Most importantly, the synergistic neurodegeneration of miR GBA+AS correlated with augmented AS accumulation and extracellular release in the striatum. Collectively, our results suggest that GCase downregulation alone is not sufficient to recapitulate key pathological features of PD in vivo, but its synergistic interplay with AS, via increased AS levels and release, drives nigrostriatal neurodegeneration. Furthermore, we demonstrate a novel double-hit model that can be used to identify putative mechanisms driving PD pathophysiology and can be subsequently used to test novel therapeutic approaches.


title: Profiling Myxococcus xanthus swarming phenotypes through mutation and environmental variation

abstract:
Myxococcus xanthus is a bacterium that lives on surfaces as a predatory biofilm called a swarm.As a growing swarm feeds on prey and expands, it displays dynamic multicellular patterns such as traveling waves called ripples and branching protrusions called flares. The rate at which a swarm expands across a surface, and the emergence of the coexisting patterns, are all controlled through coordinated cell movement. M. xanthus cells move using two motility systems known as Adventurous (A) and Social (S). Both are involved in swarm expansion and pattern formation. In this study, we describe a set of M. xanthus swarming genotype-to-phenotype associations that include both genetic and environmental perturbations. We identified new features of the swarming phenotype; recorded and measured swarm expansion using time-lapse microscopy; and compared the impact of mutation on different surfaces. These observations and analyses have increased our ability to discriminate between swarming phenotypes and provided context that allowed us to identify some phenotypes as improbable  outliers within the M. xanthus swarming phenome.

ImportanceMyxococcus xanthus grows on surfaces as a predatory biofilm called a swarm. A feeding swarm expands by moving over and consuming prey bacteria. In the laboratory, a swarm is created by spotting cell suspension onto nutrient agar in lieu of prey. The cells quickly settle on the surface and the new swarm then expands radially. An assay that measures the expansion rate of a swarm of mutant cells is the first, and sometimes only, measurement used to decide whether a particular mutation impacts swarm motility. We have broadened the scope of this assay by increasing the accuracy of measurements and reintroducing prey, resulting in new identifiable and quantifiable features that can be used to improve genotype-to-phenotype associations.

category: microbiology


title: Transfer rate of enveloped and non-enveloped viruses between fingerpads and surfaces

abstract:
Fomites can represent a reservoir for pathogens, which may be subsequently transferred from surfaces to skin. In this study we aim to understand how different factors (including virus type, surface type, time since last handwash, and direction of transfer) affect virus transfer rates, defined as the fraction of virus transferred, between fingerpads and fomites. To determine this, 360 transfer events were performed with 20 volunteers using Phi6 (a surrogate for enveloped viruses) and MS2 (a surrogate for non-enveloped viruses), and three clean surfaces (stainless steel, painted wood, and plastic). Considering all transfer events (all surfaces and both transfer directions combined), the mean transfer rates of Phi6 and MS2 were 0.17 and 0.26, respectively. Transfer of MS2 was significantly higher than Phi6 (P<0.05). Surface type was a significant factor that affected the transfer rate of Phi6: Phi6 is more easily transferred to and from stainless steel and plastic than to and from painted wood. Direction of transfer was a significant factor affecting MS2 transfer rates: MS2 is more easily transferred from surfaces to fingerpads than from fingerpads to surfaces. Data from these virus transfer events, and subsequent transfer rate distributions, provide information which can be used to refine quantitative microbial risk assessments. This study is the first to provide a large-scale data set of transfer events with a surrogate for enveloped viruses, which extends the reach of the study to the role of fomites in the transmission of human enveloped viruses like influenza and SARS-CoV-2.

ImportanceThis study created the first large-scale data set for the transfer of enveloped viruses between skin and surfaces. The data set produced by this study provides information on modelling the distribution of enveloped and non-enveloped virus transfer rates, which can aid in the implementation of risk assessment models in the future. Additionally, enveloped and non enveloped viruses were applied to experimental surfaces in an equivalent matrix to avoid matrix effects, so results between different viral species can be directly compared without confounding effects of different matrices. Our results indicating how virus type, surface type, time since last handwash, and direction of transfer affect virus transfer rates can be used in decision-making processes to lower the risk of viral infection from transmission through fomites.


title: Evidence for a bacterial Lands cycle phospholipase A: Structural and mechanistic insights into membrane phospholipid remodeling

abstract:
Cells steadily adapt their membrane glycerophospholipid (GPL) composition to changing environmental and developmental conditions. While the regulation of membrane homeostasis via GPL synthesis in bacteria has been studied in detail, the mechanisms underlying the controlled degradation of endogenous GPLs remain unknown. Thus far, the function of intracellular phospholipases A (PLAs) in GPL remodeling (Lands cycle) in bacteria is not clearly established. Here, we identified the first cytoplasmic membrane-bound phospholipase A1 (PlaF) from Pseudomonas aeruginosa involved in the Lands cycle. PlaF is an important virulence factor, as the P. aeruginosa {Delta}plaF mutant showed strongly attenuated virulence in Galleria mellonella and macrophages. We present a 2.0-[A]-resolution crystal structure of PlaF, the first structure that reveals homodimerization of a single-pass transmembrane (TM) full-length protein. PlaF dimerization, mediated solely through the intermolecular interactions of TM and juxtamembrane regions, inhibits its activity. A dimerization site and the catalytic sites are linked by an intricate ligand-mediated interaction network which likely explains the product (fatty acid) feedback inhibition observed with the purified PlaF protein. We used molecular dynamics simulations and configurational free energy computations to suggest a model of PlaF activation through a coupled monomerization and tilting of the monomer in the membrane, which constrains the active site cavity into contact with the GPL substrates. Thus, these data show the importance of the GPL remodeling pathway for virulence and pave the way for the development of a novel therapeutic class of antibiotics targeting PlaF-mediated membrane GPL remodeling.

SynopsisMembrane homeostasis can be regulated by phospholipase-controlled deacylation of endogenous glycerophospholipids (GPLs) followed by reacylation of products, known as the Lands cycle in eukaryotes. Here we show that the human pathogen Pseudomonas aeruginosa uses intracellular phospholipase A1 (PlaF) to modulate membrane GPL composition, which is the first example in bacteria. This newly identified PLA1 indirectly regulates the bacterial virulence properties by hydrolyzing a specific set of membrane GPLs. The crystal structure of full-length PlaF dimers bound to natural ligands, MD simulations, and biochemical approaches provide insights into the molecular mechanism of dimerization-mediated inactivation of this single-pass transmembrane PLA1. Our findings shed light on a mechanism by which bacterial intracellular PLAs might regulate membrane homeostasis what showcases these enzymes as a promising target for a new class of antibiotics.

category: molecular biology


title: Neuro-molecular characterization of fish cleaning interactions

abstract: Coral reef fish exhibit a large variety of behaviours crucial for fitness and survival. The cleaner wrasse Labroides dimidiatus displays cognitive abilities during interspecific interactions by providing services of ectoparasite cleaning, thus serving as a good example to understand the processes of complex social behaviour. However, little is known about the molecular underpinnings of cooperative behaviour between L. dimidiatus and a potential client fish (Acanthurus leucosternon). Therefore, we investigated the molecular mechanisms in three regions of the brain (fore-, mid-, and hindbrain) during the interaction of these fishes. Here we show, using transcriptomics, that most of the transcriptional response in both species was regulated in the hindbrain and forebrain regions and that the interacting behaviour responses of L. dimidiatus involved immediate early gene alteration, dopaminergic and glutamatergic pathways, the expression of neurohormones (such as isotocin) and steroids (e.g. progesterone and estrogen). In contrast, in the client, fewer molecular alterations were found, mostly involving pituitary hormone responses. The particular pathways found suggested learning and memory processes in the cleaner wrasse, while the client indicated stress relief.

category: animal behavior and cognition


title: Computational and experimental assessment of backbone templates for computational protein design

abstract: Computational protein design has taken big strides over the recent years, however, the tools available are still not at a state where a sequence can be designed to fold into a given protein structure at will and with high probability. We have here applied a recent release of Rosetta Design to redesign a set of structurally very similar proteins belonging to the Thioredoxin fold. We determined design success using a combination of a genetic screening tool to assay folding/stability in E. coli and selecting the best hits from this for further biochemical characterization. We have previously used this set of template proteins for redesign and found that success was highly dependent on template structure, a trait which was also found in this study. Nevertheless, state of the art design software is now able to predict the best template, most likely due to the introduction of the cart_bonded energy term. The template that led to the greatest fraction of successful designs was the same (a Thioredoxin from spinach) as that identified in our previous study. Our previously described redesign of Thioredoxin, which also used the spinach protein as template, however also performed well. In the present study, both these templates yielded proteins with compact folded structures, and enforces the conclusion that any design project must carefully consider different design templates. Fortunately, selecting designs using the cart_bonded energy term appears to correctly identify such templates.

category: biochemistry


title: CHARMM-GUI Membrane Builder for Lipid Nanoparticles with Ionizable Cationic Lipids and PEGylated Lipids

abstract: A lipid nanoparticle (LNP) formulation is a state-of-the-art delivery system for genetic drugs such as DNA, mRNA, and siRNA, which is successfully applied to COVID-19 vaccines and gains tremendous interest in therapeutic applications. Despite its importance, a molecular-level understanding of the LNP structures and dynamics is still lacking, which makes a rational LNP design almost impossible. In this work, we present an extension of CHARMM-GUI Membrane Builder to model and simulate all-atom LNPs with various (ionizable) cationic lipids and PEGylated lipids (PEG-lipids). These new lipid types can be mixed with any existing lipid types with or without a biomolecule of interest, and the generated systems can be simulated using various molecular dynamics engines. As a first illustration, we considered model LNP membranes with DLin-KC2-DMA (KC2) or DLin-MC3-DMA (MC3) without PEG-lipids. The results from these model membranes are consistent with those from the two previous studies albeit with mild accumulation of neutral MC3 in the bilayer center. To demonstrate Membrane Builders capability of building a realistic LNP patch, we generated KC2- or MC3-containing LNP membranes with high concentrations of cholesterol and ionizable cationic lipids together with 2 mol% PEG-lipids. We observe that PEG-chains are flexible, which can be more preferentially extended laterally in the presence of cationic lipids due to the attractive interactions between their head groups and PEG oxygen. The presence of PEG-lipids also relaxes the lateral packing in LNP membranes, and the area compressibility modulus (KA) of LNP membranes with cationic lipids fit into typical KA of fluid-phase membranes. Interestingly, the interactions between PEG oxygen and head group of ionizable cationic lipids induce a negative curvature. We hope that this LNP capability in Membrane Builder can be useful to better characterize various LNPs with or without genetic drugs for a rational LNP design.

category: biophysics


title: A multiscale approach for computing gated ligand binding from molecular dynamics and Brownian dynamics simulations

abstract:
We develop an approach to characterise the effects of gating by a multi-conformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterise them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multi-state gating theory from which we derive a gating factor{gamma} that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open and wide-open) and two inaccessible (closed and a newly identified,  parted) macrostate conformations. The  parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because{gamma} =0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease - that have a wide range of experimental association constants (~104 - 1010 M-1s-1) - yields gated rate constants in the range ~0.5-5.7 x 108 M-1s-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.

Graphical TOC Entry O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY


title: The Influence of Public Health Faculty on College and University Plans during the COVID-19 Pandemic

abstract: The purpose of this study was to determine whether the institutional presence of public health faculty within colleges and universities influenced operational plans for the fall semester of 2020. Using cross-sectional data collected by the College Crisis Initiative of Davidson College, six levels of instructional modalities (ranked from least to most restrictive) were compared between Council on Education of Public Health (CEPH)-accredited and non-CEPH-accredited 4-year institutions. Institutions with CEPH-accredited schools and programs were more likely to select some restrictive teaching modalities: 63.8% more likely to use hybrid/hyflex or more restrictive and 66.9% more likely to be primarily online (with some in person) or more restrictive. However, having CEPH-accredited programs did not push institutions to the most restrictive modalities. COVID-19 cases in county, enrollment, and political affiliation of the state governor were also found to influence instructional modality selection. While any ecological study has certain limitations, this study demonstrates that college and university fall plans appear to have been influenced by the presence of CEPH-accredited schools and programs of public health, and/or the input of their faculty. The influence of relevant faculty expertise on institutional decision-making can help inform college and university responses to future crises.

category: scientific communication and education


title: A Knowledge-guided Mechanistic Model of Synthetic Lethality in the HCT116 Vorinostat-resistant Colon Cancer Xenograft Model Cell-line

abstract: With an overall lifetime risk of about 4.3% and 4.0%, in men and women respectively, colorectal cancer remains the third leading cause of cancer-related deaths in the United States. In persons aged 55 and below, its rate increased at 1% per year in the years 2008 to 2017 despite the steady decline associated with improved screening, early diagnosis and treatment in the general population. Besides standardized therapeutic regimen, many trials continue to evaluate the potential benefits of vorinostat, mostly in combination with other anti-neoplastic agents for its treatment. Vorinostat, an FDA approved anti-cancer drug known as suberoylanilide hydroxamic acid (SAHA), an histone deacylase (HDAC) inhibitor, through many mechanisms, causes cancer cell arrest and death. However, like many other anti-neoplastic agents, resistance and or failures have been observed. In the HCT116 colon cancer cell line xenograft model, exploiting potential lethal molecular interactions by additional gene knockouts restored vorinotat sensitivity. This phenomenon, known as synthetic lethality, offers a promise to selectively target cancer cells. Although without clearly delineated understanding of underlying molecular processes, it has been demonstrated as an effective cancer-killing mechanism. In this study, we aimed to elucidate mechanistic interactions in multiple perturbations of identified synthetically lethal experiments, particularly in the vorinostat-resistant HCT116 (colon cancer xenograft model) cell line. Given that previous studies showed that knocking down GLI1, a downstream transcription factor involved in the Sonic Hedgehog pathway - an embryonal gene regulatory process, resulted in restoration of vorinostat sensitivity in the HCT116 colorectal cancer cell line, we hypothesized that vorinostat resistance is a result of upregulation of embryonal cellular differentiation processes; we hypothesized that elucidated regulatory mechanism would include crosstalks that regulate this biological process. We employed a knowledege-guided fuzzy logic regulatory inference method to elucidate mechanistic relationships. We validated inferred regulatory models in independent datasets. In addition, we evaluated the biomedical significance of key regulatory network genes in an independent clinically annotated dataset. We found no significant evidence that vorinostat resistance is due to an upregulation of embryonal gene regulatory pathways. Our observation rather support a topological rewiring of canonical oncogenic pathways around the PIK3CA, AKT1, RAS/BRAF etc. regulatory pathways. Reasoning that significant regulatory network genes are likely implicated in the clinical course of colorectal cancer, we show that the identified key regulatory network genes expression profile are able to predict short- to medium-term survival in colorectal cancer patients - providing a rationale basis for prognostification and potentially effective combination of therapeutics that target these genes along with vorinostat in the treatment of colorectal cancer.

category: systems biology


title: Molecular Features Underlying Shp1/Shp2 Discrimination by Immune Checkpoint Receptors

abstract: Numerous inhibitory immunoreceptors operate by recruiting phosphatase effectors Shp1 and Shp2 through conserved motifs ITIM and ITSM. Despite the similarity, these receptors exhibit distinct effector binding specificities, as exemplified by PD-1 and BTLA, which preferentially recruit Shp2 and Shp1 respectively. The molecular basis of Shp1/Shp2 discrimination is unclear. Here, we provide evidence that optimal PD-1 and BTLA binding to both Shp1 and Shp2 occurs via a bivalent, parallel mode that involves both SH2 domains of Shp1/Shp2. Moreover, PD-1 mainly uses its ITSM to discriminate Shp2 from Shp1 via their C-terminal SH2 domains. Supportive of this model, swapping the Shp1-cSH2 with Shp2-cSH2 enabled PD-1:Shp1 association in T cells. In contrast, BTLA primarily utilizes its ITIM to discriminate Shp1 from Shp2 via their N-terminal SH2 domains. Substitution of glycine at pY+1 position of the PD-1-ITIM with alanine, a residue conserved in several Shp1-recruiting receptors, was sufficient to induce PD-1:Shp1 interaction in T cells. Finally, mutagenesis screening shows that Shp1 recruitment exhibits a bell-shaped dependence on the side chain volume of the pY+1 residue of ITIM. Collectively, we provide a molecular interpretation of the Shp1/Shp2-binding specificities of PD-1 and BTLA, with general implications for the mechanism of effector discrimination by inhibitory receptors.

category: immunology


title: Fusobacterium nucleatum metabolically integrates commensals and pathogens in oral biofilms

abstract:
Fusobacterium nucleatum is a common constituent of the oral microbiota in both periodontal health and disease. Previously, we discovered ornithine cross-feeding between F. nucleatum and Streptococcus gordonii, where S. gordonii secretes ornithine via an arginine-ornithine antiporter (ArcD), which in turn supports the growth and biofilm development of F. nucleatum; however, broader metabolic aspects of F. nucleatum within polymicrobial communities and their impact on periodontal pathogenesis have not been addressed. Here, we show that when co-cultured with S. gordonii, F. nucleatum increased amino acid availability to enhance the production of butyrate and putrescine, a polyamine produced by ornithine decarboxylation. Co-culture with Veillonella parvula, another common inhabitant of the oral microbiota, also increased lysine availability, promoting cadaverine production by F. nucleatum. We confirmed that ArcD-dependent ornithine excretion by S. gordonii results in synergistic putrescine production, and mass spectrometry imaging revealed this metabolic capability creates a putrescine-rich microenvironment inside F. nucleatum biofilms. We further demonstrated that polyamines caused significant changes in the biofilm phenotype of a periodontal pathogen, Porphyromonas gingivalis, with putrescine being a potent stimulator of biofilm development and dispersal, and confirmed that F. nucleatum-mediated conversion of ornithine to putrescine enhances biofilm formation by P. gingivalis. Lastly, analysis of plaque samples revealed cooccurrence of P. gingivalis with genetic modules for putrescine production by S. gordonii and F. nucleatum. Overall, our results highlight the ability of F. nucleatum to induce synergistic polyamine production within multi-species consortia, and provide insight into how the trophic web in oral biofilm ecosystems can eventually shape disease-associated communities.

Significance StatementPeriodontitis is caused by the pathogenic transition of subgingival microbiota ecosystems, which is accompanied by alterations to microbiome functions including metabolic systems and the establishment of metabolic cross-feeding. While Fusobacterium nucleatum is a major constituent of the periodontal microbiota, its metabolic integration within polymicrobial communities and the impact on periodontal pathogenesis are poorly understood. Here, we report that amino acids supplied by other commensal bacteria induce polyamine production by F. nucleatum, creating polyamine-rich microenvironments. We further show that this trophic web results in enhancement of biofilm formation and dispersal of a periodontal pathogen, Porphyromonas gingivalis. This work provides mechanistic insight into how cooperative metabolism within oral biofilms can tip the balance toward periodontitis.


title: Bifurcation drives the evolution of assembly-line biosynthesis

abstract: Reprogramming biosynthetic assembly-lines is a topic of intense interest. This is unsurprising as the scaffolds of most antibiotics in current clinical use are produced by such pathways. The modular nature of assembly-lines provides a direct relationship between the sequence of enzymatic domains and the chemical structure of the product, but rational reprogramming efforts have been met with limited success. To gain greater insight into the design process, we wanted to examine how Nature creates assembly-lines and searched for biosynthetic pathways that might represent evolutionary transitions. By examining the biosynthesis of the anti-tubercular wollamides, we show how whole gene duplication and neofunctionalization can result in pathway bifurcation. Importantly, we show that neofunctionalization occurs primarily through intragenomic recombination. This pathway bifurcation leads to redundancy, providing the genetic robustness required to enable large structural changes during the evolution of antibiotic structures. Should the new product be none-functional, gene loss can restore the original genotype. However, if the new product confers an advantage, depreciation and eventual loss of the original gene creates a new linear pathway. This provides the blind watchmaker equivalent to the  design, build, test cycle of synthetic biology.


title: Impact of dose, duration and host immune status on ultrashort telacebec treatment in a mouse model of Buruli ulcer

abstract: Telacebec (Q203) is a new anti-tuberculosis drug in clinical development with extremely potent activity against Mycobacterium ulcerans, the causative agent of Buruli ulcer (BU). The potency of Q203 has prompted investigation of its potential role in ultra-short, even single-dose, treatment regimens for BU in mouse models. However, the relationships of Q203 dose and duration and host immune status to treatment outcomes remain unclear, as does the risk of emergence of drug resistance with Q203 monotherapy. In the present study, immunocompetent BALB/c and immunocompromised SCID-beige mice were infected in both hind footpads and treated eight weeks later. In both mouse strains, controls received rifampin-clarithromycin; others received Q203 at 0.5 or 2 mg/kg/d for 5 or 10 days. Additionally, BALB/c mice received a single dose of 2.5 or 10 mg/kg or 3.3 mg/kg/d for 3 days. Treatment response was based on changes in footpad swelling and CFU counts at the end of treatment as well as 4 and 13 weeks after stopping treatment. Efficacy depended on total dose more than duration. Total doses of 5-20 mg/kg rendered nearly all BALB/c mice culture-negative 13 weeks post-treatment without selection of Q203-resistant bacteria. Although less potent in SCID-beige mice, Q203 still rendered the majority of footpads culture-negative at total doses of 10-20 mg/kg. However, Q203 resistance was identified in relapse isolates from some SCID-beige mice. Overall, these results support the potential of Q203 monotherapy for single-dose or other ultra-short therapy for BU, although highly immunocompromised hosts may require higher doses or durations and/or combination therapy.


title: Parameter identifiability and model selection for sigmoid population growth models

abstract: Sigmoid growth models, such as the logistic, Gompertz and Richards models, are widely used to study population dynamics ranging from microscopic populations of cancer cells, to continental-scale human populations. Fundamental questions about model selection and parameter estimation are critical if these models are to be used to make practical inferences. However, the question of parameter identifiability - whether a data set contains sufficient information to give unique or sufficiently precise parameter estimates - is often overlooked. We use a profile-likelihood approach to explore practical parameter identifiability using data describing the re-growth of hard coral. With this approach, we explore the relationship between parameter identifiability and model misspecification, finding that the logistic growth model does not suffer identifiability issues for the type of data we consider whereas the Gompertz and Richards models encounter practical non-identifiability issues. This analysis of parameter identifiability and model selection is important because different growth models are used within areas of the biological modelling literature without necessarily considering whether parameters are identifiable, or checking statistical assumptions underlying model adequacy. Standard practices that do not consider parameter identifiability can lead to unreliable or imprecise parameter estimates and potentially misleading mechanistic interpretations. While tools developed here focus on three standard sigmoid growth models only, our theoretical developments are applicable to any sigmoid growth model and any relevant data set. MATLAB implementations of all software are available on GitHub.

category: ecology


title: Deep learning assessment of cultural ecosystem services from social media images

abstract: Crowdsourced social media data has become popular in the assessment of cultural ecosystem services (CES). Advances in deep learning show great potential for the timely assessment of CES at large scales. Here, we describe a procedure for automating the assessment of image elements pertaining to CES from social media. We focus on a binary (natural, human) and a multiclass (posing, species, nature, landscape, human activities, human structures) classification of those elements using two Convolutional Neural Networks (CNNs; VGG16 and ResNet152) with the weights from two large datasets - Places365 and ImageNet -, and our own dataset. We train those CNNs over Flickr and Wikiloc images from the Peneda-Geres region (Portugal) and evaluate their transferability to wider areas, using Sierra Nevada (Spain) as test. CNNs trained for Peneda-Geres performed well, with results for the binary classification (F1-score > 80%) exceeding those for the multiclass classification (> 60%). CNNs pre-trained with Places365 and ImageNet data performed significantly better than with our data. Model performance decreased when transferred to Sierra Nevada, but their performances were satisfactory (> 60%). The combination of manual annotations, freely available CNNs and pre-trained local datasets thereby show great relevance to support automated CES assessments from social media.


title: Bacterial response to spatial gradients of algal-derived nutrients in a porous microplate

abstract: AO_SCPLOWBSTRACTC_SCPLOWPhotosynthetic microalgae are responsible for 50% of the global atmospheric CO2 fixation into organic matter and hold potential as a renewable bioenergy source. Their metabolic interactions with the surrounding microbial community (the algal microbiome) play critical roles in carbon cycling, but due to methodological limitations, it has been challenging to examine how community is developed by spatial proximity to their algal host. Here we introduce a hydrogel-based porous microplate to co-culture algae and bacteria, where metabolites are constantly exchanged between the microorganisms while maintaining physical separation. In the microplate we found that the diatom Phaeodactylum tricornutum accumulated to cell abundances ~20 folds higher than under normal batch conditions due to constant replenishment of nutrients through the hydrogel. We also demonstrate that algal-associated bacteria, both single isolates and complex communities, responded to inorganic nutrients away from their host as well as organic nutrients originating from the algae in a spatially predictable manner. These experimental findings coupled with a mathematical model suggest that host proximity and algal culture growth phase impact bacterial community development in a taxon-specific manner through organic and inorganic nutrient availability. Our novel system presents a useful tool to investigate universal metabolic interactions between microbes in aquatic ecosystems.


title: De-novo draft genome sequence of Crocus Sativus L, Saffron, a golden condiment.

abstract: Crocus sativus L, saffron is the highest priced but low yielding plant of medicinal and culinary importance. Despite its economic status, the omic information on this plant is very scarce, with only a couple of transcriptomics and epigenetic studies. In the present study, the draft genome sequence of C. sativus has been assembled using Illumina sequencing. In fact, this is the first genome sequence from any member of family Iridaceae. Genome size of C. sativus was estimated to be 3.5 Gb and the draft genome is 3.01 Gb long with 84.24% genome coverage. In total, 8,62,275 repeats and 9,64,231 SSR markers have been identified. A total of 53,546 functional genes were annotated, out of which, 43,649 proteins were associated with GO annotation. 5726 proteins were identified as transcription factors, with MYB & MYB related family proteins being more abundant. Orthology analysis of C. sativus with 3 different monocot species of the same plant order and rice (model monocot plant) revealed 7328 proteins clusters to be conserved in all the five plant species, whereas 2510 proteins cluster were unique to C. sativus only. 10,912 unigenes of C. sativus were mapped to 387 KEGG pathways of monocot. The genes involved in the pathway of apocarotenoids biosynthesis (crocin, crocetin, picrocrocin, and safranal) were present in the draft genome.

category: genomics


title: Regulation of GTPase function by autophosphorylation

abstract: A unifying feature of the RAS superfamily is a functionally conserved GTPase cycle that proteins use to transition between active and inactive states. Here, we demonstrate that active site autophosphorylation of some small GTPases is an intrinsic regulatory mechanism that reduces nucleotide hydrolysis and enhances nucleotide exchange, thus altering the on/off switch that forms the basis for their signaling functions. Using x-ray crystallography, nuclear magnetic resonance spectroscopy, biolayer interferometry binding assays, and molecular dynamics on autophosphorylated mutants of H-RAS and K-RAS, we show that phosphoryl transfer from GTP requires dynamic movement of the switch II domain and that autophosphorylation promotes nucleotide exchange by opening of the active site and extraction of the stabilizing Mg. Finally, we demonstrate that autophosphorylated K-RAS exhibits altered effector interactions, including a reduced affinity for RAF proteins. Thus, autophosphorylation leads to altered active site dynamics and effector interaction properties, creating a pool of GTPases that are functionally distinct from the non-phosphorylated counterpart.


title: Cross-platform Bayesian optimization system for autonomous biological assay development

abstract: Current high-throughput screening assay optimization is often a manual and time-consuming process, even when utilizing design-of-experiment approaches. A cross-platform, Cloud-based Bayesian optimization-based algorithm was developed as part of the NCATS ASPIRE Initiative to accelerate preclinical drug discovery. A cell-free assay for papain enzymatic activity was used as proof-of-concept for biological assay development. Compared to a brute force approach that sequentially tested all 294 assay conditions to find the global optimum, the Bayesian optimization algorithm could find suitable conditions for optimal assay performance by testing only 21 assay conditions on average, with up to 20 conditions being tested simultaneously. The algorithm could achieve a seven-fold reduction in costs for lab supplies and high-throughput experimentation run-time, all while being controlled from a remote site through a secure connection. Based on this proof-of-concept, this technology is expected to be applied to more complex biological assays and automated chemistry reaction screening at NCATS, and should be transferable to other institutions.


title: INTRINSIC MOTOR NEURONE EXCITABILITY IS REDUCED IN SOLEUS AND TIBIALIS ANTERIOR OF OLDER ADULTS

abstract: Age-related deterioration within both motor neurones and monoaminergic systems should theoretically reduce neuromodulation by weakening motor neuronal persistent inward current (PIC) amplitude. However, this assumption remains untested. Surface electromyographic signals were collected using two 32-channel electrode matrices placed on soleus and tibialis anterior of 25 older adults (70{+/-}4years) and 17 young adults (29{+/-}5 years) to investigate motor unit discharge behaviours. Participants performed triangular-shaped plantar and dorsiflexion contractions to 20% of maximum torque at a rise-decline rate of 2%/s of each participants maximal torque. Pairwise and composite paired-motor unit analyses were adopted to calculate delta frequency ({Delta}F), which has been used to differentiate between the effects of synaptic excitation and intrinsic motor neuronal properties and is assumed to be proportional to PIC amplitude. Soleus and tibialis anterior motor units in older adults had lower {Delta}Fs calculated with either the pairwise [-0.99 and -1.46 pps; -35.4 and - 33.5%, respectively] or composite (-1.18 and -2.28 pps; -32.1 and -45.2%, respectively) methods. Their motor units also had lower peak discharge rates (-2.14 and -2.03 pps; -19.7 and -13.9%, respectively) and recruitment thresholds (-1.50 and -2.06% of maximum, respectively) than young adults. These results demonstrate reduced intrinsic motor neurone excitability during low-force contractions in older adults, likely mediated by decreases in the amplitude of persistent inward currents. Our findings might be explained by deterioration in the motor neurones or monoaminergic systems and could contribute to the decline in motor function during ageing; these assumptions should be explicitly tested in future investigations.

category: physiology


title: Spatial organization of face part representations within face-selective areas revealed by high-field fMRI

abstract: Regions sensitive to specific object categories as well as organized spatial patterns sensitive to different features have been found across the whole ventral temporal cortex (VTC). However, it is unclear that within each object category region, how specific feature representations are organized to support object identification. Would object features, such as object parts, be represented in fine-scale spatial organization within object category-specific regions? Here we used high-field 7T fMRI to examine the spatial organization of neural tuning to different face parts within each face-selective region. Our results show consistent spatial organization across individuals that within right posterior fusiform face area (pFFA) and right occipital face area (OFA), the posterior portion of each region was biased to eyes, while the anterior portion was biased to mouth and chin stimuli. Our results demonstrate that within the occipital and fusiform face processing regions, there exist systematic spatial organizations of neural tuning to different face parts that support further computation combining them.


title: Multi-dynamic Modelling Reveals Strongly Time-varying Resting fMRI Correlations

abstract:
The activity of functional brain networks is responsible for the emergence of time-varying cognition and behaviour. Accordingly, time-varying correlations (Functional Connectivity) in resting fMRI have been shown to be predictive of behavioural traits, and psychiatric and neurological conditions. Typically, methods that measure time varying Functional Connectivity (FC), such as sliding windows approaches, do not separately model when changes occur in the mean activity levels from when changes occur in the FC, therefore conflating these two distinct types of modulation. We show that this can bias the estimation of time-varying FC to appear more stable over time than it actually is. Here, we propose an alternative approach that models changes in the mean brain activity and in the FC as being able to occur at different times to each other. We refer to this method as the Multi-dynamic Adversarial Generator Encoder (MAGE) model, which includes a model of the network dynamics that captures long-range time dependencies, and is estimated on fMRI data using principles of Generative Adversarial Networks. We evaluated the approach across several simulation studies and resting fMRI data from the Human Connectome Project (1003 subjects), as well as from UK Biobank (13301 subjects). Importantly, we find that separating fluctuations in the mean activity levels from those in the FC reveals much stronger changes in FC over time, and is a better predictor of individual behavioural variability

Statement of SignificanceO_LIMAGE is multi-dynamic in that it models temporal fluctuations in FC independently from fluctuations in the mean of the activity.
C_LIO_LIMAGE reveals stronger changes in FC over time than single-dynamic approaches, such as sliding window correlations.
C_LIO_LIMulti-dynamic modelling provides an explanation and a solution as to why resting fMRI FC has previously looked so stable.
C_LIO_LIMAGE models fMRI data as a set of reoccurring brain states, and importantly, these states do not have to be binary and mutually exclusive (e.g., multiple states can be active at one time-point).
C_LIO_LIMAGE estimated time-varying FC is a better predictor of behavioural variability in the resting-state fMRI data than established methods.
C_LI


title: Phylodynamic assessment of control measures for highly pathogenic avian influenza epidemics in France

abstract: Phylodynamic methods have successfully been used to describe viral spread history but their applications for assessing specific control measures are rare. In 2016-17, France experienced a devastating epidemic of a highly pathogenic avian influenza virus (H5N8 clade 2.3.4.4b). Using 196 viral genomes, we conducted a phylodynamic analysis combined with generalised linear model and showed that the large-scale preventive culling of ducks significantly reduced the viral spread between departements (French administrative division). We also found that the virus likely spread more frequently between departements that shared borders, but the spread was not linked to duck transport between departements. Duck transport within departements increased the within-departement transmission intensity, although the association was weak. Together, these results indicated that the virus spread in short-distances, either between adjacent departements or within departements. Results also suggested that the restrictions on duck transport within departements might not have stopped the viral spread completely. Overall, by testing specific hypothesis related to different control measures, we demonstrated that phylodynamics methods are capable of investigating the impacts of control measures on viral spread.

category: evolutionary biology


title: A selective adenylyl cyclase 1 inhibitor relieves pain without causing tolerance

abstract: Adenylyl cyclases (ACs) catalyze the production of the second messenger cyclic adenosine monophosphate from adenosine triphosphate. Among the ten different AC isoforms, studies with knockout animals indicate that inhibition of AC1 can relieve pain and reduce behaviors linked to opioid dependence. We previously identified ST034307 as a selective inhibitor of AC1. The development of an AC1-selective inhibitor now provides the opportunity to further study the therapeutic potential of inhibiting this protein in pre-clinical animal models of pain and related adverse reactions. In the present study we have shown that ST034307 relieves pain in mouse models of formalin-induced inflammatory pain, acid-induced visceral pain, and acid-depressed nesting. In addition, ST034307 did not cause analgesic tolerance after chronic dosing. We also show that the compound is restricted to the periphery following subcutaneous injections and report the predicted molecular interaction between ST034307 and AC1. Our results indicate that AC1 inhibitors represent a promising new class of analgesic agents that treat pain and appear to produce less adverse effects than currently-used opioids.

category: pharmacology and toxicology


title: Generalising state transition structure in the human medial temporal lobe

abstract: Effective planning involves knowing where different actions will take us. However natural environments are rich and complex, leading to an exponential increase in memory demand as a plan grows in depth. One potential solution to this problem is to share the neural state transition functions used for planning between similar contexts. Here, we asked human participants to perform a sequential decision making task designed so that knowledge could be shared between some contexts but not others. Computational modelling showed that participants shared a model of state transitions between contexts where appropriate. fMRI data identified the medial temporal lobe as a locus for learning of state transitions, and within the same region, correlated BOLD patterns were observed in contexts where state transition information were shared. Finally, we show that the transition model is updated more strongly following the receipt of positive compared to negative outcomes, a finding that challenges conventional theories of planning which assume knowledge about our environment is updated independently of outcomes received. Together, these findings propose a computational and neural account of how information relevant for planning can be shared between contexts.


title: Model sharing in the human medial temporal lobe


title: SARS-CoV-2 neutralising antibodies in Dogs and Cats in the United Kingdom

abstract: 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.


title: Escherichia coli recombinant expression of SARS-CoV-2 protein fragments.

abstract: We have developed a method for the inexpensive, high-level expression of antigenic protein fragments of SARS-CoV-2 proteins in Escherichia coli. Our approach used the thermophilic family 9 carbohydrate-binding module (CBM9) as an N-terminal carrier protein and affinity tag. The CBM9 module was joined to SARS-CoV-2 protein fragments via a flexible proline-threonine linker, which proved to be resistant to E. coli proteases. Two CBM9-spike protein fragment fusion proteins and one CBM9-nucleocapsid fragment fusion protein largely resisted protease degradation, while most of the CBM9 fusion proteins were degraded at some site in the SARS-CoV-2 protein fragment. All fusion proteins were expressed in E. coli at about 0.1 g/L, and could be purified with a single affinity binding step using inexpensive cellulose powder. Three purified CBM9-SARS-CoV-2 fusion proteins were tested and found to bind antibody directed to the appropriate SARS-CoV-2 antigenic region. The largest intact CBM9 fusion protein incorporates spike protein amino acids 540-588, which is a conserved region immediately C-terminal to the receptor binding domain that is widely recognized by human convalescent sera and contains a putative protective epitope.


title: Increased lung cell entry of B.1.617.2 and evasion of antibodies induced by infection and BNT162b2 vaccination

abstract: The delta variant of SARS-CoV-2, B.1.617.2, emerged in India and has subsequently spread to over 80 countries. B.1.617.2 rapidly replaced B.1.1.7 as the dominant virus in the United Kingdom, resulting in a steep increase in new infections, and a similar development is expected for other countries. Effective countermeasures require information on susceptibility of B.1.617.2 to control by antibodies elicited by vaccines and used for COVID-19 therapy. We show, using pseudotyping, that B.1.617.2 evades control by antibodies induced upon infection and BNT162b2 vaccination, although with lower efficiency as compared to B.1.351. Further, we found that B.1.617.2 is resistant against Bamlanivimab, a monoclonal antibody with emergency use authorization for COVID-19 therapy. Finally, we show increased Calu-3-lung cell entry and enhanced cell-to-cell fusion of B.1.617.2, which may contribute to augmented transmissibility and pathogenicity of this variant. These results identify B.1.617.2 as an immune evasion variant with increased capacity to enter and fuse lung cells.


title: Local temporal Rac1-GTP nadirs and peaks restrict cell protrusions and retractions

abstract: 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.

category: cell biology


title: Unsupervised classification of SARS-CoV-2 genomic sequences uncovers hidden genetic diversity and suggests an efficient strategy for genomic surveillance

abstract:
Accurate and timely monitoring of emerging genomic diversity is crucial for limiting the spread of potentially more transmissible/pathogenic strains of SARS-CoV-2. At the time of writing, over 1.8M distinct viral genome sequences have been made publicly available, and a sophisticated nomenclature system based on phylogenetic evidence and expert manual curation has allowed the relatively rapid classification of emerging lineages of potential concern.

Here, we propose a complementary approach that integrates fine-grained spatiotemporal estimates of allele frequency with unsupervised clustering of viral haplotypes, and demonstrate that multiple highly frequent genetic variants, arising within large and/or rapidly expanding SARS-CoV-2 lineages, have highly biased geographic distributions and are not adequately captured by current SARS-CoV-2 nomenclature standards.

Our results advocate a partial revision of current methods used to track SARS-CoV-2 genomic diversity and highlight the importance of the application of strategies based on the systematic analysis and integration of regional data. Here we provide a complementary, completely automated and reproducible framework for the mapping of genetic diversity in time and across different geographic regions, and for the prioritization of virus variants of potential concern. We believe that the approach outlined in this study will contribute to relevant advances to current genomic surveillance methods.


title: Aliquoting of isobaric labeling reagents for low concentration and single cell proteomics samples

abstract:
The introduction of isobaric tagging reagents enabled more accurate, high-throughput quantitative proteomics by enabling multiple samples to be multiplexed. One drawback of these workflows is the relative expense of the proprietary isobaric reagents, which is often only second to the expense of the instruments themselves. These highly reactive chemical tags are only commercially available in relatively large aliquots compared to the typical amounts of peptides analyzed in proteomic workflows today. Excess reagents are typically disposed of following a single labeling experiment or those performed within a few days of opening a new kit. We present a simple procedure to aliquot commercial isobaric tagging reagents and demonstrate the successful and high efficiency labeling of multiple samples over a period of six months. The samples presented herein were selected as the most diverse ones labeled by prepared aliquots from a single labeling reagent kit over this period. We observe comparable labeling efficiency from 100 microgram to 100 picograms of peptide when labeling samples from both human digest standards, cancer cell lines prepared in-house and from cells directly obtained from human organ donors, despite differences in cell type, lysis, and digestion procedures. No labeling experiment of whole human proteomics samples achieved less than 92% labeling efficiency over this period. When preparing phosphoproteomic samples 6 months after the date of the aliquoting procedure, we observed a decrease in labeling efficiency to approximately 86%, indicating the end of the useful lifetime of aliquots prepared in this manner. Over this period, we have effectively reduced the reagent costs of each experiment to less than 10% of the predicted costs when following the manufacturer instructions for use and disposal. While aliquoting of reagents can be performed by hand, we provide a complete template for automatic aliquoting using an affordable liquid handling robot, including plans for 3D printing of two parts we have found useful for streamlining this procedure.

Abstract Graphic

O_FIG O_LINKSMALLFIG WIDTH=197 HEIGHT=200 SRC="FIGDIR/small/449560v1_ufig1.gif" ALT="Figure 1">
View larger version (51K):
org.highwire.dtl.DTLVardef@18cd7eeorg.highwire.dtl.DTLVardef@1b48065org.highwire.dtl.DTLVardef@159761forg.highwire.dtl.DTLVardef@5a0896_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Genetically encoded green fluorescent biosensors for monitoring UDP-GlcNAc in live cells

abstract: Uridine diphosphate N-acetylglucosamine (UDP-GlcNAc) is a nucleotide sugar used by glycosyltransferases to synthesize glycoproteins, glycosaminoglycans, glycolipids, and glycoRNA. UDP-GlcNAc also serves as the donor substrate for the formation of O-GlcNAc, a dynamic intracellular protein modification involved in diverse signaling and disease processes. UDP-GlcNAc is thus a central metabolite connecting nutrition, metabolism, signaling, and disease. There is a great interest in monitoring UDP-GlcNAc in biological systems. Here, we present the first genetically encoded, green fluorescent UDP-GlcNAc sensor (UGAcS), an optimized insertion of a circularly permuted green fluorescent protein (cpGFP) into an inactive mutant of an E. coli UDP-GlcNAc transferase, for ratiometric monitoring of UDP-GlcNAc dynamics in live mammalian cells. Although UGAcS responds to UDP-GlcNAc quite selectively among various nucleotide sugars, UDP and UTP interfere with the response. We thus developed another biosensor named UXPS, which is responsive to UDP and UTP but not UDP-GlcNAc. We demonstrated the use of the biosensors to follow UDP-GlcNAc levels in cultured mammalian cells perturbed with nutritional changes, pharmacological inhibition, and knockdown or overexpression of key enzymes in the UDP-GlcNAc synthesis pathway. We further utilized the biosensors to monitor UDP-GlcNAc concentrations in pancreatic MIN6 {beta}-cells under various culture conditions.


title: Structural insight into molecular inhibitory mechanism of InsP6 on African Swine Fever Virus mRNA-decapping enzyme g5Rp

abstract: Removal of 5' cap on cellular mRNAs by the African Swine Fever Virus (ASFV) decapping enzyme g5R protein (g5Rp) is beneficial to viral gene expression during the early stages of infection. As the only nucleoside diphosphate linked moiety X (Nudix) decapping enzyme encoded in the ASFV genome, g5Rp works in both the degradation of cellular mRNA and hydrolyzation of the diphosphoinositol polyphosphates. Here, we report the structures of dimeric g5Rp and its complex with inositol hexakisphosphate (InsP6). The two g5Rp protomers interact head-to-head to form a dimer, and the dimeric interface is formed by extensive polar and nonpolar interactions. Each protomer composed a unique N-terminal helical domain and C-terminal classic Nudix domain. As a mRNA decapping enzyme, we identified key residues, including K8, K94, K95, K98, K175, R221, and K243 located on the substrate RNA binding interfaces of g5Rp, are important to RNA binding and decapping enzyme activity. Furthermore, we identified that the g5Rp-mediated mRNA decapping was inhibited by the InsP6. The g5Rp-InsP6 complex structure showed that the InsP6 molecules occupy the same regions that primarily mediate g5Rp-RNA interaction, elucidating the roles of InsP6 in the regulation of the viral decapping activity of g5Rp in mRNA degradation. Collectively, these results provide the structural basis of interaction between RNA and g5Rp and highlight the inhibitory mechanism of InsP6 on mRNA decapping by g5Rp.


title: Comparative Structural Analyses of Selected Spike Protein-RBD Mutations in SARS-CoV-2 Lineages

abstract: The severity of the covid 19 has been observed throughout the world as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had spread globally claiming more than 2 million lives and left a devastating impact on peoples life. Recently several virulent mutant strains of this virus, such as the B.1.1.7, B.1.351, and P1 lineages have emerged. These strains are predominantly observed in UK, South Africa and Brazil. Another extremely pathogenic B.1.617 lineage and its sub-lineages, first detected in India, are now affecting some countries at notably stronger spread-rates. This paper computationally examines the time-based structures of B.1.1.7, B.1.351, P1 lineages with selected spike protein mutations. Additionally, the mutations in the more recently found B.1.617 lineage and some of its sub-lineages are explored, and the implications for multiple point mutations of the spike proteins receptor-binding domain (RBD) are described.


title: Designing BH3-mimetic Peptide Inhibitors for the Viral Bcl-2 Homologs A179L and BHRF1: Importance of long-range electrostatic interactions

abstract: 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.


title: Accurate model of liquid-liquid phase behaviour of intrinsically-disordered proteins from data-driven optimization of single-chain properties

abstract: Many intrinsically disordered proteins (IDPs) may undergo liquidliquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalisation of intracellular biochemical reactions. The phase behaviour of IDPs is sequence-dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intra- and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.


title: Dissemination and progression of pulmonary Mycobacterium avium infection in mouse model are associated with type 2 macrophage activation

abstract: Pulmonary infections caused by the group of nontuberculosis mycobacteria (NTM), Mycobacterium avium complex (MAC), are increasing worldwide and a growing public health concern. Pulmonary granulomas are the hallmark of MAC lung infection, yet reliable correlates of granuloma progression and susceptibility in immunocompetent hosts are poorly defined. The development of mouse models that recapitulate the diversity of granulomas seen in MAC pulmonary disease in humans is crucial to study mechanisms of susceptibility in humans and for preclinical evaluation of therapeutics. Unlike widely used inbred mouse strains, mice that carry the mutant allele at the genetic locus sst1 develop human-like pulmonary tuberculosis featuring well-organized caseating granulomas. These mice became instrumental in pre-clinical testing of novel interventions. In this study we tested whether the B6.Sst1S that carries the sst1 mutant allele on standard B6 background develop more advanced pulmonary infection with NTM M. avium spp. hominissuis (M.av). To assess pulmonary disease progression, we utilized traditional semi-quantitative histomorphological evaluation and fluorescent multiplex immunohistochemistry (fmIHC) in combination with whole slide imaging and digital image analysis. After infection with the laboratory M.av strain 101, the B6.Sst1S pulmonary lesions progressed 12 - 20 weeks post infection, although we did not observe the formation of necrotic granulomas during this interval. Using fmIHC, we determined that the disease progression was associated with a steadily increasing proportion of mycobacteria infected Arg1+ and double positive iNOS+/Arg1+ macrophages. The B6.Sst1S granulomas had a greater proportion of Arg1+ and double positive iNOS+/Arg1+ macrophages, and decreased T cell density, as compared to wild type B6 mice. Thus, the genetic composition of the B6.Sst1S mice renders them more susceptible to pulmonary M.av infection. In combination with more virulent clinical isolates of M.av these mice could provide an improved mouse model that recapitulates more severe pulmonary disease in humans. The Arg1 macrophage expression in this model combined with automated fmIHC could serve as a sensitive biomarker for the unbiased assessment of medical countermeasures against NTM infection.

category: pathology


title: Sustained TREM2 stabilization accelerates microglia heterogeneity and Abeta pathology in a mouse model of Alzheimer s disease

abstract: TREM2 is a transmembrane protein expressed exclusively in microglia in the brain that regulates inflammatory responses to pathological conditions. Proteolytic cleavage of membrane TREM2 affects microglial function and is associated with Alzheimers disease, but the consequence of reduced TREM2 proteolytic cleavage has not been determined. We generated a transgenic mouse model of reduced TREM2 shedding (Trem2-IPD) through amino acid substitution of ADAM-protease recognition site. We found that Trem2-IPD mice displayed increased TREM2 cell surface receptor load, survival and function in myeloid cells. Using single cell transcriptomic profiling of mouse cortex we show that sustained TREM2 stabilization induces a shift of fate in microglial maturation and accelerates microglial responses to A{beta} pathology in a mouse model of Alzheimers disease. Our data indicate that reduction of TREM2 proteolytic cleavage aggravates neuroinflammation during the course of AD pathology suggesting that TREM2 shedding is a critical regulator of microglial activity in pathological states.


title: Protein crosslinking as a therapeutic strategy for SOD1-related ALS

abstract:
Mutations in the gene encoding Cu-Zn superoxide dismutase 1 (SOD1) cause a subset of familial amyotrophic lateral sclerosis (fALS). One effect of these mutations is that SOD1, which is normally a stable dimer, dissociates into toxic monomers. Considerable research efforts have been devoted to developing compounds that stabilize the dimer of fALS SOD1 variants, but these have not yet resulted in an approved drug. We demonstrate that a cyclic thiosulfinate cross-linker can stabilize prevalent disease-causing SOD1 variants. The degree of stabilization afforded by cyclic thiosulfinates (up to 24 {degrees}C) is unprecedented. We show this compound works rapidly in vivo with a half-life of ~3 days. The efficacy, low toxicity, and pharmacodynamics of cross-linker mediated stabilization make it a promising therapeutic approach for SOD1-related fALS.

Significance statementCyclic thiosulfinate S-XL6 enables the kinetic stabilization of ALS-associated SOD1 variants, in vivo.


title: Neuronal hyperactivity in a LRRK2-G2019S cellular model of Parkinson's Disease

abstract: Monogenic Parkinsons Disease can be caused by a mutation in the leucine-rich repeat kinase 2 (LRRK2) gene, causing a late-onset autosomal dominant inherited form of Parkinsons Disease. The function of the LRRK2 gene is incompletely understood, but several in vitro studies have reported that LRRK2-G2019S mutations affect neurite branching, calcium homeostasis and mitochondrial function, but thus far, there have been no reports of effects on electrophysiological activity. We assessed the neuronal activity of induced pluripotent stem cell derived neurons from Parkinsons Disease patients with LRRK2-G2019S mutations and isogenic controls. Neuronal activity of spontaneously firing neuronal populations was recorded with a fluorescent calcium-sensitive dye (Fluo-4) and analysed with a novel image analysis pipeline that combined semi-automated neuronal segmentation and quantification of calcium transient properties. Compared with controls, LRRK2-G2019S mutants have shortened inter-spike intervals and an increased rate of spontaneous calcium transient induction.


title: Ire1α-Regulated Rate of mRNA Translation is Required for Acquisition of Identity and Polarity in Upper Layer Cortical Neurons

abstract:
Evolutionary expansion of the neocortex is associated with the increase in upper layer neurons. Here, we present Inositol-Requiring Enzyme 1, Ire1, as an essential determinant of upper layer fate, neuronal polarization and cortical lamination. We demonstrate a non-canonical function of Ire1 in the regulation of global translation rates in the developing neocortex through its dynamic interaction with the ribosome and regulation of eIF4A1 and eEF-2 expression. Inactivation of Ire1 engenders lower protein synthesis rates associated with stalled ribosomes and decreased number of translation start sites. We show unique sensitivity of upper layer fate to translation rates. Whereas eEF-2 is required for cortical lamination, eIF4A1 regulates acquisition of upper layer fate downstream of Ire1 in a mechanism of translational control dependent on 5UTR-embedded structural elements in fate determinant genes. Our data unveil developmental regulation of ribosome dynamics as post-transcriptional mechanisms orchestrating neuronal diversity establishment and assembly of cortical layers.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=93 SRC="FIGDIR/small/449563v2_ufig1.gif" ALT="Figure 1">
View larger version (19K):
org.highwire.dtl.DTLVardef@115a09aorg.highwire.dtl.DTLVardef@e30d2dorg.highwire.dtl.DTLVardef@763b27org.highwire.dtl.DTLVardef@2586e9_HPS_FORMAT_FIGEXP  M_FIG C_FIG HIGHLIGHTSO_LISmall molecule screening reveals Ire1 upstream of upper layer neuronal identity
C_LIO_LIPolarization and proper lamination of layer II/III neurons require Ire1
C_LIO_LIDevelopment of upper layers requires high translation rates driven by eIF4A1 and eEF-2 downstream of Ire1
C_LIO_LIeIF4A1-dependent Satb2 mRNA translation initiation is a mechanism of upper layer fate acquisition
C_LI


title: Ire1α-Regulated mRNA Translation Rate Controls the Identity and Polarity of Upper Layer Cortical Neurons


title: Entamoeba histolytica adaption to auranofin: a phenotypic and multi-omics characterization

abstract: 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.


title: LAP1 regulates nuclear plasticity to enable constrained migration

abstract: Metastasis involves dissemination of cancer cells away from a primary tumour and colonisation at distal sites. During this process, cancer cells must negotiate multiple physical constraints imposed by the microenvironment and tissue structure. The biophysical properties of the nucleus must be tuned since they pose a challenge to constrained migration. By analysing nuclear genes upregulated during the acquisition of metastatic potential, we discovered increased expression of the inner nuclear membrane protein LAP1 in metastatic melanoma cells and at the invasive fronts of human primary tumours and in metastases. Human cells express two LAP1 isoforms (LAP1B and LAP1C), which differ in their amino terminus. We found that whereas the longer isoform, LAP1B, binds more strongly to nuclear lamins and has restricted motility within the nuclear envelope, the shorter isoform, LAP1C, favours nuclear envelope blebbing and allows migration through constraints. We propose that LAP1 renders the nucleus plastic and enhances melanoma aggressiveness.


title: Graph-based machine learning reveals rules of spatiotemporal cell interactions in tissues

abstract: 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.


title: Probing the rules of cell coordination in live tissues by interpretable machine learning based on graph neural networks


title: A live-cell marker to visualize the dynamics of stable microtubules

abstract: 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.


title: Dystrophin regulates peripheral circadian SRF signalling

abstract: Dystrophin is a sarcolemmal protein essential for muscle contraction and maintenance, absence of which leads to the devastating muscle wasting disease Duchenne muscular dystrophy (DMD)[1, 2]. Dystrophin has an actin-binding domain [3-5], which specifically binds and stabilises filamentous (F)-actin[6], an integral component of the RhoA-actin-serum response factor (SRF)-pathway[7]. The RhoA-actin-SRF-pathway plays an essential role in circadian signalling whereby the hypothalamic suprachiasmatic nucleus, transmits systemic cues to peripheral tissues, activating SRF and transcription of clock target genes[8, 9]. Given dystrophin binds F-actin and disturbed SRF-signalling disrupts clock entrainment, we hypothesised that dystrophin loss causes circadian deficits. Here we show for the first time alterations in the RhoA-actin-SRF-signalling-pathway, in both dystrophin-deficient myotubes and dystrophic mouse models. Specifically, we demonstrate reduced F/G-actin ratios and nuclear MRTF, dysregulation of core clock and downstream target-genes, and down-regulation of key circadian genes in muscle biopsies from DMD patients harbouring an array of mutations. Further, disrupted circadian locomotor behaviour was observed in dystrophic mice indicative of disrupted SCN signalling, and indeed dystrophin protein was absent in the SCN of dystrophic animals. Dystrophin is thus a critically important component of the RhoA-actin-SRF-pathway and a novel mediator of circadian signalling in peripheral tissues, loss of which leads to circadian dysregulation.


title: Post-Insemination Selection Dominates Pre-Insemination Selection in Driving Rapid Evolution of Male Competitive Ability

abstract: Sexual reproduction is a complex process that contributes to differences between the sexes and divergence between species. From a males perspective, sexual selection can optimize reproductive success by acting on the variance in mating success (pre-insemination selection) as well as the variance in fertilization success (post-insemination selection). The balance between pre- and post-insemination selection has not yet been investigated using a strong hypothesis-testing framework that directly quantifies the effects of post-insemination selection on the evolution of reproductive success. Here we use experimental evolution of a uniquely engineered genetic system that allows sperm production to be turned off and on in obligate male-female populations of Caenorhabditis elegans. We show that enhanced post-insemination competition increases the efficacy of selection and surpasses pre-insemination sexual selection in driving a polygenic response in male reproductive success. We find that after 30 generations post-insemination selection increased male reproductive success by an average of 5- to 7-fold. Contrary to expectation, enhanced pre-insemination competition hindered selection and slowed the rate of evolution. Furthermore, we found that post-insemination selection resulted in a strong polygenic response at the whole-genome level. Our results demonstrate that post-insemination sexual selection plays a critical role in the rapid optimization of male reproductive fitness. Therefore, explicit consideration should be given to post-insemination dynamics when considering the population effects of sexual selection.


title: A TREM2-derived circRNA is upregulated in the entorhinal cortex of Alzheimer's disease patients

abstract: Circular RNAs (circRNAs) are a novel class of noncoding RNAs characterized by a covalent and stable closed loop structure. circRNAs are enriched in neural tissues, particularly at synapses, where they are involved in synaptic plasticity. Alzheimers disease (AD) is considered a synaptopathy since neurodegeneration causes loss or dysfunction of synapses. Microglia participate in synaptic pruning and also play a crucial role in developing AD. For instance, genetic variants in TREM2, a microglia-related gene, are risk factors for AD. Alterations in circRNAs expression have been described in different neurological diseases, including AD. However, no TREM2-derived circRNAs have been described so far. TREM2 has 3 linear RNA variants due to alternative splicing. We hypothesized that alternative splicing of exon 4 might be favoring circRNAs originating from TREM2 (circTREM2s), which in turn might be involved in AD pathogenesis. First, divergent primers (overlapping exons 3-4 and 4-5) were designed to amplify circRNAs by RT-PCR, which were confirmed by Sanger sequencing. Three candidate TREM2-derived circRNAs were identified on control human entorhinal samples. Then, additional primer sets were used to confirm back-splicing junctions. One of the circRNAs, circTREM2_1, was consistently amplified with all primer sets. In addition, circTREM2_1 was also present in AD entorhinal cortex samples and in HMC3 cells. We observed that circTREM2_1 is up-regulated in AD entorhinal cortex samples compared to controls, particularly at early stages of the disease, when we performed RT-qPCR. In conclusion, we have identified a novel circRNA derived from the TREM2 gene that could play a role in AD pathogenesis.

category: genetics


title: Patterns and Mechanisms of Sex Ratio Distortion in the Collaborative Cross Mouse Mapping Population

abstract:
In species with single-locus, chromosome-based mechanisms of sex determination, the laws of segregation predict an equal ratio of females to males at birth. Here, we show that departures from this Mendelian expectation are commonplace in the 8-way recombinant inbred Collaborative Cross (CC) mouse population. More than one-third of CC strains exhibit significant sex ratio distortion (SRD) at wean, with twice as many male-biased than female-biased strains. We show that these pervasive sex biases persist across multiple breeding environments, are stable over time, are not fully mediated by maternal effects, and are not explained by sex-biased neonatal mortality. SRD exhibits a heritable component, but QTL mapping analyses and targeted investigations of sex determination genes fail to nominate any large effect loci. These findings, combined with the reported absence of sex ratio biases in the CC founder strains, suggest that SRD manifests from multilocus combinations of alleles only uncovered in recombined CC genomes. We speculate that the genetic shuffling of eight diverse parental genomes during the early CC breeding generations led to the decoupling of sex-linked drivers from their co-evolved suppressors, unleashing complex, multiallelic systems of sex chromosome drive. Consistent with this interpretation, we show that several CC strains exhibit copy number imbalances at co-evolved X-and Y-linked ampliconic genes that have been previously implicated in germline genetic conflict and SRD in house mice. Overall, our findings reveal the pervasiveness of SRD in the CC population and nominate the CC as a powerful resource for investigating sex chromosome genetic conflict in action.

ARTICLE SUMMARYWe compiled breeding records from The Collaborative Cross (CC) mouse mapping population to quantify the frequency and explore potential mechanisms of sex ratio distortion. Strikingly, more than one-third of CC strains yield significantly sex-biased litters. These sex biases are not mediated by environmental effects and are moderately heritable. We conclude that the widespread sex ratio distortion in the CC manifests from multilocus permutations of selfish sex-linked elements and suppressors that are only recovered in the recombinant CC strains.


title: Parallel Evaluation of Nucleophilic and Electrophilic Chemical Probes for Sulfenic Acid: Reactivity, Selectivity and Biocompatibility

abstract: 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.


title: High quality, phased genomes of Phytophthora ramorum clonal lineages NA1 and EU1

abstract: Phytophthora ramorum is the causal agent of sudden oak death in West Coast forests and currently two clonal lineages, NA1 and EU1, cause epidemics in Oregon forests. Here, we report on two high-quality genomes of individuals belonging to the NA1 and EU1 clonal lineages respectively, using PacBio long-read sequencing. The NA1 strain Pr102, originally isolated from coast live oak in California, is the current reference genome and was previously sequenced independently using either Sanger (P. ramorum v1) or PacBio (P. ramorum v2) technology. The EU1 strain PR-15-019 was obtained from tanoak in Oregon. These new genomes have a total size of 57.5 Mb, with a contig N50 length of ~3.5-3.6 Mb and encode ~15,300 predicted protein-coding genes. Genomes were assembled into 27 and 28 scaffolds with 95% BUSCO scores and are considerably improved relative to the current JGI reference genome with 2,575 or the PacBio genomes with 1,512 scaffolds. These high-quality genomes provide a valuable resource for studying the genetics, evolution, and adaptation of these two clonal lineages.

category: plant biology


title: An F-actin mega-cable supports the migration of the sperm nucleus during the fertilization of the polarity-inverted central cell of Agave inaequidens

abstract: Asparagaceaes large embryo sacs display a central cell nucleus polarized toward the chalaza, which means the sperm nucleus that fuses it during double fertilization migrates a long distance before karyogamy. Because of the size and inverted polarity of the central cell in Asparagaceae, we hypothesize that the second fertilization process is supported by F-actin structures different from the short-range aster-like ones observed in Arabidopsis. Here, we analyzed the F-actin dynamics of Agave inaequidens, a typical Asparagaceae, before, during, and after central cell fertilization. Several parallel F-actin cables emerging from the nucleus within the central cell, enclosing the vacuole, and reaching the micropylar pole were observed. As fertilization progressed, a thick F-actin mega-cable traversing the vacuole appeared, connecting the central cell nucleus with the micropylar pole near the egg cell. This mega-cable wrapped the sperm nucleus in transit to fuse the central cell one. Once karyogamy finished, the mega-cable disassembled, but new F-actin structures formed during the endosperm development. These observations suggest that Asparagaceae, and probably other plant species with similar embryo sacs, evolved an F-actin machinery specifically adapted to support the migration of the fertilizing sperm nucleus within a large-sized and polarity-inverted central cell.