title: An ancient viral epidemic involving host coronavirus interacting genes more than 20,000 years ago in East Asia

abstract: The current SARS-CoV-2 pandemic has emphasized the vulnerability of human populations to novel viral pressures, despite the vast array of epidemiological and biomedical tools now available. Notably, modern human genomes contain evolutionary information tracing back tens of thousands of years, which may help identify the viruses that have impacted our ancestors - pointing to which viruses have future pandemic potential. Here, we apply evolutionary analyses to human genomic datasets to recover selection events involving tens of human genes that interact with coronaviruses, including SARS-CoV-2, that likely started more than 20,000 years ago. These adaptive events were limited to the population ancestral to East Asian populations. Multiple lines of functional evidence support an ancient viral selective pressure, and East Asia is the geographical origin of several modern coronavirus epidemics. An arms race with an ancient coronavirus, or with a different virus that happened to use similar interactions as coronaviruses with human hosts, may thus have taken place in ancestral East Asian populations. By learning more about our ancient viral foes, our study highlights the promise of evolutionary information to better predict the pandemics of the future. Importantly, adaptation to ancient viral epidemics in specific human populations does not necessarily imply any difference in genetic susceptibility between different human populations, and the current evidence points toward an overwhelming impact of socioeconomic factors in the case of COVID-19.


title: Distinct and Common Neural Coding of Semantic and Non-semantic Control Demands

abstract: The flexible retrieval of knowledge is critical in everyday situations involving problem solving, reasoning and social interaction. Current theories emphasise the importance of a left-lateralised semantic control network (SCN) in supporting flexible semantic behaviour, while a bilateral multiple-demand network (MDN) is implicated in executive functions across domains. No study, however, has examined whether semantic and non-semantic demands are reflected in a common neural code within regions specifically implicated in semantic control. Using functional MRI and univariate parametric modulation analysis as well as multivariate pattern analysis, we found that semantic and non-semantic demands gave rise to both similar and distinct neural responses across control-related networks. Though activity patterns in SCN and MDN could decode the difficulty of both semantic and verbal working memory decisions, there was no shared common neural coding of cognitive demands in SCN regions. In contrast, regions in MDN showed common patterns across manipulations of semantic and working memory control demands, with successful cross-classification of difficulty across tasks. Therefore, SCN and MDN can be dissociated according to the information they maintain about cognitive demands.

category: neuroscience


title: Revisiting the notion of deleterious sweeps

abstract: It has previously been shown that, conditional on its fixation, the time to fixation of a semi-dominant deleterious autosomal mutation in a randomly mating population is the same as that of an advantageous mutation. This result implies that deleterious mutations could generate selective sweep-like effects. Although their fixation probabilities greatly differ, the much larger input of deleterious relative to beneficial mutations suggests that this phenomenon could be important. We here examine how the fixation of mildly deleterious mutations affects levels and patterns of polymorphism at linked sites - both in the presence and absence of interference amongst deleterious mutations - and how this class of sites may contribute to divergence between-populations and species. We find that, while deleterious fixations are unlikely to represent a significant proportion of outliers in polymorphism-based genomic scans within populations, minor shifts in the frequencies of deleterious mutations can influence the proportions of private variants and the value of FST after a recent population split. As sites subject to deleterious mutations are necessarily found in functional genomic regions, interpretations in terms of recurrent positive selection may require reconsideration.

category: evolutionary biology


title: Plant PhysioSpace: a robust tool to compare stress response across plant species

abstract:
Generalization of transcriptomics results can be achieved by comparison across experiments, which is based on integration of interrelated transcriptomics studies into a compendium. Both characterization of the fate of the organism as well as distinguishing between generic and specific responses can be gained in such a broader context. There are numerous methods for analyzing such data sets, most focusing on gene-wise dimension reduction to obtain marker genes and gene sets, e.g. for pathway analysis. Relying only on isolated biological modules might lead to missing of important confounders and relevant context.

We have developed a novel method called Plant PhysioSpace, which provides the ability to compute experimental conditions across species and platforms without a priori reducing the reference information to specific gene-sets. It extracts physiologically relevant signatures from a reference data set, a collection of public data sets, by integrating and transforming heterogeneous reference gene expression data into a set of physiology-specific patterns. New experimental data can be mapped to these patterns, resulting in similarity scores which provide quantitative likeness of the new experiment to the a priori compendium.

Because of its robustness against noise and platform bias, Plant PhysioSpace can be used as an inter-species or cross-platform similarity measure. We have demonstrated its success in translating stress responses between different species and platforms (including single-cell technologies).

We also report the implementation of two R packages, one software and one data package, and a shiny web application, which provides plant biologists convenient ways to access our method and precomputed models from four different species.

category: plant biology


title: Bridging brain and cognition: A multilayer network analysis of brain structural covariance and general intelligence in a developmental sample of struggling learners

abstract: Network analytic methods that are ubiquitous in other areas, such as systems neuroscience, have recently been used to test network theories in psychology, including intelligence research. The network or mutualism theory of intelligence proposes that the statistical associations among cognitive abilities (e.g., specific abilities such as vocabulary or memory) stem from causal relations among them throughout development. In this study, we used network models (specifically LASSO) of cognitive abilities and brain structural covariance (grey and white matter) to simultaneously model brain-behavior relationships essential for general intelligence in a large (behavioral, N=805; cortical volume, N=246; fractional anisotropy, N=165), developmental (ages 5-18) cohort of struggling learners (CALM). We found that mostly positive, small partial correlations pervade our cognitive, neural, and multilayer networks. Moreover, using community detection (Walktrap algorithm) and calculating node centrality (absolute strength and bridge strength), we found convergent evidence that subsets of both cognitive and neural nodes play an intermediary role  between brain and behavior. We discuss implications and possible avenues for future studies.


title: Role of BMP signaling during early development of the annelid Capitella teleta

abstract: The mechanisms regulating nervous system development are still unknown for a wide variety of taxa. In insects and vertebrates, bone morphogenetic protein (BMP) signaling plays a key role in establishing the dorsal-ventral (D-V) axis and limiting the neuroectoderm to one side of that axis, leading to speculation about the conserved evolution of centralized nervous systems. Studies outside of insects and vertebrates show a more diverse picture of what, if any role, BMP signaling plays in neural development across Bilateria. This is especially true in the morphologically diverse Spiralia (~Lophotrochozoa). Despite several studies of D-V axis formation and neural induction in spiralians, there is no consensus for how these two processes are related, or whether BMP signaling may have played an ancestral role in either process. To determine the function of BMP signaling during early development of the spiralian annelid Capitella teleta, we incubated embryos and larvae in BMP4 protein for different amounts of time. Adding exogenous BMP protein to early-cleaving C. teleta embryos had a striking effect on formation of the brain, eyes, foregut, and ventral midline in a time-dependent manner. However, adding BMP did not block brain or VNC formation or majorly disrupt the D-V axis. We identified three key time windows of BMP activity. 1) BMP treatment around birth of the 3rd-quartet micromeres caused the loss of the eyes, radialization of the brain, and a reduction of the foregut, which we interpret as a loss of A- and C-quadrant identities with a possible trans-fate switch to a D-quadrant identity. 2) Treatment after birth of micromere 4d induced formation of a third ectopic brain lobe, eye, and foregut lobe, which we interpret as a trans-fate switch of B-quadrant micromeres to a C-quadrant identity. 3) Continuous BMP treatment from late cleavage (4d + 12h) through mid-larval stages resulted in a modest expansion of Ct-chrdl expression in the dorsal ectoderm and a concomitant loss of the ventral midline (neurotroch ciliary band). Loss of the ventral midline was accompanied by a collapse of the bilaterally-symmetric ventral nerve cord, although the total amount of neural tissue did not appear to be greatly affected. Our results compared to those from other annelids and molluscs suggest that BMP signaling was not ancestrally involved in delimiting neural tissue to one region of the D-V axis. However, the effects of ectopic BMP on quadrant-identity during cleavage stages may represent a non-axial organizing signal that was present in the last common ancestor of annelids and mollusks. Furthermore, in the last common ancestor of annelids, BMP signaling may have functioned in patterning ectodermal fates along the D-V axis in the trunk. Ultimately, studies on a wider range of spiralian taxa are needed to determine the role of BMP signaling during neural induction and neural patterning in the last common ancestor of this group. Ultimately, these comparisons will give us insight into the evolutionary origins of centralized nervous systems and body plans.

category: developmental biology


title: Dynamics of vesicle reshaping and scission under osmotic shocks

abstract: We study the effects of osmotic shocks on lipid vesicles via coarse-grained molecular dynamics simulations by explicitly considering the solute in the system. We find that depending on their nature (hypo- or hypertonic) such shocks can lead to bursting events or engulfing of external material into inner compartments, among other morphology transformations. We characterize the dynamics of these processes and observe a separation of time scales between the osmotic shock absorption and the shape relaxation. Our work consequently provides an insight into the dynamics of compartmentalization in vesicular systems as a result of osmotic shocks, which can be of interest in the context of early proto-cell development and proto-cell compartmentalisation.

category: biophysics


title: N2O-respiring bacteria in biogas digestates for reduced agricultural emissions

abstract: Inoculating agricultural soils with N2O-respiring bacteria (NRB) can reduce N2O-emissions, but would be impractical as a standalone operation. Here we demonstrate that digestates obtained after biogas production are suitable substrates and vectors for NRB. We show that indigenous NRB in digestates grew to high abundance during anaerobic enrichment under N2O. Gas-kinetics and meta-omic analyses showed that these NRB's, recovered as metagenome-assembled genomes (MAGs), grew by harvesting fermentation intermediates of the methanogenic consortium. Three NRB's were isolated, one of which matched the recovered MAG of a Dechloromonas, deemed by proteomics to be the dominant producer of N2O-reductase in the enrichment. While the isolates harbored genes required for a full denitrification pathway and could thus both produce and sequester N2O, their regulatory traits predicted that they act as N2O sinks in soil, which was confirmed experimentally. The isolates were grown by aerobic respiration in digestates, and fertilization with these NRB-enriched digestates reduced N2O emissions from soil. Our use of digestates for low-cost and large-scale inoculation with NRB in soil can be taken as a blueprint for future applications of this powerful instrument to engineer the soil microbiome, be it for enhancing plant growth, bioremediation, or any other desirable function.

category: microbiology


title: Population genetics of transposable element load: a mechanistic account of observed overdispersion

abstract: In an empirical analysis of transposable element (TE) abundance within natural populations of Mimulus guttatus and Drosophila melanogaster, we found a surprisingly high variance of TE count (e.g., variance-to-mean ratio on the order of 10 to 300). To obtain insight regarding those evolutionary genetic mechanisms that are may underlie the overdispersed population distributions of TE abundance, we developed a mathematical model of TE population genetics that includes the dynamics of element proliferation and purifying selection on TE load. The modeling approach begins with a master equation for a birth-death process and extends the predictions of the classical theory of TE dynamics in several ways. In particular, moment-based analyses of time-dependent and equilibrium population distributions of TE load reveal that overdispersion is likely to arise via copy-and-paste proliferation dynamics, especially when the elementary processes of proliferation and excision are first-order and approximately balanced. Parameter studies and analytic work confirm this result and further suggest that overdispersed population distributions of TE abundance are unlikely to be a consequence of purifying selection on total element load.

category: genomics


title: Celda: A Bayesian model to perform co-clustering of genes into modules and cells into subpopulations using single-cell RNA-seq data

abstract: Single-cell RNA-seq (scRNA-seq) has emerged as a powerful technique to quantify gene expression in individual cells and elucidate the molecular and cellular building blocks of complex tissues. We developed a novel Bayesian hierarchical model called Cellular Latent Dirichlet Allocation (Celda) to perform simultaneous co-clustering of genes into transcriptional modules and cells into subpopulations. Celda can quantify the probabilistic contribution of each gene to each module, each module to each cell population, and each cell population to each sample. We used Celda to identify transcriptional modules and cell subpopulations in a publicly available peripheral blood mononuclear cell (PBMC) dataset. Celda identified a population of proliferating T cells and a single plasma cell which were missed by two other clustering methods. Celda identified transcriptional modules that highlighted unique and shared biological programs across cell types. Celda also outperformed a PCA-based approach for gene clustering on simulated data. Overall, Celda presents a novel statistically principled approach towards characterizing transcriptional programs and cellular heterogeneity in single-cell RNA-seq data.

category: bioinformatics


title: Tysserand - Fast and accurate reconstruction of spatial networks from bioimages

abstract:
SummaryNetworks provide a powerful framework to analyze spatial omics experiments. However, we lack tools that integrate several methods to easily reconstruct networks for further analyses with dedicated libraries. In addition, choosing the appropriate method and parameters can be challenging.

We propose tysserand, a Python library to reconstruct spatial networks from spatially resolved omics experiments. It is intended as a common tool to which the bioinformatics community can add new methods to reconstruct networks, choose appropriate parameters, clean resulting networks and pipe data to other libraries.

Availability and implementationtysserand software and tutorials with a Jupyter notebook to reproduce the results are available at https://github.com/VeraPancaldiLab/tysserand

Supplementary informationSupplementary data are available at Bioarxiv online.


title: The gut bacteria of an invasive amphibian respond to the dual challenges of range-expansion and parasite attack

abstract: Gut bacterial communities influence, and are influenced by, the behaviour and ecology of their hosts. Those interactions have been studied primarily in humans and model organisms, but we need field research to understand the relationship between an organisms gut bacteria and its ecological challenges, such as those imposed by rapid range expansion (as in invasive species) and the presence of host-manipulating parasites. Cane toads (Rhinella marina) provide an excellent model system in this respect, because the species ongoing colonization of Australia has enforced major changes in phenotypic traits (including behaviour), and lungworm parasites (Rhabdias pseudosphaerocephala) modify host gut function in ways that enhance the viability of lungworm larvae. We collected female toads from across the species invasive range and studied their morphology, behaviour, parasite infection status and gut bacterial community. Range-core versus range-edge toads differed in morphology, behaviour, gut bacterial composition and predicted gut bacterial function but did not differ in the occurrence of parasite infection nor in the intensity of infection. Toads infected with lungworms differed from uninfected conspecifics in gut bacterial composition and diversity. Our study demonstrates strong associations between gut bacterial community and host ecology and behaviour.


title: Comprehensive generation, visualization, and reporting of quality control metrics for single-cell RNA sequencing data

abstract: Single-cell RNA sequencing (scRNA-seq) can be used to gain insights into cellular heterogeneity within complex tissues. However, a variety of technical artifacts can be present in scRNA-seq data and need to be assessed before downstream analyses can be performed. While several algorithms and tools have been developed to perform individual quality control (QC) tasks, they are scattered in different packages across several programming environments. Comprehensive pipelines to streamline the process of generating and visualizing QC metrics are lacking. To address this need, we built the SCTK-QC pipeline within the singleCellTK R package (https://github.com/compbiomed/singleCellTK). Features in this pipeline include the ability to import data from 11 different preprocessing tools or file formats, perform empty droplet detection with 2 different algorithms, generate standard quality control metrics such as number of UMIs per cell or the percentage of mitochondrial counts, predict doublets using 6 different algorithms, and estimate ambient RNA. QC data can be exported to R and/or Python objects used in popular down-stream workflows. Results are visualized in an easy-to-read HTML report. This pipeline can also be used by non-computational users with an interactive graphical user interface developed with R/Shiny. Overall, the SCTK-QC pipeline will streamline and standardize QC analysis for scRNA-seq data across a variety of different single-cell transcriptomic platforms and preprocessing tools.


title: Perphenazine-macrocycle conjugates rapidly sequester the Aβ42 monomer and inhibit amyloid formation

abstract: Alzheimers disease is imposing a growing social and economic burden worldwide and effective therapies are required. Strategies aimed at the removal of fibrillar plaques formed by the amyloid-{beta} peptide have not proved therapeutic and the focus has shifted to approaches that target the cytotoxic oligomeric amyloid-{beta} species that are populated before fibrils are deposited. We have designed and synthesized perphenazine-cyclam conjugates that specifically and rapidly bind to the monomeric form of A{beta}42, reducing the production of both cytotoxic oligomers and amyloid fibrils. We have applied detailed kinetic analysis and NMR spectroscopy to show that the perphenazine-cyclam conjugates divert the A{beta}42 monomer into amorphous aggregates that are not toxic to differentiated SH-SY5Y cells in vitro. Unlike most other amyloid inhibitors studied to date, these conjugates inhibit oligomer and fibril assembly even in the presence of pre-formed fibrillar seeds, demonstrating that they act through a monomer sequestration mechanism. These modular, three-dimensional conjugates therefore effectively prevent monomer-dependent secondary nucleation, the autocatalytic process that generates the majority of toxic oligomers.


title: Community Detection in a Weighted Directed Hypergraph Representation of Cell-to-cell Communication Networks

abstract: Cell-to-cell communication is mainly triggered by ligand-receptor activities. Through ligandreceptor pairs, cells coordinate complex processes such as development, homeostasis, and immune response. In this work, we model the ligand-receptor-mediated cell-to-cell communication network as a weighted directed hypergraph. In this mathematical model, collaborating cell types are considered as a node community while the ligand-receptor pairs connecting them are considered a hyperedge community. We first define the community structures in a weighted directed hypergraph and develop an exact community detection method to identify these communities. We then modify approximate community detection algorithms designed for simple graphs to identify the nodes and hyperedges within each community. Application to synthetic hypergraphs with known community structure confirmed that one of the proposed approximate community identification strategies, named HyperCommunity algorithm, can effectively and precisely detect embedded communities. We then applied this strategy to two organism-wide datasets and identified putative community structures. Notably the method identifies non-overlapping edge-communities mediated by different sets of ligand-receptor pairs, however node-communities can overlap.

category: systems biology


title: Mosaic nanoparticles elicit cross-reactive immune responses to zoonotic coronaviruses in mice

abstract:
Protection against SARS-CoV-2 and SARS-related emergent zoonotic coronaviruses is urgently needed. We made homotypic nanoparticles displaying the receptor-binding domain (RBD) of SARS-CoV-2 or co-displaying SARS-CoV-2 RBD along with RBDs from animal betacoronaviruses that represent threats to humans (mosaic nanoparticles; 4-8 distinct RBDs). Mice immunized with RBD-nanoparticles, but not soluble antigen, elicited cross-reactive binding and neutralization responses. Mosaic-RBD-nanoparticles elicited antibodies with superior cross-reactive recognition of heterologous RBDs compared to sera from immunizations with homotypic SARS-CoV-2-RBD-nanoparticles or COVID-19 convalescent human plasmas. Moreover, sera from mosaic-RBD-immunized mice neutralized heterologous pseudotyped coronaviruses equivalently or better after priming than sera from homotypic SARS-CoV-2-RBD-nanoparticle immunizations, demonstrating no immunogenicity loss against particular RBDs resulting from co-display. A single immunization with mosaic-RBD-nanoparticles provides a potential strategy to simultaneously protect against SARS-CoV-2 and emerging zoonotic coronaviruses.

One sentence summaryNanoparticle strategy for pan-sarbecovirus vaccine

125-character summary for online ToCImmunizing with nanoparticles displaying diverse coronavirus RBDs elicits cross-reactive and neutralizing antibody responses.

category: immunology


title: Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck

abstract:
The evolutionary mechanisms by which SARS-CoV-2 viruses adapt to mammalian hosts and, potentially, undergo antigenic evolution depend on the ways genetic variation is generated and selected within and between individual hosts. Using domestic cats as a model, we show that SARS-CoV-2 consensus sequences remain largely unchanged over time within hosts, while dynamic sub-consensus diversity reveals processes of genetic drift and weak purifying selection. We further identify a notable variant at amino acid position 655 in Spike (H655Y), which was previously shown to confer escape from human monoclonal antibodies. This variant arises rapidly and persists at intermediate frequencies in index cats. It also becomes fixed following transmission in two of three pairs. These dynamics suggest this site may be under positive selection in this system and illustrate how a variant can quickly arise and become fixed in parallel across multiple transmission pairs. Transmission of SARS-CoV-2 in cats involved a narrow bottleneck, with new infections founded by fewer than ten viruses. In RNA virus evolution, stochastic processes like narrow transmission bottlenecks and genetic drift typically act to constrain the overall pace of adaptive evolution. Our data suggest that here, positive selection in index cats followed by a narrow transmission bottleneck may have instead accelerated the fixation of S H655Y, a potentially beneficial SARS-CoV-2 variant. Overall, our study suggests species- and context-specific adaptations are likely to continue to emerge. This underscores the importance of continued genomic surveillance for new SARS-CoV-2 variants as well as heightened scrutiny for signatures of SARS-CoV-2 positive selection in humans and mammalian model systems.

Author summaryThrough ongoing human adaptation, spill-back events from other animal intermediates, or with the distribution of vaccines and therapeutics, the landscape of SARS-CoV-2 genetic variation is certain to change. The evolutionary mechanisms by which SARS-CoV-2 will continue to adapt to mammalian hosts depend on genetic variation generated within and between hosts. Here, using domestic cats as a model, we show that within-host SARS-CoV-2 genetic variation is predominantly influenced by genetic drift and purifying selection. Transmission of SARS-CoV-2 between hosts is defined by a narrow transmission bottleneck, involving 2-5 viruses. We further identify a notable variant at amino acid position 655 in Spike (H655Y), which arises rapidly and is transmitted in cats. Spike H655Y has been previously shown to confer escape from human monoclonal antibodies and is currently found in over 1000 human sequences. Overall, our study suggests species- and context-specific adaptations are likely to continue to emerge, underscoring the importance of continued genomic surveillance in humans and non-human mammalian hosts.


title: Preclinical evaluation of Imatinib does not support its use as an antiviral drug against SARS-CoV-2

abstract: Following the emergence of SARS-CoV-2, the search for an effective and rapidly available treatment was initiated worldwide based on repurposing of available drugs. Previous reports described the antiviral activity of certain tyrosine kinase inhibitors (TKIs) targeting the Abelson kinase 2 against pathogenic coronaviruses. Imatinib, one of them, has more than twenty years of safe utilization for the treatment of hematological malignancies. In this context, Imatinib was rapidly evaluated in clinical trials against Covid-19. Here, we present the pre-clinical evaluation of Imatinib in multiple models. Our results indicated that Imatinib and another TKI, the Masitinib, exhibit an antiviral activity in VeroE6 cells. However, Imatinib was inactive in a reconstructed bronchial human airway epithelium model. In vivo, Imatinib therapy failed to impair SARS-CoV-2 replication in a golden Syrian hamster model despite high concentrations in plasma and in the lung. Overall, these results do not support the use of Imatinib and similar TKIs as antivirals in the treatment of Covid-19.


title: HIV-induced membraneless organelles orchestrate post-nuclear entry steps.

abstract: HIV integration occurs in chromatin sites that favor the release of high levels of viral progeny, alternatively the virus is also able to coexist discreetly with the host. The viral infection perturbs the cellular environment inducing the remodeling of the nuclear landscape. Indeed, HIV-1 triggers the nuclear clustering of the host factor CPSF6, but the underlying mechanism is poorly understood. Our data indicate that HIV usurps a recently discovered biological phenomenon, called liquid-liquid phase separation (LLPS), to hijack the host cell. We observed CPSF6 clusters as part of HIV-induced membraneless organelles (HIV-1 MLOs) in macrophages, which are one of the main HIV target cells. We describe that HIV-1 MLOs follow phase separation rules and represent functional biomolecular condensates. We highlight HIV-1 MLOs as hubs of nuclear reverse transcription, while the double stranded viral DNA, once formed, rapidly migrates outside these structures. Transcription-competent proviruses localize outside but near HIV-1 MLOs, in LEDGF-abundant regions, known to be active chromatin sites. Therefore, HIV-1 MLOs orchestrate viral events prior to the integration step and create a favorable environment for the viral replication. This study uncovers single functional host-proviral complexes in their nuclear landscape, which is markedly restructured by HIV-1 to favor viral replication.


title: Reproductive compensation and embryo screening drive the evolution of polyembryony

abstract:
Simple polyembryony - where one gametophyte produces multiple embryos with different sires but the same maternal haplotype - is common among vascular plants. We show that together polyembryonys two benefits - "reproductive compensation" achieved by providing a backup for inviable embryos, and the opportunity to favor the fitter of surviving embryos, can favor the evolution of polyembryony. To do so, we develop an infinite-site, forward population genetics model to test how these factors can favor the evolution of polyembryony, and how these underlying benefits of polyembryony shape the genetic load under a range of biological parameters. While these two benefits are difficult to disentangle in nature, we construct variant models of polyembryony that either only include or only exclude the opportunity for reproductive compensation. We find that reproductive compensation strongly favors the evolution of polyembryony, and that polyembryony is favored much more weekly in its absence, suggesting that the benefit of a backup embryo is the force favoring polyembryony. Remarkably we find nearly identical results in cases in which mutations impact either embryo or post-embryonic fitness (no pleiotropy), and in cases in which mutations have identical fitness effects embryo or post-embryonic fitness (extreme pleiotropy). Finally, we find that the consequences of polyembryony depends on its function - polyembryony results in a decrease in mean embryonic fitness when acting as a mechanism of embryo compensation, and ultimately increases mean embryonic fitness when we exclude this potential benefit.

Nature is, above all, profligate. Dont believe them when they tell you how economical and thrifty nature is.

- Annie Dillard 1974.


title: Reproductive compensation and selection among viable embryos drive the evolution of polyembryony


title: Rapidly self-sterilizing PPE capable of 99.9% SARS-CoV-2 deactivation in 30 seconds

abstract:
The coronavirus disease 2019 (COVID-19) has created an acute worldwide demand for sustained broadband pathogen suppression in households, hospitals, and public spaces. The US recently passed a new sad milestone of 500,000 deaths due to COVID-19, the highest rate anywhere in the world. In response, we have created a rapid-acting, self-sterilizing PPE configurations capable of killing SARS-CoV-2 and other microbes in seconds. The highly active material destroys pathogens faster than any conventional copper configuration. The material maintains its antimicrobial efficacy over sustained use and is shelf stable. We have performed rigorous testing in accordance with guidelines from U.S. governing authorities and believe that the material could offer broad spectrum, non-selective defense against most microbes via integration into masks and other protective equipment.

SummaryA novel configuration of copper offering continued fast-acting protection against viruses and bacteria, including SARS-CoV-2.


title: An Integrated Workflow for Enhanced Taxonomic and Functional Coverage of the Mouse Faecal Metaproteome

abstract: Intestinal microbiota plays a key role in shaping host homeostasis by regulating metabolism, immune responses and behaviour. Its dysregulation has been associated with metabolic, immune and neuropsychiatric disorders and is accompanied by changes in bacterial metabolic regulation. Although proteomics is well suited for analysis of individual microbes, metaproteomics of faecal samples is challenging due to the physical structure of the sample, presence of contaminating host proteins and coexistence of hundreds of microorganisms. Furthermore, there is a lack of consensus regarding preparation of faecal samples, as well as downstream bioinformatic analyses following metaproteomics data acquisition. Here we assess sample preparation and data analysis strategies applied to mouse faeces in a typical mass spectrometry-based metaproteomic experiment. We show that subtle changes in sample preparation protocols may influence interpretation of biological findings. Two-step database search strategies led to significant underestimation of false positive protein identifications. Unipept software provided the highest sensitivity and specificity in taxonomic annotation of the identified peptides of unknown origin. Comparison of matching metaproteome and metagenome data revealed a positive correlation between protein and gene abundances. Notably, nearly all functional categories of detected protein groups were differentially abundant in the metaproteome compared to what would be expected from the metagenome, highlighting the need to perform metaproteomics when studying complex microbiome samples.


title: UniBind: maps of high-confidence direct TF-DNA interactions across nine species

abstract: Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. Hence, it is critical to locate these TF-DNA interactions to understand transcriptional regulation. The availability of datasets generated by chromatin immunoprecipitation followed by sequencing (ChIP-seq) empowers our efforts to predict the specific locations of TFBSs with greater confidence than previously possible by fusing computational and experimental approaches. In this work, we processed ~10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~56 million TFBSs with experimental and computational evidence for direct TF-DNA interactions for 644 TFs in >1,000 cell lines and tissues. These TFBSs were used to predict >198,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface (https://unibind.uio.no/), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions. UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.


title: Immunogenicity and protective efficacy of one- and two-dose regimens of the Ad26.COV2.S COVID-19 vaccine candidate in adult and aged rhesus macaques

abstract:
Safe and effective coronavirus disease (COVID)-19 vaccines are urgently needed to control the ongoing pandemic. While single-dose vaccine regimens would provide multiple advantages, two doses may improve the magnitude and durability of immunity and protective efficacy. We assessed one- and two-dose regimens of the Ad26.COV2.S vaccine candidate in adult and aged non-human primates (NHP). A two-dose Ad26.COV2.S regimen induced higher peak binding and neutralizing antibody responses compared to a single dose. In one-dose regimens neutralizing antibody responses were stable for at least 14 weeks, providing an early indication of durability. Ad26.COV2.S induced humoral immunity and Th1 skewed cellular responses in aged NHP that were comparable to adult animals. Importantly, aged Ad26.COV2.S-vaccinated animals challenged 3 months post -dose 1 with a SARS-CoV-2 spike G614 variant showed near complete lower and substantial upper respiratory tract protection for both regimens. These are the first NHP data showing COVID-19 vaccine protection against the SARS-CoV-2 spike G614 variant and support ongoing clinical Ad26.COV2.S development.

SummaryCOVID-19 vaccines are urgently needed and while single-dose vaccines are preferred, two-dose regimens may improve efficacy. We show improved Ad26.COV2.S immunogenicity in non-human primates after a second vaccine dose, while both regimens protected aged animals against SARS-CoV-2 disease.


title: Robustness of phylogenetic inference to model misspecification caused by pairwise epistasis

abstract: Likelihood-based phylogenetic inference posits a probabilistic model of character state change along branches of a phylogenetic tree. These models typically assume statistical independence of sites in the sequence alignment. This is a restrictive assumption that facilitates computational tractability, but ignores how epistasis, the effect of genetic background on mutational effects, influences the evolution of functional sequences. We consider the effect of using a misspecified site-independent model on the accuracy of Bayesian phylogenetic inference in the setting of pairwise-site epistasis. Previous work has shown that as alignment length increases, tree reconstruction accuracy also increases. Here, we present a simulation study demonstrating that accuracy increases with alignment size even if the additional sites are epistatically coupled. We introduce an alignment-based test statistic that is a diagnostic for pair-wise epistasis and can be used in posterior predictive checks.


title: A large panel of chicken cells are invaded in vivo by Salmonella Typhimurium even when depleted of all known invasion factors

abstract: Poultry are the main source of human infection by Salmonella. As infected poultry are asymptomatic, the identification of infected poultry farms is difficult. Controlling animal infections is thus of primary importance. As cell tropism is known to govern disease, our aim was therefore to identify the infected host-cell types in chicks and the role of the three known invasion factors in this process (T3SS-1, Rck and PagN). Chicks were inoculated with wild-type or isogenic fluorescent Salmonella Typhimurium mutants via the intraperitoneal route. Our results show that liver, spleen, gall bladder and aortic vessels could be foci of infection and that phagocytic and non-phagocytic cells, including immune, epithelial and endothelial cells, are invaded in vivo in each organ. Moreover, a mutant defective for the T3SS-1, Rck and PagN remained able to colonize organs as the wild-type strain and invaded non-phagocytic cells in each organ studied. As the infection of gall bladder was not really described in chicks, invasion of gall bladder cells was confirmed by immunohistochemistry and infection was shown to last several weeks after inoculation of chicks. All together, these findings provide new insights into the dynamics of Salmonella spread in vivo in chicks at the organ and cellular levels.


title: High content Image Analysis to study phenotypic heterogeneity in endothelial cell monolayers

abstract:
Endothelial cells (EC) are heterogeneous across and within tissues, reflecting distinct, specialised functions. EC heterogeneity has been proposed to underpin EC plasticity independently from vessel microenvironments. However, heterogeneity driven by contact-dependent or short-range cell-cell crosstalk cannot be evaluated with single cell transcriptomic approaches as spatial and contextual information is lost. Nonetheless, quantification of EC heterogeneity and understanding of its molecular drivers is key to developing novel therapeutics for cancer, cardiovascular diseases and for revascularisation in regenerative medicine.

Here, we developed an EC profiling tool (ECPT) to examine individual cells within intact monolayers. We used ECPT to characterise different phenotypes in arterial, venous and microvascular EC populations. In line with other studies, we measured heterogeneity in terms of cell cycle, proliferation, and junction organisation. ECPT uncovered a previously under-appreciated single-cell heterogeneity in NOTCH activation. We correlated cell proliferation with different NOTCH activation states at the single cell and population levels. The positional and relational information extracted with our novel approach is key to elucidating the molecular mechanisms underpinning EC heterogeneity.

Summary statementEndothelial cells heterogeneity is key to complex collective functions and cell behaviour. We developed a novel image based endothelial cell profiling tool and quantified heterogeneity in NOTCH signalling in monolayers.

category: cell biology


title: USP28 deletion and small molecule inhibition destabilises c-Myc and elicits regression of squamous cell lung carcinoma

abstract: Lung squamous cell carcinoma (LSCC) is a considerable global health burden, with an incidence of over 600,000 cases per year. Treatment options are limited, and patient 5-year survival rate is less than 5%. The ubiquitin specific protease 28 (USP28) has been implicated in tumorigenesis through its stabilization of the oncoprotein c-MYC. Here, we show that genetic inactivation of Usp28 induced regression of established murine LSCC lung tumors. We developed a small molecule that inhibits USP28 activity in the low nanomole range. While displaying cross-reactivity against the closest homologue USP25, this inhibitor showed a high degree of selectivity over other deubiquitinases. USP28 inhibitor treatment resulted in a dramatic decrease in c-Myc proteins levels and consequently induced substantial regression of autochthonous murine LSCC tumors and human LSCC xenografts, thereby phenocopying the effect observed by genetic deletion. Thus, USP28 may represent a promising therapeutic target for the treatment of squamous cell lung carcinoma.

category: cancer biology


title: New Measures of Visual Scanning Efficiency and Cognitive Effort

abstract: Visual scanning plays an important role in sampling visual information from the surrounding environments for a lot of everyday sensorimotor tasks, such as walking and car driving. In this paper, we consider the problem of visual scanning mechanism underpinning sensorimotor tasks in 3D dynamic environments. We exploit the use of eye tracking data as a behaviometric, for indicating the visuo-motor behavioral measures in the context of virtual driving. A new metric of visual scanning efficiency (VSE), which is defined as a mathematical divergence between a fixation distribution and a distribution of optical flows induced by fixations, is proposed by making use of a widely-known information theoretic tool, namely the square root of Jensen-Shannon divergence. Based on the proposed efficiency metric, a cognitive effort measure (CEM) is developed by using the concept of quantity of information. Psychophysical eye tracking studies, in virtual reality based driving, are conducted to reveal that the new metric of visual scanning efficiency can be employed very well as a proxy evaluation for driving performance. In addition, the effectiveness of the proposed cognitive effort measure is demonstrated by a strong correlation between this measure and pupil size change. These results suggest that the exploitation of eye tracking data provides an effective behaviometric for sensorimotor activity.


title: Probing the aggregated effects of purifying selection per individual on 1,380 medical phenotypes in the UK biobank

abstract:
Understanding the relationship between natural selection and phenotypic variation has been a long-standing challenge in human population genetics. With the emergence of biobank-scale datasets, along with new statistical metrics to approximate strength of purifying selection at the variant level, it is now possible to correlate a proxy of individual relative fitness with a range of medical phenotypes. We calculated a per-individual deleterious load score by summing the total number of derived alleles per individual after incorporating a weight that approximates strength of purifying selection. We assessed four methods for the weight, including GERP, phyloP, CADD, and fitcons. By quantitatively tracking each of these scores with the site frequency spectrum, we identified phyloP as the most appropriate weight. The phyloP-weighted load score was then calculated across 15,129,142 variants in 335,161 individuals from the UK Biobank and tested for association on 1,380 medical phenotypes. After accounting for multiple test correction, we observed a strong association of the load score amongst coding sites only on 27 traits including body mass, adiposity and metabolic rate. We further observed that the association signals were driven by common variants (derived allele frequency > 5%) with high phyloP score (phyloP > 2). Finally, through permutation analyses, we showed that the load score amongst coding sites had an excess of nominally significant associations on many medical phenotypes. These results suggest a broad impact of deleterious load on medical phenotypes and highlight the deleterious load score as a tool to disentangle the complex relationship between natural selection and medical phenotypes.

Author summaryThis study aims to augment our understanding between the complex relation between natural selection and human phenotypic variation. We developed a load score to approximate the relative fitness of an individual and correlate it with a set of medical phenotypes. Association tests between the load score amongst coding sites and 1,380 phenotypes in a sample of 335,161 individuals from the UK Biobank showed a strong association with 27 traits including body mass, adiposity and metabolic rate. Furthermore, an excess of nominal associations at suggestive levels was observed between the load score amongst coding sites and medical phenotypes than would be expected under a null model. These results suggest that the aggregate effect of deleterious mutations as measured by the load score has a broad effect on human phenotypes.

category: genetics


title: Fully Convolutional Networks for pathophysiological Intracranial Pressure waveform classification

abstract: It is generally assumed that the analysis of intracranial pressure (ICP) waveforms could be used for the detection of multiple cerebral pathophysiologies. A main obstacle for the analysis of ICP waveforms is given by the large variation of their generating signal, the arterial blood pressure (ABP). Using extended principal component analysis (PCA) we show that it is possible to distinguish between ICP waveforms generated by pathological ABP waveforms, e.g. in the case of a heart failure, without loosing information about the state of the cerebral compliance. We also create a dataset for ICP pulse classification that can be used to train models to distinguish between an intact and diminished cerebral compliance, as a function of the ICP-generating ABP pulse. As a baseline for classification we build a fully convolutional network (FCN) that reaches high performance with relatively few data.

category: physiology


title: PseudotimeDE: inference of differential gene expression along cell pseudotime with well-calibrated p-values from single-cell RNA sequencing data

abstract: To investigate molecular mechanisms underlying cell state changes, a crucial analysis is to identify differentially expressed (DE) genes along the pseudotime inferred from single-cell RNA-sequencing data. However, existing methods do not account for pseudotime inference uncertainty, and they have either ill-posed p-values or restrictive models. Here we propose PseudotimeDE, a DE gene identification method that adapts to various pseudotime inference methods, accounts for pseudotime inference uncertainty, and outputs well-calibrated p-values. Comprehensive simulations and real-data applications verify that PseudotimeDE outperforms existing methods in false discovery rate control and power.


title: A generalised mechano-kinetic model for use in multiscale simulation protocols.

abstract:
Many biophysical systems and proteins undergo mesoscopic conformational changes in order to perform their biological function. However, these conformational changes often result from a cascade of atomistic interactions within a much larger overall object. For simulations of such systems, the computational cost of retaining high-resolution structural and dynamical information whilst at the same time observing large scale motions over long times is high. Furthermore, this cost is only going to increase as ever larger biological objects are observed experimentally at high resolution.

We derive a generalised mechano-kinetic simulation framework which aims to compensate for these disparate time scales, capable of dynamically exploring a defined mechanical energy landscape whilst at the same time performing kinetic transitions between discretely defined states. With insight from the theory of Markov state models and transition state theory, the framework models continuous dynamical motion at coarse-grained length scales whilst simultaneously including fine-detail, chemical reactions or conformational changes implicitly using a kinetic model. The kinetic model is coupled to the dynamic model in a highly generalised manner, such that kinetic transition rates are continuously modulated by the dynamic simulation. Further, it can be applied to any defined continuous energy landscape, and hence, any existing dynamic simulation framework. We present a series of analytical examples to validate the framework, and showcase its capabilities for studying more complex systems by simulating protein unfolding via single-molecule force spectroscopy on an atomic force microscope.

Author summaryOur intention with this work is to provide a generalised, highly coarse-grained model to allow kinetic processes (conformational changes, protein unfolding, chemical reactions etc) to occur within the context of a dynamic simulation. Performing computationally intensive dynamic simulations to obtain kinetic information can be unjustifiably costly and scientifically inefficient, and so we instead want to emphasise that experimentally available kinetic information can be used to infer the underlying dynamics they result from. We hope our work can begin a discussion on the topic of computationally efficient science, and continue the drive towards collaborative science between theory and experimentation.


title: High-dimensional association detection in large scale genomic data

abstract: Joint analyses of genomic datasets obtained in multiple different conditions are essential for understanding the biological mechanism that drives tissue-specificity and cell differentiation, but they still remain computationally challenging. To address this we introduce CLIMB (Composite LIkelihood eMpirical Bayes), a statistical methodology that learns patterns of condition-specificity present in genomic data. CLIMB provides a generic framework facilitating a host of analyses, such as clustering genomic features sharing similar condition-specific patterns and identifying which of these features are involved in cell fate commitment. We apply CLIMB to three sets of hematopoietic data, which examine CTCF ChIP-seq measured in 17 different cell populations, RNA-seq measured across constituent cell populations in three committed lineages, and DNase-seq in 38 cell populations. Our results show that CLIMB improves upon existing alternatives in statistical precision, while capturing interpretable and biologically relevant clusters in the data.


title: Everything you wanted to know about Mayaro virus but were afraid to ask: Characterization and lifecycle of Mayaro virus in vertebrate and invertebrate cellular backgrounds

abstract: Mayaro virus (MAYV) is an emerging new world alphavirus (genus Alphavirus, family Togaviridae) that causes acute multiphasic febrile illness, skin rash, polyarthritis, and occasional severe clinical phenotypes. The virus lifecycle alternates between invertebrate and vertebrate hosts. Here we characterize the replication features, cell entry, life cycle, and virus-related cell pathology of MAYV using vertebrate and invertebrate in vitro models. Electron dense clathrin-coated pits in infected cells, and reduced viral production in the presence of dynasore, ammonium chloride, and bafilomycin, indicates that viral entry occurs through pH-dependent endocytosis. Increase in FITC-dextran uptake (an indicator of macropinocytosis) in MAYV-infected cells, and dose-dependent infection inhibition by 5-(N-ethyl-N-isopropyl) amiloride (a macropinocytosis inhibitor), indicated that macropinocytosis is an additional entry mechanism of MAYV in vertebrate cells. Acutely infected vertebrate and invertebrate cells formed cytoplasmic or membrane-associated extracytoplasmic replication complexes. Mosquito cells showed modified hybrid cytoplasmic vesicles that supported virus replication, nucleocapsid production, and maturation. Mature virus particles were released from cells by both exocytosis and budding from the cell membrane. MAYV replication was cytopathic and associated with induction of apoptosis by the intrinsic pathway, and later by the extrinsic pathway in infected vertebrate cells. Given that MAYV is expanding its geographical existence as a potential public health problem, this study lays the foundation of biological understanding valuable for therapeutic and preventive interventions.


title: Structural requirements for dihydrobenzoxazepinone anthelmintics: actions against medically important and model parasites - Trichuris muris, Brugia malayi, Heligmosomoides polygyrus and Schistosoma mansoni

abstract:
Nine hundred million people are infected with the soil-transmitted helminths Ascaris lumbricoides (roundworm), hookworm, and Trichuris trichiura (whipworm). However, low single-dose cure rates of the benzimidazole drugs, the mainstay of preventative chemotherapy for whipworm, together with parasite drug resistance, mean that current approaches may not be able to eliminate morbidity from Trichuriasis. We are seeking to develop new anthelmintic drugs specifically with activity against whipworm as a priority, and previously identified a hit series of dihydrobenzoxazepinone (DHB) compounds that block motility of ex vivo Trichuris muris. Here we report a systematic investigation of the structure-activity relationship of the anthelmintic activity of DHB compounds. We synthesised 47 analogues, which allowed us to define features of the molecules essential for anthelmintic action, as well as broadening the chemotype by identification of dihydrobenzoquinolinones (DBQ) with anthelmintic activity. We investigated the activity of these compounds against other parasitic nematodes, identifying DHB compounds with activity against Brugia malayi and Heligmosomoides polygyrus. We also demonstrated activity of DHB compounds against the trematode Schistosoma mansoni, a parasite that causes schistosomiasis. These results demonstrate the potential of DHB and DBQ compounds for further development as broad-spectrum anthelmintics.

Author summaryAround a billion people are infected by the soil transmitted helminths Ascaris, hookworm and whipworm. In the case of whipworm, the benzimidazole drugs, which are distributed to school children in affected areas, have low cure rates. This means that finding an improved treatment for whipworm is a priority. We previously identified five DHB compounds in a screen for new compounds active against whipworm. Here we systematically dissect these molecules, making 47 modified versions of the compounds. This allowed us to define the features of these compounds that are important for activity against whipworm. We also demonstrate activity of DHB compounds against other parasitic nematodes, and against Schistosoma mansoni, a trematode parasite. These results show the potential for further development of DHB compounds as broad-spectrum anthelmintics.

category: pharmacology and toxicology


title: Optimizing interneuron circuits for compartment-specific feedback inhibition

abstract: Cortical circuits process information by rich recurrent interactions between excitatory neurons and inhibitory interneurons. One of the prime functions of interneurons is to stabilize the circuit by feedback inhibition, but the level of specificity on which inhibitory feedback operates is not fully resolved. We hypothesized that inhibitory circuits could enable separate feedback control loops for different synaptic input streams, by means of specific feedback inhibition to different neuronal compartments. To investigate this hypothesis, we adopted an optimization approach. Leveraging recent advances in training spiking network models, we optimized the connectivity and short-term plasticity of interneuron circuits for compartment-specific feedback inhibition onto pyramidal neurons. Over the course of the optimization, the interneurons diversified into two classes that resembled parvalbumin (PV) and somatostatin (SST) expressing interneurons. The resulting circuit can be understood as a neural decoder that inverts the nonlinear biophysical computations performed within the pyramidal cells. Our model provides a proof of concept for studying structure-function relations in cortical circuits by a combination of gradient-based optimization and biologically plausible phenomenological models.


title: Substrate traits govern the assembly of microbial community engaged in metabolic division of labor

abstract: Metabolic division of labor (MDOL) is widespread in nature, whereby a complex metabolic pathway is shared between different strains within a community for mutual benefit. However, little is known about how the mutual interactions in the microbial community engaged in MDOL are regulated. We hypothesized that when degradation of an organic compound is carried out via MDOL, the substrate traits (i.e., concentration and its toxicity) modulate the benefit allocation between the two microbial populations, thus affecting the structure of this community. We tested this hypothesis by combining mathematical modelling with experiments using engineered synthetic microbial consortia. Numerous modelling analyses suggested that the proportion of the population executing the first metabolic step can be simply estimated by Monod-like formulas governed by substrate traits. The model and the proposed formula quantitatively predicted the structure of our synthetic consortia composed of two strains degrading salicylate through MDOL. Individual-based modelling and colony pattern formation assays further indicated that our rule is also applicable to estimating community structure in spatially structured environments. Our results demonstrate that the structure of the microbial communities can be quantitatively predicted from simple environmental factors, such as substrate concentration and its toxicity, which provides novel perspectives on understanding the assembly of natural communities, as well as insights into how to manage artificial microbial systems.


title: Substrate traits shape the structure of microbial community engaged in metabolic division of labor


title: A Modular Workflow for Model Building, Analysis, and Parameter Estimation in Systems Biology and Neuroscience

abstract:
Neuroscience incorporates knowledge from a range of scales, from molecular dynamics to neural networks. Modeling is a valuable tool in understanding processes at a single scale or the interactions between two adjacent scales and researchers use a variety of different software tools in the model building and analysis process. While systems biology is among the more standardized fields, conversion between different model formats and interoperability between various tools is still somewhat problematic. To offer our take on tackling these shortcomings and by keeping in mind the FAIR (findability, accessibility, interoperability, reusability) data principles, we have developed a workflow for building and analyzing biochemical pathway models, using pre-existing tools that could be utilized for the storage and refinement of models in all phases of development. We have chosen the SBtab format which allows the storage of biochemical models and associated data in a single file and provides a human readable set of syntax rules. Next, we implemented custom-made MATLAB(R) scripts to perform parameter estimation and global sensitivity analysis used in model refinement. Additionally, we have developed a web-based application for biochemical models that allows simulations with either a network free solver or stochastic solvers and incorporating geometry. Finally, we illustrate convertibility and use of a biochemical model in a biophysically detailed single neuron model by running multiscale simulations in NEURON. Using this workflow, we can simulate the same model in three different simulators, with a smooth conversion between the different model formats, enhancing the characterization of different aspects of the model.

Information Sharing StatementBoth the source code and documentation of the Subcellular Workflow are available at https://github.com/jpgsantos/Subcellular_Workflow and licensed under GNU General Public License v3.0. The model is stored in the SBtab format (Lubitz et al. 2016). Model reduction, parameter estimation and global sensitivity analysis tools are written in MATLAB(R) (RRID:SCR_001622) and require the SimBiology(R) toolbox. Conversion script to VFGEN (Weckesser 2008), MOD and SBML (RRID:SCR_007422) is written in R (RRID:SCR_001905). Conversion to SBML requires the use of libSBML (RRID:SCR_014134). Validations are run in COPASI (RRID:SCR_014260; Hoops et al. 2006), NEURON (RRID:SCR_005393; Hines and Carnevale 1997) and with the subcellular simulation setup application (RRID:SCR_018790; available at https://subcellular.humanbrainproject.eu/model/simulations) that uses a spatial solver provided by STEPS (RRID:SCR_008742; Hepburn et al. 2012) and network-free solver NFsim (available at http://michaelsneddon.net/nfsim/). The medium spiny neuron model (Lindroos et al. 2018) used in NEURON simulations is available in ModelDB database (RRID:SCR_007271) with access code 237653. The FindSim use case model is available in https://github.com/BhallaLab/FindSim (Viswan et al. 2018).


title: The ChiS-family DNA-binding domain contains a cryptic helix-turn-helix variant

abstract:
Sequence specific DNA-binding domains (DBDs) are conserved in all domains of life. These proteins carry out a variety of cellular functions, and there are a number of distinct structural domains already described that allow for sequence-specific DNA binding, including the ubiquitous helix-turn-helix (HTH) domain. In the facultative pathogen Vibrio cholerae, the chitin sensor ChiS is a transcriptional regulator that is critical for the survival of this organism in its marine reservoir. We have recently shown that ChiS contains a cryptic DBD in its C-terminus. This domain is not homologous to any known DBD, but it is a conserved domain present in other bacterial proteins. Here, we present the crystal structure of the ChiS DBD at a resolution of 1.28 [A]. We find that the ChiS DBD contains an HTH domain that is structurally similar to those found in other DNA binding proteins, like the LacI repressor. However, one striking difference observed in the ChiS DBD is that the canonical tight "turn" of the HTH is replaced with an insertion containing a {beta}-sheet, a variant which we term the "helix-sheet-helix". Through systematic mutagenesis of all positively charged residues within the ChiS DBD, we show that residues within and proximal to the ChiS helix-sheet-helix are critical for DNA binding. Finally, through phylogenetic analyses we show that the ChiS DBD is found in diverse Proteobacterial proteins that exhibit distinct domain architectures. Together, these results suggest that the structure described here represents the prototypical member of the ChiS-family of DBDs.

ImportanceRegulating gene expression is essential in all domains of life. This process is commonly facilitated by the activity of DNA-binding transcription factors. There are diverse structural domains that allow proteins to bind to specific DNA sequences. The structural basis underlying how some proteins bind to DNA, however, remains unclear. Previously, we showed that in the major human pathogen Vibrio cholerae the transcription factor ChiS directly regulates gene expression through a cryptic DNA binding domain. This domain lacked homology to any known DNA-binding protein. In the current study, we determined the structure of the ChiS DNA binding domain (DBD) and find that the ChiS-family DBD is a cryptic variant of the ubiquitous helix-turn-helix (HTH) domain. We further demonstrate that this domain is conserved in diverse proteins that may represent a novel group of transcriptional regulators.


title: Universally available herbal teas based on sage and perilla elicit potent antiviral activity against SARS-CoV-2 variants of concern by HMOX-1 upregulation in human cells

abstract: The current SARS-CoV-2/COVID-19 pandemic wreaks medical and socioeconomic havoc. Despite the availability of vaccines, cost-effective acute treatment options preventing morbidity and mortality are urgently needed. To identify affordable, ubiquitously available, and effective treatments, we tested herbs consumed worldwide as herbal teas regarding their antiviral activity against SARS-CoV-2. Aqueous infusions prepared by boiling leaves of the Lamiaceae perilla and sage elicit potent and sustained antiviral activity against SARS-CoV-2 in therapeutic as well as prophylactic regimens. The herbal infusions exerted antiviral effects comparable to interferon-{beta} and remdesivir but outperformed convalescent sera and interferon-2 upon short-term treatment early after infection. Based on protein fractionation analyses, we identified caffeic acid, perilla aldehyde, and perillyl alcohol as antiviral compounds. Global mass spectrometry (MS) analyses performed comparatively in two different cell culture infection models revealed changes of the proteome upon treatment with herbal infusions and provided insights into the mode of action. As inferred by the MS data, induction of heme oxygenase 1 (HMOX-1) was confirmed as effector mechanism by the antiviral activity of the HMOX-1-inducing compounds sulforaphane and fraxetin. In conclusion, herbal teas based on perilla and sage exhibit antiviral activity against SARS-CoV-2 including variants of concern such as Alpha, Beta, Delta, and Omicron.


title: Cell-to-cell heterogeneity in Sox2 and Brachyury expression ratios guides progenitor destiny by controlling their motility.

abstract: Although cell-to-cell heterogeneity in gene and protein expression within cell populations has been widely documented, we know little about its potential biological functions. We addressed this issue by studying progenitors that populate the posterior region of the vertebrate embryos, a cell population known for its capacity to self-renew or to contribute to the formation of the neural tube and paraxial mesoderm tissues. Posterior progenitors are characterized by the co-expression of Sox2 and Brachyury (Bra), two transcription factors related to neural and mesodermal lineages, respectively. In this study, we show that the respective levels of Sox2 and Bra proteins display a high degree of variability among posterior progenitors of the quail embryo. By developing forced expression and downregulation approaches, we further provide evidence that the value of the Sox2-to-Bra ratio in a given progenitor directs its choice of staying in place or exit the progenitor zone to generate neural or mesodermal cells. Time-lapse imaging together with mathematical modeling then reveal that variations of the Sox2-to-Bra ratio confer these cells heterogeneous motile behaviors. While high Bra levels display high motile properties that push cells to join the mesoderm, high levels of Sox2 tend to inhibit cell movement making cells get integrated into the neural tube. Our work thus provides evidence that the spatial heterogeneity of posterior progenitors, with regards to their relative contents of Sox2 and Bra and thus to their motile properties, is fundamental to maintain a pool of resident progenitors while others segregate to contribute to tissue formation. More broadly, our work reveals that heterogeneity among a population of progenitor cells is critical to ensure robust multi-tissue morphogenesis.


title: Transgenerational Regulation of Sexual Attractiveness in C. elegans Nematodes

abstract: It is unknown whether transient transgenerational epigenetic responses to environmental challenges affect the process of evolution, which typically unfolds over many generations. Here we show that in C. elegans, inherited small RNAs shape the hard-wired genome and control genetic variation by regulating the decision of whether to self-fertilize or outcross. We found that under stressful temperatures younger hermaphrodites secrete a male-attracting pheromone. Attractiveness transmits transgenerationally to unstressed progeny via heritable small RNAs and the Argonaute Heritable-RNAi-Deficient-1. We identified an endogenous small interfering RNA pathway, enriched in endo-siRNAs which target sperm genes, that transgenerationally regulates sexual attraction, male prevalence, and outcrossing rates. Multigenerational mating competitions and mathematical simulations revealed that over generations, animals that inherit attractiveness mate more, and their alleles spread in the population. We propose that sperm serves as a "stress sensor" which, via small RNA inheritance, can promote outcrossing in challenging environments, when increasing genetic variation is advantageous.

category: molecular biology


title: ABC-finder: A containerized web server for the identification and topology prediction of ABC proteins

abstract:
In view of the multiple clinical and physiological implications of ABC transporter proteins, there is a considerable interest among researchers to characterize them functionally. However, such characterizations are based on the premise that ABC proteins are accurately identified in the proteome of an organism, and their topology is correctly predicted. With this objective, we have developed ABC-finder, i.e., a Docker-based package for the identification of ABC proteins in all organisms, and visualization of the topology of ABC proteins using a web browser. ABC-finder is built and deployed in a Linux container, making it scalable for many concurrent users on our servers and enabling users to download and run it locally. Overall, ABC-finder is a convenient, portable, and platform-independent tool for the identification and topology prediction of ABC proteins.

ABC-finder is accessible at http://abc-finder.osdd.jnu.ac.in


title: Recurrent deletions in the SARS-CoV-2 spike glycoprotein drive antibody escape

abstract: Zoonotic pandemics, like that caused by SARS-CoV-2, can follow the spillover of animal viruses into highly susceptible human populations. Their descendants have adapted to the human host and evolved to evade immune pressure. Coronaviruses acquire substitutions more slowly than other RNA viruses, due to a proofreading polymerase. In the spike glycoprotein, we find recurrent deletions overcome this slow substitution rate. Deletion variants arise in diverse genetic and geographic backgrounds, transmit efficiently, and are present in novel lineages, including those of current global concern. They frequently occupy recurrent deletion regions (RDRs), which map to defined antibody epitopes. Deletions in RDRs confer resistance to neutralizing antibodies. By altering stretches of amino acids, deletions appear to accelerate SARS-CoV-2 antigenic evolution and may, more generally, drive adaptive evolution.


title: Genetic and pharmacological manipulations of glyoxalase 1 mediate ethanol withdrawal seizure susceptibility in mice

abstract: Central nervous system (CNS) hyperexcitability is a clinically significant feature of acute ethanol withdrawal. There is evidence for a genetic contribution to withdrawal severity, but specific genetic risk factors have not be identified. The gene glyoxalase 1 (Glo1) has been previously implicated in ethanol consumption in mice, and GLO1 inhibition can attenuate drinking in mice and rats. Here, we investigated whether genetic and pharmacological manipulations of GLO1 activity can also mediate ethanol withdrawal seizure severity in mice. Mice from two transgenic lines overexpressing Glo1 on different genetic backgrounds (C57BL/6J [B6] and FVB/NJ [FVB]) were tested for handling induced convulsions (HICs) as a measure of acute ethanol withdrawal. Following an injection of 4 g/kg alcohol, both B6 and FVB mice overexpressing Glo1 showed increased HICs compared to wild type littermates, though only the FVB line showed a statistically significant difference. We also administered daily ethanol injections (2 g/kg + 9 mg/kg 4-methylpyrazole) to wild type B6 mice for 10 days and tested them for HICs on the 10th day, following treatment with either vehicle or a GLO1 inhibitor (S-bromobenzylglutathione cyclopentyl diester [pBBG]). Treatment with pBBG reduced HICs, although this effect was only statistically significant following two 10-day cycles of ethanol exposure and withdrawal. These results provide converging genetic and pharmacological evidence that GLO1 can mediate ethanol withdrawal seizure susceptibility. We discuss the possible role of GLO1s substrate, methylglyoxal, acting as an agonist at GABAA receptors.


title: Complex trait-environment relationships underlie the structure of Wisconsin forest plant communities

abstract: Plant species shift in abundance as environmental conditions change because traits adapt species to particular conditions. As a result, trait values shift along environmental gradients--the so-called trait-environment relationships. These relationships are often assessed by regressing community-weighted mean (CWM) traits on environmental gradients. Such regressions (CWMr) assume that local communities exhibit centered optimum trait-abundance relationships and that traits are not independent from one another. However, the shape of trait-abundance relationships can vary widely along environmental gradients--reflecting the interaction between traits and gradients--and traits are usually interrelated. Accounting for these complexities should improve our ability to accurately describe trait-environment relationships. We tested these ideas by analyzing how abundances of 185 herbaceous understory species distributed among 189 forested sites in Wisconsin, USA, varied in response to four functional traits (vegetative height-VH, leaf size-LS, leaf mass per area-LMA, and leaf carbon content) and six soil and climate variables. A generalized linear mixed model (GLMM) allowed us to assess how the shape of trait-abundance relationships changed along environmental gradients for the 24 trait-environment combinations simultaneously. We then compared the resulting trait-environment relationships to those estimated via CWMr. The GLMM identified five significant trait-environment relationships that together explained [~]40% of variation in species abundances across sites. Temperature played important roles with warmer and more seasonal sites favoring taller plants. Soil texture and temperature seasonality affected LS and LMA more modestly; these seasonality effects declined at more seasonal sites. Only some traits under certain conditions showed centered optimum trait- abundance relationships. Concomitantly, CWMr identified 17 significant trait-environment relationships including effects of temperature, precipitation, and soil on LMA as often reported in other studies. Despite this overidentification, CWMr failed to detect significant temperature-seasonality effects found in the GLMM. Modeling the complexity of how traits and environments interact to affect plant abundance allows us to identify and rank key trait- environment relationships. Although the GLMM model was more complex compared to single CWM regressions, it identified a simpler hierarchy of trait-environment relationships that accurately and reliably predicted responses of forest understory species to gradients in environmental conditions.

category: ecology


title: Human methylome variation across Infinium 450K data on the Gene Expression Omnibus

abstract: While DNA methylation (DNAm) is the most-studied epigenetic mark, few recent studies probe the breadth of publicly available DNAm array samples. We collectively analyzed 35,360 Illumina Infinium HumanMethylation450K DNAm array samples published on the Gene Expression Omnibus (GEO). We learned a controlled vocabulary of sample labels by applying regular expressions to metadata and used existing models to predict various sample properties including epigenetic age. We found approximately two-thirds of samples were from blood, one-quarter were from brain, and one-third were from cancer patients. 19% of samples failed at least one of Illuminas 17 prescribed quality assessments; signal distributions across samples suggest modifying manufacturer-recommended thresholds for failure would make these assessments more informative. We further analyzed DNAm variances in seven tissues (adipose, nasal, blood, brain, buccal, sperm, and liver) and characterized specific probes distinguishing them. Finally, we compiled DNAm array data and metadata, including our learned and predicted sample labels, into database files accessible via the recountmethylation R/Bioconductor companion package. Its vignettes walk the user through some analyses contained in this paper.


title: Universal annotation of the human genome through integration of over a thousand epigenomic datasets

abstract:
BackgroundGenome-wide maps of chromatin marks such as histone modifications and open chromatin sites provide valuable information for annotating the non-coding genome, including identifying regulatory elements. Computational approaches such as ChromHMM have been applied to discover and annotate chromatin states defined by combinatorial and spatial patterns of chromatin marks within the same cell type. An alternative  stacked modeling approach was previously suggested, where chromatin states are defined jointly from datasets of multiple cell types to produce a single universal genome annotation based on all datasets. Despite its potential benefits for applications that are not specific to one cell type, such an approach was previously applied only for small-scale specialized purposes. Large-scale applications of stacked modeling have previously posed scalability challenges.

ResultsUsing a version of ChromHMM enhanced for large-scale applications, we applied the stacked modeling approach to produce a universal chromatin state annotation of the human genome using over 1000 datasets from more than 100 cell types, with the learned model denoted as the full-stack model. The full-stack model states show distinct enrichments for external genomic annotations, which we used in characterizing each state. Compared to per-cell-type annotations, the full-stack annotations directly differentiate constitutive from cell type specific activity and is more predictive of locations of external genomic annotations.

ConclusionsThe full-stack ChromHMM model provides a universal chromatin state annotation of the genome and a unified global view of over 1000 datasets. We expect this to be a useful resource that complements existing per-cell-type annotations for studying the non-coding human genome.


title: Topology-based Sparsification of Graph Annotations

abstract: Since the amount of published biological sequencing data is growing exponentially, efficient methods for storing and indexing this data are more needed than ever to truly benefit from this invaluable resource for biomedical research. Labeled de Bruijn graphs are a frequently-used approach for representing large sets of sequencing data. While significant progress has been made to succinctly represent the graph itself, efficient methods for storing labels on such graphs are still rapidly evolving. In this paper, we present RowDiff, a new technique for compacting graph labels by leveraging expected similarities in annotations of vertices adjacent in the graph. RowDiff can be constructed in linear time relative to the number of vertices and labels in the graph, and in space proportional to the graph size. In addition, construction can be efficiently parallelized and distributed, making the technique applicable to graphs with trillions of nodes. RowDiff can be viewed as an intermediary sparsification step of the original annotation matrix and can thus naturally be combined with existing generic schemes for compressed binary matrices. Experiments on 10,000 RNA-seq datasets show that RowDiff combined with Multi-BRWT results in a 30% reduction in annotation footprint over Mantis-MST, the previously known most compact annotation representation. Experiments on the sparser Fungi subset of the RefSeq collection show that applying RowDiff sparsification reduces the size of individual annotation columns stored as compressed bit vectors by an average factor of 42. When combining RowDiff with a Multi-BRWT representation, the resulting annotation is 26 times smaller than Mantis-MST.


title: The Mutational Signature Comprehensive Analysis Toolkit (musicatk) for the discovery, prediction, and exploration of mutational signatures

abstract: Mutational signatures are patterns of somatic alterations in the genome caused by carcinogenic exposures or aberrant cellular processes. To provide a comprehensive workflow for preprocessing, analysis, and visualization of mutational signatures we created the Mutational Signature Comprehensive Analysis Toolkit (musicatk) package. musicatk enables users to select different schemas for counting mutation types and easily combine count tables from different schemas. Multiple distinct methods are available to deconvolute signatures and exposures or to predict exposures in individual samples given a pre-existing set of signatures. Additional exploratory features include the ability to compare signatures to the COSMIC database, embed tumors in two dimensions with UMAP, cluster tumors into subgroups based on exposure frequencies, identify differentially active exposures between tumor subgroups and plot exposure distributions across user-defined annotations such as tumor type. Overall, musicatk will enable users to gain novel insights into the patterns of mutational signature observed in cancer cohorts.


title: Raman microscopy-based quantification of the physical properties of intracellular lipids

abstract: The physical properties of lipids, such as viscosity, are homeostatically maintained in cells and are intimately involved in physiological roles. Measurement of the physical properties of plasma membranes has been achieved primarily through chemical or genetically encoded fluorescent probes. However, as most probes target plasma membranes, physical properties of lipids in intracellular organelles, including lipid droplets (LDs) are yet to be analyzed. Here, we present a novel Raman microscopy-based approach for quantifying the physical properties of intracellular lipids under deuterium-labeled fatty acid treatment conditions. Focusing on the fact that Raman spectra of carbon-deuterium vibration are altered depending on the surrounding lipid species, we quantitatively represented the physical properties of lipids as the gauche/trans conformational ratio of the introduced labeled fatty acids, which can be used as an indicator of viscosity. Intracellular Raman imaging revealed that the gauche/trans ratio of cytosolic regions was robustly preserved against perturbations attempting to alter the lipid composition. This was likely due to LDs functioning as a buffer against excess gauche/trans ratio. Our novel approach enables the observation of the physical properties of organelle lipids, which is difficult to perform with conventional probes, and is useful for quantitative assessment of the subcellular lipid environment.


title: Temporal dynamics of implicit memory underlying serial dependence

abstract: Serial dependence is the effect in which the immediately preceding trial influences participants' responses to the current stimulus. But for how long does this bias last in the absence of interference from other stimuli? Here, we had 20 healthy young adult participants (12 women) perform a coincident timing task using different inter-trial intervals to characterize the serial dependence effect as the time between trials increases. Our results show that serial dependence abruptly decreases from 0.1 s to 1 s inter-trial interval, but it remains pronounced after that for up to 8 s. In addition, participants' response variability slightly decreases over longer intervals. We discuss these results in light of recent models suggesting that serial dependence might rely on a short-term memory trace kept through changes in synaptic weights, which might explain its long duration and apparent stability over time.

category: neuroscience


title: Multisensory integration-attention trade-off in cochlear-implanted deaf individuals

abstract:
AO_SCPLOWBSTRACTC_SCPLOWThe cochlear implant (CI) allows profoundly deaf individuals to partially recover hearing. Still, due to the coarse acoustic information provided by the implant, CI users have considerable difficulties in recognizing speech, especially in noisy environments. CI users therefore rely heavily on visual cues to augment speech comprehension, more so than normal-hearing individuals. However, it is unknown how attention to one (focused) or both (divided) modalities plays a role in multisensory speech recognition. Here we show that unisensory speech listening and reading were negatively impacted in divided-attention tasks for CI users - but not for normal-hearing individuals. Our psychophysical experiments revealed that, as expected, listening thresholds were consistently better for the normal-hearing, while lipreading thresholds were largely similar for the two groups. Moreover, audiovisual speech recognition for normal-hearing individuals could be described well by probabilistic summation of auditory and visual speech recognition, while CI users were better integrators than expected from statistical facilitation alone. Our results suggest that this benefit in integration comes at a cost. Unisensory speech recognition is degraded for CI users when attention needs to be divided across modalities. We conjecture that CI users exhibit an integration-attention trade-off. They focus solely on a single modality during focused-attention tasks, but need to divide their limited attentional resources in situations with uncertainty about the upcoming stimulus modality. We argue that in order to determine the benefit of a CI for speech comprehension, situational factors need to be discounted by presenting speech in realistic or complex audiovisual environments.

SO_SCPLOWIGNIFICANCEC_SCPLOWO_SCPCAP C_SCPCAPO_SCPLOWSTATEMENTC_SCPLOWDeaf individuals using a cochlear implant require significant amounts of effort to listen in noisy environments due to their impoverished hearing. Lipreading can benefit them and reduce the burden of listening by providing an additional source of information. Here we show that the improved speech recognition for audiovisual stimulation comes at a cost, however, as the cochlear-implant users now need to listen and speech-read simultaneously, paying attention to both modalities. The data suggests that cochlear-implant users run into the limits of their attentional resources, and we argue that they, unlike normal-hearing individuals, always need to consider whether a multisensory benefit outweighs the unisensory cost in everyday environments.


title: Selective vulnerability of supragranular layer neurons in schizophrenia

abstract: Schizophrenia is one of the most wide-spread mental brain disorders with complex and largely unknown etiology. To characterize the impact of schizophrenia at a cellular level, we performed single nucleus RNA sequencing of >190,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls (7 vs 11, respectively). In addition, to correlate data with cortical anatomy, >100,000 neurons were analyzed topographically by immunohistochemistry in an extended cohort of cases with schizophrenia and controls (10 vs 10). Compositional analysis of RNA sequencing data revealed reduction in relative abundance across all families of GABAergic neurons and a concomitant increase in principal neurons, which was most pronounced for supragranular subtypes (layers 2-3). Moreover, supragranular subtypes of GABAergic interneurons showed most dramatic transcriptomic changes. These results were substantiated by histological analysis, which revealed a reduction in the density of calretinin, calbindin and parvalbumin GABAergic interneurons particularly in layer 2. Common effect of schizophrenia on supragranular neuronal networks was underlined by downregulation of protein processing genes and upregulation of neuronal development/plasticity genes across supragranular subtypes of principal neurons and GABAergic interneurons. In situ hybridization and spatial transcriptomics further confirmed supragranular layer neuron vulnerability, revealing complexity of schizophrenia-affected cortical circuits. These point towards general network impairment within supragranular layers being a core substrate associated with schizophrenia symptomatology.


title: Cis regulation within a cluster of viral microRNAs

abstract: MicroRNAs (miRNAs) are small regulatory RNAs involved in virtually all biological processes. Although many of them are co-expressed from clusters, little is known regarding the impact of this organization on the regulation of their accumulation. In this study, we set to decipher a regulatory mechanism controlling the expression of the ten clustered pre-miRNAs from Kaposis sarcoma associated herpesvirus (KSHV). We measured in vitro the efficiency of cleavage of each individual pre-miRNA by the Microprocessor and found that pre-miR-K1 and -K3 were the most efficiently cleaved pre-miRNAs. A mutational analysis showed that, in addition to producing mature miRNAs, they are also important for the optimal expression of the whole set of miRNAs. We showed that this feature depends on the presence of a canonical pre-miRNA at this location since we could functionally replace pre-miR-K1 by a heterologous pre-miRNA. Further in vitro processing analysis suggests that the two stem-loops act in cis and that the cluster is cleaved in a sequential manner. Finally, we exploited this characteristic of the cluster to inhibit the expression of the whole set of miRNAs by targeting the pre-miR-K1 with LNA-based antisense oligonucleotides in cells either expressing a synthetic construct or latently infected with KSHV.


title: Ultrafast pore-loop dynamics in a AAA+ machine point to a Brownian-ratchet mechanism for protein translocation

abstract: AAA+ ring-shaped machines, such as ClpB and Hsp104, mediate substrate translocation through their central channel by a set of pore loops. Recent structural studies suggested a universal hand-over-hand translocation mechanism, in which pore loops are moving rigidly in tandem with their corresponding subunits. However, functional and biophysical studies are in discord with this model. Here, we directly measure the real-time dynamics of the pore loops of ClpB and their response to substrate binding, using single-molecule FRET spectroscopy. All pore loops undergo large-amplitude fluctuations on the microsecond timescale, and change their conformation upon interaction with substrate proteins. Pore-loop conformational dynamics are modulated by nucleotides and strongly correlate with disaggregation activity. The differential behavior of the pore loops along the axial channel points to a fast Brownian-ratchet translocation mechanism, which likely acts in parallel to the much slower hand-over-hand process.


title: Long-term stability of avalanche scaling and integrative network organization in prefrontal and premotor cortex

abstract: Ongoing neuronal activity in the cortex establishes functional networks of synchronization that reflect normal and pathological brain function. The reconstruction of these networks typically suffers from the use of indirect measures of neuronal activity at low spatiotemporal resolution and a lack of longitudinal tracking. Accordingly, the precise nature of the underlying synchronization dynamics and its translation into robust graph theoretical markers are not well characterized. Here, we studied the stability of cortical dynamics and reconstructed functional networks over many weeks in prefrontal and premotor cortex of awake nonhuman primates. We monitored neuronal population activity directly in the ongoing local field potential (LFP) at high spatial and temporal resolution using chronically implanted high-density microelectrode arrays. Ongoing activity was composed of neuronal avalanches exhibiting stable, inverted parabolic profiles with the collapse exponent of 2 in line with a critical branching process. Avalanche-based functional networks, reconstructed using a Normalized Count estimator, revealed robust integrative properties characterized by high neighborhood overlap between strongly connected nodes and robustness to weak-link pruning. "Entropy of mixing" analysis demonstrated progressive link reorganization over weeks. The long-term stability of avalanche scaling and integrative network organization should support the development of robust biomarkers to characterize normal and abnormal brain function.


title: The road to sorghum domestication: evidence from nucleotide diversity and gene expression patterns

abstract: Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and 9 wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a genetic bottleneck in the domestication history of sorghum, although several genes showed also patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyses revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totalizing 2449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that domestication and improvement do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome during the domestication event and the improvement process. Overall, these analyses pave the way towards the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

category: evolutionary biology