title: Stepwise evolution of a butterfly supergene via duplication and inversion

abstract: Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning colouration in the African Monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of >1 million bp underwent several segmental duplications at least 7.5 million years ago. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree, and suggest that structural polymorphisms have persisted for at least 4.1 million years. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes.


title: Summary statistics knockoff inference empowers identification of putative causal variants in genome-wide association studies

abstract: Recent advances in genome sequencing and imputation technologies provide an exciting opportunity to comprehensively study the contribution of genetic variants to complex phenotypes. However, our ability to translate genetic discoveries into mechanistic insights remains limited at this point. In this paper, we propose an efficient knockoff-based method, GhostKnockoff, for genome-wide association studies (GWAS) that leads to improved power and ability to prioritize putative causal variants relative to conventional GWAS approaches. The method requires only Z-scores from conventional GWAS and hence can be easily applied to enhance existing and future studies. The method can also be applied to meta-analysis of multiple GWAS allowing for arbitrary sample overlap. We demonstrate its performance using empirical simulations and two applications: (1) analysis of 1,403 binary phenotypes from the UK Biobank data in 408,961 samples of European ancestry, and (2) a meta-analysis for Alzheimers disease (AD) comprising nine overlapping large-scale GWAS, whole-exome and whole-genome sequencing studies. The UK Biobank analysis demonstrates superior performance of the proposed method compared to conventional GWAS in both statistical power (2.05-fold more discoveries) and localization of putative causal variants at each locus (46% less proxy variants due to linkage disequilibrium). The AD meta-analysis identified 55 risk loci (including 31 new loci) with ~70% of the proximal genes at these loci showing suggestive signal in downstream single-cell transcriptomic analyses. Our results demonstrate that GhostKnockoff can identify putatively functional variants with weaker statistical effects that are missed by conventional association tests.

category: genetics


title: Frontal plane ankle stiffness increases with axial load independent of muscle activity

abstract: Ankle sprains are the most common musculoskeletal injury, typically resulting from excessive inversion of the ankle. One way to prevent excessive inversion and maintain ankle stability is to generate a stiffness that is sufficient to resist externally imposed rotations. Frontal-plane ankle stiffness increases as participants place more weight on their ankle, but whether this effect is due to muscle activation or axial loading of the ankle is unknown. Identifying whether and to what extent axial loading affects ankle stiffness is important in understanding what role the passive mechanics of the ankle joint play in maintaining its stability. The objective of this study was to determine the effect of passive axial load on frontal-plane ankle stiffness. We had subjects seated in a chair as an axial load was applied to the ankle ranging from 10% to 50% body weight. Small rotational perturbations were applied to the ankle in the frontal plane to estimate stiffness. We found a significant, linear, 3-fold increase in ankle stiffness with axial load from the range of 0% bodyweight to 50% bodyweight. This increase could not be due to muscle activity as we observed no significant axial-load-dependent change in any of the recorded muscle activations. These results demonstrate that axial loading is a significant contributor to maintaining frontal-plane ankle stability, and that disruptions to the mechanism mediating this sensitivity of stiffness to axial loading may result in pathological cases of ankle instability.

category: bioengineering


title: Through-plane super-resolution with autoencoders in diffusion magnetic resonance imaging of the developing human brain

abstract: Fetal brain diffusion magnetic resonance images are often acquired with a lower through-plane than in-plane resolution. This anisotropy is often overcome by classical upsampling methods such as linear or cubic interpolation. In this work, we employ an unsupervised learning algorithm using an autoencoder neural network to enhance the through-plane resolution by leveraging a large amount of data. Our framework, which can also be used for slice outliers replacement, overperformed conventional interpolations quantitatively and qualitatively on pre-term newborns of the developing Human Connectome Project. The evaluation was performed on both the original diffusion-weighted signal and on the estimated diffusion tensor maps. A byproduct of our autoencoder was its ability to act as a denoiser. The network was able to generalize to fetal data with different levels of motion and we qualitatively showed its consistency, hence supporting the relevance of pre-term datasets to improve the processing of fetal brain images.

category: bioinformatics


title: TargetGeneReg 2.0: a comprehensive web-atlas for p53, p63, and cell cycle-dependent gene regulation

abstract: In recent years, our web-atlas at www.TargetGeneReg.org has enabled many researchers to uncover new biological insights and to identify novel regulatory mechanisms that affect p53 and the cell cycle - signaling pathways that are frequently dysregulated in diseases like cancer. Here, we provide a substantial upgrade of the database that comprises an extension to include non-coding genes and the transcription factors {Delta}Np63 and RFX7. TargetGeneReg 2.0 combines gene expression profiling and transcription factor DNA binding data to determine, for each gene, the response to p53, {Delta}Np63, and cell cycle signaling. It can be used to dissect common, cell type, and treatment-specific effects, identify the most promising candidates, and validate findings. We demonstrate the increased power and more intuitive layout of the resource using realistic examples.


title: Graphsite: ligand-binding site classification using deep graph neural network

abstract: Binding sites are concave surfaces on proteins that bind to small molecules called ligands. Types of molecules that bind to the protein determine its biological function. Meanwhile, the binding process between small molecules and the protein is also crucial to various biological functionalities. Therefore, identifying and classifying such binding sites would enormously contribute to biomedical applications such as drug repurposing. Deep learning is a modern artificial intelligence technology. It utilizes deep neural networks to handle complex tasks such as image classification and language translation. Previous work has proven the capability of deep learning models handle binding sites wherein the binding sites are represented as pixels or voxels. Graph neural networks (GNNs) are deep learning models that operate on graphs. GNNs are promising for handling binding sites related tasks - provided there is an adequate graph representation to model the binding sties. In this communication, we describe a GNN-based computational method, GraphSite, that utilizes a novel graph representation of ligand-binding sites. A state-of-the-art GNN model is trained to capture the intrinsic characteristics of these binding sites and classify them. Our model generalizes well to unseen data and achieves test accuracy of 81.28% on classifying 14 binding site classes.


title: Control of Leaf Vein Patterning by Regulated Plasmodesma Aperture

abstract: To form tissue networks, animal cells migrate and interact through proteins protruding from their plasma membranes. Plant cells can do neither, yet plants form vein networks. How plants do so is unclear, but veins are thought to form by the coordinated action of the polar transport and signal transduction of the plant hormone auxin. However, plants inhibited in both pathways still form veins. Patterning of vascular cells into veins is instead prevented in mutants lacking the function of the GNOM (GN) regulator of auxin transport and signaling, suggesting the existence of at least one more GN-dependent vein-patterning pathway. Here we show that pathway depends on the movement of an auxin signal through plasmodesmata (PDs) intercellular channels. PD permeability is high where veins are forming, lowers between veins and nonvascular tissues, but remains high between vein cells. Impaired ability to regulate PD aperture leads to defects in auxin transport and signaling, ultimately leading to vein patterning defects that are enhanced by inhibition of auxin transport or signaling. GN controls PD aperture regulation, and simultaneous inhibition of auxin signaling, auxin transport, and regulated PD aperture phenocopies null gn mutants. Therefore, veins are patterned by the coordinated action of three GN-dependent pathways: auxin signaling, polar auxin transport, and movement of an auxin signal through PDs. We have identified all the key vein-patterning pathways in plants and an unprecedented mechanism of tissue network formation in multicellular organisms.

category: plant biology


title: Phosphorylation dynamics in a flg22 induced, heterotrimeric G protein dependent signaling network in Arabidopsis thaliana reveals a candidate PP2A phosphatase involved in AtRGS1 trafficking

abstract: flg22 is a 22-amino peptide released from bacterial flagellin, a Microbe-Associated Molecular Pattern (MAMP) that is recognized by the plant cell as a signal indicating that bacteria are present. On its own, flg22 initiates a rapid increase in cytoplasmic calcium, extracellular reactive oxygen species, and activation of a Mitogen Activated Protein Kinase (MAPK) cascade, all of which are activated within 15 minutes after the cell perceives flg22. Here we show a massive change in protein abundance and phosphorylation state of the Arabidopsis root cell proteome within this 15-minute duration in wildtype and a mutant deficient in G-protein coupled signaling. Integration of phosphoproteome with protein-protein interactome data followed by network topology analyses discovered that many of the flg22-induced phosphoproteome changes fall on proteins that comprise the G protein interactome and on the most highly populated hubs of the immunity network. Approximately 95% of the phosphorylation changes in the G-protein interactome depend on a functional heterotrimeric G protein complex, some occur on proteins that interact directly with components of G-coupled signal transduction. One of these is ATB, a substrate-recognition sub-unit of the PP2A Ser/Thr phosphatase and an interactor to Arabidopsis thaliana REGULATOR OF G SIGNALING 1 protein (AtRGS1), a 7-transmembrane spanning modulator of the nucleotide-binding state of the core G protein complex. AtRGS1 is phosphorylated by BAK1, a component of the flg22 receptor, to initiate AtRGS1 endocytosis. A null mutation of ATB confers high basal endocytosis of AtRGS1, suggesting sustained phosphorylated status. Loss of ATB confers traits associated with loss of AtRGS1. Because the basal level of AtRGS1 is lower in the atb null mutant in a proteasome-dependent manner, we propose that phosphorylation-dependent endocytosis of AtRGS1 is part of a mechanism to degrade AtRGS1 which then sustains activation of the G protein complex. Thus, the role of ATB is now established as a central component of phosphorylation-dependent regulation of system dynamics in innate immunity.

category: systems biology


title: Neuronal population activity dynamics reveal a low-dimensional signature of operant learning

abstract: Learning engages a high-dimensional neuronal population space spanning multiple brain regions. We identified a low-dimensional signature associated with operant conditioning, a ubiquitous form of learning in which animals learn from the consequences of behavior. Using single-neuron resolution voltage imaging, we identified two low-dimensional motor modules in the neuronal population underlying Aplysia feeding. Our findings point to a temporal shift in module recruitment as the primary signature of operant learning.

category: neuroscience


title: Visual imagery vividness declines across the lifespan

abstract: The capacity to elicit vivid visual mental images varies within an extensive range across individuals between hyper- and aphantasia. It is not clear, however, whether imagery vividness is constant across the lifespan or changes during development and later in life. Without enforcing the constraints of strict experimental procedures and representativity across the entire population, our purpose was to take a first look at the self-reported level of imagery vividness and determine the relative proportion of aphantasic/hyperphantasic participants in different age-groups. Relying on the frequently used Vividness of Visual Mental Imagery questionnaire, we collected data on a random sample of 2252 participants between the ages of 12 to 60 years. We found a novel developmental pattern describing a declining ability to elicit vivid visual mental images in the group averages of different age-groups from adolescence to middle age. This effect involves both a decreasing proportion of individuals with very vivid imagery and an increasing proportion of individuals with weak imagery as maturation (assessed by boneage estimations in adolescents) and aging progresses. This finding may help to shed light on yet unknown developmental mechanisms of our internal, stimulus-independent processes, and might also help to determine genetic, maturational, and age-dependent factors in the extreme cases of hyper- and aphantasia.


title: Transactive response DNA-binding Protein (TDP-43) regulates early HIV-1 entry and infection.

abstract: The transactive response DNA-binding protein (TDP-43) is an important regulator of mRNA, being reported to stabilize the anti-HIV factor, histone deacetylase 6 (HDAC6). However, little is known about the role of TDP-43 in HIV infection. In this work, we seek for the TDP-43 function on regulating CD4+ T cell permissibility to HIV infection. We observed that over-expression of wt-TDP-43 in CD4+ T cells stabilized HDAC6, increasing mRNA and the protein levels of this antiviral enzyme. Under this experimental condition, HIV-1 infection was impaired, independently of the viral envelope glycoprotein (Env) complex tropism. The results obtained by using an HIV-1 Env-mediated cell-to-cell fusion model, under the same experimental conditions, suggest that the increase in TDP-43 levels negatively affects the viral Env fusion capacity. Moreover, the specific siRNA silencing of endogenous TDP-43 in target cells lead to a significant decrease in the levels of HDAC6 which consistently induces an increase in the fusogenic and infection activities of the HIV-1 Env. These observations were confirmed by using primary viral Envs from HIV+ individuals with different clinical phenotypes. An increase in the level of expression of wt-TDP-43 strongly reduced the Envs infection activity of viremic non-progressors (VNP) and rapid progressors (RP) HIV+ individuals down to the levels of the inefficient HIV-1 Envs from long-term non-progressor elite controllers (LTNP-EC) individuals. On the contrary, low levels of endogenous TDP-43, obtained after specific siRNA-TDP-43 knocking-down, significantly favors the infection activity of primary HIV-1 Envs of VNP and RP individuals, leading to an increase in the infection ability of the primary HIV-1/LTNP-EC Envs. Based on this evidence, we interpret that TDP-43 conditions cell permissibility to HIV infection by affecting viral Env fusion and infection capacities, at least by altering the cellular levels of the antiviral enzyme HDAC6.

category: immunology


title: EB3-informed dynamics of the microtubule stabilizing cap during stalled growth

abstract: Microtubule stability is known to be governed by a stabilizing GTP/GDP-Pi cap, but the exact relation between growth velocity, GTP hydrolysis and catastrophes remains unclear. We investigate the dynamics of the stabilizing cap through in vitro reconstitution of microtubule dynamics in contact with micro-fabricated barriers, using the plus-end binding protein GFP-EB3 as a marker for the nucleotide state of the tip. The interaction of growing microtubules with steric objects is known to slow down microtubule growth and accelerate catastrophes. We show that the lifetime distributions of stalled microtubules, as well as the corresponding lifetime distributions of freely growing microtubules, can be fully described with a simple phenomenological 1D model based on noisy microtubule growth and a single EB3-dependent hydrolysis rate. This same model is furthermore capable of explaining both the previously reported mild catastrophe dependence on microtubule growth rates and the catastrophe statistics during tubulin washout experiments.

category: cell biology


title: Slit-Robo Signalling Establishes a Sphingosine-1-Phosphate Gradient to Polarise Fin Mesenchyme and Establish Fin Morphology

abstract: Immigration of mesenchymal cells into the growing fin and limb buds drives distal outgrowth, with subsequent tensile forces between these cells essential for fin and limb morphogenesis. Morphogens derived from the apical domain of the fin, orientate limb mesenchyme cell polarity, migration, division and adhesion. The zebrafish mutant stomp displays defects in fin morphogenesis including blister formation and associated loss of orientation and adhesion of immigrating fin mesenchyme cells. Positional cloning of stomp identified a mutation in the gene encoding the axon guidance ligand, Slit3. We provide evidence that Slit ligands derived from immigrating mesenchyme act via Robo receptors at the Apical Ectodermal Ridge (AER) to promote release of sphingosine-1-phosphate (S1P). S1P subsequently diffuses back to the mesenchyme to promote their polarisation, orientation, positioning and adhesion to the interstitial matrix of the fin fold. We thus demonstrate coordination of the Slit-Robo and S1P signalling pathways in fin fold morphogenesis. Our work introduces a mechanism regulating the orientation, positioning and adhesion of its constituent cells.


title: The war of corals: patterns, drivers, and implications of changing coral competitive performances across reef environments

abstract: Amidst global environmental changes, predicting species responses to future environments is a critical challenge for preserving biodiversity and associated human benefits. We explored the original idea that coral competitive performances, the ability of corals to preempt ecological space on the reef through territorial warfare, serve as indicators of species ecological niches and environmental windows, and therefore, responses to future environments. Our surveys indicated that coral performances varied with taxonomic-identity, size, and position along environmental gradients, highlighting complex interplays between life-history, warfare-strategy, and niche segregation. Our results forewarn that growing alterations of coastal environments may trigger shifts in coral dominance, with decline of major reef-building taxa like acroporids, and underscore the importance of restraining human impacts for coastal resilience. Our empirical approach untangles the complexity of species battle-like interactions and can help identify winners and losers in various communities caught in the interplay between ecological niches, environmental windows, and global changes.

category: ecology


title: Using simulated temperature regimes to test growth and development of an invasive forest insect under climate change

abstract: Temperature and its impact on fitness are fundamental for understanding range shifts and population dynamics under climate change. Geographic climate heterogeneity, behavioural and physiological plasticity, and thermal adaptation to local climates makes predicting the responses of species to climate change complex. Using larvae from seven geographically distinct wild populations in the eastern United States of the non-native forest pest Lymantria dispar dispar (L.), we conducted a simulated reciprocal transplant experiment in environmental chambers using six custom temperature regimes representing contemporary conditions near the southern and northern extremes of the US invasion front and projections under two climate change scenarios for the year 2050. Larval growth rates increased with climate warming compared to current thermal regimes and responses differed by population. A significant population-by-treatment interaction indicated that growth rates increased more when a source population experienced the warming scenarios for their region, especially for southern populations. Our study demonstrates the utility of simulating thermal regimes under climate change in environmental chambers and emphasizes how the impacts from future increases in temperature can be heterogeneous due to geographic differences in climate-related performance among populations.


title: Experimental calibration of trapping methods for addressing bias in arthropod biodiversity monitoring

abstract: 1.All approaches to biodiversity monitoring have inherent biases in the taxa captured, yet, as environments, sampling goals, and conventions vary, it is not uncommon for sampling approaches to be customized to reflect the study objectives, optimizing findings to be locally relevant but at the cost of transferability. Here, we developed a calibration study to directly examine how researcher trap choice affects observations made in insect biodiversity sampling. Sampling efficiency of four types of traps: pitfalls, yellow ramp traps, a novel jar ramp trap, and yellow sticky cards, were compared with respect to an array of biodiversity metrics associated with the arthropods they captured. We found that trapping efficiency and functional groups of arthropods (flying versus ground-crawling) varied by trap type. Pitfalls and jar ramp traps performed similarly for all biodiversity metrics measured, suggesting that jar ramp traps provide a more comparable measurement of ground-crawling arthropod communities to pitfall sampling than the yellow ramp traps and should be considered when pitfall sampling cannot be used. This study illustrates the implications for biodiversity sampling of arthropods in environments with physical constraints on trapping, and the importance of directly comparing adapted methods to established sampling protocol. Future biodiversity monitoring schemes should conduct calibration experiments to provide important information on performance and potential limitations of sampling methodology.


title: A convolutional neural network highlights mutations relevant to antimicrobial resistance in Mycobacterium tuberculosis

abstract: Long diagnostic wait times hinder international efforts to address multi-drug resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited in part by a lack of interpretability and verifiability, especially in deep learning methods. Here, we present a deep convolutional neural network (CNN) that predicts the antibiotic resistance phenotypes of M. tuberculosis isolates. The CNN performs with state-of-the-art levels of predictive accuracy. Evaluation of salient sequence features permits biologically meaningful interpretation and validation of the CNNs predictions, with promising repercussions for functional variant discovery, clinical applicability, and translation to phenotype prediction in other organisms.

category: genomics


title: Fast off-rate CD229 chimeric antigen receptor T cells efficiently target multiple myeloma, spare T cells, and exhibit reduced trogocytosis

abstract:
T cells expressing chimeric antigen receptors have shown remarkable therapeutic activity against different types of cancer. However, their wider use has been hampered by the potential for life-threatening toxicities due to the unintended targeting of healthy cells expressing low levels of the targeted antigen. We have now developed an affinity-tuning approach for the generation of minimally modified, low-affinity antibody variants derived from existing high-affinity antibodies. Using this approach, we engineered low affinity variants of the fully human CD229-specific antibody 2D3. Parental 2D3 originally efficiently targeted multiple myeloma cells but also healthy T cells expressing low levels of CD229. We demonstrate that CAR T cells based on a low affinity variant of 2D3 maintain the parental antibodys anti-tumor activity, but lack its targeting of healthy T cells. In addition, variant CD229 CAR T cells show reduced trogocytosis potentially augmenting CAR T cell persistence. The fast off-rate CAR produced using our affinity tuning approach eliminates a key liability of CD229 CAR T cells and paves the way for the effective and safe treatment of patients with multiple myeloma.

One sentence summaryAffinity tuning approach yields low affinity CD229 CAR binding domain maintaining the parental clones anti-tumor activity while eliminating killing of healthy T cells, increasing CAR T cell expansion, and decreasing trogocytosis.


title: A method to build extended sequence context models of point mutations and indels

abstract: The mutation rate of a specific position in the human genome depends on the sequence context surrounding it. Modeling the mutation rate by estimating a rate for each possible k-mer, however, only works for small values of k since the data becomes too sparse for larger values of k. Here we propose a new method that solves this problem by grouping similar k-mers using IUPAC patterns. We refer to the method as k-mer pattern partition and have implemented it in a software package called kmerPaPa. We use a large set of human de novo mutations to show that this new method leads to improved prediction of mutation rates and makes it possible to create models using wider sequence contexts than previous studies. Revealing that for some mutation types, the mutation rate of a position is significantly affected by nucleotides that are up to four base pairs away. As the first method of its kind, it does not only predict rates for point mutations but also indels. We have additionally created a software package called Genovo that, given a k-mer pattern partition model, predicts the expected number of synonymous, missense, and other functional mutation types for each gene. Using this software, we show that the created mutation rate models increase the statistical power to detect genes containing disease-causing variants and to identify genes under strong constraint, e.g. haploinsufficient genes.


title: CMash: fast, multi-resolution estimation of k-mer-based Jaccard and containment indices

abstract: K-mer based methods are used ubiquitously in the field of computational biology. However, determining the optimal value of k for a specific application often remains heuristic. Simply reconstructing a new k-mer set with another k-mer size is computationally expensive, especially in metagenomic analysis where data sets are large. Here, we introduce a hashing-based technique that leverages a kind of bottom-m sketch as well as a k-mer ternary search tree (KTST) to obtain k-mer based similarity estimates for a range of k values. By truncating k-mers stored in a pre-built KTST with a large k = kmax value, we can simultaneously obtain k-mer based estimates for all k values up to kmax. This truncation approach circumvents the reconstruction of new k-mer sets when changing k values, making analysis more time and space-efficient. For example, we show that when using a KTST to estimate the containment index between a RefSeq-based microbial reference database and simulated metagenome data for 10 values of k, the running time is close to 10x faster compared to a classic MinHash approach while using less than one-fifth the space to store the data structure. A python implementation of this method, CMash, is available at https://github.com/dkoslicki/CMash. The reproduction of all experiments presented herein can be accessed via https://github.com/KoslickiLab/CMASH-reproducibles.


title: Structure-based design of a chemical probe set for the 5-HT5A serotonin receptor

abstract:
The 5-HT5A receptor (5-HT5AR), for which no selective agonists and only a few antagonists exist, remains the least understood serotonin (5-HT) receptor. A single commercial antagonist (SB-699551) has been widely used to investigate central nervous system (CNS) 5-HT5AR function in neurological disorders, including pain. However, because SB-699551 has affinity for many 5-HTRs, lacks inactive property-matched controls, and has assay interference concerns, it has liabilities as a chemical probe. To better illuminate 5-HT5AR function, we developed a probe set through iterative rounds of molecular docking, pharmacological testing, and optimization. Docking over six million lead-like molecules against a 5-HT5AR homology model identified five mid-M ligand starting points with unique scaffolds. Over multiple rounds of structure-based design and testing, a new quinoline scaffold with high affinity and enhanced selectivity for the 5-HT5AR was developed, leading to UCSF678, a 42 nM arrestin-biased partial agonist at the 5-HT5AR with a much more restricted off-target profile and decreased assay liabilities vs. SB-699551. Site-directed mutagenesis supported the docked pose of UCSF678, which was also consistent with recent published 5-HTR structures. Surprisingly, property-matched analogs of UCSF678 that were either inactive across 5-HTRs or retained affinity for UCSF678s off-targets revealed that 5-HT5AR engagement is nonessential for alleviating pain in a mouse model, contrary to previous studies using less-selective ligands. Relative to SB-699551, these molecules constitute a well-characterized and more selective probe set with which to study the function of the 5-HT5A receptor.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/471449v1_ufig1.gif" ALT="Figure 1">
View larger version (28K):
org.highwire.dtl.DTLVardef@8fb24dorg.highwire.dtl.DTLVardef@2dd84org.highwire.dtl.DTLVardef@1d3346org.highwire.dtl.DTLVardef@1117331_HPS_FORMAT_FIGEXP  M_FIG Table of Contents Graphic C_FIG

category: biophysics


title: Processing body dynamics drive non-genetic MEK inhibitors tolerance by fine-tuning KRAS and NRAS translation

abstract: Overactivation of the Mitogen-activated protein kinase (MAPK) pathway is a critical driver of many human cancers. However, therapies targeting this pathway have proven effective in only a few cancers, as cancers inevitably develop resistance. Puzzling observations suggest that MAPK targeting fails in tumors due to early compensatory RAS overexpression, albeit by unexplained mechanisms. We identified a novel mechanism of drug tolerance to MEK inhibitors (MEKi) that involves Processing Bodies (PBs), a membraneless organelle (MLO). MEKi promoted translation of the oncogenes KRAS and NRAS, which in turn triggered BRAF phosphorylation. This overexpression, which occurred in the absence of neotranscription, depended on PB dissolution as the source of the RAS mRNA. Moreover, in response to MEKi removal, the process was dynamic as PBs rapidly reformed and reduced MAPK signaling. These results highlight a dynamic spatiotemporal negative feedback loop of MAPK signaling via RAS mRNA sequestration. Furthermore, we observed a phenotype with a low number of PBs along with strong KRAS and NRAS induction capacities. Overall, we describe a new intricate mechanism involving PBs in the translational regulation of essential cellular signaling pathways like MAPKs, paving the way for future therapies altering MLO and thereby improving targeted cancer therapies.

category: cancer biology


title: A remote lecture series roadmap to equity, diversity, and inclusion in STEM

abstract: Within a year of the shutdowns caused by the COVID-19 pandemic, virtual meetings transformed from an auxiliary service to an essential work platform for hundreds of millions of people worldwide. Universities rapidly accelerated adoption of virtual platforms for remote conferences, classes, and seminars amidst a second crisis that started in May 2020 and has tested institutional commitment to the principles of equity, diversity, and inclusion. To address these concerns, we began the Diversity and Science Lecture series (DASL), a cross-institutional national platform where junior life scientists present personal stories, professional progress, and advice for their peers. We first demonstrate benchmarking of speaker diversity using public data and review scientific talk content by text mining of talk titles from DASL. We then review insight into equity, diversity, and inclusion in STEM presented by DASL speakers. Notably, nearly every speaker highlighted the importance of social or interpersonal support to their career progression. Three-fifths of speakers commented on race or ethnicity. Other identities received considerable attention from a minority of speakers: immigrant identity, gender identity, mental illness, sexual minorities, disability, and rural or agricultural background. In addition, we report that ancillary small group sessions with speakers achieve greater engagement when held before formal seminars versus after. To facilitate future programming, we provide exhaustive guides that frame the goals of virtual programming in terms of benefits for executive team members, speakers, and attendees. In total, we have curated over 75 resources for personal and professional development cited by trainee life scientists, attracted over 4,000 unique website visitors, and exceeded 3,500 hours of engagement from attendees at weekly seminars. Our success with DASL confirms that community building can flourish under a remote lecture framework.

category: scientific communication and education


title: BDNF-dependent modulation of axonal transport is selectively impaired in ALS

abstract: Axonal transport ensures long-range delivery of essential cargoes between proximal and distal compartments of neurons, and is needed for neuronal development, function, and survival. Deficits in axonal transport have been detected at pre-symptomatic stages in mouse models of amyotrophic lateral sclerosis (ALS), suggesting that impairments are fundamental for disease pathogenesis. However, the precise mechanisms responsible for the transport deficits and whether they preferentially affect -motor neuron (MN) subtypes remain unresolved. Here, we report that stimulation of wild-type neurons with brain-derived neurotrophic factor (BDNF) enhances trafficking of signalling endosomes specifically in fast MNs (FMNs). In early symptomatic SOD1G93A mice, FMNs display selective impairment of axonal transport and develop an insensitivity to BDNF stimulation, with pathology upregulating classical non-pro-survival receptors in muscles and sciatic nerves. Altogether, these data indicate that cell- and non-cell autonomous BDNF signalling is impaired in vulnerable SOD1G93A MNs, thus identifying a new key deficit in ALS amenable for future therapeutic interventions.


title: SARS-CoV-2 infection of olfactory epithelial cells and neurons drives acute lung injury and lethal COVID-19 in mice

abstract: Lethal COVID-19 is associated with respiratory failure that is thought to be caused by acute respiratory distress syndrome (ARDS) secondary to pulmonary infection. To date, the cellular pathogenesis has been inferred from studies describing the expression of ACE2, a transmembrane protein required for SARS-CoV-2 infection, and detection of viral RNA or protein in infected humans, model animals, and cultured cells. To functionally test the cellular mechanisms of COVID-19, we generated hACE2fl animals in which human ACE2 (hACE2) is expressed from the mouse Ace2 locus in a manner that permits cell-specific, Cre-mediated loss of function. hACE2fl animals developed lethal weight loss and hypoxemia within 7 days of exposure to SARS-CoV-2 that was associated with pulmonary infiltrates, intravascular thrombosis and patchy viral infection of lung epithelial cells. Deletion of hACE2 in lung epithelial cells prevented viral infection of the lung, but not weight loss, hypoxemia or death. Inhalation of SARS-CoV-2 by hACE2fl animals resulted in early infection of sustentacular cells with subsequent infection of neurons in the neighboring olfactory bulb and cerebral cortex-- events that did not require lung epithelial cell infection. Pharmacologic ablation of the olfactory epithelium or Foxg1Cre mediated deletion of hACE2 in olfactory epithelial cells and neurons prevented lethality and neuronal infection following SARS-CoV-2 infection. Conversely, transgenic expression of hACE2 specifically in olfactory epithelial cells and neurons in Foxg1Cre; LSL-hACE2 mice was sufficient to confer neuronal infection associated with respiratory failure and death. These studies establish mouse loss and gain of function genetic models with which to genetically dissect viral-host interactions and demonstrate that lethal disease due to respiratory failure may arise from extrapulmonary infection of the olfactory epithelium and brain. Future therapeutic efforts focused on preventing olfactory epithelial infection may be an effective means of protecting against severe COVID-19.

category: microbiology


title: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: the case of the spike A222V mutation

abstract: The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has now reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

category: molecular biology


title: Omicron: A heavily mutated SARS-CoV-2 variant exhibits stronger binding to ACE2 and potently escape approved COVID-19 therapeutic antibodies

abstract: The new SARS-CoV-2 variant of concern "Omicron" was recently (Nov. 24th. 2021) spotted in South Africa and already spread around the world due to its enhanced transmissibility. The variant became conspicuous as it harbors more than thirty mutations in the spike protein with 15 mutations in the RBD region alone, potentially dampening the potency of therapeutic antibodies and enhancing the ACE2 binding. More worrying, Omicron infections have been reported in individuals who have received vaccines jabs in South Africa and Hong Kong. Here, we investigated the binding strength of Omicron with ACE2 and seven monoclonal antibodies that are either approved by FDA for COVID-19 therapy or undergoing phase III clinical trials. Computational mutagenesis and binding free energies could confirm that Omicron Spike binds ACE2 stronger than prototype SARS-CoV-2. Notably, three substitutions, i.e., T478K, Q493K, and Q498R, significantly contribute to the binding energies and doubled electrostatic potential of the RBDOmic-ACE2 complex. Instead of E484K substitution that helped neutralization escape of Beta, Gamma, and Mu variants, Omicron harbors E484A substitution. Together, T478K, Q493K, Q498R, and E484A substitutions contribute to a significant drop in the electrostatic potential energies between RBDOmic-mAbs, particularly in Etesevimab, Bamlanivimab, and CT-p59. CDR diversification could help regain the neutralization strength of these antibodies; however, we could not conduct this analysis to this end. Conclusively, our findings suggest that Omicron binds ACE2 with greater affinity, enhancing its infectivity and transmissibility. Mutations in the Spike are prudently devised by the virus that enhances the receptor binding and weakens the mAbs binding to escape the immune response.


title: The High Transmission of SARS-CoV-2 Omicron (B.1.1.529) Variant is Not Only Due to Its hACE2 binding: A Free Energy of Perturbation Study

abstract:
The mutations in the spike protein of SARS-CoV-2 Omicron variant (B.1.1.529 lineage) gave rise to questions, but the data on the mechanism of action at the molecular level is limited. In this study, we present the Free energy of perturbation (FEP) data about the RBD-hACE2 binding of this new variant.

We identified two groups of mutations located close to the most contributing substitutions Q498R and Q493R, which altered significantly the RBD-hACE2 interactions. The Q498R, Y505H and G496S mutations, in addition to N501Y, highly increased the binding to hACE2. They enhanced the binding by 98, 14 and 13 folds, respectively, which transforms the S1-RBD to a picomolar binder. However, in contrast to the case in mice the Q493R/K mutations, in a combination with K417N and T478K, dramatically reduced the S1 RBD binding by over 100 folds. The N440K, G446S and T478K substitutions had lesser contribution. Thus, the total effect of these nine mutations located on the interaction surface of RBD-hACE2 turns out to be similar to that observed in the Alpha variant. In a special circumstances it could be further altered by the E484A and S477N mutations and even lower binding capacity is likely to be detected. Finally, we provide a structural basis of the observed changes in the interactions.

These data may explain only partially the observed in South Africa extremely high Omicron spread and is in support to the hypothesis for multiple mechanisms of actions involved in the transmission.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=109 SRC="FIGDIR/small/471246v1_ufig1.gif" ALT="Figure 1">
View larger version (64K):
org.highwire.dtl.DTLVardef@144d901org.highwire.dtl.DTLVardef@10310e7org.highwire.dtl.DTLVardef@4ac7dborg.highwire.dtl.DTLVardef@1870231_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: biochemistry


title: Imaging-Guided Bioreactor for De-Epithelialization and Long-Term Cultivation of Ex Vivo Rat Trachea

abstract: Recent synergistic advances in organ-on-chip and tissue engineering technologies offer opportunities to create in vitro-grown tissue or organ constructs that can faithfully recapitulate their in vivo counterparts. Such in vitro tissue or organ constructs can be utilized in multiple applications, including rapid drug screening, high-fidelity disease modeling, and precision medicine. Here, we report an imaging-guided bioreactor that allows in situ monitoring of the lumen of ex vivo airway tissues during controlled in vitro tissue manipulation and cultivation of isolated rat trachea. Using this platform, we demonstrated selective removal of the rat tracheal epithelium (i.e., de-epithelialization) without disrupting the underlying subepithelial cells and extracellular matrix. Through different tissue evaluation assays, such as immunofluorescent staining, DNA/protein quantification, and electron beam microscopy, we showed that the epithelium of the tracheal lumen can be effectively removed with negligible disruption in the underlying tissue layers, such as cartilage and blood vessel. Notably, using a custom-built micro-optical imaging device integrated with the bioreactor, the trachea lumen was visualized at the cellular level in real time, and removal of the endogenous epithelium and distribution of locally delivered exogenous cells were demonstrated in situ. Moreover, the de-epithelialized trachea supported on the bioreactor allowed attachment and growth of exogenous cells seeded topically on its denuded tissue surface. Collectively, the results suggest that our imaging-enabled rat trachea bioreactor and selective cell replacement method can facilitate creating of bioengineered in vitro airway tissue that can be used in different biomedical applications.


title: Constructing a multiple-layer interactome for SARS-CoV-2 in the context of lung disease: Linking the virus with human genes and co-infecting microbes

abstract: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has caused millions of deaths worldwide. Many efforts have focused on unraveling the mechanism of the viral infection to develop effective strategies for treatment and prevention. Previous studies have provided some clarity on the protein-protein interaction linkages occurring during the life cycle of viral infection; however, we lack a complete understanding of the full interactome, comprising human miRNAs and protein-coding genes and co-infecting microbes. To comprehensively determine this, we developed a statistical modeling method using latent Dirichlet allocation (called MLCrosstalk, for multiple-layer crosstalk) to fuse many types of data to construct the full interactome of SARS-CoV-2. Specifically, MLCrosstalk is able to integrate samples with multiple layers of information (e.g., miRNA and microbes), enforce a consistent topic distribution on all data types, and infer individual-level linkages (i.e., differing between patients). We also implement a secondary refinement with network propagation to allow our microbe-gene linkages to address larger network structures (e.g., pathways). Using MLCrosstalk, we generated a list of genes and microbes linked to SARS-CoV-2. Interestingly, we found that two of the identified microbes, Rothia mucilaginosa and Prevotella melaninogenica, show distinct patterns representing synergistic and antagonistic relationships with the virus, respectively. We also identified several SARS-COV-2-associated pathways, including the VEGFA-VEGFR2 and immune response pathways, which may provide potential targets for drug design.


title: A zebrafish model of COVID-19-associated cytokine storm syndrome reveals differential proinflammatory activities of Spike proteins of SARS-CoV-2 variants of concern.

abstract:
The sudden and unexpected appearance of the COVID-19 pandemic turned the whole world upside down in a very short time. One of the main challenges faced has been to understand COVID-19 patient heterogeneity, as a minority develop life-threatening hyperinflammation, the so-called cytokine storm syndrome (CSS). Using the unique advantages of the zebrafish model we report here the proinflammatory role of Spike (S) proteins from different SARS-CoV-2 variants of concern after injection into the hindbrain ventricle, a cavity filled with cerebrospinal fluid to which immune cells can be easily recruited and that mimics the alveolar environment of the human lung. We found that wild type/Wuhan variant S1 (S1WT) protein promoted neutrophil and macrophage recruitment, local and systemic hyperinflammation, emergency myelopoiesis, and hemorrhages. In addition, S1{gamma} protein was more proinflammatory and S1{delta} was less proinflammatory than S1WT and, strikingly, S1{beta} promoted delayed and long-lasting inflammation. Pharmacological inhibition of the canonical inflammasome robustly alleviated S1 protein-induced inflammation and emergency myelopoiesis. In contrast, genetic inhibition of angiotensin-converting enzyme 2 strengthened the proinflammatory activity of S1, and the administration of angiopoietin (1-7) fully rescued S1-induced hyperinflammation and hemorrhages. These results shed light into the mechanisms orchestrating the COVID-19-associated CSS and the host immune response to different SARS-CoV-2 S protein variants.

HighlightsO_LIS proteins of SARS-CoV-2 promote hyperinflammation, neutrophilia, monocytosis and hemorrhages in zebrafish.
C_LIO_LIS protein effects in zebrafish are mediated via the canonical inflammasome and the Ace2/Angiopoietin (1-7) axis.
C_LIO_LIDelta S1 is less proinflammatory than wild type S1 and fails to induce emergency myelopoiesis in zebrafish.
C_LIO_LINaive and primed human white blood cells are unable to respond to S proteins.
C_LI


title: The Omicron Variant Increases the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2

abstract: SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. Currently, the most common Omicron variants are the original BA.1 Omicron strain and the BA.2 variant, which became more prevalent since it first appeared. To investigate how these mutations affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ~34% and ~51% stronger binding free energies for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.


title: Omicron BA.1 and BA.2 Variants Increase the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE2


title: Receptor binding and escape from Beta antibody responses drive Omicron-B.1.1.529 evolution

abstract: On the 24th November 2021 the sequence of a new SARS CoV-2 viral isolate spreading rapidly in Southern Africa was announced, containing far more mutations in Spike (S) than previously reported variants. Neutralization titres of Omicron by sera from vaccinees and convalescent subjects infected with early pandemic as well as Alpha, Beta, Gamma, Delta are substantially reduced or fail to neutralize. Titres against Omicron are boosted by third vaccine doses and are high in cases both vaccinated and infected by Delta. Mutations in Omicron knock out or substantially reduce neutralization by most of a large panel of potent monoclonal antibodies and antibodies under commercial development. Omicron S has structural changes from earlier viruses, combining mutations conferring tight binding to ACE2 to unleash evolution driven by immune escape, leading to a large number of mutations in the ACE2 binding site which rebalance receptor affinity to that of early pandemic viruses.


title: Omicron-B.1.1.529 leads to widespread escape from neutralizing antibody responses


title: A sarbecovirus found in Russian bats uses human ACE2

abstract:
Spillover of sarbecoviruses from animals to humans has resulted in outbreaks of severe acute respiratory syndrome SARS-CoVs and the ongoing COVID-19 pandemic. Efforts to identify the origins of SARS-CoV-1 and -2 has resulted in the discovery of numerous animal sarbecoviruses - the majority of which are only distantly related to known human pathogens and do not infect human cells. The receptor binding domain (RBD) on sarbecoviruses engages receptor molecules on the host cell and mediates cell invasion. Here, we tested the receptor tropism and serological cross reactivity for RBDs from two sarbecoviruses found in Russian horseshoe bats. While these two viruses are in a viral lineage distinct from SARS-CoV-1 and -2, one virus, Khosta-2, was capable of using human ACE2 to facilitate cell entry. Viral pseudotypes with a recombinant, SARS-CoV-2 spike encoding for the Khosta 2 RBD were resistant to both SARS-CoV-2 monoclonal antibodies and serum from individuals vaccinated for SARS-CoV-2. Our findings further demonstrate that sarbecoviruses circulating in wildlife outside of Asia also pose a threat to global health and ongoing vaccine campaigns against SARS-CoV-2

ONE SENTENCE SUMMARYEuropean bat coronaviruses that are only distantly related to SARS-CoV-2 but use the same cell entry route, escape the immune response against SARS-CoV-2 vaccines, driving the need for broader vaccines.


title: An ACE2-dependent Sarbecovirus in Russian bats is resistant to SARS-CoV-2 vaccines


title: An antibody-escape calculator for mutations to the SARS-CoV-2 receptor-binding domain

abstract: A key goal of SARS-CoV-2 surveillance is to rapidly identify viral variants with mutations that reduce neutralization by polyclonal antibodies elicited by vaccination or infection. Unfortunately, direct experimental characterization of new viral variants lags their sequence-based identification. Here we help address this challenge by aggregating deep mutational scanning data into an "escape calculator" that estimates the antigenic effects of arbitrary combinations of mutations to the viruss spike receptor-binding domain (RBD). The calculator can be used to intuitively visualize how mutations impact polyclonal antibody recognition, and score the expected antigenic effect of combinations of mutations. These scores correlate with neutralization assays performed on SARS-CoV-2 variants, and emphasize the ominous antigenic properties of the recently described Omicron variant. An interactive version of the calculator is at https://jbloomlab.github.io/SARS2_RBD_Ab_escape_maps/escape-calc/, and we provide a Python module for batch processing.


title: Decoupling SARS-CoV-2 ORF6 localization and interferon antagonism

abstract: Like many pathogenic viruses, SARS-CoV-2 must overcome interferon (IFN)-mediated host defenses for infection establishment. To achieve this, SARS-CoV-2 deploys overlapping mechanisms to antagonize IFN production and signaling. The strongest IFN antagonist is the accessory protein ORF6, which localizes to multiple membranous compartments, including the nuclear envelope, where it directly binds the nuclear pore components Nup98-Rae1 to inhibit nuclear translocation of activated STAT1/IRF3 transcription factors. However, a direct cause-and-effect relationship between ORF6 localization and IFN antagonism has yet to be explored experimentally. Here, we use extensive mutagenesis studies to define the structural determinants required for steady-state localization and demonstrate that mis-localized ORF6 variants can still potently inhibit nuclear trafficking and IFN signaling. Additionally, expression of a peptide that mimics the ORF6/Nup98 interaction domain robustly inhibited nuclear trafficking. Furthermore, pharmacologic and mutational approaches combined to suggest that ORF6 is likely a peripheral-membrane protein, opposed to being a transmembrane protein as previously speculated. Thus, ORF6 localization and IFN antagonism are independent activities, which raises the possibility that ORF6 may have additional functions within membrane networks to enhance virus replication.


title: Phosphorylation dynamics in a flg22-induced, G-protein dependent network in Arabidopsis thaliana

abstract: flg22, a 22-amino peptide released from bacterial flagellin, is a Microbe-Associated Molecular Pattern (MAMP) that is recognized by the plant cell as a signal indicating that bacteria are present. On its own, flg22 initiates a rapid increase in cytoplasmic calcium, extracellular reactive oxygen species, and activation of a Mitogen Activated Protein Kinase (MAPK) cascade, all of which are activated within 15 minutes after the cell perceives flg22. Here we show a massive change in protein abundance and phosphorylation state of the Arabidopsis root cell proteome within this 15-minute duration in wildtype and a mutant deficient in G-protein coupled signaling. Integration of phosphoproteome with protein-protein interactome data followed by network topology analyses discovered that many of the flg22-induced phosphoproteome changes fall on proteins that comprise the G protein interactome and on the most highly populated hubs of the immunity network. Approximately 95% of the phosphorylation changes in the G-protein interactome depend on a functional heterotrimeric G protein complex, some occur on proteins that interact directly with components of G-coupled signal transduction. One of these is ATB, a substrate-recognition sub-unit of the PP2A Ser/Thr phosphatase and an interactor to Arabidopsis thaliana REGULATOR OF G SIGNALING 1 protein (AtRGS1), a 7-transmembrane spanning modulator of the nucleotide-binding state of the core G protein complex. AtRGS1 is phosphorylated by BAK1, a component of the flg22 receptor, to initiate AtRGS1 endocytosis. A null mutation of ATB confers high basal endocytosis of AtRGS1, suggesting sustained phosphorylated status. Loss of ATB confers traits associated with loss of AtRGS1. Because the basal level of AtRGS1 is lower in the atb null mutant in a proteasome-dependent manner, we propose that phosphorylation-dependent endocytosis of AtRGS1 is part of the mechanism to degrade AtRGS1 which then sustains activation of the G protein complex required for phosphorylation-dependent regulation of system dynamics in innate immunity.


title: Elucidating design principles for engineering cell-derived vesicles to inhibit SARS-CoV-2 infection

abstract: The ability of pathogens to develop drug resistance is a global health challenge. The SARS-CoV-2 virus presents an urgent need wherein several variants of concern resist neutralization by monoclonal antibody therapies and vaccine-induced sera. Decoy nanoparticles--cell-mimicking particles that bind and inhibit virions--are an emerging class of therapeutics that may overcome such drug resistance challenges. To date, we lack quantitative understanding as to how design features impact performance of these therapeutics. To address this gap, here we perform a systematic, comparative evaluation of various biologically-derived nanoscale vesicles, which may be particularly well-suited to sustained or repeated administration in the clinic due to low toxicity, and investigate their potential to inhibit multiple classes of model SARS-CoV-2 virions. A key finding is that such particles exhibit potent antiviral efficacy across multiple manufacturing methods, vesicle subclasses, and virus-decoy binding affinities. In addition, these cell-mimicking vesicles effectively inhibit model SARS-CoV-2 variants that evade monoclonal antibodies and recombinant protein-based decoy inhibitors. This study provides a foundation of knowledge that may guide the design of decoy nanoparticle inhibitors for SARS-CoV-2 and other viral infections.


title: Insertions in the SARS-CoV-2 Spike N-Terminal Domain May Aid COVID-19 Transmission

abstract: Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 is an ongoing pandemic that causes significant health/socioeconomic burden. Variants of concern (VOCs) have emerged which may affect transmissibility, disease severity and re-infection risk. Most studies focus on the receptor-binding domain (RBD) of the Spike protein. However, some studies suggest that the Spike N-terminal domain (NTD) may have a role in facilitating virus entry via sialic-acid receptor binding. Furthermore, most VOCs include novel NTD variants. Recent analyses demonstrated that NTD insertions in VOCs tend to lie close to loop regions likely to be involved in binding sialic acids. We extended the structural characterisation of these putative sugar binding pockets and explored whether variants could enhance the binding to sialic acids and therefore to the host membrane, thereby contributing to increased transmissibility. We found that recent NTD insertions in VOCs (i.e., Gamma, Delta and Omicron variants) and emerging variants of interest (VOIs) (i.e., Iota, Lambda, Theta variants) frequently lie close to known and putative sugar-binding pockets. For some variants, including the recent Omicron VOC, we find increases in predicted sialic acid binding energy, compared to the original SARS-CoV-2, which may contribute to increased transmission. We examined the similarity of NTD across a range of related Betacoronaviruses to determine whether the putative sugar-binding pockets are sufficiently similar to be exploited in drug design. Despite global sequence and structure similarity, most sialic-acid binding pockets of NTD vary across related coronaviruses. Typically, SARS-CoV-2 possesses additional loops in these pockets that increase contact with polysaccharides. Our work suggests ongoing evolutionary tuning of the sugar-binding pockets in the virus. Whilst three of the pockets are too structurally variable to be amenable to pan Betacoronavirus drug design, we detected a fourth pocket that is highly structurally conserved and could therefore be investigated in pursuit of a generic drug. Our structure-based analyses help rationalise the effects of VOCs and provide hypotheses for experiments. For example, the Omicron variant, which has increased binding to sialic acids in pocket 3, has a rather unique insertion near pocket 3. Our work suggests a strong need for experimental monitoring of VOC changes in NTD.


title: Sensory nerves enhance triple-negative breast cancer migration and metastasis via the axon guidance molecule PlexinB3

abstract:
In breast cancer, nerve presence has been correlated with more invasive disease and worse prognosis, yet the mechanisms by which different types of peripheral nerves drive tumor progression remain poorly understood. In this study, we identified sensory nerves as more abundant in human triple-negative breast cancer (TNBC) tumors. Coinjection of sensory neurons isolated from the dorsal root ganglia (DRG) of adult female mice with human TNBC cells in immunocompromised mice increased the number of lung metastases. Direct in vitro co-culture of human TNBC cells with the dorsal root ganglia (DRG) of adult female mice revealed that TNBC cells adhere to sensory neuron fibers leading to an increase in migration speed. Species-specific RNA sequencing revealed that co-culture of TNBC cells with sensory nerves upregulates the expression of genes associated with cell migration and adhesion in cancer cells. We demonstrate that the axon guidance molecule Plexin B3 mediates cancer cell adhesion to and migration on sensory nerves. Together, our results identify a novel mechanism by which nerves contribute to breast cancer migration and metastasis by inducing a shift in TNBC cell gene expression and support the rationale for disrupting neuron-cancer cell interactions to target metastasis.

SignificanceThe presence of nerves in breast tumors has been associated with poor outcome. Understanding the mechanisms by which nerves contribute to tumor progression could help identify novel strategies to target metastatic disease.


title: Active communities and growth of soil microorganisms are framed by mean annual precipitation in three California annual grasslands

abstract: Earth system models project altered precipitation regimes across much of the globe; in California, the winter wet season is predicted to extend into spring, and the summer dry period to lengthen. How these precipitation trends will affect microbial traits and soil carbon (C) cycling is a key knowledge gap. Specifically, we do not have a mechanistic understanding of the linkages between soil moisture legacy effects, microbial population dynamics and soil C persistence. Using quantitative stable isotope probing (qSIP), we compared total and growing soil microbial communities across three California annual grasslands that span a rainfall gradient yet developed on similar parent material. We also assessed multiple edaphic variables, including the radiocarbon (14C) age of soil C, and found soil C turnover time increased with annual precipitation, but that soil microbes respired recently-fixed C regardless of site rainfall history. Samples were assayed in the wet season, when we expected environmental conditions would be most similar across sites. We hypothesized that growing communities would be more compositionally similar across the gradient than the total background microbiome. We also predicted that the long-term legacy effect of soil water limitation would be reflected in a lower community growth capacity at the driest site. We found that the proportion of the total community that was detected as growing was 28%, 48% and 58% at the wet, intermediate and dry sites, respectively. The composition of growing communities strongly resembled that of total communities, and growing communities were no more similar across the gradient than total communities, indicating a strong effect of climate of the structure of growing microbial communities. Members of three phyla, Acidobacteria, Actinobacteria, and Proteobacteria, were responsible for [~]79% of the cumulative 18O assimilation and 80% of all taxa that we defined as  growers. Bacterial growth rates were low at the driest site relative to the intermediate and wettest sites. Reduced growth at the driest site was observed across major phyla, including the Actinobacteria, Acidobacteria, Bacteroidetes, Gemmatimonadetes and Proteobacteria. Microbial communities at the driest site displayed phylogenetic clustering, suggesting that climate history impacts microbial growth through environmental filtering for slow growing taxa. Taxonomic identity was a strong predictor of growth, such that the growth rates of a taxon at one site predicted its growth rates at the others. This cross-site coherence in growth is likely a consequence of genetically determined physiological traits, and is consistent with the idea that evolutionary history influences growth rate.


title: Gap genes are involved in inviability in hybrids between Drosophila melanogaster and D. santomea

abstract: Evolved changes within species lead to the inevitable loss of viability in hybrids. Inviability is also a convenient phenotype to genetically map and validate functionally divergent genes and pathways differentiating closely related species. Here we identify the Drosophila melanogaster form of the highly conserved essential gap gene giant (gt) as a key genetic determinant of hybrid inviability in crosses with D. santomea. We show that the coding region of this allele in D. melanogaster/D. santomea hybrids is sufficient to cause embryonic inviability not seen in either pure species. Further genetic analysis indicates that tailless (tll), another gap gene, is also involved in the hybrid defects. giant and tll are both members of the gap gene network of transcription factors that participate in establishing anterior-posterior specification of the dipteran embryo, a highly conserved developmental process. Genes whose outputs in this process are functionally conserved nevertheless evolve over short timescales to cause inviability in hybrids.

category: evolutionary biology


title: Cohesin modulates DNA replication to preserve genome integrity

abstract: Cohesin participates in loop formation by extruding DNA fibers from its ring-shaped structure. Cohesin dysfunction eliminates chromatin loops but only causes modest transcription perturbation, which cannot fully explain the frequently observed mutations of cohesin in various cancers. Here, we found that DNA replication initiates at more than one thousand extra dormant origins after acute depletion of RAD21, a core subunit of cohesin, resulting in earlier replicating timing at approximately 30% of the human genomic regions. In contrast, CTCF is dispensable for suppressing the early firing of dormant origins that are distributed away from the loop boundaries. Furthermore, greatly elevated levels of gross DNA breaks and genome-wide chromosomal translocations arise in RAD21-depleted cells, accompanied by dysregulated replication timing at dozens of hotspot genes. Thus, we conclude that cohesin coordinates DNA replication initiation to ensure proper replication timing and safeguards genome integrity.

category: genetics


title: The genetic origin and architecture of sex determination

abstract: Here, I demonstrate that sex determination and sexual dimorphism across tree of life are deeply related to polyamine biochemistry in cells, especially to the synteny of genes: [SAT1-NR0B1], [SAT2-SHBG] and DMRT1. This synteny was found to be most distinct in mammals. Further, the common protein domain of SAT1 and SAT2 - PF00583 was shown to be present in the genome of the last universal common ancestor (LUCA). Protein domain-domain interaction analysis of LUCAs genes suggests the possibility that LUCA had developed an immune defence against viruses. This domain-domain interaction analysis is the first scientific evidence indicating that viruses existed at least 3.5 billions years ago and probably co-existed with LUCA on early Hadean Earth.


title: Mutational cascade of SARS-CoV-2 leading to evolution and emergence of omicron variant

abstract:
BackgroundEmergence of new variant of SARS-CoV-2, namely omicron, has posed a global concern because of its high rate of transmissibility and mutations in its genome. Researchers worldwide are trying to understand the evolution and emergence of such variants to understand the mutational cascade events.

MethodsWe have considered all omicron genomes (n = 302 genomes) available till 2nd December 2021 in the public repository of GISAID along with representatives of variants of concern (VOC), i.e., alpha, beta, gamma, delta, and omicron; variant of interest (VOI) mu and lambda; and variant under monitoring (VUM). Whole genome-based phylogeny and mutational analysis were performed to understand the evolution of SARS CoV-2 leading to emergence of omicron variant.

ResultsWhole genome-based phylogeny depicted two phylogroups (PG-I and PG-II) forming variant specific clades except for gamma and VUM GH. Mutational analysis detected 18,261 mutations in the omicron variant, majority of which were non-synonymous mutations in spike (A67, T547K, D614G, H655Y, N679K, P681H, D796Y, N856K, Q954H), followed by RNA dependent RNA polymerase (rdrp) (A1892T, I189V, P314L, K38R, T492I, V57V), ORF6 (M19M) and nucleocapsid protein (RG203KR).

ConclusionDelta and omicron have evolutionary diverged into distinct phylogroups and do not share a common ancestry. While, omicron shares common ancestry with VOI lambda and its evolution is mainly derived by the non-synonymous mutations.


title: Protein-protein interactions on membrane surfaces analysed using SUPER template pull-downs

abstract: Discovery-based proteomics workflows that identify novel interactors rely on immunoprecipitations or pull-downs with genetically-tagged bait proteins immobilized on appropriate matrices. But strategies to analyse protein interactions on a diffusible membrane surface combined with the practical ease of pull-downs remain unavailable. Such strategies are important to analyse protein complexes that mature in composition and stability because of diffusion-based encounter between participant proteins. Here, we describe a generic pull-down strategy to analyse such complexes using chelating lipid-containing supported bilayers formed on silica beads. These templates can display desired His-tagged bait proteins on a diffusible membrane surface. Using clathrin-mediated endocytosis as a paradigm, we find that the clathrin-binding adaptor protein epsin1 displayed as bait on these templates pulls down significantly higher amounts of clathrin from brain lysates than when immobilized on conventional matrices. Together, our results establish the potential of such templates as superior matrices for analysing protein-protein interactions and resultant complexes formed on membrane surfaces.


title: Protein-protein interactions on membrane surfaces analysed using pull-downs with supported bilayers on silica beads


title: LEXAS: a web application for life science experiment search and suggestion

abstract:
MotivationIn cellular biology, researchers design wet experiments by reading the relevant articles and considering the described experiments and results. Today, researchers spend a long time exploring the literature in order to plan experiments.

ResultsTo accelerate experiment planning, we have developed a web application named LEXAS (Life-science EXperiment seArch and Suggestion). LEXAS curates the description of biomedical experiments and suggests the experiments on genes that could be performed next. To develop LEXAS, we first retrieved the descriptions of experiments from full-text biomedical articles archived in PubMed Central. Using these retrieved experiments and biomedical knowledgebases and databases, we trained a machine learning model that suggests the next experiments. This model can suggest not only reasonable genes but also novel genes as targets for the next experiment as long as they share some critical features with the gene of interest.

Availability and implementationLEXAS is available at https://lexas.f.u-tokyo.ac.jp/ and provides users with two interfaces: search and suggestion. The search interface allows users to find a comprehensive list of experiment descriptions, and the suggestion interface allows users to find a list of genes that could be analyzed along with possible experiment methods. The source code is available at https://github.com/lexas-f-utokyo/lexas.

Contactito-delightfully-kei@g.ecc.u-tokyo.ac.jp

Supplementary informationSupplementary data are available at Bioinformatics online.


title: CLIPreg: Constructing translational regulatory networks from CLIP-, Ribo- and RNA-seq

abstract:
MotivationThe creation and analysis of gene regulatory networks have been the focus of bioinformatic research and underpins much of what is known about gene regulation. However, as a result of a bias in the availability of data-types that are collected, the vast majority of gene regulatory network resources and tools have focused on either transcriptional regulation or protein-protein interactions. This has left other areas of regulation, for instance translational regulation, vastly underrepresented despite them having been shown to play a critical role in both health and disease.

ResultsIn order to address this we have developed CLIPreg, a package that integrates RNA, Ribo and CLIP-sequencing data in order to construct translational regulatory networks coordinated by RNA-binding proteins. This is the first tool of its type to be created, allowing for detailed investigation into a previously unseen layer of regulation.

Availability and implementationCLIPreg is available at https://github.com/SGDDNB/CLIPreg.

Contactowen.rackham@duke-nus.edu.sg or sonia.chothani@u.duke.nus.edu


title: Revised International Staging System (R-ISS) stage-dependent analysis uncovers oncogenes and potential immunotherapeutic targets in multiple myeloma

abstract: Multiple myeloma (MM), characterized by high intratumour heterogeneity, accounts for [~]10% of all haematologic malignancies. Stratified by the Revised International Staging System (R-ISS), little is known about R-ISS-related plasma cell (PC) heterogeneity, gene expression modules in cytotoxic T/NK cells and immunoregulatory ligands and receptors. Herein, we constructed a single-cell transcriptome atlas of bone marrow in normal and R-ISS-staged MM patients. Focusing on PCs, we identified and validated a subset of GZMA+ cytotoxic PCs. In addition, a malignant PC population with high proliferation capability (proliferating PCs) was associated with unfavourable prognosis and EBV infection in our collected samples. Ribonucleotide Reductase Regulatory Subunit M2 (RRM2), a specific marker of proliferating PCs, was shown to induce MM cell line proliferation and serve as a detrimental marker in MM. Subsequently, three R-ISS-dependent gene modules in cytotoxic CD8+ T and NKT cells were identified and functionally analysed. Finally, cell-cell communication between neutrophils and proliferating PCs with cytotoxic CD8+ T and NKT cells was investigated, which identified intercellular ligand receptors and potential immunotargets such as SIRPA-CD47 and TIGIT-NECTIN3. Collectively, this study provides an R-ISS-related single-cell MM atlas and reveals the clinical significance of two PC clusters, as well as potential immunotargets in MM progression.


title: Assessment of exposure to Di (2-ethylhexyl) phthalate (DEHP) metabolites and Bisphenol A (BPA) and its importance for the prevention of cardiometabolic diseases

abstract: Exposomics analyses have highlighted the importance of biomonitoring of human exposure to pollutants, even non-persistent, for the prevention of non-communicable diseases like obesity, diabetes, non-alcoholic fatty liver disease, atherosclerosis and cardiovascular diseases. Phthalates and bisphenol A (BPA) are endocrine disrupting chemicals (EDCs) widely used in industry and in a large range of daily life products that increase the risk of endocrine and cardiometabolic diseases especially if the exposure starts during childhood. Thus, it is important the biomonitoring of exposure to these compounds not only in adulthood but also in childhood. This was the goal of the LIFE-PERSUADED project that measured the exposure to phthalates (DEHP metabolites, MEHP, MEHHP, MEOHP) and BPA in Italian mother-children couples of different ages. In this paper we describe the method that was set up for the LIFE PERSUADED project and validated during in the proficiency test (ICI/EQUAS) showing that accurate determination of urinary phthalates and BPA can be achieved starting from small sample size (0.5ml) using two MS techniques applied in cascade on the same deconjugated matrix.

category: pharmacology and toxicology


title: Vertical transmission of tissue microbiota in 	Caenorhabditis elegans

abstract: The past forty-five years has witnessed Caenorhabditis elegans as the most significant model animal in life science since its discovery seventy years ago1,2, as it introduced principles of gene regulated organ development, and RNA interference into biology3-5. Meanwhile, it has become one of the lab animals in gut microbiota studies as these symbionts contribute significantly to many aspects in host biology6,7. Meanwhile, the origin of gut microbiota remains debatable in human8- 11, and has not been investigated in other model animals. Here we show that the symbiont bacteria in C. elegans not only vertically transmit from the parent generation to the next, but also distributes in the worm tissues parallel with its development. We found that bacteria can enter into the embryos of C. elegans, a step associated with vitellogenin, and passed to the next generation. These vertically transmitted bacteria share global similarity, and bacterial distribution in worm tissues changes as they grow at different life stages. Antibiotic treatment of worms increased their vulnerability against pathogenic bacteria, and replenishment of tissue microbiota restored their immunity. These results not only offered a molecular basis of vertical transmission of bacteria in C. elegans, but also signal a new era for the mixed tissue cell-bacteria multi-species organism study.

category: physiology


title: Neuron-like function of the nephron central command

abstract: Interoceptive neurons that sense and regulate our internal milieu have been identified in several organs except in the kidney cortex despite its major importance in maintaining body homeostasis. Here we report that the chief kidney cell type of the macula densa (MD) forms coordinated neural networks in each nephron that resemble peripheral ganglia. A combined in vivo single-cell 4D physiology (sc4DP) and scRNA sequencing approach identified the MD mechanisms of neuronal differentiation, heterogeneity (pacemaker MD cells), sensing of the local and systemic environment via multi-organ crosstalk, and regulation of organ functions by acting as the nephron central command. Consistent with their neuron-like nature, MD cells express the molecular fingerprint of neurodegeneration. Here we put forth the single-cell MD model and concept of local neural networks that control organ and body functions via interoception in normal physiological state and use an integrated mechanism of neurodegeneration in disease.


title: The oscillation of mitotic kinase governs cell cycle latches

abstract: In order to transmit a eukaryotic cells genome accurately from mother cell to daughter cells, it is essential that the basic events of the cell division cycle (DNA synthesis and mitosis) occur once and only once per cycle, i.e., that a cell progresses irreversibly from G1 to S to G2 to M and back to G1. Irreversible progression through the cell cycle is assured by a sequence of  latching molecular switches, based on molecular interactions among cyclin-dependent kinases and their auxiliary partners. Positive feedback loops (++ or --) create bistable switches with latching properties, and negative feedback loops drive progression from one stage to the next. In budding yeast (Saccharomyces cerevisiae) these events are coordinated by double-negative feedback loops between Clb-dependent kinases (Clb1-6) and their antagonists (APC:Cdh1 and Sic1). If the coordinating signal is compromised, either by deletion of Clb1-5 proteins or expression of non-degradable Clb2, then irreversibility is lost and yeast cells exhibit multiple rounds of DNA replication or mitotic exit events (Cdc14 endocycles). Using mathematical modelling of a stripped-down control network, we show how endocycles arise because the switches fail to latch, and the gates swing back and forth by the action of the negative feedback loops.


title: Quantitative and temporal measurement of dynamic autophagy rates

abstract: Autophagy is a multistep degradative process that is essential for maintaining cellular homeostasis. Systematically quantifying flux through this pathway is critical for gaining fundamental insights and effectively modulating this process that is dysregulated during many diseases. Established methods to quantify flux use steady state measurements, which provide limited information about the perturbation and the cellular response. We present a theoretical and experimental framework to measure autophagic steps in the form of rates under non-steady state conditions. We use this approach to measure temporal responses to rapamycin and wortmannin treatments, two commonly used autophagy modulators. We quantified changes in autophagy rates in as little as 10 minutes, which can establish direct mechanisms for autophagy perturbation before feedback begins. We identified concentration-dependent effects of rapamycin on the initial and temporal progression of autophagy rates. We also found variable recovery time from wortmannins inhibition of autophagy, which is further accelerated by rapamycin. In summary, this new approach enables the quantification of autophagy flux with high sensitivity and temporal resolution and facilitates a comprehensive understanding of this process.


title: AP-1 transcription factor network explains diverse patterns of cellular plasticity in melanoma

abstract: Cellular plasticity associated with fluctuations in transcriptional programs allows individual cells in a tumor to adopt heterogeneous differentiation states and switch phenotype during their adaptive responses to therapies. Despite increasing knowledge of such transcriptional programs, the molecular basis of cellular plasticity remains poorly understood. Here, we combine multiplexed transcriptional and protein measurements at population and single-cell levels with multivariate statistical modeling to show that the state of AP-1 transcription factor network plays a unifying role in explaining cellular plasticity in melanoma. We find that a tightly regulated balance between AP-1 factors cJUN, FRA2, FRA1 and cFOS determines the intrinsic diversity of differentiation states and adaptive responses to MAPK inhibitors in individual melanoma cells. Perturbing this balance through genetic depletion of specific AP-1 proteins, or by MAPK inhibitors, shifts cellular heterogeneity in a predictable fashion. Thus, AP-1 may serve as a critical node for manipulating cellular plasticity with potential therapeutic implications.


title: Modulation of certainty and valence in cardiac autonomic regulation mediated by high depressive risk

abstract: Recent studies highlight interoception as essential to processing healthy emotion and to pathology of major depressive disorder (MDD). It has been reported that interoceptive responses are impaired in MDD and in healthy individuals with a high depressive risk (HDR). However, it is unclear how individual differences in HDR relate to neurophysiological underpinnings for interoceptive and emotional reactions under different degrees of certainty. We examined whether HDR mediates the relationship between a neurophysiological marker for interoception, heartbeat-evoked potential (HEP), and an index for cardiac reactivity, heart rate (HR) or heart rate variability (HRV). In a concurrent EEG-ECG experiment, 26 healthy participants underwent an emotion-evoking picture-evaluation task. Each of differential auditory cues was associated with certainty of the upcoming pleasant or unpleasant picture. The results showed attenuated HRV activity for certain cues in HDR individuals. Neural and physiological reactions to uncertain, unpleasant pictures were enhanced by depressive risk. These results suggest reduced responses in prediction processing and enlarged precision-weighted prediction error for unpleasant pictures. Finally, HDR significantly mediated the HEP and HR relationship for unexpected unpleasant stimuli. This study provides evidence that interoceptive predictive coding in the cardiac domain is altered by HDR and suggests that investigating brain-heart interactions related to predictive coding will offer insights into how to estimate depressive risk.


title: Brain-heart dynamics in emotional regulation under uncertainty indexes depressive risk


title: Analysis of visual processing capabilities and neural coding strategies of a detailed model for laminar cortical microcircuits in mouse V1

abstract: Neural networks of the brain that process visual information have structural properties that differ significantly from those of neural networks which are commonly used for visual processing in AI, such as Convolutional Neural Networks (CNNs). But it has remained unknown how these structural differences are related to network function. We analyze visual processing capabilities of a large-scale model for area V1 that arguably provides the most comprehensive accumulation of anatomical and neurophysiological data that is currently available. Its network structure turns out to induce a number of characteristic visual processing capabilities of the brain, in particular the capability to multiplex different visual processing tasks, also on temporally dispersed visual information, with remarkable robustness to noise. This V1 model also exhibits a number of characteristic neural coding properties of the brain, which provide explanations for its superior noise robustness. Since visual processing in the brain is substantially more energy-efficient than implementations of CNNs in common computer hardware, such brain-like neural network models are likely to have also an impact on technology: As blueprints for visual processing in more energy-efficient neuromorphic hardware.


title: Impaired LAIR-1-mediated immune control due to collagen degradation in fibrosis

abstract:
Tissue repair is disturbed in fibrotic diseases like systemic sclerosis (SSc), where the deposition of large amounts of extracellular matrix components such as collagen interferes with organ function. LAIR-1 is an inhibitory collagen receptor highly expressed on tissue immune cells. We questioned whether in SSc, impaired LAIR-1-collagen interaction is contributing to the ongoing inflammation and fibrosis.

We found that SSc patients do not have an intrinsic defect in LAIR-1 expression or function. Instead, fibroblasts from SSc patients deposit disorganized collagen products in vitro, which are dysfunctional LAIR-1 ligands. This can be mimicked in healthy fibroblast stimulated by soluble factors that drive inflammation and fibrosis in SSc and is dependent of matrix metalloproteinases and platelet-derived growth factor receptor signaling.

In support of a non-redundant role of LAIR-1 in the control of fibrosis, we found that LAIR-1-deficient mice have increased skin fibrosis in the bleomycin mouse model for SSc. Thus, LAIR-1 represents an essential control mechanism for tissue repair. In fibrotic disease, excessive collagen degradation may lead to a disturbed feedback loop. The presence of functional LAIR-1 in patients provides a therapeutic opportunity to reactivate this intrinsic negative feedback mechanism in fibrotic diseases.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=93 SRC="FIGDIR/small/471103v1_ufig1.gif" ALT="Figure 1">
View larger version (19K):
org.highwire.dtl.DTLVardef@9590beorg.highwire.dtl.DTLVardef@151c1d6org.highwire.dtl.DTLVardef@8cf877org.highwire.dtl.DTLVardef@10c5e59_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Distinct within-host bacterial populations ensure function, colonization and transmission in leaf symbiosis

abstract:
Hereditary symbioses have the potential to drive transgenerational effects, yet the mechanisms responsible for transmission of heritable plant symbionts are still poorly understood. The leaf symbiosis between Dioscorea sansibarensis and the bacterium Orrella dioscoreae offers an appealing model system to study how heritable bacteria are transmitted to the next generation. Here, we demonstrate that inoculation of apical buds with a bacterial suspension is sufficient to colonize newly-formed leaves and propagules, and to ensure transmission to the next plant generation. Flagellar motility is not required for movement inside the plant, but is important for the colonization of new hosts. Further, stringent tissue-specific regulation of putative symbiotic functions highlight the presence of two distinct subpopulations of bacteria in the leaf gland and at the shoot meristem. We propose that bacteria in the leaf gland dedicate resources to symbiotic functions, while dividing bacteria in the shoot tip ensure successful colonization of meristematic tissue, glands and propagules. Compartmentalization of intra-host populations, together with tissue-specific regulation may serve as a robust mechanism for the maintenance of mutualism in leaf symbiosis.

ImportanceSeveral plant species form associations with bacteria in their leaves, called leaf symbiosis. These associations are highly specific, but the mechanisms responsible for symbiont transmission are poorly understood. Using the association between the yam species Dioscorea sansibarensis and Orrella dioscoreae as a model leaf symbiosis, we provide experimental evidence that bacteria are transmitted vertically and distributed to specific leaf structures via association with shoot meristems. Flagellar motility is required for initial infection, but does not contribute to spread within host tissue. We also provide evidence that bacterial subpopulations at the meristem or in the symbiotic leaf gland differentially express key symbiotic genes. We argue that this separation of functional symbiont populations, coupled to tight control over bacterial infection and transmission, explain the evolutionary robustness of leaf symbiosis. These findings may provide insights into how plants may recruit and maintain beneficial symbionts at the leaf surface.


title: Meet me in the middle: median temperatures impact cyanobacteria and photoautotrophy in eruptive Yellowstone hot springs

abstract:
Geographic isolation can be a main driver of microbial evolution in hot springs while temperature plays a role on local scales. For example, cyanobacteria, particularly high temperature Synechococcus spp., have undergone ecological diversification along temperature gradients in hot spring outflow channels. While water flow, and thus temperature, is largely stable in many hot springs, flow can vary in geysing/eruptive hot springs resulting in large temperature fluctuations (sometimes more than 40{degrees}C). However, the role of large temperature fluctuations in driving diversification of cyanobacteria in eruptive hot springs has not been explored. Here, we examined phototroph community composition and potential photoautotrophic activity in two alkaline eruptive hot springs with similar geochemistry in the Lower Geyser Basin in Yellowstone National Park, WY. We observed distinct cyanobacterial amplicon sequencing variants (ASVs) consistent with allopatry and levels of light-dependent inorganic carbon uptake rates similar to other hot springs, despite large temperature fluctuations. Our data suggests median temperatures may drive phototroph fitness in eruptive hot springs while future studies are necessary to determine the evolutionary consequences of thriving under continuously fluctuating temperatures. We propose that large temperature swings in eruptive hot springs offer unique environments to examine the role of allopatry vs. physical and chemical characteristics of ecosystems in driving cyanobacteria evolution and add to the debate regarding the ecology of thermal adaptation and the potential for narrowing niche breadth with increasing temperature.

ImportanceHot spring cyanobacteria have long been model systems for examining ecological diversification as well as characterizing microbial adaptation and evolution to extreme environments. These studies have reported cyanobacterial diversification in hot spring outflow channels that can be defined by distinct temperature ranges. Our study builds on these previous studies by examining cyanobacteria in geysing hot springs. Geysing hot springs result in outflow channel that experience regular and large temperature fluctuations. While community composition is similar between geysing and nongeysing hot spring outflow channels, our data suggests median, rather than high temperature, drive the fitness of cyanobacteria in geysing hot springs. We propose that large temperature swings may result in patterns of ecological diversification that are distinct from more stable outflows.


title: Stem cell specific interferon stimulated gene expression is regulated by the formative pluripotency network through IRF1

abstract: Stem cells intrinsically express a subset of genes which are normally associated with interferon stimulation, and thus the innate immunity response. Expression of these interferon stimulated genes (ISG) in stem cells is independent of external stimuli such as viral infection. Here we show that the interferon regulatory factor 1, Irf1, is directly controlled by the murine formative pluripotency gene regulatory network and therefore upregulated in the transition from naive to formative pluripotency. IRF1 in turn binds to regulatory regions of a conserved set of ISGs and is required for their faithful expression in formative pluripotent cells. IRF1 also binds to an enhancer of the formative pluripotency transcription factor Oct6 and is partially required for upregulation of Oct6. IRF1 therefore acts as a link between the formative pluripotency network and the regulation of innate immunity genes in formative pluripotency.

category: developmental biology


title: Adventive larval parasitoids reconstruct their close association with spotted-wing drosophila in the invaded North American range

abstract: Two species of larval parasitoids of the globally invasive fruit pest, Drosophila suzukii (Diptera: Drosophilidae), Leptopilina japonica and Ganaspis brasiliensis (both Hymenoptera: Figitidae), were detected in British Columbia, Canada in 2019. Both are presumed to have been unintentionally introduced from Asia, however; the extent of their establishment across different habitats with diverse host plants used by D. suzukii was unclear. In addition, there was no knowledge of the temporal dynamics of parasitism of D. suzukii by these two parasitoids. We repeatedly sampled the fruits of known host plants of D. suzukii over the entire 2020 growing season in British Columbia. We documented the presence of L. japonica and G. brasiliensis and estimated the apparent percentage of D. suzukii parasitized. Across a large region of southwestern British Columbia, both L. japonica and G. brasiliensis were found to be very common across a variety of mostly unmanaged habitats over the entire course of the season (May-October) in the fruits of most host plants known to host D. suzukii larvae. The two parasitoids were responsible for more than 98% of D. suzukii larval parasitism and usually co-existed. Parasitism of D. suzukii was variable among hosts plants and sites (0-66% percent parasitism) and appeared to be time-structured. Our study demonstrates that the close association between the two larval parasitoids and D. suzukii that exists in Asia has evidently been reconstructed in North America, resulting in the highest parasitism levels of D. suzukii yet recorded outside of its area of origin.


title: A single nucleotide variant in the PPARgamma-homolog Eip75B affects fecundity in Drosophila

abstract: Still little is understood about the nucleotide changes that underlie variation in complex phenotypes. Variation in the PPAR{gamma}-homolog Eip75B has previously been suggested to be associated with longevity and life-history differences in the fruit fly Drosophila melanogaster. Using RNAi knockdown, we first demonstrate that reduced expression of Eip75B in adults affects lifespan, egg laying rate and egg volume. To then test the effect of a naturally occurring SNP variant within a cis-regulatory domain of Eip75B, we screened wildtype lines with alternative alleles and conducted precise genome editing using CRISPR/Cas9. These experiments revealed that this natural polymorphism has a significant effect on fecundity and egg-to-adult viability, but not on longevity or other life-history traits. These results provide a rare functional validation for the role of a natural allelic variant in adaptation of life-history traits directly linked to fitness at the single nucleotide level.